| | |
|---|
| | | centermass(vesicle); |
|---|
| | | cell_occupation(vesicle); |
|---|
| | | vesicle_volume(vesicle); //needed for constant volume at this moment |
|---|
| | | vesicle_area(vesicle); //needed for constant area at this moment |
|---|
| | | vesicle_area(vesicle); //needed for constant area and stretching energy at this moment |
|---|
| | | if(V0<0.000001) |
|---|
| | | V0=vesicle->volume; |
|---|
| | | ts_fprintf(stdout,"Setting volume V0=%.17f\n",V0); |
|---|
| | |
|---|
| | | cell_occupation(vesicle); |
|---|
| | | dump_state(vesicle,i); |
|---|
| | | if(i>=inititer){ |
|---|
| | | vesicle_volume(vesicle); //calculates just volume. |
|---|
| | | vesicle_area(vesicle); //calculates area. |
|---|
| | | write_vertex_xml_file(vesicle,i-inititer,NULL); |
|---|
| | | write_master_xml_file(command_line_args.output_fullfilename); |
|---|
| | | epochtime=get_epoch(); |
|---|
| | | gyration_eigen(vesicle, &l1, &l2, &l3); |
|---|
| | | vesicle_volume(vesicle); //calculates just volume. |
|---|
| | | vesicle_area(vesicle); //calculates area. |
|---|
| | | r0=getR0(vesicle); |
|---|
| | | if(vesicle->sphHarmonics!=NULL){ |
|---|
| | | preparationSh(vesicle,r0); |
|---|
