Homec4science
Diffusion CADD_mesher (master) Browse

rCADDMESH/test/block_test.cc
2380c45f99c4master

Last Change
Raw File

block_test.cc
View Options

#include "cadd_mesher.hh"
#include <math.h>
int main(int argc, char *argv[])
{
Real tol = 1e-4;
Real point1 [3] = {0, 1, 0};
Real point2 [3] = {0, -1, 0};
Real point3 [3] = {0, 3.2, -1};
Real min_max [6] = {0,1,-.2,3.2, -.9, 1};//{0,1./sqrt(3),0,1/sqrt(3),0,1/sqrt(3)};
Block<2> block( min_max, "affe");
std::cout << block << std::endl;
int ret_val = 0;
if (block.is_inside(point1, tol) > 0) {
ret_val += 1;}
if (block.is_inside(point2, tol) <= 0) {
ret_val += 10;}
if (block.is_inside(point3, tol) > 0) {
ret_val += 100;}
const Uint num_dirs = 4;
Real directions [num_dirs][3] = {{1,0,0}, {0,1,0}, {0,0,1}, {1,1,1}};
Block<3> block3( min_max, "3 dimensinal");
std::cout << block3 << std::endl;
for (Uint i = 0; i < num_dirs; ++i) {
std::cout << "len in angs = \t";
std::cout.width (13);
std::cout << block3.size_in_direction(directions[i], 1)
<< ", direction = (" << directions[i][0] << ", "
<< directions[i][1] << ", " << directions[i][2] << ")"
<< std::endl;
std::cout << "len in latt = \t";
std::cout.width (13);
std::cout << block3.size_in_direction(directions[i], 4.04)
<< ", direction = (" << directions[i][0] << ", "
<< directions[i][1] << ", " << directions[i][2] << ")"
<< std::endl;
std::cout << "min in angs = \t";
std::cout.width (13);
std::cout << block3.min_in_direction(directions[i], 1)
<< ", direction = (" << directions[i][0] << ", "
<< directions[i][1] << ", " << directions[i][2] << ")"
<< std::endl;
std::cout << "min in latt = \t";
std::cout.width (13);
std::cout << block3.min_in_direction(directions[i], 4.04)
<< ", direction = (" << directions[i][0] << ", "
<< directions[i][1] << ", " << directions[i][2] << ")"
<< std::endl;
std::cout << "max in angs = \t";
std::cout.width (13);
std::cout << block3.max_in_direction(directions[i], 1)
<< ", direction = (" << directions[i][0] << ", "
<< directions[i][1] << ", " << directions[i][2] << ")"
<< std::endl;
std::cout << "max in latt = \t";
std::cout.width (13);
std::cout << block3.max_in_direction(directions[i], 4.04)
<< ", direction = (" << directions[i][0] << ", "
<< directions[i][1] << ", " << directions[i][2] << ")"
<< std::endl;
}
if (block3.is_inside(point1, tol) > 0) {
ret_val +=1000;}
if (block3.is_inside(point2, tol) <= 0) {
ret_val += 10000;}
if (block3.is_inside(point3, tol) <= 0) {
ret_val += 100000;}
if (ret_val == 0) {
std::cout << "Test passed." << std::endl;
} else {
std::cout << "Test failed!!! code is "<<ret_val<< std::endl;
}
return ret_val;
}
c4science · Help