Homec4science
Diffusion akantu (bugfixes/46-bug-in-strain-energy-split-and-missing-constraint-in-phasefield) Browse

rAKA/test/test_surface_extraction/test_surface_extraction_tetrahedron_4.cc
72afa32339fabugfixes/46-bug-in…

Last Change
Raw File

test_surface_extraction_tetrahedron_4.cc
View Options

/**
* @file test_surface_extraction_tetrahedron_4.cc
*
* @author Leonardo Snozzi <leonardo.snozzi@epfl.ch>
* @author Marco Vocialta <marco.vocialta@epfl.ch>
*
* @date creation: Tue Oct 26 2010
* @date last modification: Thu Mar 27 2014
*
* @brief 3d surface extraction tests
*
* @section LICENSE
*
* Copyright (©) 2010-2012, 2014 EPFL (Ecole Polytechnique Fédérale de Lausanne)
* Laboratory (LSMS - Laboratoire de Simulation en Mécanique des Solides)
*
* Akantu is free software: you can redistribute it and/or modify it under the
* terms of the GNU Lesser General Public License as published by the Free
* Software Foundation, either version 3 of the License, or (at your option) any
* later version.
*
* Akantu is distributed in the hope that it will be useful, but WITHOUT ANY
* WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS FOR
* A PARTICULAR PURPOSE. See the GNU Lesser General Public License for more
* details.
*
* You should have received a copy of the GNU Lesser General Public License
* along with Akantu. If not, see <http://www.gnu.org/licenses/>.
*
*/
/* -------------------------------------------------------------------------- */
#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"
#include <cstdio>
#ifdef AKANTU_USE_IOHELPER
# include "io_helper.hh"
#endif //AKANTU_USE_IOHELPER
using namespace akantu;
int main(int argc, char *argv[])
{
int dim = 3;
akantu::initialize(argc, argv);
Mesh mesh(dim);
MeshIOMSH mesh_io;
mesh_io.read("cubes.msh", mesh);
MeshUtils::buildFacets(mesh);
mesh.createBoundaryGroupFromGeometry();
// MeshUtils::buildSurfaceID(mesh);
unsigned int nb_nodes = mesh.getNbNodes();
#ifdef AKANTU_USE_IOHELPER
iohelper::DumperParaview dumper;
dumper.setMode(iohelper::TEXT);
dumper.setPoints(mesh.getNodes().storage(), dim, nb_nodes, "test-surface-extraction");
dumper.setConnectivity((int*)mesh.getConnectivity(_tetrahedron_4).storage(),
iohelper::TETRA1, mesh.getNbElement(_tetrahedron_4), iohelper::C_MODE);
dumper.setPrefix("paraview/");
dumper.init();
dumper.dump();
iohelper::DumperParaview dumper_surface;
dumper_surface.setMode(iohelper::TEXT);
dumper_surface.setPoints(mesh.getNodes().storage(), dim, nb_nodes, "test-surface-extraction_boundary");
dumper_surface.setConnectivity((int *)mesh.getConnectivity(_triangle_3).storage(),
iohelper::TRIANGLE1, mesh.getNbElement(_triangle_3), iohelper::C_MODE);
UInt * surf_id = new UInt[mesh.getNbElement(_triangle_3)];
for(Mesh::const_element_group_iterator it(mesh.element_group_begin());
it != mesh.element_group_end(); ++it) {
const Array<UInt> & element_ids = it->second->getElements(_triangle_3);
for(UInt i(0); i << element_ids.getSize(); ++i) {
UInt elem_idx = element_ids(i);
surf_id[elem_idx] = atoi((it->first.c_str()));
}
}
// for (UInt i = 0; i < mesh.getSurfaceID(_triangle_3).getSize(); ++i)
// surf_id[i] = (double)mesh.getSurfaceID(_triangle_3).storage()[i];
dumper_surface.addElemDataField("surface_id", surf_id, 1, mesh.getNbElement(_triangle_3));
dumper_surface.setPrefix("paraview/");
dumper_surface.init();
dumper_surface.dump();
delete [] surf_id;
#endif //AKANTU_USE_IOHELPER
finalize();
return EXIT_SUCCESS;
}
c4science · Help