Homec4science
Diffusion akantu (features/40-contact-using-the-internodes-method) Browse

rAKA/test/test_mesh/test_mesh_periodic.cc
e71e540aaa48features/40-contact…

Last Change
Raw File

test_mesh_periodic.cc
View Options

/**
* @file test_mesh_periodic.cc
*
* @author Nicolas Richart
*
* @date creation Sun Feb 11 2018
*
* @brief test makePeriodic
*
* @section LICENSE
*
* Copyright (©) 2010-2011 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 "mesh.hh"
/* -------------------------------------------------------------------------- */
#include "dumpable_inline_impl.hh"
#include "dumper_element_partition.hh"
#include "dumper_iohelper_paraview.hh"
/* -------------------------------------------------------------------------- */
using namespace akantu;
int main(int argc, char ** argv) {
initialize(argc, argv);
constexpr UInt dim = 3;
Mesh mesh(dim);
if (Communicator::getStaticCommunicator().whoAmI() == 0) {
mesh.read("cube_periodic.msh");
}
mesh.distribute();
auto * dumper = new DumperParaview("periodic", "./paraview");
mesh.registerExternalDumper(*dumper, "periodic", true);
mesh.addDumpMesh(mesh);
if (mesh.isDistributed())
mesh.addDumpFieldExternalToDumper(
"periodic", "node_type", const_cast<const Mesh &>(mesh).getNodesFlags());
mesh.dump();
mesh.makePeriodic(_x);
mesh.makePeriodic(_y);
mesh.makePeriodic(_z);
Array<Real> periodic(mesh.getNbNodes(), 1, 0.);
UInt prev_node = -1;
UInt value = 0;
const auto & periodic_ms = mesh.getPeriodicMasterSlaves();
for (auto & pair : periodic_ms) {
if (prev_node != pair.first) {
++value;
}
prev_node = pair.first;
periodic(pair.first) = value;
periodic(pair.second) = value;
auto it = periodic_ms.find(pair.second);
if (it != periodic_ms.end()) {
std::cout << pair.first << " " << pair.second << std::endl;
std::cout << it->first << " " << it->second << std::endl;
AKANTU_EXCEPTION("Shit meeeen!");
}
}
mesh.addDumpFieldExternalToDumper("periodic", "periodic", periodic);
mesh.dump();
}
c4science · Help