Re: [petsc-users] Extending petsc4py

I am sorry. I meant to ask about petsc4py only, so please replace every mpi4py in my question with petsc4py. The "demo/wrap-swig" has solved my problems. Thanks Ingo 2018-07-11 11:26 GMT+02:00 Lisandro Dalcin : > > > On Tue, 10 Jul 2018 at 20:17, Ingo Gaertner > wrote

[petsc-users] Extending petsc4py

Hello, can I find somewhere an example how to wrap my own PETSc-based code so that it can be used together with petsc4py? Let's assume I have a C function, which I want to access from python: Vec createMyVec(); I have written a SWIG interface file to generate a python wrapper for this function.

Re: [petsc-users] dmplex face normals orientation

Thank you, Matt, this answers my question. Ingo 2017-04-18 22:21 GMT+02:00 Matthew Knepley <knep...@gmail.com>: > On Tue, Apr 18, 2017 at 2:33 PM, Ingo Gaertner <ingogaertner@gmail.com > > wrote: > >> The part that does not make sense is: >> The code calcu

Re: [petsc-users] dmplex face normals orientation

l other faces the normal direction and their orientation sign agree with what you said. Thanks Ingo 2017-04-14 11:28 GMT+02:00 Ingo Gaertner <ingogaertner@gmail.com>: > Thank you, Matt, > as you say, the face orientations do change sign when switching between > the two adjacent c

Re: [petsc-users] Off-diagonal matrix-vector product y=(A-diag(A))x

sig-email?utm_medium=email_source=link_campaign=sig-email_content=webmail> Virenfrei. www.avast.com <https://www.avast.com/sig-email?utm_medium=email_source=link_campaign=sig-email_content=webmail> <#DAB4FAD8-2DD7-40BB-A1B8-4E2AA1F9FDF2> 2017-04-14 16:50 GMT+02:00 Jed Brown <j...@j

Re: [petsc-users] dmplex face normals orientation

plemented.\n",dim); exit(1); } } ierr = PetscFinalize(); } 2017-04-14 11:00 GMT+02:00 Matthew Knepley <knep...@gmail.com>: > On Wed, Apr 12, 2017 at 10:52 AM, Ingo Gaertner < > ingogaertner@gmail.com> wrote: > >> Hello, >> I have problem

[petsc-users] Off-diagonal matrix-vector product y=(A-diag(A))x

Does PETSc include an efficient implementation for the operation y=(A-diag(A))x or y_i=\sum_{j!=i}A_{ij}x_j on a sparse matrix A? In words, I need a matrix-vector product after the matrix diagonal has been set to zero. For efficiency reasons I can't copy and modify the matrix or first calculate

[petsc-users] dmplex face normals orientation

Hello, I have problems determining the orientation of the face normals of a DMPlex. I create a DMPlex, for example with DMPlexCreateHexBoxMesh(). Next, I get the face normals using DMPlexComputeGeometryFVM(DM dm, Vec *cellgeom, Vec *facegeom). facegeom gives the correct normals, but I don't know

Re: [petsc-users] examples of DMPlex*FVM methods

2017-04-05 19:56 GMT+02:00 Jed Brown <j...@jedbrown.org>: > Ingo Gaertner <ingogaertner@gmail.com> writes: > > > Hi Matt, > > I don't care if FV is suboptimal to solve the Poisson equation. I only > want > > to better understand the method by getting

Re: [petsc-users] examples of DMPlex*FVM methods

, because it seems like a matter of taste, and our professor of numerical methods for CFD seems to strongly prefer FVM without a solid basis to justify his preference. Thanks Ingo 2017-04-05 18:56 GMT+02:00 Matthew Knepley <knep...@gmail.com>: > On Wed, Apr 5, 2017 at 11:50 AM, Ingo

Re: [petsc-users] examples of DMPlex*FVM methods

Hi Jed, thank you for your reply. Two followup questions below: 2017-04-04 22:18 GMT+02:00 Jed Brown <j...@jedbrown.org>: > Ingo Gaertner <ingogaertner@gmail.com> writes: > > > We have never talked about Riemann solvers in our CFD course, and I don't > > unde

Re: [petsc-users] examples of DMPlex*FVM methods

2017-04-03 23:58 GMT+02:00 Matthew Knepley : > There are no tutorials, and almost no documentation. > Uhh, I'm not sure whether it makes sense for me to use PETSc then. > The best thing to look at is TS ex11. This solves a bunch of different > equations > (advection, shallow

[petsc-users] examples of DMPlex*FVM methods

Dear all, as part of my studies I would like to implement a simple finite volume CFD solver (following the textbook by Ferziger) on an unstructured, distributed mesh. It seems like the DMPlex class with its DMPlex*FVM methods has prepared much of what is needed for such a CFD solver. Unfortunately