On 20/07/15 22:51, Matthew Knepley wrote:

It looks a bit like DMPlexConstructGhostCells() might be messing with the DM's sfNatural property? or does this not work the way I think it should?

Thats a bug. CreateGhostCells() makes an entirely new DM, and the SetNatural is a piece of state which did not get carried over. I will fix it.

Was there a fix for this? I'm just starting to look at this stuff again, and it still seems to be having the same problem. Cheers, Adrian -- Dr Adrian Croucher Senior Research Fellow Department of Engineering Science University of Auckland, New Zealand email: [email protected] tel: +64 (0)9 923 84611