Dynamics.hpp

Go to the documentation of this file.

 
#pragma once
 
#include <math.h>
#include <cmath>
#include <vector>
#include <iostream>
#include <fstream>
#include <sstream>
#include <omp.h>
#include <random>
#include <algorithm>
// model setup
#include "mesh/Mesh.hpp"
#include "mesh/Face.hpp"
#include "mesh/Vertex.hpp"
#include "mesh/Gauss_quadrature.hpp"
#include "energy_force/Energy.hpp"
#include "energy_force/Force.hpp"
#include "model/Model.hpp"
 
// matrix math
#include "linalg/Linear_algebra.hpp"
// parameters
#include "Parameters.hpp"
 
class DynamicMesh : public Mesh
{
public:
    Matrix mesh2surface;
    Matrix surface2mesh;
    Matrix matMesh;
    Matrix matSurface;
 
    DynamicMesh(Param &srcParam);
 
    void setup_flat();
 
    void assign_mesh2surface();
 
    void update_vertices_mat_with_vector();
 
    void update_vertices_vector_with_mat();
 
    void postprocess_ghost_periodic();
 
protected:
    int get_relative_pt_periodic(int i, int n, int m);
};
 
class DynamicModel : public Model
{
public:
    DynamicMesh &mesh; 
    double randScaleConst;
    double forceScaleConst;
    unsigned int randomSeed;
    std::mt19937 gen;
    std::normal_distribution<> normal_dist;
 
    DynamicModel(DynamicMesh &mesh_, Record &record_);
 
    void next_step();
};
 
/*
//boundary condition
int get_relatve_pt_periodic(int i, int n, int m);
void postprocess_ghost_periodic(Param& param, gsl_matrix *verticesOnMesh, std::vector<Vertex>& vertex);
*/