From 620ba72765eda244ca3f3806b8ad709929e74b3c Mon Sep 17 00:00:00 2001
From: Samo Penic <samo.penic@gmail.com>
Date: Thu, 21 Jun 2018 10:50:23 +0000
Subject: [PATCH] Isak's proposal for stretching energy calculation, take 1
---
src/io.c | 44 ++++++++++++++++----------------------------
1 files changed, 16 insertions(+), 28 deletions(-)
diff --git a/src/io.c b/src/io.c
index 83bb26d..3ea3662 100644
--- a/src/io.c
+++ b/src/io.c
@@ -826,7 +826,7 @@
ts_bond_list *blist=vesicle->blist;
ts_vertex **vtx=vlist->vtx;
ts_uint i,j;
- ts_double senergy=0.0;
+ //ts_double senergy=0.0;
char filename[10000];
char just_name[255];
FILE *fh;
@@ -974,39 +974,27 @@
}
}
fprintf(fh,"</DataArray>\n");
+
+
+ fprintf(fh,"</PointData>\n<CellData>\n");
+
if(vesicle->tape->stretchswitch==1){
fprintf(fh,"<DataArray type=\"Float64\" Name=\"stretching_energy\" format=\"ascii\">");
- for(i=0;i<vlist->n;i++){
- senergy=0.0;
- for(j=0;j<vtx[i]->tristar_no;j++){
- senergy+=vtx[i]->tristar[j]->energy;
- }
- fprintf(fh,"%.17e ",senergy);
+ for(i=0;i<blist->n;i++){
+ fprintf(fh, "0.0 ");
}
- //polymeres
- if(poly){
- poly_idx=vlist->n;
- for(i=0;i<vesicle->poly_list->n;i++){
- for(j=0;j<vesicle->poly_list->poly[i]->vlist->n;j++,poly_idx++){
- fprintf(fh,"%.17e ", 0.0);
- }
- }
- }
- //filaments
- if(fil){
- poly_idx=vlist->n+monono*polyno;
- for(i=0;i<vesicle->filament_list->n;i++){
- for(j=0;j<vesicle->filament_list->poly[i]->vlist->n;j++,poly_idx++){
- // fprintf(stderr,"was here\n");
- fprintf(fh,"%.17e ", 0.0);
- }
- }
- }
+ for(i=0;i<monono*polyno+filno*(fonono-1);i++){
+ fprintf(fh,"0.0 ");
+ }
+ for(i=0;i<vesicle->tlist->n;i++){
+ fprintf(fh,"%.17e ",vesicle->tlist->tria[i]->energy);
+ }
fprintf(fh,"</DataArray>\n");
}
-
- fprintf(fh,"</PointData>\n<CellData>\n</CellData>\n<Points>\n<DataArray type=\"Float64\" Name=\"Koordinate tock\" NumberOfComponents=\"3\" format=\"ascii\">\n");
+
+
+ fprintf(fh,"</CellData>\n<Points>\n<DataArray type=\"Float64\" Name=\"Koordinate tock\" NumberOfComponents=\"3\" format=\"ascii\">\n");
for(i=0;i<vlist->n;i++){
fprintf(fh,"%.17e %.17e %.17e\n",vtx[i]->x,vtx[i]->y, vtx[i]->z);
}
--
Gitblit v1.9.3