src/math/Matrix3x3d.hpp

/*
 * Copyright (C) 2000-2004  Object Oriented Parallel Simulation Engine (OOPSE) project
 * 
 * Contact: oopse@oopse.org
 * 
 * This program is free software; you can redistribute it and/or
 * modify it under the terms of the GNU Lesser General Public License
 * as published by the Free Software Foundation; either version 2.1
 * of the License, or (at your option) any later version.
 * All we ask is that proper credit is given for our work, which includes
 * - but is not limited to - adding the above copyright notice to the beginning
 * of your source code files, and to any copyright notice that you may distribute
 * with programs based on this work.
 * 
 * This program is distributed in the hope that it will be useful,
 * but WITHOUT ANY WARRANTY; without even the implied warranty of
 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
 * GNU Lesser General Public License for more details.
 * 
 * You should have received a copy of the GNU Lesser General Public License
 * along with this program; if not, write to the Free Software
 * Foundation, Inc., 59 Temple Place - Suite 330, Boston, MA  02111-1307, USA.
 *
 */

/**
 * @file Matrix3x3d.hpp
 * @author Teng Lin
 * @date 10/11/2004
 * @version 1.0
 */
#ifndef MATH_MATRIX3X3D_HPP 
#define MATH_MATRIX3X3D_HPP 

#include "Vector3d.hpp"

namespace oopse {

    /**
     * @class Matrix3x3d Matrix3x3d.hpp "math/Matrix3x3d.hpp"
     * @brief A 3x3 matrix class
     */
    class Matrix3x3d{
        public:

        /** default constructor */
        Matrix3x3d();

        /** Constructs and initializes every element of this matrix to a scalar */ 
        Matrix3x3d(double s);
        
        /**
         * Return the reference of a single element of this matrix.
         * @return the reference of a single element of this matrix 
         * @param i row index
         * @param j colum index
         */
        double& operator()(unsigned int i, unsigned int j) {  return data_[i][j];  }

        /**
         * Return the value of a single element of this matrix.
         * @return the value of a single element of this matrix 
         * @param i row index
         * @param j colum index
         */        
        double operator()(unsigned int i, unsigned int j) const  {  return data_[i][j];  }

        /**
         * Returns a row of  this matrix as a vector.
         * @return a row of  this matrix as a vector 
         * @param i the row index
         */                
        Vector3d getRow(unsigned int i);

        /**
         * Sets a row of  this matrix
         * @param i the row index
         * @param v the vector to be set
         */                
         void setRow(unsigned int i, const Vector3d& v);

        /**
         * Returns a column of  this matrix as a vector.
         * @return a column of  this matrix as a vector 
         * @param i the column index
         */                
        Vector3d getColum(unsigned int i);

        /**
         * Sets a column of  this matrix
         * @param i the column index
         * @param v the vector to be set
         */                
         void setColum(unsigned int i, const Vector3d& v);

        /**
         *  Adds a scalar into every element of this matrix.
         * @param s the scalar value
         */
        Matrix3x3d& operator +=(const double s);

        /**
         *  Subtracts a scalar from every element of this matrix.
         * @param s the scalar value
         */
        Matrix3x3d& operator -=(const double s);

        /**
         *  Multiples a scalar into every element of this matrix.
         * @param s the scalar value
         */
        Matrix3x3d& operator *=(const double s);

        Matrix3x3d& operator /=(const double s);

        /**
         * Sets the value of this matrix to the sum of the other matrix and itself (*this += m).
         * @param m the other matrix
         */
        Matrix3x3d& operator += (const Matrix3x3d& m);

        /**
         * Sets the value of this matrix to the differerence of itself and the other matrix (*this -= m) 
         * @param m the other matrix
         */
        Matrix3x3d& operator -= (const Matrix3x3d& m);

        /** Returns the inverse of this matrix. */
        Matrix3x3d inverse();

        /** Returns the transpose of this matrix. */
        Matrix3x3d transpose();

        /** Returns the determinant of this matrix. */
        double determinant() const;

        /** Returns the trace of this matrix. */
        double trace() const;

        /** Tests if this matrix is symmetrix. */            
        bool isSymmetric() const;

        /** Tests if this matrix is orthogona. */            
        bool isOrthogonal() const;

        /** Tests if this matrix is diagonal. */
        bool isDiagonal() const;

        /** Tests if this matrix is the unit matrix. */
        bool isUnitMatrix() const;
        
        /** Negate the value of every element of this matrix. */
        friend inline Matrix3x3d operator -(const Matrix3x3d& m);
        
        /**
        * Return the sum of two matrixes  (m1 + m2). 
        * @return the sum of two matrixes
        * @param m1 the first matrix
        * @param m2 the second matrix
        */ 
        friend inline Matrix3x3d operator + (const Matrix3x3d& m1, const Matrix3x3d& m2);

        /**
        * Return the difference of two matrixes  (m1 - m2). 
        * @return the sum of two matrixes
        * @param m1 the first matrix
        * @param m2 the second matrix
        */
        friend inline Matrix3x3d operator - (const Matrix3x3d& m1, const Matrix3x3d& m2);

        /**
        * Return the multiplication of two matrixes  (m1 * m2). 
        * @return the multiplication of two matrixes
        * @param m1 the first matrix
        * @param m2 the second matrix
        */
        friend inline Matrix3x3d operator * (const Matrix3x3d& m1, const Matrix3x3d& m2);

        /**
        * Return the multiplication of  matrixes and vector (m1 * v1). 
        * @return the multiplication of matrixes and vector
        * @param m1 the matrix
        * @param v1 the vector
        */
        friend inline Matrix3x3d operator * (const Matrix3x3d& m1, const Matrix3x3d& v1);
        
        protected:
            double data_[3][3]; /**< matrix element */            

    }
}
#endif //MATH_MATRIX3X3D_HPP 
Revision:	1561
Committed:	Tue Oct 12 23:57:40 2004 UTC (21 years, 1 month ago) by tim
File size:	6379 byte(s)
Log Message:	add Vector3d.hpp and Matrix3x3d.hpp
#	Content
1	/*
2	* Copyright (C) 2000-2004 Object Oriented Parallel Simulation Engine (OOPSE) project
3	*
4	* Contact: oopse@oopse.org
5	*
6	* This program is free software; you can redistribute it and/or
7	* modify it under the terms of the GNU Lesser General Public License
8	* as published by the Free Software Foundation; either version 2.1
9	* of the License, or (at your option) any later version.
10	* All we ask is that proper credit is given for our work, which includes
11	* - but is not limited to - adding the above copyright notice to the beginning
12	* of your source code files, and to any copyright notice that you may distribute
13	* with programs based on this work.
14	*
15	* This program is distributed in the hope that it will be useful,
16	* but WITHOUT ANY WARRANTY; without even the implied warranty of
17	* MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
18	* GNU Lesser General Public License for more details.
19	*
20	* You should have received a copy of the GNU Lesser General Public License
21	* along with this program; if not, write to the Free Software
22	* Foundation, Inc., 59 Temple Place - Suite 330, Boston, MA 02111-1307, USA.
23	*
24	*/
25
26	/**
27	* @file Matrix3x3d.hpp
28	* @author Teng Lin
29	* @date 10/11/2004
30	* @version 1.0
31	*/
32	#ifndef MATH_MATRIX3X3D_HPP
33	#define MATH_MATRIX3X3D_HPP
34
35	#include "Vector3d.hpp"
36
37	namespace oopse {
38
39	/**
40	* @class Matrix3x3d Matrix3x3d.hpp "math/Matrix3x3d.hpp"
41	* @brief A 3x3 matrix class
42	*/
43	class Matrix3x3d{
44	public:
45
46	/** default constructor */
47	Matrix3x3d();
48
49	/** Constructs and initializes every element of this matrix to a scalar */
50	Matrix3x3d(double s);
51
52	/**
53	* Return the reference of a single element of this matrix.
54	* @return the reference of a single element of this matrix
55	* @param i row index
56	* @param j colum index
57	*/
58	double& operator()(unsigned int i, unsigned int j) { return data_[i][j]; }
59
60	/**
61	* Return the value of a single element of this matrix.
62	* @return the value of a single element of this matrix
63	* @param i row index
64	* @param j colum index
65	*/
66	double operator()(unsigned int i, unsigned int j) const { return data_[i][j]; }
67
68	/**
69	* Returns a row of this matrix as a vector.
70	* @return a row of this matrix as a vector
71	* @param i the row index
72	*/
73	Vector3d getRow(unsigned int i);
74
75	/**
76	* Sets a row of this matrix
77	* @param i the row index
78	* @param v the vector to be set
79	*/
80	void setRow(unsigned int i, const Vector3d& v);
81
82	/**
83	* Returns a column of this matrix as a vector.
84	* @return a column of this matrix as a vector
85	* @param i the column index
86	*/
87	Vector3d getColum(unsigned int i);
88
89	/**
90	* Sets a column of this matrix
91	* @param i the column index
92	* @param v the vector to be set
93	*/
94	void setColum(unsigned int i, const Vector3d& v);
95
96	/**
97	* Adds a scalar into every element of this matrix.
98	* @param s the scalar value
99	*/
100	Matrix3x3d& operator +=(const double s);
101
102	/**
103	* Subtracts a scalar from every element of this matrix.
104	* @param s the scalar value
105	*/
106	Matrix3x3d& operator -=(const double s);
107
108	/**
109	* Multiples a scalar into every element of this matrix.
110	* @param s the scalar value
111	*/
112	Matrix3x3d& operator *=(const double s);
113
114	Matrix3x3d& operator /=(const double s);
115
116	/**
117	* Sets the value of this matrix to the sum of the other matrix and itself (*this += m).
118	* @param m the other matrix
119	*/
120	Matrix3x3d& operator += (const Matrix3x3d& m);
121
122	/**
123	* Sets the value of this matrix to the differerence of itself and the other matrix (*this -= m)
124	* @param m the other matrix
125	*/
126	Matrix3x3d& operator -= (const Matrix3x3d& m);
127
128	/** Returns the inverse of this matrix. */
129	Matrix3x3d inverse();
130
131	/** Returns the transpose of this matrix. */
132	Matrix3x3d transpose();
133
134	/** Returns the determinant of this matrix. */
135	double determinant() const;
136
137	/** Returns the trace of this matrix. */
138	double trace() const;
139
140	/** Tests if this matrix is symmetrix. */
141	bool isSymmetric() const;
142
143	/** Tests if this matrix is orthogona. */
144	bool isOrthogonal() const;
145
146	/** Tests if this matrix is diagonal. */
147	bool isDiagonal() const;
148
149	/** Tests if this matrix is the unit matrix. */
150	bool isUnitMatrix() const;
151
152	/** Negate the value of every element of this matrix. */
153	friend inline Matrix3x3d operator -(const Matrix3x3d& m);
154
155	/**
156	* Return the sum of two matrixes (m1 + m2).
157	* @return the sum of two matrixes
158	* @param m1 the first matrix
159	* @param m2 the second matrix
160	*/
161	friend inline Matrix3x3d operator + (const Matrix3x3d& m1, const Matrix3x3d& m2);
162
163	/**
164	* Return the difference of two matrixes (m1 - m2).
165	* @return the sum of two matrixes
166	* @param m1 the first matrix
167	* @param m2 the second matrix
168	*/
169	friend inline Matrix3x3d operator - (const Matrix3x3d& m1, const Matrix3x3d& m2);
170
171	/**
172	* Return the multiplication of two matrixes (m1 * m2).
173	* @return the multiplication of two matrixes
174	* @param m1 the first matrix
175	* @param m2 the second matrix
176	*/
177	friend inline Matrix3x3d operator * (const Matrix3x3d& m1, const Matrix3x3d& m2);
178
179	/**
180	* Return the multiplication of matrixes and vector (m1 * v1).
181	* @return the multiplication of matrixes and vector
182	* @param m1 the matrix
183	* @param v1 the vector
184	*/
185	friend inline Matrix3x3d operator * (const Matrix3x3d& m1, const Matrix3x3d& v1);
186
187	protected:
188	double data_[3][3]; /*< matrix element /
189
190	}
191	}
192	#endif //MATH_MATRIX3X3D_HPP