diff --git a/src/mhd/DirectInduction/mhd.h b/src/mhd/DirectInduction/mhd.h index ef83e645851d100fe2c43ea108d5236a8bf43c7b..170e670928ea60a9b37fd4e881c5b7bb144cdcf6 100644 --- a/src/mhd/DirectInduction/mhd.h +++ b/src/mhd/DirectInduction/mhd.h @@ -371,9 +371,9 @@ __attribute__((always_inline)) INLINE static void mhd_end_gradient( Cos_Ind_Diff = (Adv_B_source[0]*Delta_B[0]+Adv_B_source[1]*Delta_B[1]+Adv_B_source[2]*Delta_B[2]); - p->mhd_data.eta_OWAR += 1.0f/OW * ( 0.5f * p->h * p->h / (absB+FLT_MIN)) * 0.5f * fmaxf(0.0f,1.0f-Cos_Ind_Diff) * Abs_Adv_B_source * p->rho * (p->mass * kernel_root); + p->mhd_data.eta_OWAR += 1.0f/OW * ( 0.5f * p->h * p->h / (absB+FLT_MIN)) * 0.5f * fmaxf(0.0f,1.0f-Cos_Ind_Diff) * Abs_Adv_B_source * p->rho; /** * (p->mass * kernel_root); - p->mhd_data.eta_OWAR *= pow_dimension(1.0f/p->h) / p->rho; + p->mhd_data.eta_OWAR *= pow_dimension(1.0f/p->h) / p->rho; **/ if (p->mhd_data.eta_OWAR<0.0f){ error( diff --git a/src/mhd/DirectInduction/mhd_iact.h b/src/mhd/DirectInduction/mhd_iact.h index 5f9f65cebe20d0e1234f883555099b417e04bb79..2ab6195e2b10307e1aa7ba4a53f4122e167f52d9 100644 --- a/src/mhd/DirectInduction/mhd_iact.h +++ b/src/mhd/DirectInduction/mhd_iact.h @@ -158,7 +158,7 @@ __attribute__((always_inline)) INLINE static void runner_iact_mhd_gradient( } /* Calculate OWAR */ - +/** float OW; OW = 1.0f; @@ -216,7 +216,7 @@ __attribute__((always_inline)) INLINE static void runner_iact_mhd_gradient( "Error: incorrect OWAR " ); } - +**/ } @@ -307,7 +307,7 @@ runner_iact_nonsym_mhd_gradient(const float r2, const float dx[3], } /* Calculate OWAR */ - +/** float OW; OW = 1.0f; @@ -340,7 +340,7 @@ runner_iact_nonsym_mhd_gradient(const float r2, const float dx[3], "Error: incorrect OWAR " ); } - +**/ } /**