Changes in src/unittests/MatrixUnittest.cpp [407782:5630bd]
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
src/unittests/MatrixUnittest.cpp
r407782 r5630bd 39 39 } 40 40 perm1 = new Matrix(); 41 perm1->column(0) = unitVec[0];42 perm1->column(1) = unitVec[2];43 perm1->column(2) = unitVec[1];41 perm1->column(0) = e1; 42 perm1->column(1) = e3; 43 perm1->column(2) = e2; 44 44 45 45 46 46 perm2 = new Matrix(); 47 perm2->column(0) = unitVec[1];48 perm2->column(1) = unitVec[0];49 perm2->column(2) = unitVec[2];47 perm2->column(0) = e2; 48 perm2->column(1) = e1; 49 perm2->column(2) = e3; 50 50 51 51 perm3 = new Matrix(); 52 perm3->column(0) = unitVec[1];53 perm3->column(1) = unitVec[2];54 perm3->column(2) = unitVec[0];52 perm3->column(0) = e2; 53 perm3->column(1) = e3; 54 perm3->column(2) = e1; 55 55 56 56 perm4 = new Matrix(); 57 perm4->column(0) = unitVec[2];58 perm4->column(1) = unitVec[1];59 perm4->column(2) = unitVec[0];57 perm4->column(0) = e3; 58 perm4->column(1) = e2; 59 perm4->column(2) = e1; 60 60 61 61 perm5 = new Matrix(); 62 perm5->column(0) = unitVec[2];63 perm5->column(1) = unitVec[0];64 perm5->column(2) = unitVec[1];62 perm5->column(0) = e3; 63 perm5->column(1) = e1; 64 perm5->column(2) = e2; 65 65 66 66 } … … 108 108 109 109 mat.one(); 110 CPPUNIT_ASSERT_EQUAL(mat.row(0), unitVec[0]);111 CPPUNIT_ASSERT_EQUAL(mat.row(1), unitVec[1]);112 CPPUNIT_ASSERT_EQUAL(mat.row(2), unitVec[2]);113 CPPUNIT_ASSERT_EQUAL(mat.column(0), unitVec[0]);114 CPPUNIT_ASSERT_EQUAL(mat.column(1), unitVec[1]);115 CPPUNIT_ASSERT_EQUAL(mat.column(2), unitVec[2]);110 CPPUNIT_ASSERT_EQUAL(mat.row(0),e1); 111 CPPUNIT_ASSERT_EQUAL(mat.row(1),e2); 112 CPPUNIT_ASSERT_EQUAL(mat.row(2),e3); 113 CPPUNIT_ASSERT_EQUAL(mat.column(0),e1); 114 CPPUNIT_ASSERT_EQUAL(mat.column(1),e2); 115 CPPUNIT_ASSERT_EQUAL(mat.column(2),e3); 116 116 117 117 Vector t1=Vector(1.,1.,1.); … … 204 204 205 205 res = (*diagonal)*(*perm1); 206 CPPUNIT_ASSERT_EQUAL(res.column(0), unitVec[0]);207 CPPUNIT_ASSERT_EQUAL(res.column(1),3* unitVec[2]);208 CPPUNIT_ASSERT_EQUAL(res.column(2),2* unitVec[1]);206 CPPUNIT_ASSERT_EQUAL(res.column(0),e1); 207 CPPUNIT_ASSERT_EQUAL(res.column(1),3*e3); 208 CPPUNIT_ASSERT_EQUAL(res.column(2),2*e2); 209 209 res = (*diagonal)*(*perm2); 210 CPPUNIT_ASSERT_EQUAL(res.column(0),2* unitVec[1]);211 CPPUNIT_ASSERT_EQUAL(res.column(1), unitVec[0]);212 CPPUNIT_ASSERT_EQUAL(res.column(2),3* unitVec[2]);210 CPPUNIT_ASSERT_EQUAL(res.column(0),2*e2); 211 CPPUNIT_ASSERT_EQUAL(res.column(1),e1); 212 CPPUNIT_ASSERT_EQUAL(res.column(2),3*e3); 213 213 res = (*diagonal)*(*perm3); 214 CPPUNIT_ASSERT_EQUAL(res.column(0),2* unitVec[1]);215 CPPUNIT_ASSERT_EQUAL(res.column(1),3* unitVec[2]);216 CPPUNIT_ASSERT_EQUAL(res.column(2), unitVec[0]);214 CPPUNIT_ASSERT_EQUAL(res.column(0),2*e2); 215 CPPUNIT_ASSERT_EQUAL(res.column(1),3*e3); 216 CPPUNIT_ASSERT_EQUAL(res.column(2),e1); 217 217 res = (*diagonal)*(*perm4); 218 CPPUNIT_ASSERT_EQUAL(res.column(0),3* unitVec[2]);219 CPPUNIT_ASSERT_EQUAL(res.column(1),2* unitVec[1]);220 CPPUNIT_ASSERT_EQUAL(res.column(2), unitVec[0]);218 CPPUNIT_ASSERT_EQUAL(res.column(0),3*e3); 219 CPPUNIT_ASSERT_EQUAL(res.column(1),2*e2); 220 CPPUNIT_ASSERT_EQUAL(res.column(2),e1); 221 221 res = (*diagonal)*(*perm5); 222 CPPUNIT_ASSERT_EQUAL(res.column(0),3* unitVec[2]);223 CPPUNIT_ASSERT_EQUAL(res.column(1), unitVec[0]);224 CPPUNIT_ASSERT_EQUAL(res.column(2),2* unitVec[1]);222 CPPUNIT_ASSERT_EQUAL(res.column(0),3*e3); 223 CPPUNIT_ASSERT_EQUAL(res.column(1),e1); 224 CPPUNIT_ASSERT_EQUAL(res.column(2),2*e2); 225 225 } 226 226 … … 261 261 262 262 void MatrixUnittest::VecMultTest(){ 263 CPPUNIT_ASSERT_EQUAL((*zero)* unitVec[0],zeroVec);264 CPPUNIT_ASSERT_EQUAL((*zero)* unitVec[1],zeroVec);265 CPPUNIT_ASSERT_EQUAL((*zero)* unitVec[2],zeroVec);263 CPPUNIT_ASSERT_EQUAL((*zero)*e1,zeroVec); 264 CPPUNIT_ASSERT_EQUAL((*zero)*e2,zeroVec); 265 CPPUNIT_ASSERT_EQUAL((*zero)*e3,zeroVec); 266 266 CPPUNIT_ASSERT_EQUAL((*zero)*zeroVec,zeroVec); 267 267 268 CPPUNIT_ASSERT_EQUAL((*one)* unitVec[0],unitVec[0]);269 CPPUNIT_ASSERT_EQUAL((*one)* unitVec[1],unitVec[1]);270 CPPUNIT_ASSERT_EQUAL((*one)* unitVec[2],unitVec[2]);268 CPPUNIT_ASSERT_EQUAL((*one)*e1,e1); 269 CPPUNIT_ASSERT_EQUAL((*one)*e2,e2); 270 CPPUNIT_ASSERT_EQUAL((*one)*e3,e3); 271 271 CPPUNIT_ASSERT_EQUAL((*one)*zeroVec,zeroVec); 272 272 273 CPPUNIT_ASSERT_EQUAL((*diagonal)* unitVec[0],unitVec[0]);274 CPPUNIT_ASSERT_EQUAL((*diagonal)* unitVec[1],2*unitVec[1]);275 CPPUNIT_ASSERT_EQUAL((*diagonal)* unitVec[2],3*unitVec[2]);273 CPPUNIT_ASSERT_EQUAL((*diagonal)*e1,e1); 274 CPPUNIT_ASSERT_EQUAL((*diagonal)*e2,2*e2); 275 CPPUNIT_ASSERT_EQUAL((*diagonal)*e3,3*e3); 276 276 CPPUNIT_ASSERT_EQUAL((*diagonal)*zeroVec,zeroVec); 277 277 278 CPPUNIT_ASSERT_EQUAL((*perm1)* unitVec[0],unitVec[0]);279 CPPUNIT_ASSERT_EQUAL((*perm1)* unitVec[1],unitVec[2]);280 CPPUNIT_ASSERT_EQUAL((*perm1)* unitVec[2],unitVec[1]);278 CPPUNIT_ASSERT_EQUAL((*perm1)*e1,e1); 279 CPPUNIT_ASSERT_EQUAL((*perm1)*e2,e3); 280 CPPUNIT_ASSERT_EQUAL((*perm1)*e3,e2); 281 281 CPPUNIT_ASSERT_EQUAL((*perm1)*zeroVec,zeroVec); 282 282 283 CPPUNIT_ASSERT_EQUAL((*perm2)* unitVec[0],unitVec[1]);284 CPPUNIT_ASSERT_EQUAL((*perm2)* unitVec[1],unitVec[0]);285 CPPUNIT_ASSERT_EQUAL((*perm2)* unitVec[2],unitVec[2]);283 CPPUNIT_ASSERT_EQUAL((*perm2)*e1,e2); 284 CPPUNIT_ASSERT_EQUAL((*perm2)*e2,e1); 285 CPPUNIT_ASSERT_EQUAL((*perm2)*e3,e3); 286 286 CPPUNIT_ASSERT_EQUAL((*perm2)*zeroVec,zeroVec); 287 287 288 CPPUNIT_ASSERT_EQUAL((*perm3)* unitVec[0],unitVec[1]);289 CPPUNIT_ASSERT_EQUAL((*perm3)* unitVec[1],unitVec[2]);290 CPPUNIT_ASSERT_EQUAL((*perm3)* unitVec[2],unitVec[0]);288 CPPUNIT_ASSERT_EQUAL((*perm3)*e1,e2); 289 CPPUNIT_ASSERT_EQUAL((*perm3)*e2,e3); 290 CPPUNIT_ASSERT_EQUAL((*perm3)*e3,e1); 291 291 CPPUNIT_ASSERT_EQUAL((*perm3)*zeroVec,zeroVec); 292 292 293 CPPUNIT_ASSERT_EQUAL((*perm4)* unitVec[0],unitVec[2]);294 CPPUNIT_ASSERT_EQUAL((*perm4)* unitVec[1],unitVec[1]);295 CPPUNIT_ASSERT_EQUAL((*perm4)* unitVec[2],unitVec[0]);293 CPPUNIT_ASSERT_EQUAL((*perm4)*e1,e3); 294 CPPUNIT_ASSERT_EQUAL((*perm4)*e2,e2); 295 CPPUNIT_ASSERT_EQUAL((*perm4)*e3,e1); 296 296 CPPUNIT_ASSERT_EQUAL((*perm4)*zeroVec,zeroVec); 297 297 298 CPPUNIT_ASSERT_EQUAL((*perm5)* unitVec[0],unitVec[2]);299 CPPUNIT_ASSERT_EQUAL((*perm5)* unitVec[1],unitVec[0]);300 CPPUNIT_ASSERT_EQUAL((*perm5)* unitVec[2],unitVec[1]);298 CPPUNIT_ASSERT_EQUAL((*perm5)*e1,e3); 299 CPPUNIT_ASSERT_EQUAL((*perm5)*e2,e1); 300 CPPUNIT_ASSERT_EQUAL((*perm5)*e3,e2); 301 301 CPPUNIT_ASSERT_EQUAL((*perm5)*zeroVec,zeroVec); 302 302
Note:
See TracChangeset
for help on using the changeset viewer.