[petsc-dev] mesh face areas with DMPlexComputeCellGeometryFVM()

[Adrian Croucher]

On 13/04/14 16:08, Matthew Knepley wrote:
>
> Thanks for the test. It was a bug in a special case:
>
> -        coords[3] = PetscRealPart(coords[3*dim+0] - coords[0*dim+0]);
> +        coords[3] = PetscRealPart(coords[3*dim+1] - coords[0*dim+1]);
>
> I will push it to next

Thanks very much- that's fixed it.

The problems I was seeing last week with VecView() writing to VTK seem 
to have gone now too. Wunderbar!

Cheers, Adrian

-- 
Dr Adrian Croucher
Department of Engineering Science
University of Auckland
New Zealand
tel 64-9-373-7599 ext 84611

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20140414/d4c76346/attachment.html>