Trisurf Monte Carlo simulator
Samo Penic
2013-11-24 84af719f37b54185ddd7daf0627421143874905c
commit | author | age
4e56fe 1 #include<stdio.h>
SP 2 #include<math.h>
3 #include "general.h"
4 #include "vertex.h"
5 #include "bond.h"
6 #include "triangle.h"
7 #include "cell.h"
8 #include "vesicle.h"
9 #include "io.h"
10 #include "initial_distribution.h"
11 #include "frame.h"
12 #include "timestep.h"
13
14
15 int main(int argv, char *argc[]){
16 ts_vesicle *vesicle;
17 ts_uint inititer,mcsweeps, iterations;
18 vesicle=parsetape(&mcsweeps, &inititer, &iterations);
19
20 int i,j,k;
21 /* TAINTING TEST */
22 for(k=0;k<2;k++){
23
24 printf("*** Starting test No. %d\n",k+1);
25 for(i=0;i<vesicle->vlist->n;i++){
26     vertex_taint(vesicle->vlist->vtx[i],1);
27     vertex_untaint(vesicle->vlist->vtx[i],1);
c85b51 28     printf("%d taint amount %u\n",i,*(vesicle->vlist->vtx[i]->locked));
4e56fe 29     if(vertex_tainted(vesicle->vlist->vtx[i], 0, 1)) printf("Ooops, Main  %d is tainted!\n",i);
SP 30     for(j=0;j<vesicle->vlist->vtx[i]->neigh_no;j++){
31         if(vertex_tainted(vesicle->vlist->vtx[i]->neigh[j], 0, 1)){
32             printf("Ooops, The neighbour %d of %d is tainted!\n", j, i);
33         }
34     }
35     
36 }
37 }
38
39
50f10f 40 printf("At the end you should see no Ooopses and all taint amounts should be 0!\n");
4e56fe 41 vesicle_free(vesicle);
SP 42 return 0;
43 }