test/math/RectMatrixTestCase.cpp

#include "RectMatrixTestCase.hpp"

// Registers the fixture into the 'registry'
CPPUNIT_TEST_SUITE_REGISTRATION( RectMatrixTestCase );

void RectMatrixTestCase::setUp(){

    m1(0, 0) = 1.0;
    m1(0, 1) = 1.0;
    m1(1, 0) = 3.0;
    m1(1, 1) = 1.0;
    
    m2(0, 0) = 1.0;
    m2(0, 1) = 0.0;
    m2(1, 0) = 0.0;
    m2(1, 1) = 1.0;
 
    
    m3(0, 0) = 1.0;
    m3(0, 1) = 0.0;
    m3(1, 0) = 0.0;
    m3(1, 1) = 1.0;

    m4(0, 0) = -1.0;
    m4(0, 1) = -1.0;
    m4(1, 0) = -3.0;
    m4(1, 1) = -1.0;

    zero(0, 0) = 0.0;
    zero(0, 1) = 0.0;
    zero(1, 0) = 0.0;
    zero(1, 1) = 0.0;    

    one(0, 0) = 1.0;
    one(0, 1) = 1.0;
    one(1, 0) = 1.0;
    one(1, 1) = 1.0;    

    two(0, 0) = 2.0;
    two(0, 1) = 2.0;
    two(1, 0) = 2.0;
    two(1, 1) = 2.0;  
    
    a(0, 0) = 1.0;
    a(0, 1) = 0.0;
    a(0, 2) = 0.0;
    a(1, 0) = 0.0;
    a(1, 1) = 1.0;
    a(1, 2) = 0.0;

    b(0, 0) = 1.0;
    b(0, 1) = 0.0;
    b(1, 0) = 0.0;
    b(1, 1) = 1.0;
    b(2, 0) = 0.0;
    b(2, 1) = 0.0;

    c(0, 0) = 1.0;
    c(0, 1) = 0.0;
    c(1, 0) = 0.0;
    c(1, 1) = 1.0;       

    d(0, 0) = 1.0;
    d(0, 1) = 0.0;
    d(0, 2) = 0.0;
    d(1, 0) = 0.0;
    d(1, 1) = 0.0;
    d(1, 2) = 1.0;    
    d(2, 0) = 0.0;
    d(2, 1) = 1.0;
    d(2, 2) = 0.0;  

    e(0, 0) = 2.0;
    e(0, 1) = 4.0;
    e(0, 2) = 1.0;
    e(1, 0) = 0.0;
    e(1, 1) = 0.0;
    e(1, 2) = 3.0;    
    e(2, 0) = 0.0;
    e(2, 1) = 6.0;
    e(2, 2) = 5.0;  

    f(0, 0) = 2.0;
    f(0, 1) = 4.0;
    f(0, 2) = 1.0;
    f(1, 0) = 0.0;
    f(1, 1) = 6.0;
    f(1, 2) = 5.0;    
    f(2, 0) = 0.0;
    f(2, 1) = 0.0;
    f(2, 2) = 3.0;  

    f(0, 0) = 2.0;
    f(0, 1) = 4.0;
    f(0, 2) = 1.0;
    f(1, 0) = 0.0;
    f(1, 1) = 6.0;
    f(1, 2) = 5.0;    
    f(2, 0) = 0.0;
    f(2, 1) = 0.0;
    f(2, 2) = 3.0;  

    g(0, 0) = 1.0;
    g(0, 1) = 0.0;
    g(0, 2) = 0.0;
    g(1, 0) = -2.0;
    g(1, 1) = 1.0;
    g(1, 2) = 0.0;    
    g(2, 0) = 0.0;
    g(2, 1) = 0.0;
    g(2, 2) = 1.0;  

    h(0, 0) = 2.0;
    h(0, 1) = 4.0;
    h(0, 2) = -2.0;
    h(0, 3) = 2.0;
    h(1, 0) = 4.0;
    h(1, 1) = 9.0;
    h(1, 2) = -3.0;
    h(1, 3) = 8.0;
    h(2, 0) = -2.0;
    h(2, 1) = -3.0;
    h(2, 2) = 7.0;
    h(2, 3) = 10.0;

    i(0, 0) = 2.0;
    i(0, 1) = 4.0;
    i(0, 2) = -2.0;
    i(0, 3) = 2.0;
    i(1, 0) = 0.0;
    i(1, 1) = 1.0;
    i(1, 2) = 1.0;
    i(1, 3) = 4.0;
    i(2, 0) = -2.0;
    i(2, 1) = -3.0;
    i(2, 2) = 7.0;
    i(2, 3) = 10.0;    


    v1(0) = 2.0;
    v1(1) = 4.0;
    v1(2) = -2.0;

    v2(0) = 2.0;
    v2(1) = 0.0;
    v2(2) = -2.0;

}

void RectMatrixTestCase::testConstructor(){

    //test default constructor
    RMat2x2 tmp1;
    CPPUNIT_ASSERT(tmp1 == zero);

    //test RectMatrix(Real s)
    RMat2x2 tmp2;
    CPPUNIT_ASSERT(tmp2 == zero);

    //test copy constructor
    RMat2x2 tmp3(m1);
    CPPUNIT_ASSERT(tmp3 == m1);

    //test copy assignment
    RMat2x2 tmp4 = m2;
    tmp4 = tmp4;
    CPPUNIT_ASSERT(tmp4 == m2);
    
    
}

void RectMatrixTestCase::testEqual() {
    CPPUNIT_ASSERT(m2 == m3);
    CPPUNIT_ASSERT(m2 != m3);
    
}

void RectMatrixTestCase::testNegate() {

    CPPUNIT_ASSERT(m1 == -m4);
}

void RectMatrixTestCase::testAdd() {
    RMat2x2 tmp;
    
    CPPUNIT_ASSERT(m1 + m4 == zero);

    tmp.add(m1, m1);

    CPPUNIT_ASSERT(m1 * 2.0 == tmp);

    tmp = one;
    tmp *= 2.0;

    CPPUNIT_ASSERT(tmp == two);
    
}

void RectMatrixTestCase::testSub() {
    
    RMat2x2 tmp(m2);
    tmp.sub(m2);
    CPPUNIT_ASSERT(tmp == zero);

    tmp = m2;
    tmp -= m2;
    CPPUNIT_ASSERT(tmp == zero);    

    tmp = m1;
    tmp.sub(m1 , m4);
    CPPUNIT_ASSERT(tmp == m1 * 2.0);
    
    CPPUNIT_ASSERT(m1 -m4 == m1 * 2.0);

    
}

void RectMatrixTestCase::testMul() {

    CPPUNIT_ASSERT(m1 * 1.0 == m1);
    CPPUNIT_ASSERT(m1 * 0.0 == zero);
    CPPUNIT_ASSERT(2.0 *m1 == m1 + m1);

    //test matrix multiplication
    CPPUNIT_ASSERT(a * b == c);
    CPPUNIT_ASSERT(d  * e == f);
    CPPUNIT_ASSERT(g  * h == i);

    //test matrix vector multiplication
    CPPUNIT_ASSERT(g  * v1 == v2);

}

void RectMatrixTestCase::testDiv() {

    CPPUNIT_ASSERT(m1 / 1.0 == m1);

    CPPUNIT_ASSERT(m1 / 2.0 * 2.0 == m1);


}

void RectMatrixTestCase::testAccessEntries(){
    CPPUNIT_ASSERT_DOUBLES_EQUAL(m1(1, 0), 3.0, oopse::epsilon);
}

void RectMatrixTestCase::testRowColOperations() {
    Vec3 row;
    Vec3 col;
    RMat3x3 m;
    
    //test getRow
    row = e.getRow(0);
    CPPUNIT_ASSERT_DOUBLES_EQUAL(row[0], 2.0, oopse::epsilon);
    CPPUNIT_ASSERT_DOUBLES_EQUAL(row[1], 4.0, oopse::epsilon);
    CPPUNIT_ASSERT_DOUBLES_EQUAL(row[2], 1.0, oopse::epsilon);
    //test setRow
    row[0] = 2.0;    
    row[1] = 4.0;    
    row[2] = 1.0;    
    m.setRow(0, row);
    row[0] = 0.0;    
    row[1] = 0.0;    
    row[2] = 3.0;    
    m.setRow(1, row);
    row[0] = 0.0;    
    row[1] = 6.0;    
    row[2] = 5.0;    
    m.setRow(2, row);
    CPPUNIT_ASSERT(m == e);
    
    //test getCol
    col = e.getColum(1);
    CPPUNIT_ASSERT_DOUBLES_EQUAL(col[0], 4.0, oopse::epsilon);    
    CPPUNIT_ASSERT_DOUBLES_EQUAL(col[1], 0.0, oopse::epsilon);
    CPPUNIT_ASSERT_DOUBLES_EQUAL(col[2], 6.0, oopse::epsilon);
    //test setCol
    col[0] = 2.0;    
    col[1] = 0.0;    
    col[2] = 0.0;    
    m.setColum(0, col);
    col[0] = 4.0;    
    col[1] = 0.0;    
    col[2] = 6.0;    
    m.setColum(1, col);
    col[0] = 1.0;    
    col[1] = 3.0;    
    col[2] = 5.0;    
    m.setColum(2, col);
    CPPUNIT_ASSERT(m == e);

    //test swapRow
    RMat2x3 r;
    r(0, 0) = 0.0;
    r(0, 1) = 1.0;
    r(0, 2) = 0.0;
    r(1, 0) = 1.0;
    r(1, 1) = 0.0;
    r(1, 2) = 0.0;
    r.swapRow(0, 1);
    CPPUNIT_ASSERT(r == a);

    //test swapCol
    RMat3x3 s;
    s(0, 0) = 4.0;
    s(0, 1) = 2.0;
    s(0, 2) = 1.0;
    s(1, 0) = 0.0;    
    s(1, 1) = 0.0;
    s(1, 2) = 3.0;    
    s(2, 0) = 6.0;
    s(2, 1) = 0.0;
    s(2, 2) = 5.0;

    s.swapColum(0, 1);
    CPPUNIT_ASSERT(s == e);

    double* p = s.getArrayPointer();

    p[0] = 2.0;
    p[1] = 4.0;
    p[2] = 1.0;
    p[3] = 0.0;
    p[4] = 6.0;
    p[5] = 5.0;    
    p[6] = 0.0;
    p[7] = 0.0;
    p[8] = 3.0;  

    CPPUNIT_ASSERT(s == f);    
}    

void RectMatrixTestCase::testOtherMemberFunctions(){
    //test transpose
    CPPUNIT_ASSERT((a.transpose()).transpose() == a);
    
    CPPUNIT_ASSERT(a.transpose() == b);
}
Revision:	1630
Committed:	Thu Oct 21 21:31:39 2004 UTC (19 years, 8 months ago) by tim
File size:	6357 byte(s)
Log Message:	Snapshot and SnapshotManager in progress
#	User	Rev	Content
1	tim	1565	#include "RectMatrixTestCase.hpp"
2
3			// Registers the fixture into the 'registry'
4			CPPUNIT_TEST_SUITE_REGISTRATION( RectMatrixTestCase );
5
6			void RectMatrixTestCase::setUp(){
7
8			m1(0, 0) = 1.0;
9			m1(0, 1) = 1.0;
10			m1(1, 0) = 3.0;
11			m1(1, 1) = 1.0;
12
13			m2(0, 0) = 1.0;
14			m2(0, 1) = 0.0;
15			m2(1, 0) = 0.0;
16			m2(1, 1) = 1.0;
17
18
19			m3(0, 0) = 1.0;
20			m3(0, 1) = 0.0;
21			m3(1, 0) = 0.0;
22			m3(1, 1) = 1.0;
23
24			m4(0, 0) = -1.0;
25			m4(0, 1) = -1.0;
26			m4(1, 0) = -3.0;
27			m4(1, 1) = -1.0;
28
29			zero(0, 0) = 0.0;
30			zero(0, 1) = 0.0;
31			zero(1, 0) = 0.0;
32			zero(1, 1) = 0.0;
33
34			one(0, 0) = 1.0;
35			one(0, 1) = 1.0;
36			one(1, 0) = 1.0;
37			one(1, 1) = 1.0;
38
39			two(0, 0) = 2.0;
40			two(0, 1) = 2.0;
41			two(1, 0) = 2.0;
42			two(1, 1) = 2.0;
43
44			a(0, 0) = 1.0;
45			a(0, 1) = 0.0;
46			a(0, 2) = 0.0;
47			a(1, 0) = 0.0;
48			a(1, 1) = 1.0;
49			a(1, 2) = 0.0;
50
51			b(0, 0) = 1.0;
52			b(0, 1) = 0.0;
53			b(1, 0) = 0.0;
54			b(1, 1) = 1.0;
55			b(2, 0) = 0.0;
56			b(2, 1) = 0.0;
57
58	tim	1571	c(0, 0) = 1.0;
59			c(0, 1) = 0.0;
60			c(1, 0) = 0.0;
61			c(1, 1) = 1.0;
62
63			d(0, 0) = 1.0;
64			d(0, 1) = 0.0;
65			d(0, 2) = 0.0;
66			d(1, 0) = 0.0;
67			d(1, 1) = 0.0;
68			d(1, 2) = 1.0;
69			d(2, 0) = 0.0;
70			d(2, 1) = 1.0;
71			d(2, 2) = 0.0;
72
73			e(0, 0) = 2.0;
74			e(0, 1) = 4.0;
75			e(0, 2) = 1.0;
76			e(1, 0) = 0.0;
77			e(1, 1) = 0.0;
78			e(1, 2) = 3.0;
79			e(2, 0) = 0.0;
80			e(2, 1) = 6.0;
81			e(2, 2) = 5.0;
82
83			f(0, 0) = 2.0;
84			f(0, 1) = 4.0;
85			f(0, 2) = 1.0;
86			f(1, 0) = 0.0;
87			f(1, 1) = 6.0;
88			f(1, 2) = 5.0;
89			f(2, 0) = 0.0;
90			f(2, 1) = 0.0;
91			f(2, 2) = 3.0;
92
93			f(0, 0) = 2.0;
94			f(0, 1) = 4.0;
95			f(0, 2) = 1.0;
96			f(1, 0) = 0.0;
97			f(1, 1) = 6.0;
98			f(1, 2) = 5.0;
99			f(2, 0) = 0.0;
100			f(2, 1) = 0.0;
101			f(2, 2) = 3.0;
102
103			g(0, 0) = 1.0;
104			g(0, 1) = 0.0;
105			g(0, 2) = 0.0;
106			g(1, 0) = -2.0;
107			g(1, 1) = 1.0;
108			g(1, 2) = 0.0;
109			g(2, 0) = 0.0;
110			g(2, 1) = 0.0;
111			g(2, 2) = 1.0;
112
113			h(0, 0) = 2.0;
114			h(0, 1) = 4.0;
115			h(0, 2) = -2.0;
116			h(0, 3) = 2.0;
117			h(1, 0) = 4.0;
118			h(1, 1) = 9.0;
119			h(1, 2) = -3.0;
120			h(1, 3) = 8.0;
121			h(2, 0) = -2.0;
122			h(2, 1) = -3.0;
123			h(2, 2) = 7.0;
124			h(2, 3) = 10.0;
125
126			i(0, 0) = 2.0;
127			i(0, 1) = 4.0;
128			i(0, 2) = -2.0;
129			i(0, 3) = 2.0;
130			i(1, 0) = 0.0;
131			i(1, 1) = 1.0;
132			i(1, 2) = 1.0;
133			i(1, 3) = 4.0;
134			i(2, 0) = -2.0;
135			i(2, 1) = -3.0;
136			i(2, 2) = 7.0;
137			i(2, 3) = 10.0;
138
139
140			v1(0) = 2.0;
141			v1(1) = 4.0;
142			v1(2) = -2.0;
143
144			v2(0) = 2.0;
145			v2(1) = 0.0;
146			v2(2) = -2.0;
147
148	tim	1565	}
149
150			void RectMatrixTestCase::testConstructor(){
151
152			//test default constructor
153			RMat2x2 tmp1;
154			CPPUNIT_ASSERT(tmp1 == zero);
155
156			//test RectMatrix(Real s)
157			RMat2x2 tmp2;
158			CPPUNIT_ASSERT(tmp2 == zero);
159
160			//test copy constructor
161			RMat2x2 tmp3(m1);
162			CPPUNIT_ASSERT(tmp3 == m1);
163
164			//test copy assignment
165			RMat2x2 tmp4 = m2;
166			tmp4 = tmp4;
167			CPPUNIT_ASSERT(tmp4 == m2);
168
169
170			}
171
172			void RectMatrixTestCase::testEqual() {
173			CPPUNIT_ASSERT(m2 == m3);
174	tim	1571	CPPUNIT_ASSERT(m2 != m3);
175
176	tim	1565	}
177
178			void RectMatrixTestCase::testNegate() {
179
180			CPPUNIT_ASSERT(m1 == -m4);
181			}
182
183			void RectMatrixTestCase::testAdd() {
184			RMat2x2 tmp;
185
186			CPPUNIT_ASSERT(m1 + m4 == zero);
187
188			tmp.add(m1, m1);
189
190			CPPUNIT_ASSERT(m1 * 2.0 == tmp);
191
192			tmp = one;
193			tmp *= 2.0;
194
195			CPPUNIT_ASSERT(tmp == two);
196
197			}
198
199			void RectMatrixTestCase::testSub() {
200
201			RMat2x2 tmp(m2);
202			tmp.sub(m2);
203			CPPUNIT_ASSERT(tmp == zero);
204
205			tmp = m2;
206			tmp -= m2;
207			CPPUNIT_ASSERT(tmp == zero);
208
209			tmp = m1;
210			tmp.sub(m1 , m4);
211			CPPUNIT_ASSERT(tmp == m1 * 2.0);
212
213			CPPUNIT_ASSERT(m1 -m4 == m1 * 2.0);
214
215
216			}
217
218			void RectMatrixTestCase::testMul() {
219
220			CPPUNIT_ASSERT(m1 * 1.0 == m1);
221			CPPUNIT_ASSERT(m1 * 0.0 == zero);
222			CPPUNIT_ASSERT(2.0 *m1 == m1 + m1);
223	tim	1571
224			//test matrix multiplication
225			CPPUNIT_ASSERT(a * b == c);
226			CPPUNIT_ASSERT(d * e == f);
227			CPPUNIT_ASSERT(g * h == i);
228
229			//test matrix vector multiplication
230			CPPUNIT_ASSERT(g * v1 == v2);
231
232	tim	1565	}
233
234			void RectMatrixTestCase::testDiv() {
235
236			CPPUNIT_ASSERT(m1 / 1.0 == m1);
237
238			CPPUNIT_ASSERT(m1 / 2.0 * 2.0 == m1);
239
240	tim	1603
241	tim	1565	}
242
243			void RectMatrixTestCase::testAccessEntries(){
244	tim	1571	CPPUNIT_ASSERT_DOUBLES_EQUAL(m1(1, 0), 3.0, oopse::epsilon);
245	tim	1565	}
246
247	tim	1595	void RectMatrixTestCase::testRowColOperations() {
248			Vec3 row;
249			Vec3 col;
250			RMat3x3 m;
251
252			//test getRow
253			row = e.getRow(0);
254	tim	1603	CPPUNIT_ASSERT_DOUBLES_EQUAL(row[0], 2.0, oopse::epsilon);
255			CPPUNIT_ASSERT_DOUBLES_EQUAL(row[1], 4.0, oopse::epsilon);
256			CPPUNIT_ASSERT_DOUBLES_EQUAL(row[2], 1.0, oopse::epsilon);
257	tim	1595	//test setRow
258			row[0] = 2.0;
259			row[1] = 4.0;
260			row[2] = 1.0;
261			m.setRow(0, row);
262			row[0] = 0.0;
263			row[1] = 0.0;
264			row[2] = 3.0;
265			m.setRow(1, row);
266			row[0] = 0.0;
267			row[1] = 6.0;
268			row[2] = 5.0;
269			m.setRow(2, row);
270			CPPUNIT_ASSERT(m == e);
271
272			//test getCol
273	tim	1597	col = e.getColum(1);
274	tim	1603	CPPUNIT_ASSERT_DOUBLES_EQUAL(col[0], 4.0, oopse::epsilon);
275			CPPUNIT_ASSERT_DOUBLES_EQUAL(col[1], 0.0, oopse::epsilon);
276			CPPUNIT_ASSERT_DOUBLES_EQUAL(col[2], 6.0, oopse::epsilon);
277	tim	1595	//test setCol
278			col[0] = 2.0;
279			col[1] = 0.0;
280			col[2] = 0.0;
281	tim	1596	m.setColum(0, col);
282	tim	1595	col[0] = 4.0;
283			col[1] = 0.0;
284			col[2] = 6.0;
285	tim	1596	m.setColum(1, col);
286	tim	1595	col[0] = 1.0;
287			col[1] = 3.0;
288			col[2] = 5.0;
289	tim	1596	m.setColum(2, col);
290	tim	1595	CPPUNIT_ASSERT(m == e);
291
292			//test swapRow
293			RMat2x3 r;
294			r(0, 0) = 0.0;
295			r(0, 1) = 1.0;
296			r(0, 2) = 0.0;
297			r(1, 0) = 1.0;
298			r(1, 1) = 0.0;
299			r(1, 2) = 0.0;
300			r.swapRow(0, 1);
301			CPPUNIT_ASSERT(r == a);
302
303			//test swapCol
304			RMat3x3 s;
305			s(0, 0) = 4.0;
306			s(0, 1) = 2.0;
307			s(0, 2) = 1.0;
308			s(1, 0) = 0.0;
309			s(1, 1) = 0.0;
310			s(1, 2) = 3.0;
311			s(2, 0) = 6.0;
312			s(2, 1) = 0.0;
313			s(2, 2) = 5.0;
314
315	tim	1596	s.swapColum(0, 1);
316	tim	1595	CPPUNIT_ASSERT(s == e);
317	tim	1630
318			double* p = s.getArrayPointer();
319
320			p[0] = 2.0;
321			p[1] = 4.0;
322			p[2] = 1.0;
323			p[3] = 0.0;
324			p[4] = 6.0;
325			p[5] = 5.0;
326			p[6] = 0.0;
327			p[7] = 0.0;
328			p[8] = 3.0;
329
330			CPPUNIT_ASSERT(s == f);
331	tim	1595	}
332
333			void RectMatrixTestCase::testOtherMemberFunctions(){
334			//test transpose
335	tim	1565	CPPUNIT_ASSERT((a.transpose()).transpose() == a);
336
337			CPPUNIT_ASSERT(a.transpose() == b);
338			}