2DNodeMesh.cpp File Reference

A simple check of the TwoD_Node_Mesh container that stores nodal data over a given mesh. More...

#include <Generic_bundle.h>
#include <TwoD_Node_Mesh.h>

Go to the source code of this file.

Functions
int	main ()

Detailed Description

A simple check of the TwoD_Node_Mesh container that stores nodal data over a given mesh.

This simply does a remeshing check.

Definition in file 2DNodeMesh.cpp.

Function Documentation

main()

int main

(

)

Definition at line 14 of file 2DNodeMesh.cpp.

{
  cout.precision( 10 );
  cout << "\n";
  cout << "=== TwoD_Node_Mesh: very basic read/write test ======\n";
  cout << "\n";
 
  std::size_t nx( 11 );
  std::size_t ny( 21 );
  DenseVector<double> x = Utility::uniform_node_vector( -1.0, 1.0, nx );
  DenseVector<double> y = Utility::uniform_node_vector(  0.0, 1.0, ny );
  std::size_t nx2( 31 );
  std::size_t ny2( 41 );
  DenseVector<double> x2 = Utility::uniform_node_vector( -1.0, 1.0, nx2 );
  DenseVector<double> y2 = Utility::uniform_node_vector(  0.0, 1.0, ny2 );
 
 
  TwoD_Node_Mesh<double> mesh( x, y, 1 );
  // write
  for ( std::size_t i = 0; i < nx; ++i )
  {
    for ( std::size_t j = 0; j < ny; ++j )
    {
      mesh( i, j, 0 ) = cos( M_PI * x[ i ] ) * sin( M_PI * y[ j ] );
    }
  }
 
  //mesh.dump_gnu( "./DATA/mesh1.dat" );
  mesh.remesh1( x2, y2 );
  //mesh.dump_gnu( "./DATA/mesh2.dat" );
  mesh.remesh1( x, y );
  //mesh.dump_gnu( "./DATA/mesh3.dat" );
 
  // read
  for ( std::size_t j = 0; j < ny; ++j )
  {
    for ( std::size_t i = 0; i < nx; ++i )
    {
      mesh( i, j, 0 ) -= cos( M_PI * x[ i ] ) * sin( M_PI * y[ j ] );
    }
  }
 
  // check diff
  const double diff = mesh.get_var_as_matrix( 0 ).inf_norm();
  if ( diff > 1.e-14 )
  {
    cout << "\033[1;31;48m  * FAILED \033[0m\n";
    cout.precision( 10 );
    cout << "Difference following remeshing is " << diff << "\n";
    return 1;
  }
  else
  {
    cout << "\033[1;32;48m  * PASSED \033[0m\n";
    return 0;
  }
 
}

References CppNoddy::TwoD_Node_Mesh< _Type >::get_var_as_matrix(), CppNoddy::TwoD_Node_Mesh< _Type >::remesh1(), and CppNoddy::Utility::uniform_node_vector().