PARFEAP requests

[Colin McAuliffe]

Hello all, here are a few requests for PARFEAP:

1. Use petsc pcfieldsplit as a parallel implementation of the PARTition command

2. Write petsc ksp view and memory usage information to the O files

Both or either of these would be hugely useful!

All the best,

Colin

[FEAP_Admin]

Do you know of 'easy' ways to either of these?

[Colin McAuliffe]

For the petsc monitor/memory logging it is possible to make a user defined monitor routine which can make any relevant data available for writing to a file. Here is an example:

http://www.mcs.anl.gov/petsc/petsc-dev/src/ksp/ksp/examples/tutorials/ex2f.F.html


PCFieldSplit would be more complicated, I briefly looked into how to this in August but did not have much time to spend on it and didn't really make progress. The idea with the new Field Split stuff in petsc is that the user specifies the data layout defining the fields in a PetscSection, which is then set as the default section in a DM. Then, the solution methods for the fields or any combination of fields can be specified from the command line. This is described in this paper:

http://www.mcs.anl.gov/uploads/cels/papers/P2017-0112.pdf

Also this email thread might help

http://lists.mcs.anl.gov/pipermail/petsc-users/2012-August/014611.html

I am still not entirely clear on how this can be done for unblocked format though.

Thanks
Colin