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(row[0] == 2.0 && row[1] == 4.0 && row[2] == 1.0)

    //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.getCol(1);
    CPPUNIT_ASSERT(col[0] == 4.0 && col[1] == 0.0 && col[2] == 6.0)    

    //test setCol
    col[0] = 2.0;    
    col[1] = 0.0;    
    col[2] = 0.0;    
    m.setCol(0, col);
    col[0] = 4.0;    
    col[1] = 0.0;    
    col[2] = 6.0;    
    m.setCol(1, col);
    col[0] = 1.0;    
    col[1] = 3.0;    
    col[2] = 5.0;    
    m.setCol(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.swapCol(0, 1);
    CPPUNIT_ASSERT(s == e);
    
}    

void RectMatrixTestCase::testOtherMemberFunctions(){
    //test transpose
    CPPUNIT_ASSERT((a.transpose()).transpose() == a);
    
    CPPUNIT_ASSERT(a.transpose() == b);
}
Revision:	102
Committed:	Tue Oct 19 04:21:07 2004 UTC (21 years ago) by tim
Original Path:	*trunk/test/math/RectMatrixTestCase.cpp*
File size:	5889 byte(s)
Log Message:	implement more test cases
#	Content
1	#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	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	}
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	CPPUNIT_ASSERT(m2 != m3);
175
176	}
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
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	}
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	}
241
242	void RectMatrixTestCase::testAccessEntries(){
243	CPPUNIT_ASSERT_DOUBLES_EQUAL(m1(1, 0), 3.0, oopse::epsilon);
244	}
245
246	void RectMatrixTestCase::testRowColOperations() {
247	Vec3 row;
248	Vec3 col;
249	RMat3x3 m;
250
251	//test getRow
252	row = e.getRow(0);
253	CPPUNIT_ASSERT(row[0] == 2.0 && row[1] == 4.0 && row[2] == 1.0)
254
255	//test setRow
256	row[0] = 2.0;
257	row[1] = 4.0;
258	row[2] = 1.0;
259	m.setRow(0, row);
260	row[0] = 0.0;
261	row[1] = 0.0;
262	row[2] = 3.0;
263	m.setRow(1, row);
264	row[0] = 0.0;
265	row[1] = 6.0;
266	row[2] = 5.0;
267	m.setRow(2, row);
268	CPPUNIT_ASSERT(m == e);
269
270	//test getCol
271	col = e.getCol(1);
272	CPPUNIT_ASSERT(col[0] == 4.0 && col[1] == 0.0 && col[2] == 6.0)
273
274	//test setCol
275	col[0] = 2.0;
276	col[1] = 0.0;
277	col[2] = 0.0;
278	m.setCol(0, col);
279	col[0] = 4.0;
280	col[1] = 0.0;
281	col[2] = 6.0;
282	m.setCol(1, col);
283	col[0] = 1.0;
284	col[1] = 3.0;
285	col[2] = 5.0;
286	m.setCol(2, col);
287	CPPUNIT_ASSERT(m == e);
288
289	//test swapRow
290	RMat2x3 r;
291	r(0, 0) = 0.0;
292	r(0, 1) = 1.0;
293	r(0, 2) = 0.0;
294	r(1, 0) = 1.0;
295	r(1, 1) = 0.0;
296	r(1, 2) = 0.0;
297	r.swapRow(0, 1);
298	CPPUNIT_ASSERT(r == a);
299
300	//test swapCol
301	RMat3x3 s;
302	s(0, 0) = 4.0;
303	s(0, 1) = 2.0;
304	s(0, 2) = 1.0;
305	s(1, 0) = 0.0;
306	s(1, 1) = 0.0;
307	s(1, 2) = 3.0;
308	s(2, 0) = 6.0;
309	s(2, 1) = 0.0;
310	s(2, 2) = 5.0;
311
312	s.swapCol(0, 1);
313	CPPUNIT_ASSERT(s == e);
314
315	}
316
317	void RectMatrixTestCase::testOtherMemberFunctions(){
318	//test transpose
319	CPPUNIT_ASSERT((a.transpose()).transpose() == a);
320
321	CPPUNIT_ASSERT(a.transpose() == b);
322	}