From ad4e7069f10b80446be8d13eba28a03c2c642a3b Mon Sep 17 00:00:00 2001 From: Samo Penic <samo.penic@gmail.com> Date: Fri, 13 Jun 2014 20:37:10 +0000 Subject: [PATCH] Now it works after remake! However, some crossection segments are in the list multiple times. --- src/main.c | 22 +++++++++++++++++++++- 1 files changed, 21 insertions(+), 1 deletions(-) diff --git a/src/main.c b/src/main.c index 0343378..f300719 100644 --- a/src/main.c +++ b/src/main.c @@ -1,5 +1,6 @@ #include<stdio.h> #include<math.h> +#include<stdlib.h> #include "general.h" #include "vertex.h" #include "bond.h" @@ -11,6 +12,8 @@ #include "frame.h" #include "timestep.h" #include "poly.h" +#include "sh.h" +#include "shcomplex.h" /** Entrance function to the program * @param argv is a number of parameters used in program call (including the program name @@ -38,12 +41,29 @@ ts_fprintf(stdout,"**********************************************************************\n\n"); tape=parsetape("tape"); vesicle=restore_state(&start_iteration); + if(vesicle==NULL){ + ts_fprintf(stderr, "Dump file does not exist or is not a regular file! Did you mean to invoke trisurf with --force-from-tape option?\n\n"); + return 1; + } // nove vrednosti iz tapea... vesicle->bending_rigidity=tape->xk0; vtx_set_global_values(vesicle); + vesicle->pswitch =tape->pswitch; vesicle->pressure=tape->pressure; + vesicle->dmax=tape->dmax*tape->dmax; + poly_assign_filament_xi(vesicle,tape); + free(vesicle->tape); + vesicle->tape=tape; + vesicle->clist->dmin_interspecies = tape->dmin_interspecies*tape->dmin_interspecies; + /* spherical harmonics */ + if(tape->shc>0){ + vesicle->sphHarmonics=complex_sph_init(vesicle->vlist,tape->shc); + } + else { + vesicle->sphHarmonics=NULL; + } - if(command_line_args.reset_iteration_count) start_iteration=tape->inititer+1; + if(command_line_args.reset_iteration_count) start_iteration=tape->inititer; else start_iteration++; if(start_iteration>=tape->iterations){ -- Gitblit v1.9.3