Trisurf Monte Carlo simulator
mihaf
2014-03-25 fe24d29b3c8684f08dccd01f5785aa48b7137322
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);
ea1cce 51         vesicle->clist->dmin_interspecies = tape->dmin_interspecies*tape->dmin_interspecies;
d9cb01 52
144784 53         if(command_line_args.reset_iteration_count) start_iteration=tape->inititer;
8a6614 54         else start_iteration++;
083e03 55
8a6614 56         if(start_iteration>=tape->iterations){
SP 57             ts_fprintf(stdout, "Simulation already completed. if you want to rerun it try with --force-from-tape or --reset-iteration-count\n\n");
58             return 0;
59         }
60     }
61
62     run_simulation(vesicle, tape->mcsweeps, tape->inititer, tape->iterations, start_iteration);
63     write_master_xml_file("test.pvd");
64     write_dout_fcompat_file(vesicle,"dout");
65     vesicle_free(vesicle);
66     tape_free(tape);
67     return 0; //program finished perfectly ok. We return 0.
d7639a 68