Yes good call, I pulled master and got the patch.
I don't have errors anymore.
Do this mean that when a dm is active, the dm is assumed to create hold
the information about the vectors (residual, solution and internal work
vectors) for the problem, whereas when the dm is not active or is not
used, the ksp is assumed to hold the information about vectors format?
Best,
Luc
On 11/21/2014 12:08 PM, Lawrence Mitchell wrote:
On 21 Nov 2014, at 16:55, Luc Berger-Vergiat <[email protected]
<mailto:[email protected]>> wrote:
Sure, so with gamg (which was my first choice too), I get the following:
[0]PETSC ERROR: --------------------- Error Message
--------------------------------------------------------------
[0]PETSC ERROR:
[0]PETSC ERROR: Must call DMShellSetGlobalVector() or
DMShellSetCreateGlobalVector()
[0]PETSC ERROR: See
http://www.mcs.anl.gov/petsc/documentation/faq.html for trouble shooting.
[0]PETSC ERROR: Petsc Release Version 3.5.2, Sep, 08, 2014
[0]PETSC ERROR: /home/luc/research/feap_repo/ShearBands/parfeap/feap
on a arch-opt named euler by
I wonder if this is the same problem described here
https://bitbucket.org/petsc/petsc/pull-request/221
Where the ksp contains an inactive dmshell that doesn't know how to
create vectors
Lawrence