escript/test/DataAlgorithmAdapterTestCase.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/DataExpanded.h"
#include "escript/DataArrayView.h"
#include "escript/DataAlgorithm.h"
#include "DataAlgorithmAdapterTestCase.h"

#include <iostream>
#include <algorithm>
#include <math.h>
#include <limits>

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

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

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

void DataAlgorithmAdapterTestCase::testAll() {

  cout << endl;

  cout << "\tTesting FMax." << endl;

  FMax fmax;
  assert(std::abs(fmax(5,6)-6)<=REL_TOL*6);
  assert(std::abs(fmax(5,-6)-5)<=REL_TOL*5);
  assert(std::abs(fmax(0,0)-0)<=REL_TOL*0);
  assert(std::abs(fmax(15,-96)-15)<=REL_TOL*15);

  DataAlgorithmAdapter<FMax> sup(numeric_limits<double>::max()*-1);
  sup.resetResult();
  sup(-1);
  sup(-2);
  sup(-14);
  sup(3);
  assert(std::abs(sup.getResult()-3)<=REL_TOL*3);

  cout << "\tTesting AbsMax." << endl;

  AbsMax absmax;
  assert(std::abs(absmax(5,6)-6)<=REL_TOL*6);
  assert(std::abs(absmax(5,-6)-6)<=REL_TOL*6);
  assert(std::abs(absmax(0,0)-0)<=REL_TOL*6);
  assert(std::abs(absmax(15,-96)-96)<=REL_TOL*6);

  DataAlgorithmAdapter<AbsMax> Lsup(0);
  Lsup.resetResult();
  Lsup(-2);
  Lsup(2);
  Lsup(5);
  Lsup(-10);
  assert(std::abs(Lsup.getResult()-10)<=REL_TOL*10);

  cout << "\tTesting FMin." << endl;

  FMin fmin;
  assert(std::abs(fmin(5,6)-5)<=REL_TOL*5);
  assert(std::abs(fmin(5,-6)-(-6))<=REL_TOL*6);
  assert(std::abs(fmin(0,0)-0)<=REL_TOL*0);
  assert(std::abs(fmin(15,-96)-(-96))<=REL_TOL*96);

  DataAlgorithmAdapter<FMin> inf(numeric_limits<double>::max());
  inf.resetResult();
  inf(1);
  inf(12);
  inf(2);
  inf(99);
  assert(std::abs(inf.getResult()-1)<=REL_TOL*1);

  cout << "\tTesting Length." << endl;

  Length lngth;
  assert(std::abs(lngth(5,6)-std::sqrt(61.0))<=REL_TOL*std::sqrt(61.0));
  assert(std::abs(lngth(5,-6)-std::sqrt(61.0))<=REL_TOL*std::sqrt(61.0));
  assert(std::abs(lngth(0,0)-std::sqrt(0.0))<=REL_TOL*std::sqrt(61.0));
  assert(std::abs(lngth(15,-96)-std::sqrt(9441.0))<=REL_TOL*std::sqrt(61.0));

  DataAlgorithmAdapter<Length> length(0);
  length.resetResult();
  length(2);
  length(4);
  length(6);
  length(8);
  assert(std::abs(length.getResult()-std::sqrt(120.0))<=REL_TOL*std::sqrt(120.0));
  length.resetResult();
  length(1.5);
  length(2.5);
  length(3.5);
  length(4.5);
  assert(std::abs(length.getResult()-std::sqrt(41.0))<=REL_TOL*std::sqrt(41.0));

  cout << "\tTesting Trace." << endl;

  Trace trce;
  assert(std::abs(trce(5,6)-11)<=REL_TOL*11);
  assert(std::abs(trce(5,-6)-(-1))<=REL_TOL*1);
  assert(std::abs(trce(0,0)-0)<=REL_TOL*0);
  assert(std::abs(trce(15,-96)-(-81))<=REL_TOL*81);

  DataAlgorithmAdapter<Trace> trace(0);
  trace.resetResult();
  trace(1);
  trace(2);
  trace(3);
  trace(4);
  trace(5);
  assert(std::abs(trace.getResult()-15)<=REL_TOL*15);
  trace.resetResult();
  trace(1.5);
  trace(2.5);
  trace(3.5);
  trace(4.5);
  trace(5.5);
  assert(std::abs(trace.getResult()-17.5)<=REL_TOL*17.5);

}

void DataAlgorithmAdapterTestCase::testAlgorithm() {

  cout << endl;

  {

    cout << "\tTest algorithm on Data objects with a single rank 2 data-point." << endl;

    // define the shape for the DataArrayView
    DataArrayView::ShapeType shape;
    shape.push_back(2);
    shape.push_back(3);

    // allocate the data for the DataArrayView
    DataArrayView::ValueType dataArray(DataArrayView::noValues(shape),0);

    // construct DataArrayView
    DataArrayView dataView(dataArray,shape);

    // assign values to the data point
    for (int i=0;i<shape[0];i++) {
      for (int j=0;j<shape[1];j++) {
        dataView(i,j)=dataView.index(i,j);
      }
    }

    // create a few Data objects from the created DataArrayView
    DataExpanded dataExp(dataView,FunctionSpace());
    DataConstant dataCon(dataView,FunctionSpace());
    DataTagged   dataTag(dataCon);

    // test algorithm on DataExpanded
    FMin fmin_func;
    assert(std::abs(escript::algorithm(dataExp,fmin_func,numeric_limits<double>::max())-0)<=REL_TOL*0);
    FMax fmax_func;
    assert(std::abs(escript::algorithm(dataExp,fmax_func,numeric_limits<double>::max()*-1)-5)<=REL_TOL*5);

    // test algorithm on DataTagged
    assert(std::abs(escript::algorithm(dataTag,fmin_func,numeric_limits<double>::max())-0)<=REL_TOL*0);
    assert(std::abs(escript::algorithm(dataTag,fmax_func,numeric_limits<double>::max()*-1)-5)<=REL_TOL*5);

    // test algorithm on DataConstant
    assert(std::abs(escript::algorithm(dataCon,fmin_func,numeric_limits<double>::max())-0)<=REL_TOL*0);
    assert(std::abs(escript::algorithm(dataCon,fmax_func,numeric_limits<double>::max()*-1)-5)<=REL_TOL*5);

  }

}

void DataAlgorithmAdapterTestCase::testDpAlgorithm() {

  cout << endl;

  {

    cout << "\tTest dp_algorithm on Data objects with a single rank 2 data-point." << endl;

    // define the shapes for the DataArrayViews
    DataArrayView::ShapeType shape;
    shape.push_back(2);
    shape.push_back(3);
    DataArrayView::ShapeType shape2;

    // allocate the data for the DataArrayViews
    DataArrayView::ValueType dataArray(DataArrayView::noValues(shape),0);
    DataArrayView::ValueType dataArray2(DataArrayView::noValues(shape2),0);

    // construct DataArrayViews
    DataArrayView dataView(dataArray,shape);
    DataArrayView dataView2(dataArray2,shape2);

    // assign values to the data point
    for (int i=0;i<shape[0];i++) {
      for (int j=0;j<shape[1];j++) {
        dataView(i,j)=dataView.index(i,j);
      }
    }

    // create a few Data objects from the created DataArrayViews
    DataExpanded dataExp(dataView,FunctionSpace());
    DataConstant dataCon(dataView,FunctionSpace());
    DataTagged   dataTag(dataCon);

    // and create Data objects to receive the results of the dp_algorithm calls
    DataExpanded dataExp2(dataView2,FunctionSpace());
    DataConstant dataCon2(dataView2,FunctionSpace());
    DataTagged   dataTag2(dataCon2);

    // test dp_algorithm on DataExpanded
    FMin fmin_func;
    escript::dp_algorithm(dataExp,dataExp2,fmin_func,numeric_limits<double>::max());
    assert(std::abs(dataExp2.getDataPoint(0,0)()-0)<=REL_TOL*0);
    FMax fmax_func;
    escript::dp_algorithm(dataExp,dataExp2,fmax_func,numeric_limits<double>::max()*-1);
    assert(std::abs(dataExp2.getDataPoint(0,0)()-5)<=REL_TOL*5);

    // test dp_algorithm on DataTagged
    escript::dp_algorithm(dataTag,dataTag2,fmin_func,numeric_limits<double>::max());
    assert(std::abs(dataTag2.getDataPoint(0,0)()-0)<=REL_TOL*0);
    escript::dp_algorithm(dataTag,dataTag2,fmax_func,numeric_limits<double>::max()*-1);
    assert(std::abs(dataTag2.getDataPoint(0,0)()-5)<=REL_TOL*5);

    // test dp_algorithm on DataConstant
    escript::dp_algorithm(dataCon,dataCon2,fmin_func,numeric_limits<double>::max());
    assert(std::abs(dataCon2.getDataPoint(0,0)()-0)<=REL_TOL*0);
    escript::dp_algorithm(dataCon,dataCon2,fmax_func,numeric_limits<double>::max()*-1);
    assert(std::abs(dataCon2.getDataPoint(0,0)()-5)<=REL_TOL*5);

  }

}

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

  testSuite->addTest (new TestCaller< DataAlgorithmAdapterTestCase>("testAll",&DataAlgorithmAdapterTestCase::testAll));
  testSuite->addTest (new TestCaller<DataAlgorithmAdapterTestCase>("testAlgorithm",&DataAlgorithmAdapterTestCase::testAlgorithm));
  testSuite->addTest (new TestCaller<DataAlgorithmAdapterTestCase>("testDpAlgorithm",&DataAlgorithmAdapterTestCase::testDpAlgorithm));
  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
16	#include "escript/DataExpanded.h"
17	#include "escript/DataArrayView.h"
18	#include "escript/DataAlgorithm.h"
19	#include "DataAlgorithmAdapterTestCase.h"
20
21	#include <iostream>
22	#include <algorithm>
23	#include <math.h>
24	#include <limits>
25
26	using namespace CppUnitTest;
27	using namespace std;
28	using namespace escript;
29
30	void DataAlgorithmAdapterTestCase::setUp() {
31	//
32	// This is called before each test is run
33
34	}
35
36	void DataAlgorithmAdapterTestCase::tearDown() {
37	//
38	// This is called after each test has been run
39
40	}
41
42	void DataAlgorithmAdapterTestCase::testAll() {
43
44	cout << endl;
45
46	cout << "\tTesting FMax." << endl;
47
48	FMax fmax;
49	assert(std::abs(fmax(5,6)-6)<=REL_TOL*6);
50	assert(std::abs(fmax(5,-6)-5)<=REL_TOL*5);
51	assert(std::abs(fmax(0,0)-0)<=REL_TOL*0);
52	assert(std::abs(fmax(15,-96)-15)<=REL_TOL*15);
53
54	DataAlgorithmAdapter<FMax> sup(numeric_limits<double>::max()*-1);
55	sup.resetResult();
56	sup(-1);
57	sup(-2);
58	sup(-14);
59	sup(3);
60	assert(std::abs(sup.getResult()-3)<=REL_TOL*3);
61
62	cout << "\tTesting AbsMax." << endl;
63
64	AbsMax absmax;
65	assert(std::abs(absmax(5,6)-6)<=REL_TOL*6);
66	assert(std::abs(absmax(5,-6)-6)<=REL_TOL*6);
67	assert(std::abs(absmax(0,0)-0)<=REL_TOL*6);
68	assert(std::abs(absmax(15,-96)-96)<=REL_TOL*6);
69
70	DataAlgorithmAdapter<AbsMax> Lsup(0);
71	Lsup.resetResult();
72	Lsup(-2);
73	Lsup(2);
74	Lsup(5);
75	Lsup(-10);
76	assert(std::abs(Lsup.getResult()-10)<=REL_TOL*10);
77
78	cout << "\tTesting FMin." << endl;
79
80	FMin fmin;
81	assert(std::abs(fmin(5,6)-5)<=REL_TOL*5);
82	assert(std::abs(fmin(5,-6)-(-6))<=REL_TOL*6);
83	assert(std::abs(fmin(0,0)-0)<=REL_TOL*0);
84	assert(std::abs(fmin(15,-96)-(-96))<=REL_TOL*96);
85
86	DataAlgorithmAdapter<FMin> inf(numeric_limits<double>::max());
87	inf.resetResult();
88	inf(1);
89	inf(12);
90	inf(2);
91	inf(99);
92	assert(std::abs(inf.getResult()-1)<=REL_TOL*1);
93
94	cout << "\tTesting Length." << endl;
95
96	Length lngth;
97	assert(std::abs(lngth(5,6)-std::sqrt(61.0))<=REL_TOL*std::sqrt(61.0));
98	assert(std::abs(lngth(5,-6)-std::sqrt(61.0))<=REL_TOL*std::sqrt(61.0));
99	assert(std::abs(lngth(0,0)-std::sqrt(0.0))<=REL_TOL*std::sqrt(61.0));
100	assert(std::abs(lngth(15,-96)-std::sqrt(9441.0))<=REL_TOL*std::sqrt(61.0));
101
102	DataAlgorithmAdapter<Length> length(0);
103	length.resetResult();
104	length(2);
105	length(4);
106	length(6);
107	length(8);
108	assert(std::abs(length.getResult()-std::sqrt(120.0))<=REL_TOL*std::sqrt(120.0));
109	length.resetResult();
110	length(1.5);
111	length(2.5);
112	length(3.5);
113	length(4.5);
114	assert(std::abs(length.getResult()-std::sqrt(41.0))<=REL_TOL*std::sqrt(41.0));
115
116	cout << "\tTesting Trace." << endl;
117
118	Trace trce;
119	assert(std::abs(trce(5,6)-11)<=REL_TOL*11);
120	assert(std::abs(trce(5,-6)-(-1))<=REL_TOL*1);
121	assert(std::abs(trce(0,0)-0)<=REL_TOL*0);
122	assert(std::abs(trce(15,-96)-(-81))<=REL_TOL*81);
123
124	DataAlgorithmAdapter<Trace> trace(0);
125	trace.resetResult();
126	trace(1);
127	trace(2);
128	trace(3);
129	trace(4);
130	trace(5);
131	assert(std::abs(trace.getResult()-15)<=REL_TOL*15);
132	trace.resetResult();
133	trace(1.5);
134	trace(2.5);
135	trace(3.5);
136	trace(4.5);
137	trace(5.5);
138	assert(std::abs(trace.getResult()-17.5)<=REL_TOL*17.5);
139
140	}
141
142	void DataAlgorithmAdapterTestCase::testAlgorithm() {
143
144	cout << endl;
145
146	{
147
148	cout << "\tTest algorithm on Data objects with a single rank 2 data-point." << endl;
149
150	// define the shape for the DataArrayView
151	DataArrayView::ShapeType shape;
152	shape.push_back(2);
153	shape.push_back(3);
154
155	// allocate the data for the DataArrayView
156	DataArrayView::ValueType dataArray(DataArrayView::noValues(shape),0);
157
158	// construct DataArrayView
159	DataArrayView dataView(dataArray,shape);
160
161	// assign values to the data point
162	for (int i=0;i<shape[0];i++) {
163	for (int j=0;j<shape[1];j++) {
164	dataView(i,j)=dataView.index(i,j);
165	}
166	}
167
168	// create a few Data objects from the created DataArrayView
169	DataExpanded dataExp(dataView,FunctionSpace());
170	DataConstant dataCon(dataView,FunctionSpace());
171	DataTagged dataTag(dataCon);
172
173	// test algorithm on DataExpanded
174	FMin fmin_func;
175	assert(std::abs(escript::algorithm(dataExp,fmin_func,numeric_limits<double>::max())-0)<=REL_TOL*0);
176	FMax fmax_func;
177	assert(std::abs(escript::algorithm(dataExp,fmax_func,numeric_limits<double>::max()-1)-5)<=REL_TOL5);
178
179	// test algorithm on DataTagged
180	assert(std::abs(escript::algorithm(dataTag,fmin_func,numeric_limits<double>::max())-0)<=REL_TOL*0);
181	assert(std::abs(escript::algorithm(dataTag,fmax_func,numeric_limits<double>::max()-1)-5)<=REL_TOL5);
182
183	// test algorithm on DataConstant
184	assert(std::abs(escript::algorithm(dataCon,fmin_func,numeric_limits<double>::max())-0)<=REL_TOL*0);
185	assert(std::abs(escript::algorithm(dataCon,fmax_func,numeric_limits<double>::max()-1)-5)<=REL_TOL5);
186
187	}
188
189	}
190
191	void DataAlgorithmAdapterTestCase::testDpAlgorithm() {
192
193	cout << endl;
194
195	{
196
197	cout << "\tTest dp_algorithm on Data objects with a single rank 2 data-point." << endl;
198
199	// define the shapes for the DataArrayViews
200	DataArrayView::ShapeType shape;
201	shape.push_back(2);
202	shape.push_back(3);
203	DataArrayView::ShapeType shape2;
204
205	// allocate the data for the DataArrayViews
206	DataArrayView::ValueType dataArray(DataArrayView::noValues(shape),0);
207	DataArrayView::ValueType dataArray2(DataArrayView::noValues(shape2),0);
208
209	// construct DataArrayViews
210	DataArrayView dataView(dataArray,shape);
211	DataArrayView dataView2(dataArray2,shape2);
212
213	// assign values to the data point
214	for (int i=0;i<shape[0];i++) {
215	for (int j=0;j<shape[1];j++) {
216	dataView(i,j)=dataView.index(i,j);
217	}
218	}
219
220	// create a few Data objects from the created DataArrayViews
221	DataExpanded dataExp(dataView,FunctionSpace());
222	DataConstant dataCon(dataView,FunctionSpace());
223	DataTagged dataTag(dataCon);
224
225	// and create Data objects to receive the results of the dp_algorithm calls
226	DataExpanded dataExp2(dataView2,FunctionSpace());
227	DataConstant dataCon2(dataView2,FunctionSpace());
228	DataTagged dataTag2(dataCon2);
229
230	// test dp_algorithm on DataExpanded
231	FMin fmin_func;
232	escript::dp_algorithm(dataExp,dataExp2,fmin_func,numeric_limits<double>::max());
233	assert(std::abs(dataExp2.getDataPoint(0,0)()-0)<=REL_TOL*0);
234	FMax fmax_func;
235	escript::dp_algorithm(dataExp,dataExp2,fmax_func,numeric_limits<double>::max()*-1);
236	assert(std::abs(dataExp2.getDataPoint(0,0)()-5)<=REL_TOL*5);
237
238	// test dp_algorithm on DataTagged
239	escript::dp_algorithm(dataTag,dataTag2,fmin_func,numeric_limits<double>::max());
240	assert(std::abs(dataTag2.getDataPoint(0,0)()-0)<=REL_TOL*0);
241	escript::dp_algorithm(dataTag,dataTag2,fmax_func,numeric_limits<double>::max()*-1);
242	assert(std::abs(dataTag2.getDataPoint(0,0)()-5)<=REL_TOL*5);
243
244	// test dp_algorithm on DataConstant
245	escript::dp_algorithm(dataCon,dataCon2,fmin_func,numeric_limits<double>::max());
246	assert(std::abs(dataCon2.getDataPoint(0,0)()-0)<=REL_TOL*0);
247	escript::dp_algorithm(dataCon,dataCon2,fmax_func,numeric_limits<double>::max()*-1);
248	assert(std::abs(dataCon2.getDataPoint(0,0)()-5)<=REL_TOL*5);
249
250	}
251
252	}
253
254	TestSuite* DataAlgorithmAdapterTestCase::suite ()
255	{
256	//
257	// create the suite of tests to perform.
258	TestSuite *testSuite = new TestSuite ("DataAlgorithmAdapterTestCase");
259
260	testSuite->addTest (new TestCaller< DataAlgorithmAdapterTestCase>("testAll",&DataAlgorithmAdapterTestCase::testAll));
261	testSuite->addTest (new TestCaller<DataAlgorithmAdapterTestCase>("testAlgorithm",&DataAlgorithmAdapterTestCase::testAlgorithm));
262	testSuite->addTest (new TestCaller<DataAlgorithmAdapterTestCase>("testDpAlgorithm",&DataAlgorithmAdapterTestCase::testDpAlgorithm));
263	return testSuite;
264	}
Name	Value
svn:eol-style	native
svn:keywords	Author Date Id Revision