22 #include <LbmLib/include/solver/BioSolver/tutorial_01_BioSolverMembraneTension.hpp>
23 #include <LbmLib/include/geometry/GeometryHandler.hpp>
24 #include <LbmLib/include/solver/ForceSolver.hpp>
25 #include <LbmLib/include/GlobalSimulationParameters.hpp>
34 const unsigned int FREQUENCY = 10;
35 const double FORCE = 0.01;
39 tutorial_01_BioSolverMembraneTension::tutorial_01_BioSolverMembraneTension() : BioBaseSolver()
46 std::stringstream forcedescriptor;
49 if (Parameters.getCurrentIteration()%FREQUENCY != 0) {
57 #pragma omp parallel for schedule(static) private(forcedescriptor)
59 j<tempconnectionmap.size();
61 forcedescriptor.str( std::string() );
62 forcedescriptor.clear();
63 forcedescriptor <<
"6\t"
64 << tempconnectionmap[j]->getGeometryNodes().first->getId() <<
"\t"
65 << tempconnectionmap[j]->getGeometryNodes().second->getId() <<
"\t"
67 forcesolver.
addForce(&forcedescriptor);
69 forcedescriptor.str( std::string() );
70 forcedescriptor.clear();
71 forcedescriptor <<
"6\t"
72 << tempconnectionmap[j]->getGeometryNodes().second->getId() <<
"\t"
73 << tempconnectionmap[j]->getGeometryNodes().first->getId() <<
"\t"
75 forcesolver.
addForce(&forcedescriptor);
79 const std::string tutorial_01_BioSolverMembraneTension::name =
"tutorial_01_BioSolverMembraneTension";
void deleteForceType(const unsigned int forcetype)
reset all forces
const Geometry & getGeometry() const
getter for the geometry
virtual void applyBioProcess(geometry::GeometryHandler &geometryhandler, solver::ForceSolver &forcesolver)
Applies biological processes.
const std::vector< std::shared_ptr< Connection > > & getConnections() const
getConnections Getter for connections
void addForce(std::stringstream *const forcedescriptor)
add a force
class responsible for generating the internal geometry representation