stardis

commit 02bc98ff940b7eebde778c6e200041c2bf43743a
parent 7f9079b9c9e7ac53eb71d1e26d7c0a26fed818e0
Author: christophe coustet <christophe.coustet@meso-star.com>
Date:   Tue, 23 Oct 2018 17:36:40 +0200

Adapt to stardis 0.5 API changes

Diffstat:
M
src/stardis-compute.c
 | 
9
+++++----

1 file changed, 5 insertions(+), 4 deletions(-)
diff --git a/src/stardis-compute.c b/src/stardis-compute.c
@@ -397,7 +397,7 @@ stardis_compute(struct stardis* stardis, enum stardis_mode mode)
 
   int* bound2fluid = NULL;
   double pos[3] = {0,0,0};
-  double time = 0;
+  double time[2] = { 0, 0};
   size_t nfailures;
   int i = 0;
 
@@ -555,7 +555,7 @@ stardis_compute(struct stardis* stardis, enum stardis_mode mode)
     pos[0] = stardis->probe[0];
     pos[1] = stardis->probe[1];
     pos[2] = stardis->probe[2];
-    time   = stardis->probe[3];
+    time[0] = time[1]  = stardis->probe[3];
     {
       double uv[2] = {0,0};
       size_t iprim = -1;
@@ -570,7 +570,8 @@ stardis_compute(struct stardis* stardis, enum stardis_mode mode)
       if(iprim == -1){
         SDIS(solve_probe(scn,
               stardis->N, 
-              pos,time, 
+              pos,
+              time, 
               stardis->scale_factor,
               stardis->radiative_temp[0], 
               stardis->radiative_temp[1], 
@@ -595,7 +596,7 @@ stardis_compute(struct stardis* stardis, enum stardis_mode mode)
 
     /* Print the results */
     printf("Temperature at [%g, %g, %g, %g] = %g +/- %g\n",
-        pos[0], pos[1], pos[2], time,
+        pos[0], pos[1], pos[2], time[0],
         temperature.E, /* Expected value */
         temperature.SE); /* Standard error */
     printf("#failures: %lu/%lu\n",