Changes in src/unittests/LineUnittest.cpp [407782:3dcb1f]

File:

• src/unittests/LineUnittest.cpp (modified) (9 diffs)

src/unittests/LineUnittest.cpp

@@ -29 +29 @@
 void LineUnittest::setUp(){
   // three lines along the axes
-  la1 = new Line(zeroVec,unitVec[0]);
-  la2 = new Line(zeroVec,unitVec[1]);
-  la3 = new Line(zeroVec,unitVec[2]);
+  la1 = new Line(zeroVec,e1);
+  la2 = new Line(zeroVec,e2);
+  la3 = new Line(zeroVec,e3);
 
   // the lines along the planes defined by two coordinate axes
-  lp1 = new Line(unitVec[0],unitVec[0]-unitVec[1]);
-  lp2 = new Line(unitVec[1],unitVec[1]-unitVec[2]);
-  lp3 = new Line(unitVec[2],unitVec[2]-unitVec[0]);
+  lp1 = new Line(e1,e1-e2);
+  lp2 = new Line(e2,e2-e3);
+  lp3 = new Line(e3,e3-e1);
 }
 void LineUnittest::tearDown(){
@@ -52 +52 @@
 
   // direction+origin should never fail
-  CPPUNIT_ASSERT_NO_THROW(Line(zeroVec,unitVec[0]));
-  CPPUNIT_ASSERT_NO_THROW(Line(zeroVec,unitVec[1]));
-  CPPUNIT_ASSERT_NO_THROW(Line(zeroVec,unitVec[2]));
+  CPPUNIT_ASSERT_NO_THROW(Line(zeroVec,e1));
+  CPPUNIT_ASSERT_NO_THROW(Line(zeroVec,e2));
+  CPPUNIT_ASSERT_NO_THROW(Line(zeroVec,e3));
 
   // two points fails if both points are the same
-  CPPUNIT_ASSERT_NO_THROW(makeLineThrough(unitVec[0],unitVec[1]));
-  CPPUNIT_ASSERT_NO_THROW(makeLineThrough(unitVec[1],unitVec[2]));
-  CPPUNIT_ASSERT_NO_THROW(makeLineThrough(unitVec[2],unitVec[0]));
+  CPPUNIT_ASSERT_NO_THROW(makeLineThrough(e1,e2));
+  CPPUNIT_ASSERT_NO_THROW(makeLineThrough(e2,e3));
+  CPPUNIT_ASSERT_NO_THROW(makeLineThrough(e3,e1));
   // for zerovectors
-  CPPUNIT_ASSERT_NO_THROW(makeLineThrough(unitVec[0],zeroVec));
-  CPPUNIT_ASSERT_NO_THROW(makeLineThrough(unitVec[1],zeroVec));
-  CPPUNIT_ASSERT_NO_THROW(makeLineThrough(unitVec[2],zeroVec));
+  CPPUNIT_ASSERT_NO_THROW(makeLineThrough(e1,zeroVec));
+  CPPUNIT_ASSERT_NO_THROW(makeLineThrough(e2,zeroVec));
+  CPPUNIT_ASSERT_NO_THROW(makeLineThrough(e3,zeroVec));
   // now we pass two times the same point
   CPPUNIT_ASSERT_THROW(makeLineThrough(zeroVec,zeroVec),LinearDependenceException);
-  CPPUNIT_ASSERT_THROW(makeLineThrough(unitVec[0],unitVec[0]),LinearDependenceException);
-  CPPUNIT_ASSERT_THROW(makeLineThrough(unitVec[1],unitVec[1]),LinearDependenceException);
-  CPPUNIT_ASSERT_THROW(makeLineThrough(unitVec[2],unitVec[2]),LinearDependenceException);
+  CPPUNIT_ASSERT_THROW(makeLineThrough(e1,e1),LinearDependenceException);
+  CPPUNIT_ASSERT_THROW(makeLineThrough(e2,e2),LinearDependenceException);
+  CPPUNIT_ASSERT_THROW(makeLineThrough(e3,e3),LinearDependenceException);
 
 }
@@ -78 +78 @@
 void LineUnittest::constructionResultTest(){
   // test all directions
-  CPPUNIT_ASSERT(testDirection(la1->getDirection(),unitVec[0]));
-  CPPUNIT_ASSERT(testDirection(la2->getDirection(),unitVec[1]));
-  CPPUNIT_ASSERT(testDirection(la3->getDirection(),unitVec[2]));
+  CPPUNIT_ASSERT(testDirection(la1->getDirection(),e1));
+  CPPUNIT_ASSERT(testDirection(la2->getDirection(),e2));
+  CPPUNIT_ASSERT(testDirection(la3->getDirection(),e3));
 
   // test origins
@@ -92 +92 @@
   CPPUNIT_ASSERT(la3->isContained(zeroVec));
 
-  CPPUNIT_ASSERT(la1->isContained(unitVec[0]));
-  CPPUNIT_ASSERT(la2->isContained(unitVec[1]));
-  CPPUNIT_ASSERT(la3->isContained(unitVec[2]));
-
-  CPPUNIT_ASSERT(lp1->isContained(unitVec[0]));
-  CPPUNIT_ASSERT(lp2->isContained(unitVec[1]));
-  CPPUNIT_ASSERT(lp3->isContained(unitVec[2]));
-
-  CPPUNIT_ASSERT(lp1->isContained(unitVec[1]));
-  CPPUNIT_ASSERT(lp2->isContained(unitVec[2]));
-  CPPUNIT_ASSERT(lp3->isContained(unitVec[0]));
+  CPPUNIT_ASSERT(la1->isContained(e1));
+  CPPUNIT_ASSERT(la2->isContained(e2));
+  CPPUNIT_ASSERT(la3->isContained(e3));
+
+  CPPUNIT_ASSERT(lp1->isContained(e1));
+  CPPUNIT_ASSERT(lp2->isContained(e2));
+  CPPUNIT_ASSERT(lp3->isContained(e3));
+
+  CPPUNIT_ASSERT(lp1->isContained(e2));
+  CPPUNIT_ASSERT(lp2->isContained(e3));
+  CPPUNIT_ASSERT(lp3->isContained(e1));
 }
 
@@ -112 +112 @@
 
   // multiples of the second support vector
-  CPPUNIT_ASSERT(la1->isContained(unitVec[0]));
-  CPPUNIT_ASSERT(la2->isContained(unitVec[1]));
-  CPPUNIT_ASSERT(la3->isContained(unitVec[2]));
-
-  CPPUNIT_ASSERT(la1->isContained(2*unitVec[0]));
-  CPPUNIT_ASSERT(la2->isContained(2*unitVec[1]));
-  CPPUNIT_ASSERT(la3->isContained(2*unitVec[2]));
-
-  CPPUNIT_ASSERT(la1->isContained(3*unitVec[0]));
-  CPPUNIT_ASSERT(la2->isContained(3*unitVec[1]));
-  CPPUNIT_ASSERT(la3->isContained(3*unitVec[2]));
+  CPPUNIT_ASSERT(la1->isContained(e1));
+  CPPUNIT_ASSERT(la2->isContained(e2));
+  CPPUNIT_ASSERT(la3->isContained(e3));
+
+  CPPUNIT_ASSERT(la1->isContained(2*e1));
+  CPPUNIT_ASSERT(la2->isContained(2*e2));
+  CPPUNIT_ASSERT(la3->isContained(2*e3));
+
+  CPPUNIT_ASSERT(la1->isContained(3*e1));
+  CPPUNIT_ASSERT(la2->isContained(3*e2));
+  CPPUNIT_ASSERT(la3->isContained(3*e3));
 
   // negative multiples
-  CPPUNIT_ASSERT(la1->isContained(-1*unitVec[0]));
-  CPPUNIT_ASSERT(la2->isContained(-1*unitVec[1]));
-  CPPUNIT_ASSERT(la3->isContained(-1*unitVec[2]));
-
-  CPPUNIT_ASSERT(la1->isContained(-2*unitVec[0]));
-  CPPUNIT_ASSERT(la2->isContained(-2*unitVec[1]));
-  CPPUNIT_ASSERT(la3->isContained(-2*unitVec[2]));
+  CPPUNIT_ASSERT(la1->isContained(-1*e1));
+  CPPUNIT_ASSERT(la2->isContained(-1*e2));
+  CPPUNIT_ASSERT(la3->isContained(-1*e3));
+
+  CPPUNIT_ASSERT(la1->isContained(-2*e1));
+  CPPUNIT_ASSERT(la2->isContained(-2*e2));
+  CPPUNIT_ASSERT(la3->isContained(-2*e3));
 
   // points that should not be on the lines
-  CPPUNIT_ASSERT(!la1->isContained(unitVec[1]));
-  CPPUNIT_ASSERT(!la2->isContained(unitVec[2]));
-  CPPUNIT_ASSERT(!la3->isContained(unitVec[0]));
-
-  CPPUNIT_ASSERT(!la1->isContained(2*unitVec[1]));
-  CPPUNIT_ASSERT(!la2->isContained(2*unitVec[2]));
-  CPPUNIT_ASSERT(!la3->isContained(2*unitVec[0]));
-
-  CPPUNIT_ASSERT(!la1->isContained(-1*unitVec[1]));
-  CPPUNIT_ASSERT(!la2->isContained(-1*unitVec[2]));
-  CPPUNIT_ASSERT(!la3->isContained(-1*unitVec[0]));
+  CPPUNIT_ASSERT(!la1->isContained(e2));
+  CPPUNIT_ASSERT(!la2->isContained(e3));
+  CPPUNIT_ASSERT(!la3->isContained(e1));
+
+  CPPUNIT_ASSERT(!la1->isContained(2*e2));
+  CPPUNIT_ASSERT(!la2->isContained(2*e3));
+  CPPUNIT_ASSERT(!la3->isContained(2*e1));
+
+  CPPUNIT_ASSERT(!la1->isContained(-1*e2));
+  CPPUNIT_ASSERT(!la2->isContained(-1*e3));
+  CPPUNIT_ASSERT(!la3->isContained(-1*e1));
 
   // For the plane lines
-  CPPUNIT_ASSERT(lp1->isContained(unitVec[0]));
-  CPPUNIT_ASSERT(lp2->isContained(unitVec[1]));
-  CPPUNIT_ASSERT(lp3->isContained(unitVec[2]));
-
-  CPPUNIT_ASSERT(lp1->isContained(unitVec[1]));
-  CPPUNIT_ASSERT(lp2->isContained(unitVec[2]));
-  CPPUNIT_ASSERT(lp3->isContained(unitVec[0]));
-
-  CPPUNIT_ASSERT(lp1->isContained(unitVec[0]+2*(unitVec[0]-unitVec[1])));
-  CPPUNIT_ASSERT(lp2->isContained(unitVec[1]+2*(unitVec[1]-unitVec[2])));
-  CPPUNIT_ASSERT(lp3->isContained(unitVec[2]+2*(unitVec[2]-unitVec[0])));
-
-  CPPUNIT_ASSERT(lp1->isContained(unitVec[0]-2*(unitVec[0]-unitVec[1])));
-  CPPUNIT_ASSERT(lp2->isContained(unitVec[1]-2*(unitVec[1]-unitVec[2])));
-  CPPUNIT_ASSERT(lp3->isContained(unitVec[2]-2*(unitVec[2]-unitVec[0])));
+  CPPUNIT_ASSERT(lp1->isContained(e1));
+  CPPUNIT_ASSERT(lp2->isContained(e2));
+  CPPUNIT_ASSERT(lp3->isContained(e3));
+
+  CPPUNIT_ASSERT(lp1->isContained(e2));
+  CPPUNIT_ASSERT(lp2->isContained(e3));
+  CPPUNIT_ASSERT(lp3->isContained(e1));
+
+  CPPUNIT_ASSERT(lp1->isContained(e1+2*(e1-e2)));
+  CPPUNIT_ASSERT(lp2->isContained(e2+2*(e2-e3)));
+  CPPUNIT_ASSERT(lp3->isContained(e3+2*(e3-e1)));
+
+  CPPUNIT_ASSERT(lp1->isContained(e1-2*(e1-e2)));
+  CPPUNIT_ASSERT(lp2->isContained(e2-2*(e2-e3)));
+  CPPUNIT_ASSERT(lp3->isContained(e3-2*(e3-e1)));
 }
 
@@ -177 +177 @@
   // axes and plane lines
   fixture = la1->getIntersection(*lp1);
-  CPPUNIT_ASSERT_EQUAL(fixture,unitVec[0]);
+  CPPUNIT_ASSERT_EQUAL(fixture,e1);
   fixture = la2->getIntersection(*lp2);
-  CPPUNIT_ASSERT_EQUAL(fixture,unitVec[1]);
+  CPPUNIT_ASSERT_EQUAL(fixture,e2);
   fixture = la3->getIntersection(*lp3);
-  CPPUNIT_ASSERT_EQUAL(fixture,unitVec[2]);
+  CPPUNIT_ASSERT_EQUAL(fixture,e3);
 
   fixture = la1->getIntersection(*lp3);
-  CPPUNIT_ASSERT_EQUAL(fixture,unitVec[0]);
+  CPPUNIT_ASSERT_EQUAL(fixture,e1);
   fixture = la2->getIntersection(*lp1);
-  CPPUNIT_ASSERT_EQUAL(fixture,unitVec[1]);
+  CPPUNIT_ASSERT_EQUAL(fixture,e2);
   fixture = la3->getIntersection(*lp2);
-  CPPUNIT_ASSERT_EQUAL(fixture,unitVec[2]);
+  CPPUNIT_ASSERT_EQUAL(fixture,e3);
 
   // two plane lines
   fixture = lp1->getIntersection(*lp2);
-  CPPUNIT_ASSERT_EQUAL(fixture,unitVec[1]);
+  CPPUNIT_ASSERT_EQUAL(fixture,e2);
   fixture = lp2->getIntersection(*lp3);
-  CPPUNIT_ASSERT_EQUAL(fixture,unitVec[2]);
+  CPPUNIT_ASSERT_EQUAL(fixture,e3);
   fixture = lp3->getIntersection(*lp1);
-  CPPUNIT_ASSERT_EQUAL(fixture,unitVec[0]);
+  CPPUNIT_ASSERT_EQUAL(fixture,e1);
 
   // When we have two times the same line, we check if the point is on the line
@@ -242 +242 @@
 
   // rotate vectors on the axis around their lines
-  fixture = la1->rotateVector(unitVec[0],1.);
-  CPPUNIT_ASSERT_EQUAL(fixture,unitVec[0]);
-  fixture = la2->rotateVector(unitVec[1],1.);
-  CPPUNIT_ASSERT_EQUAL(fixture,unitVec[1]);
-  fixture = la3->rotateVector(unitVec[2],1.);
-  CPPUNIT_ASSERT_EQUAL(fixture,unitVec[2]);
-
-  fixture = la1->rotateVector(unitVec[0],2.);
-  CPPUNIT_ASSERT_EQUAL(fixture,unitVec[0]);
-  fixture = la2->rotateVector(unitVec[1],2.);
-  CPPUNIT_ASSERT_EQUAL(fixture,unitVec[1]);
-  fixture = la3->rotateVector(unitVec[2],2.);
-  CPPUNIT_ASSERT_EQUAL(fixture,unitVec[2]);
+  fixture = la1->rotateVector(e1,1.);
+  CPPUNIT_ASSERT_EQUAL(fixture,e1);
+  fixture = la2->rotateVector(e2,1.);
+  CPPUNIT_ASSERT_EQUAL(fixture,e2);
+  fixture = la3->rotateVector(e3,1.);
+  CPPUNIT_ASSERT_EQUAL(fixture,e3);
+
+  fixture = la1->rotateVector(e1,2.);
+  CPPUNIT_ASSERT_EQUAL(fixture,e1);
+  fixture = la2->rotateVector(e2,2.);
+  CPPUNIT_ASSERT_EQUAL(fixture,e2);
+  fixture = la3->rotateVector(e3,2.);
+  CPPUNIT_ASSERT_EQUAL(fixture,e3);
 
   // more vectors on the axis
-  fixture = la1->rotateVector(2*unitVec[0],1.);
-  CPPUNIT_ASSERT_EQUAL(fixture,2*unitVec[0]);
-  fixture = la2->rotateVector(2*unitVec[1],1.);
-  CPPUNIT_ASSERT_EQUAL(fixture,2*unitVec[1]);
-  fixture = la3->rotateVector(2*unitVec[2],1.);
-  CPPUNIT_ASSERT_EQUAL(fixture,2*unitVec[2]);
-
-  fixture = la1->rotateVector(2*unitVec[0],2.);
-  CPPUNIT_ASSERT_EQUAL(fixture,2*unitVec[0]);
-  fixture = la2->rotateVector(2*unitVec[1],2.);
-  CPPUNIT_ASSERT_EQUAL(fixture,2*unitVec[1]);
-  fixture = la3->rotateVector(2*unitVec[2],2.);
-  CPPUNIT_ASSERT_EQUAL(fixture,2*unitVec[2]);
+  fixture = la1->rotateVector(2*e1,1.);
+  CPPUNIT_ASSERT_EQUAL(fixture,2*e1);
+  fixture = la2->rotateVector(2*e2,1.);
+  CPPUNIT_ASSERT_EQUAL(fixture,2*e2);
+  fixture = la3->rotateVector(2*e3,1.);
+  CPPUNIT_ASSERT_EQUAL(fixture,2*e3);
+
+  fixture = la1->rotateVector(2*e1,2.);
+  CPPUNIT_ASSERT_EQUAL(fixture,2*e1);
+  fixture = la2->rotateVector(2*e2,2.);
+  CPPUNIT_ASSERT_EQUAL(fixture,2*e2);
+  fixture = la3->rotateVector(2*e3,2.);
+  CPPUNIT_ASSERT_EQUAL(fixture,2*e3);
 
   // negative factors
-  fixture = la1->rotateVector(-1*unitVec[0],1.);
-  CPPUNIT_ASSERT_EQUAL(fixture,-1*unitVec[0]);
-  fixture = la2->rotateVector(-1*unitVec[1],1.);
-  CPPUNIT_ASSERT_EQUAL(fixture,-1*unitVec[1]);
-  fixture = la3->rotateVector(-1*unitVec[2],1.);
-  CPPUNIT_ASSERT_EQUAL(fixture,-1*unitVec[2]);
-
-  fixture = la1->rotateVector(-1*unitVec[0],2.);
-  CPPUNIT_ASSERT_EQUAL(fixture,-1*unitVec[0]);
-  fixture = la2->rotateVector(-1*unitVec[1],2.);
-  CPPUNIT_ASSERT_EQUAL(fixture,-1*unitVec[1]);
-  fixture = la3->rotateVector(-1*unitVec[2],2.);
-  CPPUNIT_ASSERT_EQUAL(fixture,-1*unitVec[2]);
+  fixture = la1->rotateVector(-1*e1,1.);
+  CPPUNIT_ASSERT_EQUAL(fixture,-1*e1);
+  fixture = la2->rotateVector(-1*e2,1.);
+  CPPUNIT_ASSERT_EQUAL(fixture,-1*e2);
+  fixture = la3->rotateVector(-1*e3,1.);
+  CPPUNIT_ASSERT_EQUAL(fixture,-1*e3);
+
+  fixture = la1->rotateVector(-1*e1,2.);
+  CPPUNIT_ASSERT_EQUAL(fixture,-1*e1);
+  fixture = la2->rotateVector(-1*e2,2.);
+  CPPUNIT_ASSERT_EQUAL(fixture,-1*e2);
+  fixture = la3->rotateVector(-1*e3,2.);
+  CPPUNIT_ASSERT_EQUAL(fixture,-1*e3);
 
 
 
   // now the real rotations
-  // unitVec[1] around unitVec[0]
-  fixture = la1->rotateVector(unitVec[1],0);
-  CPPUNIT_ASSERT_EQUAL(fixture,unitVec[1]);
-  fixture = la1->rotateVector(unitVec[1],1./2.*M_PI);
-  CPPUNIT_ASSERT_EQUAL(fixture,-1*unitVec[2]);
-  fixture = la1->rotateVector(unitVec[1],M_PI);
-  CPPUNIT_ASSERT_EQUAL(fixture,-1*unitVec[1]);
-  fixture = la1->rotateVector(unitVec[1],2*M_PI);
-  CPPUNIT_ASSERT_EQUAL(fixture,unitVec[1]);
-
-  // unitVec[2] around unitVec[1]
-  fixture = la2->rotateVector(unitVec[2],0);
-  CPPUNIT_ASSERT_EQUAL(fixture,unitVec[2]);
-  fixture = la2->rotateVector(unitVec[2],1./2.*M_PI);
-  CPPUNIT_ASSERT_EQUAL(fixture,-1*unitVec[0]);
-  fixture = la2->rotateVector(unitVec[2],M_PI);
-  CPPUNIT_ASSERT_EQUAL(fixture,-1*unitVec[2]);
-  fixture = la2->rotateVector(unitVec[2],2*M_PI);
-  CPPUNIT_ASSERT_EQUAL(fixture,unitVec[2]);
-
-  // unitVec[0] around unitVec[2]
-  fixture = la3->rotateVector(unitVec[0],0);
-  CPPUNIT_ASSERT_EQUAL(fixture,unitVec[0]);
-  fixture = la3->rotateVector(unitVec[0],1./2.*M_PI);
-  CPPUNIT_ASSERT_EQUAL(fixture,-1*unitVec[1]);
-  fixture = la3->rotateVector(unitVec[0],M_PI);
-  CPPUNIT_ASSERT_EQUAL(fixture,-1*unitVec[0]);
-  fixture = la3->rotateVector(unitVec[0],2*M_PI);
-  CPPUNIT_ASSERT_EQUAL(fixture,unitVec[0]);
+  // e2 around e1
+  fixture = la1->rotateVector(e2,0);
+  CPPUNIT_ASSERT_EQUAL(fixture,e2);
+  fixture = la1->rotateVector(e2,1./2.*M_PI);
+  CPPUNIT_ASSERT_EQUAL(fixture,-1*e3);
+  fixture = la1->rotateVector(e2,M_PI);
+  CPPUNIT_ASSERT_EQUAL(fixture,-1*e2);
+  fixture = la1->rotateVector(e2,2*M_PI);
+  CPPUNIT_ASSERT_EQUAL(fixture,e2);
+
+  // e3 around e2
+  fixture = la2->rotateVector(e3,0);
+  CPPUNIT_ASSERT_EQUAL(fixture,e3);
+  fixture = la2->rotateVector(e3,1./2.*M_PI);
+  CPPUNIT_ASSERT_EQUAL(fixture,-1*e1);
+  fixture = la2->rotateVector(e3,M_PI);
+  CPPUNIT_ASSERT_EQUAL(fixture,-1*e3);
+  fixture = la2->rotateVector(e3,2*M_PI);
+  CPPUNIT_ASSERT_EQUAL(fixture,e3);
+
+  // e1 around e3
+  fixture = la3->rotateVector(e1,0);
+  CPPUNIT_ASSERT_EQUAL(fixture,e1);
+  fixture = la3->rotateVector(e1,1./2.*M_PI);
+  CPPUNIT_ASSERT_EQUAL(fixture,-1*e2);
+  fixture = la3->rotateVector(e1,M_PI);
+  CPPUNIT_ASSERT_EQUAL(fixture,-1*e1);
+  fixture = la3->rotateVector(e1,2*M_PI);
+  CPPUNIT_ASSERT_EQUAL(fixture,e1);
 
 
@@ -323 +323 @@
 
   // Vectors on the line
-  fixture = lp1->rotateVector(unitVec[0],1.);
-  CPPUNIT_ASSERT_EQUAL(fixture,unitVec[0]);
-  fixture = lp1->rotateVector(unitVec[1],1.);
-  CPPUNIT_ASSERT_EQUAL(fixture,unitVec[1]);
-
-  fixture = lp2->rotateVector(unitVec[1],1.);
-  CPPUNIT_ASSERT_EQUAL(fixture,unitVec[1]);
-  fixture = lp2->rotateVector(unitVec[2],1.);
-  CPPUNIT_ASSERT_EQUAL(fixture,unitVec[2]);
-
-  fixture = lp3->rotateVector(unitVec[2],1.);
-  CPPUNIT_ASSERT_EQUAL(fixture,unitVec[2]);
-  fixture = lp3->rotateVector(unitVec[0],1.);
-  CPPUNIT_ASSERT_EQUAL(fixture,unitVec[0]);
+  fixture = lp1->rotateVector(e1,1.);
+  CPPUNIT_ASSERT_EQUAL(fixture,e1);
+  fixture = lp1->rotateVector(e2,1.);
+  CPPUNIT_ASSERT_EQUAL(fixture,e2);
+
+  fixture = lp2->rotateVector(e2,1.);
+  CPPUNIT_ASSERT_EQUAL(fixture,e2);
+  fixture = lp2->rotateVector(e3,1.);
+  CPPUNIT_ASSERT_EQUAL(fixture,e3);
+
+  fixture = lp3->rotateVector(e3,1.);
+  CPPUNIT_ASSERT_EQUAL(fixture,e3);
+  fixture = lp3->rotateVector(e1,1.);
+  CPPUNIT_ASSERT_EQUAL(fixture,e1);
 
   // the real stuff
@@ -353 +353 @@
   CPPUNIT_ASSERT_EQUAL(fixture,zeroVec);
 }
-
-void LineUnittest::sphereIntersectionTest(){
-  {
-    std::vector<Vector> res = la1->getSphereIntersections();
-    CPPUNIT_ASSERT_EQUAL(res.size(),(size_t)2);
-    CPPUNIT_ASSERT(testDirection(res[0],unitVec[0]));
-    CPPUNIT_ASSERT(testDirection(res[1],unitVec[0]));
-    CPPUNIT_ASSERT(res[0]!=res[1]);
-  }
-
-  {
-    std::vector<Vector> res = la2->getSphereIntersections();
-    CPPUNIT_ASSERT_EQUAL(res.size(),(size_t)2);
-    CPPUNIT_ASSERT(testDirection(res[0],unitVec[1]));
-    CPPUNIT_ASSERT(testDirection(res[1],unitVec[1]));
-    CPPUNIT_ASSERT(res[0]!=res[1]);
-  }
-
-  {
-    std::vector<Vector> res = la3->getSphereIntersections();
-    CPPUNIT_ASSERT_EQUAL(res.size(),(size_t)2);
-    CPPUNIT_ASSERT(testDirection(res[0],unitVec[2]));
-    CPPUNIT_ASSERT(testDirection(res[1],unitVec[2]));
-    CPPUNIT_ASSERT(res[0]!=res[1]);
-  }
-
-  {
-    std::vector<Vector> res = lp1->getSphereIntersections();
-    CPPUNIT_ASSERT_EQUAL(res.size(),(size_t)2);
-    CPPUNIT_ASSERT((res[0]==unitVec[0]) || (res[0]==unitVec[1]));
-    CPPUNIT_ASSERT((res[1]==unitVec[0]) || (res[1]==unitVec[1]));
-    CPPUNIT_ASSERT(res[0]!=res[1]);
-  }
-
-  {
-    std::vector<Vector> res = lp2->getSphereIntersections();
-    CPPUNIT_ASSERT_EQUAL(res.size(),(size_t)2);
-    CPPUNIT_ASSERT((res[0]==unitVec[1]) || (res[0]==unitVec[2]));
-    CPPUNIT_ASSERT((res[1]==unitVec[1]) || (res[1]==unitVec[2]));
-    CPPUNIT_ASSERT(res[0]!=res[1]);
-  }
-
-  {
-    std::vector<Vector> res = lp3->getSphereIntersections();
-    CPPUNIT_ASSERT_EQUAL(res.size(),(size_t)2);
-    CPPUNIT_ASSERT((res[0]==unitVec[2]) || (res[0]==unitVec[0]));
-    CPPUNIT_ASSERT((res[1]==unitVec[2]) || (res[1]==unitVec[0]));
-    CPPUNIT_ASSERT(res[0]!=res[1]);
-  }
-}