math.hpp

Go to the documentation of this file.

 
#ifndef MATH_HPP
#define MATH_HPP
 
#include <vector>
#include "Simulation.hpp"
 
#include <eigen3/Eigen/Eigenvalues>
 
std::pair<float, float> calculate_center_of_mass_cell(const std::vector<std::pair<int, int>>& m_position);
 
std::array<float, 3> calculate_inertia_tensor_cell(const std::vector<std::pair<int, int>>& m_position);
 
std::array<std::array<float, 2>, 2> calculate_eigenvectors_inertia_tensor_cell(const std::vector<std::pair<int, int>>& m_position);
 
std::pair<float, float> calculate_projection_on_principal_axis(const std::vector<std::pair<int, int>>& m_position);
 
float get_length_cell(const std::vector<std::pair<int, int>>& m_position);
 
// PDE solver
 
void initializeGrid(std::vector<std::vector<double>>& grid, double initConc);
 
void applyBoundaryConditions(std::vector<std::vector<double>>& grid, double boundaryConc);
 
double secretionTerm(double secretionRate, int id_cell);
 
double decayTerm(double decayRate, int id_cell);
 
void solve_aStepOfDiffusionSecretionDecay(std::vector<std::vector<double>>& grid, double dx, double dy, double dt, double D, 
        double boundaryConc, double secretionRate, double decayRate, std::vector<std::vector<int>> lattice);
 
#endif