Re: [Users] w_lorentz factor in TOVSolver external.inc

On Wed, May 26, 2010 at 07:49:17AM -0400, Roland Haas wrote:

•         w_lorentz_2*w_lorentz_2 *
  

•         (-2 * v_2 / (source[i3D]/my_psi4/my_psi4 + TOV_K[0]*pow(rhonew, TOV_Gamma[0]))) *
  

•         TOV_K[0]*TOV_Gamma[0]*pow(rhonew, TOV_Gamma[0]-1);
  

Maybe it would be better to write it a bit different, along the lines:

-2 * w_lorentz_2*w_lorentz_2 * v_2 * TOV_K[0]*TOV_Gamma[0]*pow(rhonew, TOV_Gamma[0]-1) / (source[i3D]/my_psi4/my_psi4 + TOV_K[0]*pow(rhonew, TOV_Gamma[0]))

Frank