test/SystemMatrixAdapter/SystemMatrixAdapterTestCase.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 "finley/CPPAdapter/MeshAdapter.h"
#include "finley/CPPAdapter/MeshAdapterFactory.h"
#include "escript/Data/AbstractContinuousDomain.h"
#include "finley/CPPAdapter/SystemMatrixAdapter.h"
#include "finley/CPPAdapter/FinleyAdapterException.h"
#include "finley/CPPAdapter/FinleyError.h"

#include "SystemMatrixAdapterTestCase.h"

using namespace std;

using namespace CppUnitTest;

using namespace escript;
using namespace finley;

static Finley_Mesh *mesh;
static Finley_SystemMatrix *system_matrix;
static Finley_SystemMatrixPattern *matrix_pattern;

static Finley_SystemMatrixType type;

static int symmetric;

static int row_blocksize;
static int column_blocksize;

static int reduce_row_order;
static int reduce_col_order;

static FunctionSpace row_functionspace;
static FunctionSpace colum_functionspace;

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

  mesh = Finley_Mesh_alloc("foo", 2, 1);

  type = CSR;

  symmetric = 0;

  row_blocksize = 10;
  column_blocksize = 10;

  reduce_row_order = 0;
  reduce_col_order = 0;

  matrix_pattern = Finley_getPattern(mesh, reduce_row_order, reduce_col_order);
  system_matrix = Finley_SystemMatrix_alloc(type,matrix_pattern,row_blocksize, column_blocksize);

}

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

  // this fails on gcc ...
  //Finley_SystemMatrix_dealloc(system_matrix);

}

void SystemMatrixAdapterTestCase::testAll() {
  //
  // The test code may be entered here
  // There is nothing special about the function name, it may be renamed to
  // something more suitable. 
  // As many test methods as desired may be added.

  cout << endl;

  try {
    cout << "\tTest illegal default construction." << endl;
    SystemMatrixAdapter system_matrix_adapter;
    assert(false);
  }
  catch (FinleyAdapterException& e) {
    cout << "\t" << e.toString() << endl;
    assert(true);
  }

  cout << "\tTest constructor." << endl;

  SystemMatrixAdapter system_matrix_adapter(system_matrix, row_blocksize, row_functionspace, column_blocksize, colum_functionspace);
  Finley_SystemMatrix* adapter_system_matrix = system_matrix_adapter.getFinley_SystemMatrix();
  assert(adapter_system_matrix == system_matrix);

  cout << "\tExercise nullifyRowsAndCols." << endl;

  Data row_q;
  Data col_q;
  double mdv = 1.0;

  system_matrix_adapter.nullifyRowsAndCols(row_q, col_q, mdv);

}

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

  testSuite->addTest (new TestCaller< SystemMatrixAdapterTestCase>("testAll",&SystemMatrixAdapterTestCase::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 "finley/CPPAdapter/MeshAdapter.h"
16	#include "finley/CPPAdapter/MeshAdapterFactory.h"
17	#include "escript/Data/AbstractContinuousDomain.h"
18	#include "finley/CPPAdapter/SystemMatrixAdapter.h"
19	#include "finley/CPPAdapter/FinleyAdapterException.h"
20	#include "finley/CPPAdapter/FinleyError.h"
21
22	#include "SystemMatrixAdapterTestCase.h"
23
24	using namespace std;
25
26	using namespace CppUnitTest;
27
28	using namespace escript;
29	using namespace finley;
30
31	static Finley_Mesh *mesh;
32	static Finley_SystemMatrix *system_matrix;
33	static Finley_SystemMatrixPattern *matrix_pattern;
34
35	static Finley_SystemMatrixType type;
36
37	static int symmetric;
38
39	static int row_blocksize;
40	static int column_blocksize;
41
42	static int reduce_row_order;
43	static int reduce_col_order;
44
45	static FunctionSpace row_functionspace;
46	static FunctionSpace colum_functionspace;
47
48	void SystemMatrixAdapterTestCase::setUp() {
49	//
50	// This is called before each test is run
51
52	mesh = Finley_Mesh_alloc("foo", 2, 1);
53
54	type = CSR;
55
56	symmetric = 0;
57
58	row_blocksize = 10;
59	column_blocksize = 10;
60
61	reduce_row_order = 0;
62	reduce_col_order = 0;
63
64	matrix_pattern = Finley_getPattern(mesh, reduce_row_order, reduce_col_order);
65	system_matrix = Finley_SystemMatrix_alloc(type,matrix_pattern,row_blocksize, column_blocksize);
66
67	}
68
69	void SystemMatrixAdapterTestCase::tearDown() {
70	//
71	// This is called after each test has been run
72
73	// this fails on gcc ...
74	//Finley_SystemMatrix_dealloc(system_matrix);
75
76	}
77
78	void SystemMatrixAdapterTestCase::testAll() {
79	//
80	// The test code may be entered here
81	// There is nothing special about the function name, it may be renamed to
82	// something more suitable.
83	// As many test methods as desired may be added.
84
85	cout << endl;
86
87	try {
88	cout << "\tTest illegal default construction." << endl;
89	SystemMatrixAdapter system_matrix_adapter;
90	assert(false);
91	}
92	catch (FinleyAdapterException& e) {
93	cout << "\t" << e.toString() << endl;
94	assert(true);
95	}
96
97	cout << "\tTest constructor." << endl;
98
99	SystemMatrixAdapter system_matrix_adapter(system_matrix, row_blocksize, row_functionspace, column_blocksize, colum_functionspace);
100	Finley_SystemMatrix* adapter_system_matrix = system_matrix_adapter.getFinley_SystemMatrix();
101	assert(adapter_system_matrix == system_matrix);
102
103	cout << "\tExercise nullifyRowsAndCols." << endl;
104
105	Data row_q;
106	Data col_q;
107	double mdv = 1.0;
108
109	system_matrix_adapter.nullifyRowsAndCols(row_q, col_q, mdv);
110
111	}
112
113	TestSuite* SystemMatrixAdapterTestCase::suite ()
114	{
115	//
116	// create the suite of tests to perform.
117	TestSuite *testSuite = new TestSuite ("SystemMatrixAdapterTestCase");
118
119	testSuite->addTest (new TestCaller< SystemMatrixAdapterTestCase>("testAll",&SystemMatrixAdapterTestCase::testAll));
120	return testSuite;
121	}
Name	Value
svn:eol-style	native
svn:keywords	Author Date Id Revision