Changes between Version 22 and Version 23 of Piernik

Timestamp:

03/11/13 00:25:57 (12 years ago)

Author:

mmamonski

Comment:

--

Piernik

@@ -206 +206 @@
 In the original code, before every MC step the rho_gas and m_gas multidimensional arrays were recomputed, in the MUSCLE variant we replace the code with MUSCLE_Receive calls
 {{{
-   call MUSCLE_Receive("rho_gas", rho_gas, %REF(size(rho_gas)), MUSCLE_DOUBLE) !here we assume that real is always real*8
-   call MUSCLE_Receive("m_gas", m_gas, %REF(size(m_gas)), MUSCLE_DOUBLE) !here we assume that real is always real*8
+   call MUSCLE_Receive("rho_gas", rho_gas, size(rho_gas), MUSCLE_DOUBLE) !here we assume that real is always real*8
+   call MUSCLE_Receive("m_gas", m_gas, size(m_gas), MUSCLE_DOUBLE) !here we assume that real is always real*8
 }}}
 As mentioned before the arrays are multidimensional (4 and 3 dimensional saying more precisely) but taking into consideration fact that both kernels are FORTRAN codes and that FORTRAN multidimensional arrays have continues memory layout we can safely cast it to single dimension.
@@ -213 +213 @@
 Now we need to add corresponding `MUSCLE_Send` calls to the MHD kernel after the gas density and momentum values are recomputed:
 {{{
-   call MUSCLE_Send("rho_gas", rho_gas, %REF(size(rho_gas)), MUSCLE_DOUBLE) !here we assume that real is always real*8
-   call MUSCLE_Send("m_gas", m_gas, %REF(size(m_gas)), MUSCLE_DOUBLE) !here we assume that real is always real*8
+   call MUSCLE_Send("rho_gas", rho_gas, size(rho_gas), MUSCLE_DOUBLE) !here we assume that real is always real*8
+   call MUSCLE_Send("m_gas", m_gas, size(m_gas), MUSCLE_DOUBLE) !here we assume that real is always real*8
 }}}
 
@@ -230 +230 @@
 }
 }}}
+And the corresponding code:
+{{{
+#ifdef MHD_KERNEL
+      call time_step(dt)
+      call MUSCLE_Send("dt_final", dt, dt_len, MUSCLE_DOUBLE)
+#else /*MC_KERNEL*/
+      dt = timestep_mc()
+      call MUSCLE_Send("dt_mc", dt, dt_len, MUSCLE_DOUBLE)
+      call MUSCLE_Receive("dt_final", dt, dt_len, MUSCLE_DOUBLE)
+#endif
+}}}
 == MPI enabling ==