From 8db569a42c280be13ea9edbe4c528e0041b6fd3f Mon Sep 17 00:00:00 2001 From: mihaf <miha.fosnaric@gmail.com> Date: Tue, 03 Dec 2013 13:05:49 +0000 Subject: [PATCH] Samo&Miha poly_list_free added. --- src/main.c | 38 ++++++++++++++++++++++++-------------- 1 files changed, 24 insertions(+), 14 deletions(-) diff --git a/src/main.c b/src/main.c index 96473ba..b756320 100644 --- a/src/main.c +++ b/src/main.c @@ -9,6 +9,8 @@ #include "io.h" #include "initial_distribution.h" #include "frame.h" +#include "timestep.h" +#include "poly.h" /** Entrance function to the program * @param argv is a number of parameters used in program call (including the program name @@ -17,37 +19,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_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]); @@ -58,12 +63,17 @@ vtx_list_free(vlist1); printf("Tests complete.\n"); -vesicle=initial_distribution_dipyramid(17,60,60,60,0.15); +*/ +vesicle=parsetape(&mcsweeps, &inititer, &iterations); -centermass(vesicle); -cell_occupation(vesicle); +/*Testing */ +vesicle->poly_list=init_poly_list(1400,20,vesicle->vlist); -write_vertex_xml_file(vesicle,0); +poly_list_free(vesicle->poly_list); +/*End testing*/ + +run_simulation(vesicle, mcsweeps, inititer, iterations); + write_master_xml_file("test.pvd"); write_dout_fcompat_file(vesicle,"dout"); vesicle_free(vesicle); -- Gitblit v1.9.3