math.cpp File Reference

#include <eigen3/Eigen/Eigenvalues>
#include <fstream>
#include <vector>
#include <iostream>
#include "math.hpp"

Include dependency graph for math.cpp:

Functions
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)
	Calculate the inertia tensor of a cell.
std::array< std::array< float, 2 >, 2 >	calculate_eigenvectors_inertia_tensor_cell (const std::vector< std::pair< int, int > > &m_position)
	Calculate the eigenvectors of the inertia tensor of a cell.
std::pair< float, float >	calculate_projection_on_principal_axis (const std::vector< std::pair< int, int > > &m_position)
	Calculate the projection of a cell on its principal axis.
std::pair< float, float >	calculate_projection_std_on_principal_axis (const std::vector< std::pair< int, int > > &m_position)
	Calculate the standard deviation of the projection of a cell on its principal axis. This way, we can have an idea of the cell length and not the absolute length.
float	get_length_cell (const std::vector< std::pair< int, int > > &m_position)
	Calculate the length of a cell.
void	testGrid ()
	Function to test the grid.
void	initializeGrid (std::vector< std::vector< double > > &grid, double initConc)
	Function to initialize the grid with a given initial concentration.
void	applyBoundaryConditions (std::vector< std::vector< double > > &grid, double boundaryConc)
	Function to apply the boundary conditions to the grid.
double	secretionTerm (double secretionRate, int id_cell)
	Function to define the secretion term S(x, y, t)
double	decayTerm (double decayRate, int id_cell)
	Function to define the decay term D(x, y, t)
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)
	Solve a step of the diffusion-secretion-decay PDE.

Detailed Description

Math functions

Date

28/11/2024

Version

1.0

Author

Lucas BOISTAY

Function Documentation

applyBoundaryConditions()

void applyBoundaryConditions	(	std::vector< std::vector< double > > &	grid,
		double	boundaryConc
	)

Function to apply the boundary conditions to the grid.

Parameters

grid	the grid to apply the boundary conditions
boundaryConc	the boundary concentration value

calculate_eigenvectors_inertia_tensor_cell()

std::array< std::array< float, 2 >, 2 > calculate_eigenvectors_inertia_tensor_cell

(

const std::vector< std::pair< int, int > > &

m_position

)

Calculate the eigenvectors of the inertia tensor of a cell.

Parameters

m_position

All the positions of the cell

Returns

the eigenvectors of the inertia tensor of the cell

calculate_inertia_tensor_cell()

std::array< float, 3 > calculate_inertia_tensor_cell

(

const std::vector< std::pair< int, int > > &

m_position

)

Calculate the inertia tensor of a cell.

Parameters

m_position

All the positions of the cell

Returns

the inertia tensor of the cell

calculate_projection_on_principal_axis()

std::pair< float, float > calculate_projection_on_principal_axis

(

const std::vector< std::pair< int, int > > &

m_position

)

Calculate the projection of a cell on its principal axis.

Parameters

m_position

All the positions of the cell

Returns

the projection of the cell on its 2 principal axis

calculate_projection_std_on_principal_axis()

std::pair< float, float > calculate_projection_std_on_principal_axis

(

const std::vector< std::pair< int, int > > &

m_position

)

Calculate the standard deviation of the projection of a cell on its principal axis. This way, we can have an idea of the cell length and not the absolute length.

Parameters

m_position

All the positions of the cell

Returns

the standard deviation of the projection of the cell on its principal axis

decayTerm()

double decayTerm	(	double	decayRate,
		int	id_cell
	)

Function to define the decay term D(x, y, t)

Parameters

decayRate	the decay rate
id_cell	the id of the cell

Returns

double the value of the decay term

get_length_cell()

float get_length_cell

(

const std::vector< std::pair< int, int > > &

m_position

)

Calculate the length of a cell.

Parameters

m_position

All the positions of the cell

Returns

the max length of the cell

initializeGrid()

void initializeGrid	(	std::vector< std::vector< double > > &	grid,
		double	initConc
	)

Function to initialize the grid with a given initial concentration.

Parameters

grid	the grid to initialize
initConc	the initial concentration value

secretionTerm()

double secretionTerm	(	double	secretionRate,
		int	id_cell
	)

Function to define the secretion term S(x, y, t)

Parameters

secretionRate	the secretion rate
id_cell	the id of the cell

Returns

double the value of the secretion term

solve_aStepOfDiffusionSecretionDecay()

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
	)

Solve a step of the diffusion-secretion-decay PDE.

Parameters

grid	the grid to solve
dx	the x step
dy	the y step (usually equal to dx)
dt	the time step
D	the diffusion coefficient
boundaryConc	the boundary concentration value
secretionRate	the secretion rate
decayRate	the decay rate
lattice	the lattice of the cells