Trisurf Monte Carlo simulator
Samo Penic
2016-02-15 383644169457b45c9c6bc047c469c8e1ca1d2562
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
#include<stdio.h>
#include<general.h>
#include<snapshot.h>
ts_bool xml_trisurf_data(FILE *fh, ts_vesicle *vesicle){
    xml_trisurf_header(fh, vesicle);
    xml_trisurf_tria(fh,vesicle->tlist);
    xml_trisurf_tria_neigh(fh,vesicle->tlist);
    xml_trisurf_vtx_neigh(fh,vesicle->vlist);
    xml_trisurf_vtx_tristar(fh,vesicle->vlist);
    xml_trisurf_footer(fh);
    return TS_SUCCESS;
}
 
ts_bool xml_trisurf_header(FILE *fh, ts_vesicle *vesicle){
    fprintf(fh, "<trisurf nvtx=\"%u\" npoly=\"%u\" nfono=\"%u\">\n", vesicle->vlist->n, vesicle->poly_list->n, vesicle->poly_list->poly[0]->vlist->n);
    return TS_SUCCESS;
}
 
ts_bool xml_trisurf_footer(FILE *fh){
    fprintf(fh, "</trisurf>\n");
    return TS_SUCCESS;
}
 
ts_bool xml_trisurf_tria(FILE *fh, ts_triangle_list *tlist){
    ts_uint i;
    fprintf(fh,"<tria>\n");
    for(i=0; i<tlist->n;i++){
        fprintf(fh,"%u %u %u\n",tlist->tria[i]->vertex[0]->idx, tlist->tria[i]->vertex[1]->idx, tlist->tria[i]->vertex[2]->idx);
    }
    fprintf(fh,"</tria>\n");
    return TS_SUCCESS;
}
 
ts_bool xml_trisurf_tria_neigh(FILE *fh, ts_triangle_list *tlist){
    ts_uint i;
    fprintf(fh,"<trianeigh>\n");
    for(i=0; i<tlist->n;i++){
        fprintf(fh,"%u %u %u\n",tlist->tria[i]->neigh[0]->idx, tlist->tria[i]->neigh[1]->idx, tlist->tria[i]->neigh[2]->idx);
    }
    fprintf(fh,"</trianeigh>\n");
    return TS_SUCCESS;
}
 
ts_bool xml_trisurf_vtx_neigh(FILE *fh, ts_vertex_list *vlist){
    ts_uint i,j;
    for(i=0;i<vlist->n;i++){
        fprintf(fh,"<vtxn idx=\"%u\">",vlist->vtx[i]->idx);
        for(j=0;j<vlist->vtx[i]->neigh_no;j++){
            fprintf(fh,"%u ",vlist->vtx[i]->neigh[j]->idx);
        }
        fprintf(fh, "</vtxn>\n");
    }
    return TS_SUCCESS;
}
 
ts_bool xml_trisurf_vtx_tristar(FILE *fh, ts_vertex_list *vlist){
    ts_uint i,j;
    for(i=0;i<vlist->n;i++){
        fprintf(fh,"<tristar idx=\"%u\">",vlist->vtx[i]->idx);
        for(j=0;j<vlist->vtx[i]->tristar_no;j++){
            fprintf(fh,"%u ",vlist->vtx[i]->tristar[j]->idx);
        }
        fprintf(fh, "</tristar>\n");
    }
    return TS_SUCCESS;
}