test/brains/DataStorageTestCase.cpp

#include "brains/DataStorageTestCase.hpp"
// Registers the fixture into the 'registry'
CPPUNIT_TEST_SUITE_REGISTRATION( DataStorageTestCase );

void DataStorageTestCase::testDataStorage() {
    //test constructor
    DataStorage ds1(10020);
    CPPUNIT_ASSERT_EQUAL(ds1.getStorageLayout(), 0x11111111);
    CPPUNIT_ASSERT_EQUAL(ds1.getSize(), 10020);

    //test resize
    ds1.resize(56);
    CPPUNIT_ASSERT_EQUAL(ds1.getSize(), 56);

    //test reserve
    
    DataStorage ds2(10000, DataStorage::dslForce |DataStorage::dslAmat |DataStorage::dslZAngle);

    CPPUNIT_ASSERT(!(ds2.getStorageLayout() & DataStorage::dslPosition));
    CPPUNIT_ASSERT(!(ds2.getStorageLayout() & DataStorage::dslVelocity));
    CPPUNIT_ASSERT(ds2.getStorageLayout() & DataStorage::dslAmat);
    CPPUNIT_ASSERT(!(ds2.getStorageLayout() & DataStorage::dslAngularMomentum));
    CPPUNIT_ASSERT(!(ds2.getStorageLayout() & DataStorage::dslUnitVector));
    CPPUNIT_ASSERT(ds2.getStorageLayout() & DataStorage::dslZAngle);    
    CPPUNIT_ASSERT(ds2.getStorageLayout() & DataStorage::dslForce);
    CPPUNIT_ASSERT(!(ds2.getStorageLayout() & DataStorage::dslTorque));


    CPPUNIT_ASSERT(ds2.position.size() == 0);
    CPPUNIT_ASSERT(ds2.velocity.size() == 0);
    CPPUNIT_ASSERT(ds2.aMat.size() == 10000);
    CPPUNIT_ASSERT(ds2.angularMomentum.size() == 0);
    CPPUNIT_ASSERT(ds2.unitVector.size() == 0);
    CPPUNIT_ASSERT(ds2.zAngle.size() == 10000);
    CPPUNIT_ASSERT(ds2.force.size() == 10000);
    CPPUNIT_ASSERT(ds2.torque.size() == 0);
    
    //test getArrayPointer

    double* pamat = ds2.getArrayPointer(DataStorage::dslAmat);
    double* pzangle = ds2.getArrayPointer(DataStorage::dslZAngle);
    double* pforce = ds2.getArrayPointer(DataStorage::dslForce);

    pamat[9] = 2.0993;
    CPPUNIT_ASSERT_DOUBLES_EQUAL(ds2.aMat[1](0,0), 2.0993, OpenMD::NumericConstant::epsilon);

    pzangle[9743] = 8464353.129574;
    CPPUNIT_ASSERT_DOUBLES_EQUAL(ds2.zAngle[9743], 8464353.129574, OpenMD::NumericConstant::epsilon);

    pforce[12] = 345546.3413263;
    CPPUNIT_ASSERT_DOUBLES_EQUAL(ds2.force[4][0], 345546.3413263, OpenMD::NumericConstant::epsilon);

    //test copy
    ds2.zAngle[0] = 7.0;
    ds2.zAngle[1] = 3.5;

    ds2.force[0][0]= 6;
    ds2.force[0][1]= 7.0;
    ds2.force[0][2]= 8.0;
    
    ds2.force[1][0]= 9.0;
    ds2.force[1][1]= 10.0;
    ds2.force[1][2]= 11.0;
    
    ds2.aMat[0](0,0) = 1.0;
    ds2.aMat[0](0,1) = 2.0;
    ds2.aMat[0](0,2) = 3.0;
    ds2.aMat[0](1,0) = 4.0;
    ds2.aMat[0](1,1) = 5.0;
    ds2.aMat[0](1,2) = 6.0;
    ds2.aMat[0](2,0) = 7.0;
    ds2.aMat[0](2,1) = 8.0;
    ds2.aMat[0](2,2) = 9.0;

    ds2.aMat[1](0,0) = 11.0;
    ds2.aMat[1](0,1) = 12.0;
    ds2.aMat[1](0,2) = 13.0;
    ds2.aMat[1](1,0) = 14.0;
    ds2.aMat[1](1,1) = 15.0;
    ds2.aMat[1](1,2) = 16.0;
    ds2.aMat[1](2,0) = 17.0;
    ds2.aMat[1](2,1) = 18.0;
    ds2.aMat[1](2,2) = 19.0;

    ds2.copy(0, 2, 10);

    CPPUNIT_ASSERT_DOUBLES_EQUAL(ds2.zAngle[10], 7.0, OpenMD::NumericConstant::epsilon);
    CPPUNIT_ASSERT_DOUBLES_EQUAL(ds2.zAngle[11], 3.5, OpenMD::NumericConstant::epsilon);

    CPPUNIT_ASSERT_DOUBLES_EQUAL(ds2.force[10][0], 6.0, OpenMD::NumericConstant::epsilon);
    CPPUNIT_ASSERT_DOUBLES_EQUAL(ds2.force[10][1], 7.0, OpenMD::NumericConstant::epsilon);
    CPPUNIT_ASSERT_DOUBLES_EQUAL(ds2.force[10][2], 8.0, OpenMD::NumericConstant::epsilon);

    CPPUNIT_ASSERT_DOUBLES_EQUAL(ds2.force[11][0], 9.0, OpenMD::NumericConstant::epsilon);
    CPPUNIT_ASSERT_DOUBLES_EQUAL(ds2.force[11][1], 10.0, OpenMD::NumericConstant::epsilon);
    CPPUNIT_ASSERT_DOUBLES_EQUAL(ds2.force[11][2], 11.0, OpenMD::NumericConstant::epsilon);

    CPPUNIT_ASSERT_DOUBLES_EQUAL(ds2.aMat[10](0, 0), 1.0, OpenMD::NumericConstant::epsilon);
    CPPUNIT_ASSERT_DOUBLES_EQUAL(ds2.aMat[10](0, 1), 2.0, OpenMD::NumericConstant::epsilon);
    CPPUNIT_ASSERT_DOUBLES_EQUAL(ds2.aMat[10](0, 2), 3.0, OpenMD::NumericConstant::epsilon);
    CPPUNIT_ASSERT_DOUBLES_EQUAL(ds2.aMat[10](1, 0), 4.0, OpenMD::NumericConstant::epsilon);
    CPPUNIT_ASSERT_DOUBLES_EQUAL(ds2.aMat[10](1, 1) , 5.0, OpenMD::NumericConstant::epsilon);
    CPPUNIT_ASSERT_DOUBLES_EQUAL(ds2.aMat[10](1, 2), 6.0, OpenMD::NumericConstant::epsilon);
    CPPUNIT_ASSERT_DOUBLES_EQUAL(ds2.aMat[10](2, 0), 7.0, OpenMD::NumericConstant::epsilon);
    CPPUNIT_ASSERT_DOUBLES_EQUAL(ds2.aMat[10](2, 1), 8.0, OpenMD::NumericConstant::epsilon);
    CPPUNIT_ASSERT_DOUBLES_EQUAL(ds2.aMat[10](2, 2), 9.0, OpenMD::NumericConstant::epsilon);

    CPPUNIT_ASSERT_DOUBLES_EQUAL(ds2.aMat[11](0, 0), 11.0, OpenMD::NumericConstant::epsilon);
    CPPUNIT_ASSERT_DOUBLES_EQUAL(ds2.aMat[11](0, 1), 12.0, OpenMD::NumericConstant::epsilon);
    CPPUNIT_ASSERT_DOUBLES_EQUAL(ds2.aMat[11](0, 2), 13.0, OpenMD::NumericConstant::epsilon);
    CPPUNIT_ASSERT_DOUBLES_EQUAL(ds2.aMat[11](1, 0), 14.0, OpenMD::NumericConstant::epsilon);
    CPPUNIT_ASSERT_DOUBLES_EQUAL(ds2.aMat[11](1, 1) , 15.0, OpenMD::NumericConstant::epsilon);
    CPPUNIT_ASSERT_DOUBLES_EQUAL(ds2.aMat[11](1, 2), 16.0, OpenMD::NumericConstant::epsilon);
    CPPUNIT_ASSERT_DOUBLES_EQUAL(ds2.aMat[11](2, 0), 17.0, OpenMD::NumericConstant::epsilon);
    CPPUNIT_ASSERT_DOUBLES_EQUAL(ds2.aMat[11](2, 1), 18.0, OpenMD::NumericConstant::epsilon);
    CPPUNIT_ASSERT_DOUBLES_EQUAL(ds2.aMat[11](2, 2), 19.0, OpenMD::NumericConstant::epsilon);

    
}

Revision:	1390
Committed:	Wed Nov 25 20:02:06 2009 UTC (15 years, 5 months ago) by gezelter
Original Path:	*trunk/test/brains/DataStorageTestCase.cpp*
File size:	5341 byte(s)
Log Message:	Almost all of the changes necessary to create OpenMD out of our old project (OOPSE-4)
#	User	Rev	Content
1	tim	155	#include "brains/DataStorageTestCase.hpp"
2			// Registers the fixture into the 'registry'
3			CPPUNIT_TEST_SUITE_REGISTRATION( DataStorageTestCase );
4
5			void DataStorageTestCase::testDataStorage() {
6			//test constructor
7			DataStorage ds1(10020);
8			CPPUNIT_ASSERT_EQUAL(ds1.getStorageLayout(), 0x11111111);
9			CPPUNIT_ASSERT_EQUAL(ds1.getSize(), 10020);
10
11			//test resize
12			ds1.resize(56);
13			CPPUNIT_ASSERT_EQUAL(ds1.getSize(), 56);
14
15			//test reserve
16
17			DataStorage ds2(10000, DataStorage::dslForce \|DataStorage::dslAmat \|DataStorage::dslZAngle);
18
19			CPPUNIT_ASSERT(!(ds2.getStorageLayout() & DataStorage::dslPosition));
20			CPPUNIT_ASSERT(!(ds2.getStorageLayout() & DataStorage::dslVelocity));
21			CPPUNIT_ASSERT(ds2.getStorageLayout() & DataStorage::dslAmat);
22			CPPUNIT_ASSERT(!(ds2.getStorageLayout() & DataStorage::dslAngularMomentum));
23			CPPUNIT_ASSERT(!(ds2.getStorageLayout() & DataStorage::dslUnitVector));
24			CPPUNIT_ASSERT(ds2.getStorageLayout() & DataStorage::dslZAngle);
25			CPPUNIT_ASSERT(ds2.getStorageLayout() & DataStorage::dslForce);
26			CPPUNIT_ASSERT(!(ds2.getStorageLayout() & DataStorage::dslTorque));
27
28
29			CPPUNIT_ASSERT(ds2.position.size() == 0);
30			CPPUNIT_ASSERT(ds2.velocity.size() == 0);
31			CPPUNIT_ASSERT(ds2.aMat.size() == 10000);
32			CPPUNIT_ASSERT(ds2.angularMomentum.size() == 0);
33			CPPUNIT_ASSERT(ds2.unitVector.size() == 0);
34			CPPUNIT_ASSERT(ds2.zAngle.size() == 10000);
35			CPPUNIT_ASSERT(ds2.force.size() == 10000);
36			CPPUNIT_ASSERT(ds2.torque.size() == 0);
37
38			//test getArrayPointer
39
40			double* pamat = ds2.getArrayPointer(DataStorage::dslAmat);
41			double* pzangle = ds2.getArrayPointer(DataStorage::dslZAngle);
42			double* pforce = ds2.getArrayPointer(DataStorage::dslForce);
43
44			pamat[9] = 2.0993;
45	gezelter	1390	CPPUNIT_ASSERT_DOUBLES_EQUAL(ds2.aMat[1](0,0), 2.0993, OpenMD::NumericConstant::epsilon);
46	tim	155
47			pzangle[9743] = 8464353.129574;
48	gezelter	1390	CPPUNIT_ASSERT_DOUBLES_EQUAL(ds2.zAngle[9743], 8464353.129574, OpenMD::NumericConstant::epsilon);
49	tim	155
50			pforce[12] = 345546.3413263;
51	gezelter	1390	CPPUNIT_ASSERT_DOUBLES_EQUAL(ds2.force[4][0], 345546.3413263, OpenMD::NumericConstant::epsilon);
52	tim	155
53			//test copy
54			ds2.zAngle[0] = 7.0;
55			ds2.zAngle[1] = 3.5;
56
57			ds2.force[0][0]= 6;
58			ds2.force[0][1]= 7.0;
59			ds2.force[0][2]= 8.0;
60
61			ds2.force[1][0]= 9.0;
62			ds2.force[1][1]= 10.0;
63			ds2.force[1][2]= 11.0;
64
65			ds2.aMat[0](0,0) = 1.0;
66			ds2.aMat[0](0,1) = 2.0;
67			ds2.aMat[0](0,2) = 3.0;
68			ds2.aMat[0](1,0) = 4.0;
69			ds2.aMat[0](1,1) = 5.0;
70			ds2.aMat[0](1,2) = 6.0;
71			ds2.aMat[0](2,0) = 7.0;
72			ds2.aMat[0](2,1) = 8.0;
73			ds2.aMat[0](2,2) = 9.0;
74
75			ds2.aMat[1](0,0) = 11.0;
76			ds2.aMat[1](0,1) = 12.0;
77			ds2.aMat[1](0,2) = 13.0;
78			ds2.aMat[1](1,0) = 14.0;
79			ds2.aMat[1](1,1) = 15.0;
80			ds2.aMat[1](1,2) = 16.0;
81			ds2.aMat[1](2,0) = 17.0;
82			ds2.aMat[1](2,1) = 18.0;
83			ds2.aMat[1](2,2) = 19.0;
84
85			ds2.copy(0, 2, 10);
86
87	gezelter	1390	CPPUNIT_ASSERT_DOUBLES_EQUAL(ds2.zAngle[10], 7.0, OpenMD::NumericConstant::epsilon);
88			CPPUNIT_ASSERT_DOUBLES_EQUAL(ds2.zAngle[11], 3.5, OpenMD::NumericConstant::epsilon);
89	tim	155
90	gezelter	1390	CPPUNIT_ASSERT_DOUBLES_EQUAL(ds2.force[10][0], 6.0, OpenMD::NumericConstant::epsilon);
91			CPPUNIT_ASSERT_DOUBLES_EQUAL(ds2.force[10][1], 7.0, OpenMD::NumericConstant::epsilon);
92			CPPUNIT_ASSERT_DOUBLES_EQUAL(ds2.force[10][2], 8.0, OpenMD::NumericConstant::epsilon);
93	tim	155
94	gezelter	1390	CPPUNIT_ASSERT_DOUBLES_EQUAL(ds2.force[11][0], 9.0, OpenMD::NumericConstant::epsilon);
95			CPPUNIT_ASSERT_DOUBLES_EQUAL(ds2.force[11][1], 10.0, OpenMD::NumericConstant::epsilon);
96			CPPUNIT_ASSERT_DOUBLES_EQUAL(ds2.force[11][2], 11.0, OpenMD::NumericConstant::epsilon);
97	tim	155
98	gezelter	1390	CPPUNIT_ASSERT_DOUBLES_EQUAL(ds2.aMat[10](0, 0), 1.0, OpenMD::NumericConstant::epsilon);
99			CPPUNIT_ASSERT_DOUBLES_EQUAL(ds2.aMat[10](0, 1), 2.0, OpenMD::NumericConstant::epsilon);
100			CPPUNIT_ASSERT_DOUBLES_EQUAL(ds2.aMat[10](0, 2), 3.0, OpenMD::NumericConstant::epsilon);
101			CPPUNIT_ASSERT_DOUBLES_EQUAL(ds2.aMat[10](1, 0), 4.0, OpenMD::NumericConstant::epsilon);
102			CPPUNIT_ASSERT_DOUBLES_EQUAL(ds2.aMat[10](1, 1) , 5.0, OpenMD::NumericConstant::epsilon);
103			CPPUNIT_ASSERT_DOUBLES_EQUAL(ds2.aMat[10](1, 2), 6.0, OpenMD::NumericConstant::epsilon);
104			CPPUNIT_ASSERT_DOUBLES_EQUAL(ds2.aMat[10](2, 0), 7.0, OpenMD::NumericConstant::epsilon);
105			CPPUNIT_ASSERT_DOUBLES_EQUAL(ds2.aMat[10](2, 1), 8.0, OpenMD::NumericConstant::epsilon);
106			CPPUNIT_ASSERT_DOUBLES_EQUAL(ds2.aMat[10](2, 2), 9.0, OpenMD::NumericConstant::epsilon);
107	tim	155
108	gezelter	1390	CPPUNIT_ASSERT_DOUBLES_EQUAL(ds2.aMat[11](0, 0), 11.0, OpenMD::NumericConstant::epsilon);
109			CPPUNIT_ASSERT_DOUBLES_EQUAL(ds2.aMat[11](0, 1), 12.0, OpenMD::NumericConstant::epsilon);
110			CPPUNIT_ASSERT_DOUBLES_EQUAL(ds2.aMat[11](0, 2), 13.0, OpenMD::NumericConstant::epsilon);
111			CPPUNIT_ASSERT_DOUBLES_EQUAL(ds2.aMat[11](1, 0), 14.0, OpenMD::NumericConstant::epsilon);
112			CPPUNIT_ASSERT_DOUBLES_EQUAL(ds2.aMat[11](1, 1) , 15.0, OpenMD::NumericConstant::epsilon);
113			CPPUNIT_ASSERT_DOUBLES_EQUAL(ds2.aMat[11](1, 2), 16.0, OpenMD::NumericConstant::epsilon);
114			CPPUNIT_ASSERT_DOUBLES_EQUAL(ds2.aMat[11](2, 0), 17.0, OpenMD::NumericConstant::epsilon);
115			CPPUNIT_ASSERT_DOUBLES_EQUAL(ds2.aMat[11](2, 1), 18.0, OpenMD::NumericConstant::epsilon);
116			CPPUNIT_ASSERT_DOUBLES_EQUAL(ds2.aMat[11](2, 2), 19.0, OpenMD::NumericConstant::epsilon);
117	tim	155
118
119			}
120