Trisurf Monte Carlo simulator
Samo Penic
2014-03-24 144784496550cccc927607be8edc487af7aa851f
commit | author | age
d7639a 1 #include<stdio.h>
SP 2 #include<math.h>
3 #include "general.h"
4 #include "vertex.h"
a10dd5 5 #include "bond.h"
a2a676 6 #include "triangle.h"
bb77ca 7 #include "cell.h"
7958e9 8 #include "vesicle.h"
a6b1b5 9 #include "io.h"
7958e9 10 #include "initial_distribution.h"
dac2e5 11 #include "frame.h"
aec47d 12 #include "timestep.h"
8db569 13 #include "poly.h"
d7639a 14
SP 15 /** Entrance function to the program
16   * @param argv is a number of parameters used in program call (including the program name
17   * @param argc is a pointer to strings (character arrays) which holds the arguments
18   * @returns returns 0 on success, any other number on fail.
19 */
20
21 int main(int argv, char *argc[]){
1ab449 22     ts_vesicle *vesicle;
SP 23     ts_tape *tape;
8a6614 24     ts_uint start_iteration=0;
152052 25     force_from_tape=0;
8a6614 26     parse_args(argv,argc); // sets global variable command_line_args (defined in io.h)
SP 27     ts_fprintf(stdout,"Starting program...\n\n");
152052 28     if(command_line_args.force_from_tape){
8a6614 29         ts_fprintf(stdout,"************************************************\n");
SP 30         ts_fprintf(stdout,"**** Generating initial geometry from tape *****\n");
31         ts_fprintf(stdout,"************************************************\n\n");
32         tape=parsetape("tape");
33         vesicle=create_vesicle_from_tape(tape);
34     } else {
083e03 35
8a6614 36         ts_fprintf(stdout,"**********************************************************************\n");
SP 37         ts_fprintf(stdout,"**** Recreating vesicle from dump file and continuing simulation *****\n");
38         ts_fprintf(stdout,"**********************************************************************\n\n");
39         tape=parsetape("tape");
40         vesicle=restore_state(&start_iteration);
3f5c83 41         if(vesicle==NULL){
SP 42             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");
43             return 1;
44         }
152052 45         // nove vrednosti iz tapea...
d9cb01 46         vesicle->bending_rigidity=tape->xk0;
M 47         vtx_set_global_values(vesicle);
152052 48         vesicle->pressure=tape->pressure;
2dcc65 49         vesicle->dmax=tape->dmax*tape->dmax;
bcf455 50         poly_assign_filament_xi(vesicle,tape);
d9cb01 51
144784 52         if(command_line_args.reset_iteration_count) start_iteration=tape->inititer;
8a6614 53         else start_iteration++;
083e03 54
8a6614 55         if(start_iteration>=tape->iterations){
SP 56             ts_fprintf(stdout, "Simulation already completed. if you want to rerun it try with --force-from-tape or --reset-iteration-count\n\n");
57             return 0;
58         }
59     }
60
61     run_simulation(vesicle, tape->mcsweeps, tape->inititer, tape->iterations, start_iteration);
62     write_master_xml_file("test.pvd");
63     write_dout_fcompat_file(vesicle,"dout");
64     vesicle_free(vesicle);
65     tape_free(tape);
66     return 0; //program finished perfectly ok. We return 0.
d7639a 67