From 19c1fd84d1fd75e9bfd55feb670fda87a2e02086 Mon Sep 17 00:00:00 2001 From: Samo Penic <samo.penic@gmail.com> Date: Fri, 08 Feb 2019 08:32:29 +0000 Subject: [PATCH] Merge branch 'master' of ssh://git.penic.eu:29418/trisurf-ng --- src/io.c | 24 ++++++++++++++++++++++-- 1 files changed, 22 insertions(+), 2 deletions(-) diff --git a/src/io.c b/src/io.c index c6a35cb..943d96b 100644 --- a/src/io.c +++ b/src/io.c @@ -826,6 +826,7 @@ ts_bond_list *blist=vesicle->blist; ts_vertex **vtx=vlist->vtx; ts_uint i,j; + //ts_double senergy=0.0; char filename[10000]; char just_name[255]; FILE *fh; @@ -974,9 +975,26 @@ } 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,"</PointData>\n<CellData>\n"); + + if(vesicle->tape->stretchswitch==1){ + fprintf(fh,"<DataArray type=\"Float64\" Name=\"stretching_energy\" format=\"ascii\">"); + for(i=0;i<blist->n;i++){ + fprintf(fh, "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,"</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); } @@ -1165,6 +1183,8 @@ CFG_SIMPLE_INT("constvolswitch",&tape->constvolswitch), CFG_SIMPLE_INT("constareaswitch",&tape->constareaswitch), CFG_SIMPLE_FLOAT("constvolprecision",&tape->constvolprecision), + CFG_SIMPLE_INT("stretchswitch",&tape->stretchswitch), + CFG_SIMPLE_FLOAT("xkA0",&tape->xkA0), CFG_SIMPLE_FLOAT("pressure",&tape->pressure), CFG_SIMPLE_FLOAT("k_spring",&tape->kspring), CFG_SIMPLE_FLOAT("xi",&tape->xi), -- Gitblit v1.9.3