Hi,
> Once I get to the office, I will prepare a short list of all the
sequence of calls to PETSc and settings that I use for that case. With
that in principle, you should be able to reproduce the scenario, once
you get a sparsity pattern, an actual matrix and a rhs vector.
Thanks. I hope that the options suffice in order to reproduce it with
one of our examples.
Actually, I can also send you files with the Petsc matrix, the vector and the
sparsity pattern data on, say, 2 processors so that you can try to read them in
and solve the exact same system...
I'll try to reproduce with our examples first and let you know if I need
the vector and matrix data.
Best regards,
Karli
On Jan 22, 2014, at 10:07 AM, Karl Rupp <[email protected]> wrote:
Hi,
Well, I don't usually use PETSc runtime options. I have only a subset of
parameters (PC types, KSP types, etc.) which are interfaced by my solver and
exposed to the end-user. Other few settings (matrix types for instance) are
hardcoded in some objects (which my linear system solver interface makes use
of).
Okay, which parameters are these? Apparently you set PCASM and KSPGMRES at some
point, but are there any others I need to be aware of?
I can setup a testcase, but it cannot be too simple. You'll have to install my
code to run it. I can setup a script for this. If you want to give it a try,
I'll send you some instructions tomorrow (here is night already).
I think it's easier if I try to reproduce it directly with some of our examples
first. If I fail to reproduce it, I'll ask you for a test case :-)
Best regards,
Karli