Ignore:
File:
1 edited

Legend:

Unmodified
Added
Removed
  • src/unittests/MatrixUnittest.cpp

    r5630bd r407782  
    3939  }
    4040  perm1 = new Matrix();
    41   perm1->column(0) = e1;
    42   perm1->column(1) = e3;
    43   perm1->column(2) = e2;
     41  perm1->column(0) = unitVec[0];
     42  perm1->column(1) = unitVec[2];
     43  perm1->column(2) = unitVec[1];
    4444
    4545
    4646  perm2 = new Matrix();
    47   perm2->column(0) = e2;
    48   perm2->column(1) = e1;
    49   perm2->column(2) = e3;
     47  perm2->column(0) = unitVec[1];
     48  perm2->column(1) = unitVec[0];
     49  perm2->column(2) = unitVec[2];
    5050
    5151  perm3 = new Matrix();
    52   perm3->column(0) = e2;
    53   perm3->column(1) = e3;
    54   perm3->column(2) = e1;
     52  perm3->column(0) = unitVec[1];
     53  perm3->column(1) = unitVec[2];
     54  perm3->column(2) = unitVec[0];
    5555
    5656  perm4 = new Matrix();
    57   perm4->column(0) = e3;
    58   perm4->column(1) = e2;
    59   perm4->column(2) = e1;
     57  perm4->column(0) = unitVec[2];
     58  perm4->column(1) = unitVec[1];
     59  perm4->column(2) = unitVec[0];
    6060
    6161  perm5 = new Matrix();
    62   perm5->column(0) = e3;
    63   perm5->column(1) = e1;
    64   perm5->column(2) = e2;
     62  perm5->column(0) = unitVec[2];
     63  perm5->column(1) = unitVec[0];
     64  perm5->column(2) = unitVec[1];
    6565
    6666}
     
    108108
    109109  mat.one();
    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);
     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]);
    116116
    117117  Vector t1=Vector(1.,1.,1.);
     
    204204
    205205  res = (*diagonal)*(*perm1);
    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);
     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]);
    209209  res = (*diagonal)*(*perm2);
    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);
     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]);
    213213  res = (*diagonal)*(*perm3);
    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);
     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]);
    217217  res = (*diagonal)*(*perm4);
    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);
     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]);
    221221  res = (*diagonal)*(*perm5);
    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);
     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]);
    225225}
    226226
     
    261261
    262262void MatrixUnittest::VecMultTest(){
    263   CPPUNIT_ASSERT_EQUAL((*zero)*e1,zeroVec);
    264   CPPUNIT_ASSERT_EQUAL((*zero)*e2,zeroVec);
    265   CPPUNIT_ASSERT_EQUAL((*zero)*e3,zeroVec);
     263  CPPUNIT_ASSERT_EQUAL((*zero)*unitVec[0],zeroVec);
     264  CPPUNIT_ASSERT_EQUAL((*zero)*unitVec[1],zeroVec);
     265  CPPUNIT_ASSERT_EQUAL((*zero)*unitVec[2],zeroVec);
    266266  CPPUNIT_ASSERT_EQUAL((*zero)*zeroVec,zeroVec);
    267267
    268   CPPUNIT_ASSERT_EQUAL((*one)*e1,e1);
    269   CPPUNIT_ASSERT_EQUAL((*one)*e2,e2);
    270   CPPUNIT_ASSERT_EQUAL((*one)*e3,e3);
     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]);
    271271  CPPUNIT_ASSERT_EQUAL((*one)*zeroVec,zeroVec);
    272272
    273   CPPUNIT_ASSERT_EQUAL((*diagonal)*e1,e1);
    274   CPPUNIT_ASSERT_EQUAL((*diagonal)*e2,2*e2);
    275   CPPUNIT_ASSERT_EQUAL((*diagonal)*e3,3*e3);
     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]);
    276276  CPPUNIT_ASSERT_EQUAL((*diagonal)*zeroVec,zeroVec);
    277277
    278   CPPUNIT_ASSERT_EQUAL((*perm1)*e1,e1);
    279   CPPUNIT_ASSERT_EQUAL((*perm1)*e2,e3);
    280   CPPUNIT_ASSERT_EQUAL((*perm1)*e3,e2);
     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]);
    281281  CPPUNIT_ASSERT_EQUAL((*perm1)*zeroVec,zeroVec);
    282282
    283   CPPUNIT_ASSERT_EQUAL((*perm2)*e1,e2);
    284   CPPUNIT_ASSERT_EQUAL((*perm2)*e2,e1);
    285   CPPUNIT_ASSERT_EQUAL((*perm2)*e3,e3);
     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]);
    286286  CPPUNIT_ASSERT_EQUAL((*perm2)*zeroVec,zeroVec);
    287287
    288   CPPUNIT_ASSERT_EQUAL((*perm3)*e1,e2);
    289   CPPUNIT_ASSERT_EQUAL((*perm3)*e2,e3);
    290   CPPUNIT_ASSERT_EQUAL((*perm3)*e3,e1);
     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]);
    291291  CPPUNIT_ASSERT_EQUAL((*perm3)*zeroVec,zeroVec);
    292292
    293   CPPUNIT_ASSERT_EQUAL((*perm4)*e1,e3);
    294   CPPUNIT_ASSERT_EQUAL((*perm4)*e2,e2);
    295   CPPUNIT_ASSERT_EQUAL((*perm4)*e3,e1);
     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]);
    296296  CPPUNIT_ASSERT_EQUAL((*perm4)*zeroVec,zeroVec);
    297297
    298   CPPUNIT_ASSERT_EQUAL((*perm5)*e1,e3);
    299   CPPUNIT_ASSERT_EQUAL((*perm5)*e2,e1);
    300   CPPUNIT_ASSERT_EQUAL((*perm5)*e3,e2);
     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]);
    301301  CPPUNIT_ASSERT_EQUAL((*perm5)*zeroVec,zeroVec);
    302302
Note: See TracChangeset for help on using the changeset viewer.