| 1 | /* | 
|---|
| 2 | * Project: MoleCuilder | 
|---|
| 3 | * Description: creates and alters molecular systems | 
|---|
| 4 | * Copyright (C)  2010 University of Bonn. All rights reserved. | 
|---|
| 5 | * Please see the LICENSE file or "Copyright notice" in builder.cpp for details. | 
|---|
| 6 | */ | 
|---|
| 7 |  | 
|---|
| 8 | /* | 
|---|
| 9 | * BaseShapesUnittest.cpp | 
|---|
| 10 | * | 
|---|
| 11 | *  Created on: Jun 18, 2010 | 
|---|
| 12 | *      Author: crueger | 
|---|
| 13 | */ | 
|---|
| 14 |  | 
|---|
| 15 | #include "BaseShapesUnittest.hpp" | 
|---|
| 16 |  | 
|---|
| 17 | #include <cppunit/CompilerOutputter.h> | 
|---|
| 18 | #include <cppunit/extensions/TestFactoryRegistry.h> | 
|---|
| 19 | #include <cppunit/ui/text/TestRunner.h> | 
|---|
| 20 |  | 
|---|
| 21 | #include <cmath> | 
|---|
| 22 |  | 
|---|
| 23 | #ifdef HAVE_TESTRUNNER | 
|---|
| 24 | #include "UnitTestMain.hpp" | 
|---|
| 25 | #endif /*HAVE_TESTRUNNER*/ | 
|---|
| 26 |  | 
|---|
| 27 | #include "LinearAlgebra/Vector.hpp" | 
|---|
| 28 | #include "Shapes/Shape.hpp" | 
|---|
| 29 |  | 
|---|
| 30 | #include "Shapes/BaseShapes.hpp" | 
|---|
| 31 |  | 
|---|
| 32 | // Registers the fixture into the 'registry' | 
|---|
| 33 | CPPUNIT_TEST_SUITE_REGISTRATION( BaseShapesUnittest ); | 
|---|
| 34 |  | 
|---|
| 35 | void BaseShapesUnittest::setUp() | 
|---|
| 36 | { | 
|---|
| 37 | v000 =  0*e1+0*e2+0*e3; | 
|---|
| 38 | v100 =  1*e1+0*e2+0*e3; | 
|---|
| 39 | v200 = -1*e1+0*e2+0*e3; | 
|---|
| 40 | v010 =  0*e1+1*e2+0*e3; | 
|---|
| 41 | v110 =  1*e1+1*e2+0*e3; | 
|---|
| 42 | v210 = -1*e1+1*e2+0*e3; | 
|---|
| 43 | v020 =  0*e1-1*e2+0*e3; | 
|---|
| 44 | v120 =  1*e1-1*e2+0*e3; | 
|---|
| 45 | v220 = -1*e1-1*e2+0*e3; | 
|---|
| 46 | v001 =  0*e1+0*e2+1*e3; | 
|---|
| 47 | v101 =  1*e1+0*e2+1*e3; | 
|---|
| 48 | v201 = -1*e1+0*e2+1*e3; | 
|---|
| 49 | v011 =  0*e1+1*e2+1*e3; | 
|---|
| 50 | v111 =  1*e1+1*e2+1*e3; | 
|---|
| 51 | v211 = -1*e1+1*e2+1*e3; | 
|---|
| 52 | v021 =  0*e1-1*e2+1*e3; | 
|---|
| 53 | v121 =  1*e1-1*e2+1*e3; | 
|---|
| 54 | v221 = -1*e1-1*e2+1*e3; | 
|---|
| 55 | v002 =  0*e1+0*e2-1*e3; | 
|---|
| 56 | v102 =  1*e1+0*e2-1*e3; | 
|---|
| 57 | v202 = -1*e1+0*e2-1*e3; | 
|---|
| 58 | v012 =  0*e1+1*e2-1*e3; | 
|---|
| 59 | v112 =  1*e1+1*e2-1*e3; | 
|---|
| 60 | v212 = -1*e1+1*e2-1*e3; | 
|---|
| 61 | v022 =  0*e1-1*e2-1*e3; | 
|---|
| 62 | v122 =  1*e1-1*e2-1*e3; | 
|---|
| 63 | v222 = -1*e1-1*e2-1*e3; | 
|---|
| 64 | } | 
|---|
| 65 |  | 
|---|
| 66 | void BaseShapesUnittest::tearDown() | 
|---|
| 67 | {} | 
|---|
| 68 |  | 
|---|
| 69 |  | 
|---|
| 70 | void BaseShapesUnittest::isInsideTest(){ | 
|---|
| 71 | Shape s1 = Sphere(); | 
|---|
| 72 |  | 
|---|
| 73 | CPPUNIT_ASSERT(s1.isInside(v000)); | 
|---|
| 74 | CPPUNIT_ASSERT(s1.isInside(v100)); | 
|---|
| 75 | CPPUNIT_ASSERT(s1.isInside(v200)); | 
|---|
| 76 | CPPUNIT_ASSERT(s1.isInside(v010)); | 
|---|
| 77 | CPPUNIT_ASSERT(!s1.isInside(v110)); | 
|---|
| 78 | CPPUNIT_ASSERT(!s1.isInside(v210)); | 
|---|
| 79 | CPPUNIT_ASSERT(s1.isInside(v020)); | 
|---|
| 80 | CPPUNIT_ASSERT(!s1.isInside(v120)); | 
|---|
| 81 | CPPUNIT_ASSERT(!s1.isInside(v220)); | 
|---|
| 82 | CPPUNIT_ASSERT(s1.isInside(v001)); | 
|---|
| 83 | CPPUNIT_ASSERT(!s1.isInside(v101)); | 
|---|
| 84 | CPPUNIT_ASSERT(!s1.isInside(v201)); | 
|---|
| 85 | CPPUNIT_ASSERT(!s1.isInside(v011)); | 
|---|
| 86 | CPPUNIT_ASSERT(!s1.isInside(v111)); | 
|---|
| 87 | CPPUNIT_ASSERT(!s1.isInside(v211)); | 
|---|
| 88 | CPPUNIT_ASSERT(!s1.isInside(v021)); | 
|---|
| 89 | CPPUNIT_ASSERT(!s1.isInside(v121)); | 
|---|
| 90 | CPPUNIT_ASSERT(!s1.isInside(v221)); | 
|---|
| 91 | CPPUNIT_ASSERT(s1.isInside(v002)); | 
|---|
| 92 | CPPUNIT_ASSERT(!s1.isInside(v102)); | 
|---|
| 93 | CPPUNIT_ASSERT(!s1.isInside(v202)); | 
|---|
| 94 | CPPUNIT_ASSERT(!s1.isInside(v012)); | 
|---|
| 95 | CPPUNIT_ASSERT(!s1.isInside(v112)); | 
|---|
| 96 | CPPUNIT_ASSERT(!s1.isInside(v212)); | 
|---|
| 97 | CPPUNIT_ASSERT(!s1.isInside(v022)); | 
|---|
| 98 | CPPUNIT_ASSERT(!s1.isInside(v122)); | 
|---|
| 99 | CPPUNIT_ASSERT(!s1.isInside(v222)); | 
|---|
| 100 |  | 
|---|
| 101 | Shape s2 = Cuboid(); | 
|---|
| 102 |  | 
|---|
| 103 | CPPUNIT_ASSERT(s2.isInside(v000)); | 
|---|
| 104 | CPPUNIT_ASSERT(s2.isInside(v100)); | 
|---|
| 105 | CPPUNIT_ASSERT(!s2.isInside(v200)); | 
|---|
| 106 | CPPUNIT_ASSERT(s2.isInside(v010)); | 
|---|
| 107 | CPPUNIT_ASSERT(s2.isInside(v110)); | 
|---|
| 108 | CPPUNIT_ASSERT(!s2.isInside(v210)); | 
|---|
| 109 | CPPUNIT_ASSERT(!s2.isInside(v020)); | 
|---|
| 110 | CPPUNIT_ASSERT(!s2.isInside(v120)); | 
|---|
| 111 | CPPUNIT_ASSERT(!s2.isInside(v220)); | 
|---|
| 112 | CPPUNIT_ASSERT(s2.isInside(v001)); | 
|---|
| 113 | CPPUNIT_ASSERT(s2.isInside(v101)); | 
|---|
| 114 | CPPUNIT_ASSERT(!s2.isInside(v201)); | 
|---|
| 115 | CPPUNIT_ASSERT(s2.isInside(v011)); | 
|---|
| 116 | CPPUNIT_ASSERT(s2.isInside(v111)); | 
|---|
| 117 | CPPUNIT_ASSERT(!s2.isInside(v211)); | 
|---|
| 118 | CPPUNIT_ASSERT(!s2.isInside(v021)); | 
|---|
| 119 | CPPUNIT_ASSERT(!s2.isInside(v121)); | 
|---|
| 120 | CPPUNIT_ASSERT(!s2.isInside(v221)); | 
|---|
| 121 | CPPUNIT_ASSERT(!s2.isInside(v002)); | 
|---|
| 122 | CPPUNIT_ASSERT(!s2.isInside(v102)); | 
|---|
| 123 | CPPUNIT_ASSERT(!s2.isInside(v202)); | 
|---|
| 124 | CPPUNIT_ASSERT(!s2.isInside(v012)); | 
|---|
| 125 | CPPUNIT_ASSERT(!s2.isInside(v112)); | 
|---|
| 126 | CPPUNIT_ASSERT(!s2.isInside(v212)); | 
|---|
| 127 | CPPUNIT_ASSERT(!s2.isInside(v022)); | 
|---|
| 128 | CPPUNIT_ASSERT(!s2.isInside(v122)); | 
|---|
| 129 | CPPUNIT_ASSERT(!s2.isInside(v222)); | 
|---|
| 130 |  | 
|---|
| 131 | } | 
|---|
| 132 |  | 
|---|
| 133 | void BaseShapesUnittest::PointsOnSurfaceTest(){ | 
|---|
| 134 | Shape s = Sphere(); | 
|---|
| 135 | const size_t N = 200; | 
|---|
| 136 | std::vector<Vector> PointsOnSurface = s.getHomogeneousPointsOnSurface(N); | 
|---|
| 137 | for (std::vector<Vector>::const_iterator iter = PointsOnSurface.begin(); iter != PointsOnSurface.end(); ++iter) { | 
|---|
| 138 | CPPUNIT_ASSERT(fabs(1. - (*iter).NormSquared()) < MYEPSILON); | 
|---|
| 139 | } | 
|---|
| 140 | CPPUNIT_ASSERT_EQUAL(N, PointsOnSurface.size()); | 
|---|
| 141 | } | 
|---|
| 142 |  | 
|---|