[Users] w_lorentz factor in TOVSolver external.inc

[Frank Loeffler]

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

-------------- next part --------------
A non-text attachment was scrubbed...
Name: not available
Type: application/pgp-signature
Size: 835 bytes
Desc: not available
Url : http://lists.einsteintoolkit.org/pipermail/users/attachments/20100527/44dcd604/attachment.bin