[petsc-users] Error while calling a fieldsplit preconditioner containing a mg/gamg preconditioner

Matthew Knepley knepley at gmail.com
Fri Nov 21 15:16:10 CST 2014

On Fri, Nov 21, 2014 at 3:12 PM, Luc Berger-Vergiat <lb2653 at columbia.edu>

>  Yes good call, I pulled master and got the patch.
> I don't have errors anymore.

+1 Lawrence

> 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?

Yes, exactly.

We wanted a way to have a DM associated with a problem, and yet have the
Vec/Mat passed in
define it. We were trying to force the DM organization down the throat of



> Best,
> Luc
> On 11/21/2014 12:08 PM, Lawrence Mitchell wrote:
>  On 21 Nov 2014, at 16:55, Luc Berger-Vergiat <lb2653 at columbia.edu> wrote:
>  Sure, so with gamg (which was my first choice too), I get the following:
> [0]PETSC ERROR: --------------------- Error Message
> --------------------------------------------------------------
> [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

What most experimenters take for granted before they begin their
experiments is infinitely more interesting than any results to which their
experiments lead.
-- Norbert Wiener
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20141121/3715056f/attachment.html>

More information about the petsc-users mailing list