ViewVC Help
View File | Revision Log | Show Annotations | View Changeset | Root Listing
root/OpenMD/trunk/test/math/RectMatrixTestCase.cpp
(Generate patch)

Comparing trunk/test/math/RectMatrixTestCase.cpp (file contents):
Revision 72 by tim, Wed Oct 13 22:27:58 2004 UTC vs.
Revision 137 by tim, Thu Oct 21 21:31:39 2004 UTC

# Line 55 | Line 55 | void RectMatrixTestCase::setUp(){
55      b(2, 0) = 0.0;
56      b(2, 1) = 0.0;
57  
58 <    
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(){
# Line 82 | Line 171 | void RectMatrixTestCase::testEqual() {
171  
172   void RectMatrixTestCase::testEqual() {
173      CPPUNIT_ASSERT(m2 == m3);
174 +    CPPUNIT_ASSERT(m2 != m3);
175 +    
176   }
177  
178   void RectMatrixTestCase::testNegate() {
# Line 129 | Line 220 | void RectMatrixTestCase::testMul() {
220      CPPUNIT_ASSERT(m1 * 1.0 == m1);
221      CPPUNIT_ASSERT(m1 * 0.0 == zero);
222      CPPUNIT_ASSERT(2.0 *m1 == m1 + m1);
223 <    
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() {
# Line 138 | Line 237 | void RectMatrixTestCase::testDiv() {
237  
238      CPPUNIT_ASSERT(m1 / 2.0 * 2.0 == m1);
239  
240 +
241   }
242  
243   void RectMatrixTestCase::testAccessEntries(){
244 <    CPPUNIT_ASSERT(m1(1, 0) == 3.0);
244 >    CPPUNIT_ASSERT_DOUBLES_EQUAL(m1(1, 0), 3.0, oopse::epsilon);
245   }
246  
247 < void RectMatrixTestCase::testTranspose(){
247 > void RectMatrixTestCase::testRowColOperations() {
248 >    Vec3 row;
249 >    Vec3 col;
250 >    RMat3x3 m;
251 >    
252 >    //test getRow
253 >    row = e.getRow(0);
254 >    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 >    //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 >    col = e.getColum(1);
274 >    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 >    //test setCol
278 >    col[0] = 2.0;    
279 >    col[1] = 0.0;    
280 >    col[2] = 0.0;    
281 >    m.setColum(0, col);
282 >    col[0] = 4.0;    
283 >    col[1] = 0.0;    
284 >    col[2] = 6.0;    
285 >    m.setColum(1, col);
286 >    col[0] = 1.0;    
287 >    col[1] = 3.0;    
288 >    col[2] = 5.0;    
289 >    m.setColum(2, col);
290 >    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 +    s.swapColum(0, 1);
316 +    CPPUNIT_ASSERT(s == e);
317 +
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 + }    
332 +
333 + void RectMatrixTestCase::testOtherMemberFunctions(){
334 +    //test transpose
335      CPPUNIT_ASSERT((a.transpose()).transpose() == a);
336      
337      CPPUNIT_ASSERT(a.transpose() == b);

Diff Legend

– Removed lines
+ Added lines
< Changed lines
> Changed lines