tests_vorel.cpp

Go to the documentation of this file.

#include "problem.h"
#include "gelib.h"
#include "homogenization.h"
#include "mesh.h"

using namespace mumech;


int cmp_fields (int n, const double *f1, const double *f2, double eps)
{
  int nd = 0;
  
  while (n--)
    if (fabs( (f1[n] - f2[n]) / f1[n]) > eps)
      nd++;
  
  return nd;
}

void vzor (void)
{
  // char vtkGeomFile[]    = "vzor/X3D_3incl.geom.vtk";
  // char vtkAmndFile[]    = "vzor/X3D_3incl.amnd.vtk";
  // char vtkMeshFile[]    = "vzor/X3D_3incl.mesh";
  // 
  // // regular mesh
  // double p1[3] = {-3.0, -3.0, -3.0 };
  // double p2[3] = { 3.0,  3.0,  3.0 };
  // long n[3];
  // n[0] = (p2[0] - p1[0]) / 0.2;
  // n[1] = (p2[1] - p1[1]) / 0.2;
  // n[2] = (p2[2] - p1[2]) / 0.2;
  // 
  // 
  // Problem p; // mumech problem
  // p.read_inclusions_plus_initialize_and_print (vtkGeomFile, vtkAmndFile, DT_ANALITICAL);
  // 
  // 
  // Mesh *mesh = p.give_new_mesh ();
  // mesh->generate_regular_mesh (p1, p2, n);
  // mesh->compute_node_fields ('p', true, true, true, 0, 1, PFCM_FULL);
  // mesh->print_geometry_file_vtk (vtkMeshFile, 0, 1);
  // 
  // Homogenization *h = p.give_new_homogenization (HT_Dilute);
  // h->set_mesh (mesh);
  // h->compute_fields_HT_reg_grid (0, 0, PFCM_FULL);
  
  // // kontrola
  // double stress[6] = { 1.0, 1.0, 1.0, 1.0, 1.0, 1.0 };
  // 
  // const double **field = h->give_stress();
  // if (cmp_fields (6, stress, field[0], 1e-7))  _warningg("TestVerified_triple_incl_3D_homog_regular results differs");
}


void tests_vorel(void) {
  
  vzor ();
  
} 

/*end of file*/