From 4e6285b2342038e3740ff9cea7c31364b4caa5e5 Mon Sep 17 00:00:00 2001 From: Samo Penic <samo.penic@fe.uni-lj.si> Date: Tue, 11 Feb 2014 12:03:29 +0000 Subject: [PATCH] Merge branch 'master' of bitbucket.org:samop/trisurf-ng --- src/main.c | 35 ++++++++++++++++------------------- 1 files changed, 16 insertions(+), 19 deletions(-) diff --git a/src/main.c b/src/main.c index c402ebf..7e12c1c 100644 --- a/src/main.c +++ b/src/main.c @@ -18,38 +18,40 @@ */ int main(int argv, char *argc[]){ +ts_uint inititer,mcsweeps, iterations; +ts_vesicle *vesicle; +/* THIS SHOULD GO INTO UNIT TEST ts_bool retval; -ts_uint i; ts_vertex_list *vlist=init_vertex_list(5); ts_vertex_list *vlist1; ts_bond_list *blist=init_bond_list(); ts_triangle_list *tlist=init_triangle_list(); ts_cell_list *clist=init_cell_list(3,3,3,0.3); -ts_vesicle *vesicle; retval=vtx_add_cneighbour(blist,vlist->vtx[1],vlist->vtx[0]); if(retval==TS_FAIL) printf("1. already a member or vertex is null!\n"); -retval=vtx_add_cneighbour(blist,vlist->vtx[0],vlist->vtx[1]); +retval=vtx_add_neighbour(vlist->vtx[0],vlist->vtx[1]); if(retval==TS_FAIL) printf("2. already a member or vertex is null!\n"); - -retval=vtx_remove_neighbour(vlist->vtx[0],vlist->vtx[1]); +fprintf(stderr,"Was here"); +retval=vtx_remove_neighbour(vlist->vtx[1],vlist->vtx[0]); vtx_add_neighbour(vlist->vtx[0],vlist->vtx[1]); +fprintf(stderr,"Was here too!\n"); -vlist->vtx[0]->data->x=1.0; -vlist->vtx[0]->data->x=1.1; +vlist->vtx[0]->x=1.0; +vlist->vtx[0]->x=1.1; vlist1=vertex_list_copy(vlist); bond_add(blist, vlist->vtx[1],vlist->vtx[0]); triangle_add(tlist,vlist->vtx[1],vlist->vtx[2],vlist->vtx[3]); triangle_add(tlist,vlist->vtx[1],vlist->vtx[2],vlist->vtx[3]); -printf("Cell idx=1 has vertices=%u\n",clist->cell[0]->data->nvertex); +printf("Cell idx=1 has vertices=%u\n",clist->cell[0]->nvertex); cell_add_vertex(clist->cell[0], vlist->vtx[0]); -printf("Cell idx=1 has vertices=%u\n",clist->cell[0]->data->nvertex); -printf("Cell idx=1 has vertex[0] has x coordinate=%e\n",clist->cell[0]->data->vertex[0]->data->x); +printf("Cell idx=1 has vertices=%u\n",clist->cell[0]->nvertex); +printf("Cell idx=1 has vertex[0] has x coordinate=%e\n",clist->cell[0]->vertex[0]->x); cell_list_cell_occupation_clear(clist); -printf("Cell idx=1 has vertices=%u\n",clist->cell[0]->data->nvertex); +printf("Cell idx=1 has vertices=%u\n",clist->cell[0]->nvertex); cell_add_vertex(clist->cell[0], vlist->vtx[0]); @@ -60,15 +62,10 @@ vtx_list_free(vlist1); printf("Tests complete.\n"); -vesicle=initial_distribution_dipyramid(17,60,60,60,0.15); +*/ +vesicle=parsetape(&mcsweeps, &inititer, &iterations); +run_simulation(vesicle, mcsweeps, inititer, iterations); -centermass(vesicle); -cell_occupation(vesicle); - -for(i=0;i<100;i++){ -single_timestep(vesicle); -write_vertex_xml_file(vesicle,i); -} write_master_xml_file("test.pvd"); write_dout_fcompat_file(vesicle,"dout"); vesicle_free(vesicle); -- Gitblit v1.9.3