Re: [petsc-users] Error: DM global to natural SF was not created when DMSetUseNatural has already been called

On Wed, Nov 28, 2018 at 8:58 PM Danyang Su via petsc-users < petsc-users@mcs.anl.gov> wrote: > Dear All, > > I got the following error when using DMPlexGlobalToNatural function > using 1 processor. > We do not create that mapping on 1 proc because the orderings are the same. Reordering happens

Re: [petsc-users] Do we support to use NORM_X instead of NORM_2 to check the convergence of KSP and SNES?

> On Nov 28, 2018, at 3:20 PM, Fande Kong via petsc-users > wrote: > > Hi Developers, > > I just checked into the SNES and KSP code. We always hard code the Vec Norm > as NORM_2 when computing the linear and nonlinear residuals. > > Does this mean we have to use norm_2 to check the

Re: [petsc-users] A question regarding a potential use case for DMNetwork

Markus : Can you provide a concrete example of your application? Is your model built on a network? Is so, give us a simple example of such network and the math systems over it. This would help us understand your request. Hong Dear PETSc users and developers, in particular those experienced with

[petsc-users] A question regarding a potential use case for DMNetwork

Dear PETSc users and developers, in particular those experienced with the relatively new DMNetwork object, I would like to get some advice on wether it makes sense for my application to be built using this PETSc feature or if I am equally well served if I use plain Vec and Mat objects. In the

Re: [petsc-users] DMPlex global to natural ordering

Have a look at src/dm/impls/plex/examples/tests/ex26.c It does exactly this (among other stuff) using exodus format Blaise > On Nov 28, 2018, at 1:51 PM, Danyang Su via petsc-users > wrote: > > Dear All, > > My simulation needs to pass initial condition from external file to the code. >

[petsc-users] IS Invert Not Permutation

Hey, thanks for your quick reply! As far as I understand how MatSetLocalToGlobalMapping works, there is no way to create such a mapping that can be used for MatSetValuesLocal? What I need is basically MatSetValuesLocal row/col argument = 3 / 4 / 6 maps to local row/col 0 / 1 / 2 I seems to

Re: [petsc-users] Solving complex linear sparse matrix in parallel + external library

On Wed, Nov 28, 2018 at 10:26 AM Sal Am wrote: > Thank you indeed --download-mpich and using PETSC_ARCH/bin/mpiexec seems > to work. > > Now I am wondering about the other problem namely getting the residual, is > the residual only computed when using iterative solvers? Cause using > richardson

Re: [petsc-users] Solving complex linear sparse matrix in parallel + external library

On Wed, Nov 28, 2018 at 10:27 AM Sal Am via petsc-users < petsc-users@mcs.anl.gov> wrote: > Thank you indeed --download-mpich and using PETSC_ARCH/bin/mpiexec seems > to work. > > Now I am wondering about the other problem namely getting the residual, is > the residual only computed when using