Laboratoire RDP » Biophysics Team » Vascular Networks in Growing Leaves

#include "step.h"

#include <stdio.h>

#include <math.h>

#include "minimize_BFGS.h"

#include "growth_and_chemistry.h"

#include <gsl/gsl_odeiv2.h>

#include <gsl/gsl_errno.h>

#include "compute_energy_and_forces.h"

#include "save_data.h"

#include <cstdlib>

#include "constants.h"

#include <algorithm>

#include "compute_lengths.h"

/* this function performs a time step (growth and chemistry then mechanics) */

void step(double *vertices, int vertex_number, int edge_number, int *edges_vertices, double *lengths, double *rest_lengths, double *edges_width, int *periodizer, double *width, double *height){

/* minimization of energy: move the vertices in order to make the energy as small as possible, i.e. in order to make cells

   as close as possible of their targets */

minimize_BFGS(vertex_number, vertices, edge_number, edges_vertices, lengths, rest_lengths, edges_width, periodizer, width, height);

compute_lengths(vertices, edges_vertices, edge_number, lengths, periodizer, width[0], height[0]);

struct parameters_list2 {int edge_number; double *lengths; double *edges_width;};

struct parameters_list2 parameters2;

parameters2.edge_number=edge_number;

parameters2.lengths=lengths;

parameters2.edges_width=edges_width;

void *pparameters2=&parameters2;

int system_size=edge_number;

double x[system_size];

for (int i=0;i<edge_number;i++){ x[i]=rest_lengths[i];}

gsl_odeiv2_system sys={growth_and_chemistry, NULL, system_size, pparameters2};

gsl_odeiv2_driver *d=gsl_odeiv2_driver_alloc_y_new (&sys, gsl_odeiv2_step_rk2, 1.e-8, 1.e-8, 0.);

int s; double t=0.0;

s=gsl_odeiv2_driver_apply(d, &t, 1.e-7, x);

for (int i=0;i<edge_number;i++){ rest_lengths[i]=x[i];}

gsl_odeiv2_driver_free(d);

}