Laboratoire RDP » Biophysics Team » Vascular Networks in Growing Leaves

#include "optimize_width.h"

#include "compute_energy_width.h"

#include <nlopt.h>

#include <stdio.h>

#include <stdlib.h>

#include <cmath>

void optimize_width(int vertex_number, double *vertices, int edge_number, int *edges_vertices,

                   double *lengths, double *rest_lengths, double *edges_width, int *periodizer, double *width, double *height){

nlopt_opt opt;

opt=nlopt_create(NLOPT_LN_BOBYQA,edge_number);

nlopt_set_stopval(opt, -HUGE_VAL);

nlopt_set_maxeval(opt, -1.);

nlopt_set_maxtime(opt, -1.);

nlopt_set_xtol_rel(opt, -1.);

nlopt_set_xtol_abs1(opt, 1.e-6);

nlopt_set_ftol_rel(opt, -1.);

nlopt_set_ftol_abs(opt, -1.);

struct parameters_list3 {int vertex_number; double *vertices; int edge_number; int *edges_vertices;

                        double *lengths; double *rest_lengths; int *periodizer; double *width; double *height;};

struct parameters_list3 parameters;

parameters.vertex_number=vertex_number;

parameters.vertices=vertices;

parameters.edge_number=edge_number;

parameters.edges_vertices=edges_vertices;

parameters.lengths=lengths;

parameters.rest_lengths=rest_lengths;

parameters.periodizer=periodizer;

parameters.width=width;

parameters.height=height;

void* pparameters=&parameters;

double width_optimize[edge_number]; double energy;

for (int i=0; i<edge_number;i++){ width_optimize[i]=edges_width[i];}

/// first resolution: find the reference width and height without the anisotropic stress

int reso=nlopt_set_min_objective(opt, compute_energy_width, pparameters);

nlopt_optimize(opt, width_optimize, &energy);

printf("%d\n",reso);

for (int i=0; i<edge_number;i++){edges_width[i]=width_optimize[i];}

nlopt_destroy(opt);

}