[petsc-users] Cannot use KSP if I choose local Matrix/Vector ownership

[William A. Perkins]

Hello,

I have run into a problem using KSP.  Attached is a hopefully short,
simple example program (which has nothing to do with my actual
application) demonstrating the problem.  Also attached is an output
log from several runs of the program.  The program runs fine if I
choose local ownership according to PetscSplitOwnership().  If,
however, I slightly modify that local ownership scheme, the program
will not run in parallel.  In my view, KSP appears to make its own
coefficient (or preconditioner, or whatever) matrix that has different
local ownership than the one I give to KSPSetOperators(), and so
cannot do matrix/vector multiplies.

I would like to be able to choose my own local ownership. Any help
would be appreciated. Hopefully, I'm not doing something stupid.

Thanks.

Bill

-- 
Bill Perkins
Research Engineer
Hydrology Group

Pacific Northwest National Laboratory
902 Battelle Boulevard
P.O. Box 999, MSIN K9-36
Richland, WA  99352 USA
Tel:  509-372-6131
Fax: 509-372-6089
william.perkins at pnnl.gov
www.pnnl.gov
-------------- next part --------------
An embedded and charset-unspecified text was scrubbed...
Name: versteeg.cpp
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20130926/9dacd86a/attachment-0002.ksh>
-------------- next part --------------
An embedded and charset-unspecified text was scrubbed...
Name: versteeg.log
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20130926/9dacd86a/attachment-0003.ksh>