Trisurf Monte Carlo simulator
Samo Penic
2010-12-28 dac2e5020dc34c236b741ff5c4591244e73f56f2
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"
d7639a 12
SP 13 /** Entrance function to the program
14   * @param argv is a number of parameters used in program call (including the program name
15   * @param argc is a pointer to strings (character arrays) which holds the arguments
16   * @returns returns 0 on success, any other number on fail.
17 */
18
19 int main(int argv, char *argc[]){
20 ts_bool retval;
3131dc 21     ts_vertex_list *vlist=init_vertex_list(5);
b01cc1 22 ts_vertex_list *vlist1;
a10dd5 23 ts_bond_list *blist=init_bond_list();
a2a676 24 ts_triangle_list *tlist=init_triangle_list();
bb77ca 25 ts_cell_list *clist=init_cell_list(3,3,3,0.3);
7958e9 26 ts_vesicle *vesicle;
d7639a 27
a10dd5 28 retval=vtx_add_cneighbour(blist,vlist->vtx[1],vlist->vtx[0]);
d7639a 29 if(retval==TS_FAIL) printf("1. already a member or vertex is null!\n");
a10dd5 30
SP 31 retval=vtx_add_cneighbour(blist,vlist->vtx[0],vlist->vtx[1]);
d7639a 32 if(retval==TS_FAIL) printf("2. already a member or vertex is null!\n");
a10dd5 33
9802f1 34 retval=vtx_remove_neighbour(vlist->vtx[0],vlist->vtx[1]);
SP 35 vtx_add_neighbour(vlist->vtx[0],vlist->vtx[1]);
36
a10dd5 37 vlist->vtx[0]->data->x=1.0;
SP 38 vlist->vtx[0]->data->x=1.1;
b01cc1 39 vlist1=vertex_list_copy(vlist);
a10dd5 40 bond_add(blist, vlist->vtx[1],vlist->vtx[0]);
a2a676 41 triangle_add(tlist,vlist->vtx[1],vlist->vtx[2],vlist->vtx[3]);
a10dd5 42
a2a676 43 triangle_add(tlist,vlist->vtx[1],vlist->vtx[2],vlist->vtx[3]);
SP 44
bb77ca 45 printf("Cell idx=1 has vertices=%u\n",clist->cell[0]->data->nvertex);
SP 46 cell_add_vertex(clist->cell[0], vlist->vtx[0]);
47 printf("Cell idx=1 has vertices=%u\n",clist->cell[0]->data->nvertex);
48 printf("Cell idx=1 has vertex[0] has x coordinate=%e\n",clist->cell[0]->data->vertex[0]->data->x);
49 cell_list_cell_occupation_clear(clist);
50 printf("Cell idx=1 has vertices=%u\n",clist->cell[0]->data->nvertex);
51 cell_add_vertex(clist->cell[0], vlist->vtx[0]);
52
53
a2a676 54 triangle_list_free(tlist);
a10dd5 55 bond_list_free(blist);
73f967 56 vtx_list_free(vlist);
bb77ca 57 cell_list_free(clist);
7958e9 58
b01cc1 59 vtx_list_free(vlist1);
SP 60 printf("Tests complete.\n");
dac2e5 61 vesicle=initial_distribution_dipyramid(17,60,60,60,0.15);
SP 62
63 centermass(vesicle);
64 cell_occupation(vesicle);
65
a6b1b5 66 write_vertex_xml_file(vesicle,0);
SP 67 write_master_xml_file("test.pvd");
3131dc 68 write_dout_fcompat_file(vesicle,"dout");
7958e9 69 vesicle_free(vesicle);
SP 70
d7639a 71 return 0; //program finished perfectly ok. We return 0.
SP 72