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::testTranspose(){

    CPPUNIT_ASSERT((a.transpose()).transpose() == a);
    
    CPPUNIT_ASSERT(a.transpose() == b);
}
Revision:	78
Committed:	Fri Oct 15 15:19:28 2004 UTC (20 years, 6 months ago) by tim
File size:	4437 byte(s)
Log Message:	add SquareMatrixTestCase class
#	User	Rev	Content
1	tim	72	#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	78	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	72	}
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	78	CPPUNIT_ASSERT(m2 != m3);
175
176	tim	72	}
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	78
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	72	}
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	tim	78	CPPUNIT_ASSERT_DOUBLES_EQUAL(m1(1, 0), 3.0, oopse::epsilon);
244	tim	72	}
245
246			void RectMatrixTestCase::testTranspose(){
247
248			CPPUNIT_ASSERT((a.transpose()).transpose() == a);
249
250			CPPUNIT_ASSERT(a.transpose() == b);
251			}