Laboratoire RDP » Biophysics Team » Cell Division Rule

#include "compute_areas.h"

#include "compute_cell_areas.h"

#include "compute_barycenters.h"

#include "compute_shapes.h"

#include "compute_lengths.h"

#include "initialize.h"

#include "constants.h"

#include <stdio.h>

#include "compute_cell_walls.h"

/* this function initializes the various lists used in the program */

void initialize(double *vertices, int vertices_number, int *cells_vertices, int *vertices_number_in_each_cell, int cell_number, int *wall_vertices, int wall_number, int *cell_walls,

                double *barycenters, double *shapes, double *targets, double *lengths, double *target_lengths, double *areas, double *cell_areas, double *forces, double *stress_periclinal, double *stress_anticlinal){

int i;

compute_barycenters(vertices,vertices_number,cells_vertices,vertices_number_in_each_cell,cell_number,barycenters);

compute_shapes(vertices, cells_vertices, vertices_number_in_each_cell, cell_number, barycenters, shapes);

compute_areas(vertices, cells_vertices, vertices_number_in_each_cell, cell_number, areas);

compute_cell_areas(vertices, cells_vertices, vertices_number_in_each_cell, cell_number, areas, cell_areas);

compute_cell_walls(vertices, cells_vertices, vertices_number_in_each_cell, cell_number, cell_walls, wall_vertices, wall_number);

compute_lengths(vertices, wall_vertices, wall_number, lengths);

for (i=0;i<cell_number;i++){

    targets[3*i]=initial_growth*shapes[3*i]; targets[3*i+1]=initial_growth*shapes[3*i+1]; targets[3*i+2]=initial_growth*shapes[3*i+2];

    stress_periclinal[3*i]=0.; stress_periclinal[3*i+1]=0.; stress_periclinal[3*i+2]=0.;

}

for (i=0;i<wall_number;i++){

    target_lengths[i]=lengths[i]; stress_anticlinal[3*i]=0.; stress_anticlinal[3*i+1]=0.; stress_anticlinal[3*i+2]=0.;

}

}