Re: [petsc-users] How to create a quadrature object?

On Thu, Jul 28, 2016 at 7:08 PM, Andrew Ho wrote: > > On Thu, Jul 28, 2016 at 4:35 PM, Matthew Knepley > wrote: > >> >> This is strange. First, you should only have uninitialized classids if >> you build without dynamics libraries. Did you? If you >> include

Re: [petsc-users] How to create a quadrature object?

On Thu, Jul 28, 2016 at 4:35 PM, Matthew Knepley wrote: > > This is strange. First, you should only have uninitialized classids if you > build without dynamics libraries. Did you? If you > include the entire error output, it would show us. > > This is pretty much the entire

Re: [petsc-users] How to create a quadrature object?

On Thu, Jul 28, 2016 at 3:43 PM, Andrew Ho wrote: > I am trying to create a very simple quadrature object, but for some reason > PETSc keeps giving me an "invalid argument" error. > > Relevant code: > > #include >> int main(int argc, char** argv) >> { >>

[petsc-users] How to create a quadrature object?

I am trying to create a very simple quadrature object, but for some reason PETSc keeps giving me an "invalid argument" error. Relevant code: #include > int main(int argc, char** argv) > { > CHKERRQ(PetscInitialize(, , nullptr, "quadrature testing")); > PetscQuadrature quad; >

[petsc-users] Implementing discontinuous Galerkin FEM?

I am trying to implement a discontinuous Galerkin discretization using the PETSc DM features to handle most of the topology/geometry specific functions. However, I'm not really sure which direction to approach this from since DG is kind of a middle ground between finite volume and traditional

[petsc-users] PCMG with matrix-free operators accessing DMDA

Hi everyone, this time I am using PETSc to do something that is more complicated than my usual and I want to do it at the highest possible abstraction level. To put it in a nutshell, my intent is to build a parallel multigrid to solve a linear system via DM, KSP and PCMG (I would like to use DMMG

Re: [petsc-users] Comprehensive example for KSPRegister()

On Thu, Jul 28, 2016 at 9:07 AM, JC wrote: > Hey everyone, > > I was wondering if any of you had a comprehensive example of KSPRegister() > to create our own KSP solver in fortran? > I have tried to look online but have not been able to find it. > We do not currently have

[petsc-users] Comprehensive example for KSPRegister()

Hey everyone, I was wondering if any of you had a comprehensive example of KSPRegister() to create our own KSP solver in fortran? I have tried to look online but have not been able to find it. Thanks a lot, JC

[petsc-users] block matrix without MatCreateNest

I'm trying to understand how to assemble a block matrix in a format-independent manner, so that I can switch between types mpiaij and matnest. The manual states that the key to format-independent assembly is to use MatGetLocalSubMatrix. So, in the code below, I'm using this to assemble a 3-by-3

Re: [petsc-users] Nested Fieldsplit for custom index sets

Dear Artur, On 28/07/16 02:20, Safin, Artur wrote: > Barry, Lawrence, > >> I think the SubKSPs (and therefore SubPCs) are not set up until you call >> KSPSetUp(ksp) which your code does not do explicitly and is therefore done >> in KSPSolve. > > I added KSPSetUp(), but unfortunately the issue

Re: [petsc-users] Ignore command line arguments with fortran code using PETSc

2016-07-27 21:42 GMT+02:00 Barry Smith : > > Actually there is currently no way to PetscInitialize from Fortran without > adding the command line options to the database. In the middle > of petscinitialize_() is the code fragment > > PETScParseFortranArgs_Private(,); >