solstice-pp

commit 4f3de7960d725181ae08d10776f6bcc3b2438a7a
parent 93dc656ea6c758b1260164f2ce08546782bf4227
Author: Christophe Coustet <christophe.coustet@meso-star.com>
Date:   Tue,  4 Sep 2018 17:10:21 +0200

Fix output data of solpp

Set data type to double in vtk headers as data are actually double.

Diffstat:
M
solpp.c
 | 
12
++++++------

1 file changed, 6 insertions(+), 6 deletions(-)
diff --git a/solpp.c b/solpp.c
@@ -45,7 +45,7 @@ mesh_write_vtk(FILE* output, struct mesh* msh)
   fprintf(output, "DATASET POLYDATA\n");
 
   n = BUF_SZ(msh->coords)/3;
-  fprintf(output, "POINTS %zu float\n", n);
+  fprintf(output, "POINTS %zu double\n", n);
   FOR_EACH(i, 0, n) {
     fprintf(output, "%g %g %g\n",
       BUF_AT(msh->coords, i*3+0),
@@ -99,7 +99,7 @@ mesh_write_prim_data_vtk
   fprintf(output, "CELL_DATA %zu\n", n);
 
   fprintf(output, "FIELD PrimaryData %zu\n", 2 + BUF_SZ(simul->rcvs)*6);
-  fprintf(output, "cos_factor 2 %zu float\n", n);
+  fprintf(output, "cos_factor 2 %zu double\n", n);
   FOR_EACH(iprim, 0, BUF_SZ(msh->entities)) {
     CHK(prim = find_primary_by_id(simul, BUF_AT(msh->entities, iprim)));
     FOR_EACH(icell, 0, BUF_AT(msh->ncells, iprim)) {
@@ -107,7 +107,7 @@ mesh_write_prim_data_vtk
     }
   }
 
-  fprintf(output, "shadow_loss 2 %zu float\n", n);
+  fprintf(output, "shadow_loss 2 %zu double\n", n);
   FOR_EACH(iprim, 0, BUF_SZ(msh->entities)) {
     CHK(prim = find_primary_by_id(simul, BUF_AT(msh->entities, iprim)));
     FOR_EACH(icell, 0, BUF_AT(msh->ncells, iprim)) {
@@ -116,7 +116,7 @@ mesh_write_prim_data_vtk
   }
 
   #define WRITE(Side, Name) {                                                  \
-    fprintf(output, "%s_"STR(Side)"_"STR(Name)" 2 %zu float\n", rcv->name, n); \
+    fprintf(output, "%s_"STR(Side)"_"STR(Name)" 2 %zu double\n", rcv->name, n); \
     FOR_EACH(iprim, 0, BUF_SZ(msh->entities)) {                                \
       CHK(rcvXprim = find_rcvXprim(simul,rcv->id,BUF_AT(msh->entities,iprim)));\
       FOR_EACH(icell, 0, BUF_AT(msh->ncells, iprim)) {                         \
@@ -152,7 +152,7 @@ mesh_write_rcv_data_vtk
   fprintf(output, "FIELD PrimaryData 12\n");
 
   #define WRITE(Side, Name) {                                                  \
-    fprintf(output, STR(Side)"_"STR(Name)" 2 %zu float\n", n);                 \
+    fprintf(output, STR(Side)"_"STR(Name)" 2 %zu double\n", n);                 \
     FOR_EACH(ircv, 0, BUF_SZ(msh->entities)) {                                 \
       CHK(rcv = find_receiver_by_id(simul, BUF_AT(msh->entities, ircv)));      \
       FOR_EACH(icell, 0, BUF_AT(msh->ncells, ircv)) {                          \
@@ -171,7 +171,7 @@ mesh_write_rcv_data_vtk
   #undef WRITE
 
   #define WRITE_MAP(Flux, Side) {                                              \
-    fprintf(output, STR(Side)"_"STR(Flux)"_map 2 %zu float\n", n);             \
+    fprintf(output, STR(Side)"_"STR(Flux)"_map 2 %zu double\n", n);             \
     FOR_EACH(ircv, 0, BUF_SZ(msh->entities)) {                                 \
       CHK(rcv = find_receiver_by_id(simul, BUF_AT(msh->entities, ircv)));      \
       if(!BUF_SZ(rcv->map[Flux][Side])) {                                      \