[Users] [Commits] [svn:einsteintoolkit] GRHydro/trunk/src/ (Rev. 319)

Bruno C. Mundim bcmsma at astro.rit.edu
Sat May 12 11:34:24 CDT 2012


On 05/11/2012 06:46 PM, schnetter at cct.lsu.edu wrote:
> @@ -478,7 +480,7 @@
>   	 end do
>         else  !!$ Really implement C&W, page 197; which requires stencil 4.
>   	 do i = 4, nx - 3
> -	    s=sign(1.d0, -dpress(i))
> +	    s=sign(one, -dpress(i))
>   	    flatten = max(tilde_flatten(i), tilde_flatten(i+s))
>   	    if (abs(1.d0 - flatten)>  0.d0) then
>   	       trivial_rp(i-1) = .false.
Hi Erik,

you missed the line above: if (abs(1.d0 - flatten)>  0.d0) then

Cheers,
Bruno.


More information about the Users mailing list