On 27/06/16 13:44, Barry Smith wrote:

Yes, you are correct, there is no MatFDSetType() and in the SetFromOptions we totally cheat http://www.mcs.anl.gov/petsc/petsc-current/src/mat/matfd/fdmatrix.c.html#Mat... and just stick the value directly into the data structure. For now you can call PetscOptionsSetValue(NULL,"-mat_fd_type","ds") in your code (before you start using the MatFD stuff) and we'll add the routine to the master branch.

Thanks, that workaround does the trick for now.

I didn't think anyone would care about the difference between ds and wp with MatFD. It makes a noticeable difference for the problem I'm solving right now. That may be because the problem is not non-dimensionalized (as yet anyway) and so some of the solution elements are of different magnitudes than others. If I non-dimensionalize them it may possibly not matter any more about wp or ds.

- 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