test/DataConstant/DataConstantTestCase.cpp

// $Id$
/* 
 *****************************************************************************
 *                                                                           *
 *       COPYRIGHT  ACcESS  -  All Rights Reserved                           *
 *                                                                           *
 * This software is the property of ACcESS. No part of this code             *
 * may be copied in any form or by any means without the expressed written   *
 * consent of ACcESS.  Copying, use or modification of this software         *
 * by any unauthorised person is illegal unless that person has a software   *
 * license agreement with ACcESS.                                            *
 *                                                                           *
 *****************************************************************************
*/
#include "escript/Data/DataConstant.h"
#include "escript/Data/FunctionSpace.h"
#include "esysUtils/EsysException.h"

#include "DataConstantTestCase.h"

#include <iostream>

using namespace CppUnitTest;
using namespace escript;
using namespace std;
using namespace esysUtils;

void DataConstantTestCase::setUp() {
  //
  // This is called before each test is run
}

void DataConstantTestCase::tearDown() {
  //
  // This is called after each test has been run
}

void DataConstantTestCase::testAll() {

  cout << endl;

  //
  // Create a scalar pointData
  DataArrayView::ShapeType shape;
  DataArrayView::ValueType data(DataArrayView::noValues(shape),0);
  DataArrayView pointData(data,shape);

  //
  // assign an arbitrary value
  pointData()=1.0;

  //
  // Test construction
  cout << "\tTesting default constructor." << endl;
  DataConstant testData(pointData, FunctionSpace());

  cout << "\tTest getLength." << endl;
  assert(testData.getLength()==1);

  cout << "\tTest reshape." << endl;
  shape.push_back(2);
  shape.push_back(3);
  shape.push_back(21);
  testData.reshapeDataPoint(shape);
  assert(testData.getPointDataView().getRank()==shape.size());

  cout << "\tTest getPointDataView." << endl;
  for (int k=0;k<shape[2];k++) {
    for (int j=0;j<shape[1];j++) {
      for (int i=0;i<shape[0];i++) {
    assert(testData.getPointDataView()(i,j,k)==pointData());
      }
    }
  }

  try {
    cout << "\tTest illegal reshape." << endl;
    testData.reshapeDataPoint(shape);
    assert(false);
  }
  catch (EsysException& e) {
    //cout << e.toString() << endl;
    assert(true);
  }

  cout << "\tVerify data point attributes." << endl;
  DataArrayView dataView=testData.getPointDataView();
  assert(dataView.getRank()==3);
  assert(dataView.noValues()==126);
  assert(dataView.getShape()[0]==2);
  assert(dataView.getShape()[1]==3);
  assert(dataView.getShape()[2]==21);

  cout << "\tTesting alternative constructor." << endl;
  DataArrayView::ValueType data1(DataArrayView::noValues(shape),1.0);
  DataConstant testData1(FunctionSpace(), shape, data1);

  for (int k=0;k<shape[2];k++) {
    for (int j=0;j<shape[1];j++) {
      for (int i=0;i<shape[0];i++) {
    assert(testData1.getPointDataView()(i,j,k)==1.0);
      }
    }
  }

  cout << "\tTest getLength." << endl;
  assert(testData1.getLength()==126);

  cout << "\tVerify data point attributes." << endl;
  dataView=testData1.getPointDataView();
  assert(dataView.getRank()==3);
  assert(dataView.noValues()==126);
  assert(dataView.getShape()[0]==2);
  assert(dataView.getShape()[1]==3);
  assert(dataView.getShape()[2]==21);

  cout << "\tTesting copy constructor." << endl;
  DataConstant testData2(testData);

  for (int k=0;k<shape[2];k++) {
    for (int j=0;j<shape[1];j++) {
      for (int i=0;i<shape[0];i++) {
    assert(testData2.getPointDataView()(i,j,k)==pointData());
      }
    }
  }

  cout << "\tTest getLength." << endl;
  assert(testData2.getLength()==126);

  cout << "\tVerify data point attributes." << endl;
  dataView=testData2.getPointDataView();
  assert(dataView.getRank()==3);
  assert(dataView.noValues()==126);
  assert(dataView.getShape()[0]==2);
  assert(dataView.getShape()[1]==3);
  assert(dataView.getShape()[2]==21);

  cout << "\tTest slicing (whole object)." << endl;

  DataArrayView::RegionType region;
  region.push_back(DataArrayView::RegionType::value_type(0,shape[0]));
  region.push_back(DataArrayView::RegionType::value_type(0,shape[1]));
  region.push_back(DataArrayView::RegionType::value_type(0,shape[2]));

  DataAbstract* testData3=testData2.getSlice(region);

  for (int k=0;k<shape[2];k++) {
    for (int j=0;j<shape[1];j++) {
      for (int i=0;i<shape[0];i++) {
    assert(testData3->getPointDataView()(i,j,k)==pointData());
      }
    }
  }

  assert(testData3->getLength()==126);

  cout << "\tVerify data point attributes." << endl;
  dataView=testData3->getPointDataView();
  assert(dataView.getRank()==3);
  assert(dataView.noValues()==126);
  assert(dataView.getShape()[0]==2);
  assert(dataView.getShape()[1]==3);
  assert(dataView.getShape()[2]==21);

  cout << "\tTest slicing (part object)." << endl;

  DataArrayView::RegionType region2;
  region2.push_back(DataArrayView::RegionType::value_type(0,2));
  region2.push_back(DataArrayView::RegionType::value_type(0,2));
  region2.push_back(DataArrayView::RegionType::value_type(0,2));

  DataAbstract* testData4=testData3->getSlice(region2);

  for (int k=0;k<2;k++) {
    for (int j=0;j<2;j++) {
      for (int i=0;i<2;i++) {
    assert(testData4->getPointDataView()(i,j,k)==pointData());
      }
    }
  }

  assert(testData4->getLength()==8);

  cout << "\tVerify data point attributes." << endl;
  dataView=testData4->getPointDataView();
  assert(dataView.getRank()==3);
  assert(dataView.noValues()==8);
  assert(dataView.getShape()[0]==2);
  assert(dataView.getShape()[1]==2);
  assert(dataView.getShape()[2]==2);

  cout << "\tTest slicing (part object)." << endl;

  DataArrayView::RegionType region3;
  region3.push_back(DataArrayView::RegionType::value_type(1,2));
  region3.push_back(DataArrayView::RegionType::value_type(1,3));
  region3.push_back(DataArrayView::RegionType::value_type(5,9));

  DataAbstract* testData5=testData3->getSlice(region3);

  for (int k=0;k<4;k++) {
    for (int j=0;j<2;j++) {
      for (int i=0;i<1;i++) {
    assert(testData5->getPointDataView()(i,j,k)==pointData());
      }
    }
  }

  assert(testData5->getLength()==8);

  cout << "\tVerify data point attributes." << endl;
  dataView=testData5->getPointDataView();
  assert(dataView.getRank()==3);
  assert(dataView.noValues()==8);
  assert(dataView.getShape()[0]==1);
  assert(dataView.getShape()[1]==2);
  assert(dataView.getShape()[2]==4);

  cout << "\tTest slice setting (1)." << endl;

  DataArrayView::RegionType region4;
  region4.push_back(DataArrayView::RegionType::value_type(0,1));
  region4.push_back(DataArrayView::RegionType::value_type(0,2));
  region4.push_back(DataArrayView::RegionType::value_type(0,4));

  DataAbstract* testData6=testData3->getSlice(region3);

  testData5->setSlice(testData6,region4);

  for (int k=0;k<4;k++) {
    for (int j=0;j<2;j++) {
      for (int i=0;i<1;i++) {
    assert(testData5->getPointDataView()(i,j,k)==pointData());
      }
    }
  }

  assert(testData5->getLength()==8);

  cout << "\tVerify data point attributes." << endl;
  dataView=testData5->getPointDataView();
  assert(dataView.getRank()==3);
  assert(dataView.noValues()==8);
  assert(dataView.getShape()[0]==1);
  assert(dataView.getShape()[1]==2);
  assert(dataView.getShape()[2]==4);

  //cout << "\tTest slice setting (2)." << endl;

  //DataArrayView::RegionType region5;
  //region5.push_back(DataArrayView::RegionType::value_type(1,2));
  //region5.push_back(DataArrayView::RegionType::value_type(1,3));
  //region5.push_back(DataArrayView::RegionType::value_type(14,18));

  //DataAbstract* testData7=testData3->getSlice(region);

  //testData5->setSlice(testData7,region5);

  //for (int k=0;k<4;k++) {
  //  for (int j=0;j<2;j++) {
  //    for (int i=0;i<1;i++) {
  //      assert(testData5->getPointDataView()(i,j,k)==pointData());
  //    }
  //  }
  // }

  //assert(testData5->getLength()==8);

  //cout << "\tVerify data point attributes." << endl;
  //dataView=testData5->getPointDataView();
  //assert(dataView.getRank()==3);
  //assert(dataView.noValues()==8);
  //assert(dataView.getShape()[0]==1);
  //assert(dataView.getShape()[1]==2);
  //assert(dataView.getShape()[2]==4);

}

TestSuite* DataConstantTestCase::suite ()
{
  //
  // create the suite of tests to perform.
  TestSuite *testSuite = new TestSuite ("DataConstantTestCase");

  testSuite->addTest (new TestCaller< DataConstantTestCase>("testAll",&DataConstantTestCase::testAll));
  return testSuite;
}
1	// $Id$
2	/*
3	*****************************************************************************
4	* *
5	* COPYRIGHT ACcESS - All Rights Reserved *
6	* *
7	* This software is the property of ACcESS. No part of this code *
8	* may be copied in any form or by any means without the expressed written *
9	* consent of ACcESS. Copying, use or modification of this software *
10	* by any unauthorised person is illegal unless that person has a software *
11	* license agreement with ACcESS. *
12	* *
13	*****************************************************************************
14	*/
15	#include "escript/Data/DataConstant.h"
16	#include "escript/Data/FunctionSpace.h"
17	#include "esysUtils/EsysException.h"
18
19	#include "DataConstantTestCase.h"
20
21	#include <iostream>
22
23	using namespace CppUnitTest;
24	using namespace escript;
25	using namespace std;
26	using namespace esysUtils;
27
28	void DataConstantTestCase::setUp() {
29	//
30	// This is called before each test is run
31	}
32
33	void DataConstantTestCase::tearDown() {
34	//
35	// This is called after each test has been run
36	}
37
38	void DataConstantTestCase::testAll() {
39
40	cout << endl;
41
42	//
43	// Create a scalar pointData
44	DataArrayView::ShapeType shape;
45	DataArrayView::ValueType data(DataArrayView::noValues(shape),0);
46	DataArrayView pointData(data,shape);
47
48	//
49	// assign an arbitrary value
50	pointData()=1.0;
51
52	//
53	// Test construction
54	cout << "\tTesting default constructor." << endl;
55	DataConstant testData(pointData, FunctionSpace());
56
57	cout << "\tTest getLength." << endl;
58	assert(testData.getLength()==1);
59
60	cout << "\tTest reshape." << endl;
61	shape.push_back(2);
62	shape.push_back(3);
63	shape.push_back(21);
64	testData.reshapeDataPoint(shape);
65	assert(testData.getPointDataView().getRank()==shape.size());
66
67	cout << "\tTest getPointDataView." << endl;
68	for (int k=0;k<shape[2];k++) {
69	for (int j=0;j<shape[1];j++) {
70	for (int i=0;i<shape[0];i++) {
71	assert(testData.getPointDataView()(i,j,k)==pointData());
72	}
73	}
74	}
75
76	try {
77	cout << "\tTest illegal reshape." << endl;
78	testData.reshapeDataPoint(shape);
79	assert(false);
80	}
81	catch (EsysException& e) {
82	//cout << e.toString() << endl;
83	assert(true);
84	}
85
86	cout << "\tVerify data point attributes." << endl;
87	DataArrayView dataView=testData.getPointDataView();
88	assert(dataView.getRank()==3);
89	assert(dataView.noValues()==126);
90	assert(dataView.getShape()[0]==2);
91	assert(dataView.getShape()[1]==3);
92	assert(dataView.getShape()[2]==21);
93
94	cout << "\tTesting alternative constructor." << endl;
95	DataArrayView::ValueType data1(DataArrayView::noValues(shape),1.0);
96	DataConstant testData1(FunctionSpace(), shape, data1);
97
98	for (int k=0;k<shape[2];k++) {
99	for (int j=0;j<shape[1];j++) {
100	for (int i=0;i<shape[0];i++) {
101	assert(testData1.getPointDataView()(i,j,k)==1.0);
102	}
103	}
104	}
105
106	cout << "\tTest getLength." << endl;
107	assert(testData1.getLength()==126);
108
109	cout << "\tVerify data point attributes." << endl;
110	dataView=testData1.getPointDataView();
111	assert(dataView.getRank()==3);
112	assert(dataView.noValues()==126);
113	assert(dataView.getShape()[0]==2);
114	assert(dataView.getShape()[1]==3);
115	assert(dataView.getShape()[2]==21);
116
117	cout << "\tTesting copy constructor." << endl;
118	DataConstant testData2(testData);
119
120	for (int k=0;k<shape[2];k++) {
121	for (int j=0;j<shape[1];j++) {
122	for (int i=0;i<shape[0];i++) {
123	assert(testData2.getPointDataView()(i,j,k)==pointData());
124	}
125	}
126	}
127
128	cout << "\tTest getLength." << endl;
129	assert(testData2.getLength()==126);
130
131	cout << "\tVerify data point attributes." << endl;
132	dataView=testData2.getPointDataView();
133	assert(dataView.getRank()==3);
134	assert(dataView.noValues()==126);
135	assert(dataView.getShape()[0]==2);
136	assert(dataView.getShape()[1]==3);
137	assert(dataView.getShape()[2]==21);
138
139	cout << "\tTest slicing (whole object)." << endl;
140
141	DataArrayView::RegionType region;
142	region.push_back(DataArrayView::RegionType::value_type(0,shape[0]));
143	region.push_back(DataArrayView::RegionType::value_type(0,shape[1]));
144	region.push_back(DataArrayView::RegionType::value_type(0,shape[2]));
145
146	DataAbstract* testData3=testData2.getSlice(region);
147
148	for (int k=0;k<shape[2];k++) {
149	for (int j=0;j<shape[1];j++) {
150	for (int i=0;i<shape[0];i++) {
151	assert(testData3->getPointDataView()(i,j,k)==pointData());
152	}
153	}
154	}
155
156	assert(testData3->getLength()==126);
157
158	cout << "\tVerify data point attributes." << endl;
159	dataView=testData3->getPointDataView();
160	assert(dataView.getRank()==3);
161	assert(dataView.noValues()==126);
162	assert(dataView.getShape()[0]==2);
163	assert(dataView.getShape()[1]==3);
164	assert(dataView.getShape()[2]==21);
165
166	cout << "\tTest slicing (part object)." << endl;
167
168	DataArrayView::RegionType region2;
169	region2.push_back(DataArrayView::RegionType::value_type(0,2));
170	region2.push_back(DataArrayView::RegionType::value_type(0,2));
171	region2.push_back(DataArrayView::RegionType::value_type(0,2));
172
173	DataAbstract* testData4=testData3->getSlice(region2);
174
175	for (int k=0;k<2;k++) {
176	for (int j=0;j<2;j++) {
177	for (int i=0;i<2;i++) {
178	assert(testData4->getPointDataView()(i,j,k)==pointData());
179	}
180	}
181	}
182
183	assert(testData4->getLength()==8);
184
185	cout << "\tVerify data point attributes." << endl;
186	dataView=testData4->getPointDataView();
187	assert(dataView.getRank()==3);
188	assert(dataView.noValues()==8);
189	assert(dataView.getShape()[0]==2);
190	assert(dataView.getShape()[1]==2);
191	assert(dataView.getShape()[2]==2);
192
193	cout << "\tTest slicing (part object)." << endl;
194
195	DataArrayView::RegionType region3;
196	region3.push_back(DataArrayView::RegionType::value_type(1,2));
197	region3.push_back(DataArrayView::RegionType::value_type(1,3));
198	region3.push_back(DataArrayView::RegionType::value_type(5,9));
199
200	DataAbstract* testData5=testData3->getSlice(region3);
201
202	for (int k=0;k<4;k++) {
203	for (int j=0;j<2;j++) {
204	for (int i=0;i<1;i++) {
205	assert(testData5->getPointDataView()(i,j,k)==pointData());
206	}
207	}
208	}
209
210	assert(testData5->getLength()==8);
211
212	cout << "\tVerify data point attributes." << endl;
213	dataView=testData5->getPointDataView();
214	assert(dataView.getRank()==3);
215	assert(dataView.noValues()==8);
216	assert(dataView.getShape()[0]==1);
217	assert(dataView.getShape()[1]==2);
218	assert(dataView.getShape()[2]==4);
219
220	cout << "\tTest slice setting (1)." << endl;
221
222	DataArrayView::RegionType region4;
223	region4.push_back(DataArrayView::RegionType::value_type(0,1));
224	region4.push_back(DataArrayView::RegionType::value_type(0,2));
225	region4.push_back(DataArrayView::RegionType::value_type(0,4));
226
227	DataAbstract* testData6=testData3->getSlice(region3);
228
229	testData5->setSlice(testData6,region4);
230
231	for (int k=0;k<4;k++) {
232	for (int j=0;j<2;j++) {
233	for (int i=0;i<1;i++) {
234	assert(testData5->getPointDataView()(i,j,k)==pointData());
235	}
236	}
237	}
238
239	assert(testData5->getLength()==8);
240
241	cout << "\tVerify data point attributes." << endl;
242	dataView=testData5->getPointDataView();
243	assert(dataView.getRank()==3);
244	assert(dataView.noValues()==8);
245	assert(dataView.getShape()[0]==1);
246	assert(dataView.getShape()[1]==2);
247	assert(dataView.getShape()[2]==4);
248
249	//cout << "\tTest slice setting (2)." << endl;
250
251	//DataArrayView::RegionType region5;
252	//region5.push_back(DataArrayView::RegionType::value_type(1,2));
253	//region5.push_back(DataArrayView::RegionType::value_type(1,3));
254	//region5.push_back(DataArrayView::RegionType::value_type(14,18));
255
256	//DataAbstract* testData7=testData3->getSlice(region);
257
258	//testData5->setSlice(testData7,region5);
259
260	//for (int k=0;k<4;k++) {
261	// for (int j=0;j<2;j++) {
262	// for (int i=0;i<1;i++) {
263	// assert(testData5->getPointDataView()(i,j,k)==pointData());
264	// }
265	// }
266	// }
267
268	//assert(testData5->getLength()==8);
269
270	//cout << "\tVerify data point attributes." << endl;
271	//dataView=testData5->getPointDataView();
272	//assert(dataView.getRank()==3);
273	//assert(dataView.noValues()==8);
274	//assert(dataView.getShape()[0]==1);
275	//assert(dataView.getShape()[1]==2);
276	//assert(dataView.getShape()[2]==4);
277
278	}
279
280	TestSuite* DataConstantTestCase::suite ()
281	{
282	//
283	// create the suite of tests to perform.
284	TestSuite *testSuite = new TestSuite ("DataConstantTestCase");
285
286	testSuite->addTest (new TestCaller< DataConstantTestCase>("testAll",&DataConstantTestCase::testAll));
287	return testSuite;
288	}
Name	Value
svn:eol-style	native
svn:keywords	Author Date Id Revision