Changes between Version 17 and Version 18 of Piernik
- Timestamp:
- 03/10/13 14:21:14 (12 years ago)
Legend:
- Unmodified
- Added
- Removed
- Modified
-
Piernik
v17 v18 1 1 = Introduction = 2 This page ilustrates MUSCLization of real astrophysics code [http://piernik.astri.umk.pl/ PIERNIK] ( "Gingerbread").2 This page ilustrates MUSCLization of real astrophysics code [http://piernik.astri.umk.pl/ PIERNIK] (what stands for "Gingerbread" in polish). 3 3 4 4 = The orginal code = … … 212 212 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 214 Now we need to add corresponding `MUSCLE_Send` calls to the MHD kernel after the gas density and momentum isrecomputed:214 Now we need to add corresponding `MUSCLE_Send` calls to the MHD kernel after the gas density and momentum values are recomputed: 215 215 {{{ 216 216 call MUSCLE_Send("rho_gas", rho_gas, %REF(size(rho_gas)), MUSCLE_DOUBLE) !here we assume that real is always real*8 … … 218 218 }}} 219 219 220 == Synchronizing delta timesteps == 221 Until now we silently omitted one important aspect of the PIERNIK simulaitons: the delta timestep is not constant and depends on the current state of both modules: MHD and MC. In the original code the lower of the `dt,,mc,,` and dt,,mhd,, were chosen. 222 {{{ 223 224 }}} 225 == MPI enabling == 226 227 == Results == 228 229 == Further works ==