Homec4science
Diffusion akantu (old-master) Browse

rAKA/test/test_facet_extraction_tetra1.cc
08a273bbb6beold-master

Last Change
Raw File

test_facet_extraction_tetra1.cc
View Options

/**
* @file test_facet_extraction.cc
* @author Guillaume ANCIAUX <anciaux@lsmscluster1.epfl.ch>
* @date Thu Aug 19 13:05:27 2010
*
* @brief test of internal facet extraction
*
* @section LICENSE
*
* <insert license here>
*
*/
/* -------------------------------------------------------------------------- */
/* -------------------------------------------------------------------------- */
#include "aka_common.hh"
#include "mesh.hh"
#include "mesh_io.hh"
#include "mesh_io_msh.hh"
#include "mesh_utils.hh"
#include "solid_mechanics_model.hh"
#include "material.hh"
/* -------------------------------------------------------------------------- */
#ifdef AKANTU_USE_IOHELPER
# include "io_helper.h"
#endif //AKANTU_USE_IOHELPER
int main(int argc, char *argv[])
{
akantu::UInt max_steps = 10000;
akantu::Real epot, ekin;
int dim = 3;
akantu::Mesh mesh(dim);
akantu::MeshIOMSH mesh_io;
mesh_io.read("cube.msh", mesh);
akantu::MeshUtils::buildFacets(mesh,1,1);
unsigned int nb_nodes = mesh.getNbNodes();
#ifdef AKANTU_USE_IOHELPER
DumperParaview dumper;
dumper.SetMode(TEXT);
dumper.SetPoints(mesh.getNodes().values, dim, nb_nodes, "test-facet-extraction");
dumper.SetConnectivity((int*)mesh.getConnectivity(akantu::_tetrahedra_1).values,
TETRA1, mesh.getNbElement(akantu::_tetrahedra_1), C_MODE);
dumper.SetPrefix("paraview/");
dumper.Init();
dumper.Dump();
DumperParaview dumper_facet;
dumper_facet.SetMode(TEXT);
dumper_facet.SetPoints(mesh.getNodes().values, dim, nb_nodes, "test-facet-extraction_boundary");
dumper_facet.SetConnectivity((int*)mesh.getConnectivity(akantu::_triangle_1).values,
TRIANGLE1, mesh.getNbElement(akantu::_triangle_1), C_MODE);
dumper_facet.SetPrefix("paraview/");
dumper_facet.Init();
dumper_facet.Dump();
dumper_facet.SetPoints(mesh.getNodes().values, dim, nb_nodes, "test-facet-extraction_internal");
dumper_facet.SetConnectivity((int*)mesh.getInternalFacetsMesh().getConnectivity(akantu::_triangle_1).values,
TRIANGLE1, mesh.getInternalFacetsMesh().getNbElement(akantu::_triangle_1), C_MODE);
dumper_facet.Init();
dumper_facet.Dump();
#endif //AKANTU_USE_IOHELPER
return EXIT_SUCCESS;
}
c4science · Help