Re: [petsc-users] Element to local dof map using dmplex

2016-10-25 Thread Morten Nobel-Jørgensen
Dear Matt Did you (or anyone else) find time to look at our issue? We are really looking forward to your answer :) Kind regards, Morten From: Matthew Knepley [knep...@gmail.com] Sent: Wednesday, October 12, 2016 3:41 PM To: Morten Nobel-Jørgensen Cc: petsc-users

[petsc-users] Element to local dof map using dmplex

2016-10-12 Thread Morten Nobel-Jørgensen
Dear PETSc developers / Matt Thanks for your suggestions regarding our use of dmplex in a FEM context. However, Matt's advise on using the PetscFE is not sufficient for our needs (our end goal is a topology optimization framework - not just FEM) and we must honestly admit that we do not see

Re: [petsc-users] DMPlex problem

2016-09-26 Thread Morten Nobel-Jørgensen
ind regards, Morten From: Matthew Knepley [knep...@gmail.com] Sent: Monday, September 26, 2016 2:19 PM To: Morten Nobel-Jørgensen Cc: PETSc ‎[petsc-users@mcs.anl.gov]‎ Subject: Re: [petsc-users] DMPlex problem On Mon, Sep 26, 2016 at 7:00 AM, Morten Nobel-Jørgen

Re: [petsc-users] DMPlex problem

2016-09-26 Thread Morten Nobel-Jørgensen
From: petsc-users-boun...@mcs.anl.gov [petsc-users-boun...@mcs.anl.gov] on behalf of Morten Nobel-Jørgensen [m...@dtu.dk] Sent: Sunday, September 25, 2016 11:15 AM To: Matthew Knepley Cc: PETSc ‎[petsc-users@mcs.anl.gov]‎ Subject: Re: [petsc-users] DMPlex problem Hi

Re: [petsc-users] DMPlex problem

2016-09-25 Thread Morten Nobel-Jørgensen
Nobel-Jørgensen Cc: PETSc ‎[petsc-users@mcs.anl.gov]‎ Subject: Re: [petsc-users] DMPlex problem On Fri, Sep 23, 2016 at 7:45 AM, Matthew Knepley <knep...@gmail.com<mailto:knep...@gmail.com>> wrote: On Fri, Sep 23, 2016 at 3:48 AM, Morten Nobel-Jørgensen <m...@dtu.dk<mailto:m

Re: [petsc-users] DMPlex problem

2016-09-23 Thread Morten Nobel-Jørgensen
From: Matthew Knepley [knep...@gmail.com] Sent: Friday, September 09, 2016 12:21 PM To: Morten Nobel-Jørgensen Cc: PETSc ‎[petsc-users@mcs.anl.gov]‎ Subject: Re: [petsc-users] DMPlex problem On Fri, Sep 9, 2016 at 4:04 AM, Morten Nobel-Jørgensen <m...@dtu.dk<mailto:m...@dtu.dk>> wrote

[petsc-users] DMPlex problem

2016-09-09 Thread Morten Nobel-Jørgensen
Dear PETSc developers and users, Last week we posted a question regarding an error with DMPlex and multiple dofs and have not gotten any feedback yet. This is uncharted waters for us, since we have gotten used to an extremely fast feedback from the PETSc crew. So - with the chance of sounding

Re: [petsc-users] Distribution of DMPlex for FEM

2016-08-30 Thread Morten Nobel-Jørgensen
or tips will be appreciated :) Kind regards, Morten From: Morten Nobel-Jørgensen Sent: Thursday, July 14, 2016 9:45 AM To: Matthew Knepley Cc: petsc-users@mcs.anl.gov Subject: Re: [petsc-users] Distribution of DMPlex for FEM Hi Matthew Thanks for your answer

Re: [petsc-users] Distribution of DMPlex for FEM

2016-07-14 Thread Morten Nobel-Jørgensen
sers@mcs.anl.gov<mailto:petsc-users@mcs.anl.gov>" <petsc-users@mcs.anl.gov<mailto:petsc-users@mcs.anl.gov>> Emne: Re: [petsc-users] Distribution of DMPlex for FEM On Wed, Jul 13, 2016 at 3:57 AM, Morten Nobel-Jørgensen <m...@dtu.dk<mailto:m...@dtu.dk>> wrote:

[petsc-users] Duplicate cells when exporting a distributed dmplex

2016-07-05 Thread Morten Nobel-Jørgensen
Hi all, I hope someone can help me with the following: I’m having some problems when exporting a distributed DMPlex – the cells (+cell types) seems to be duplicated. When I’m running the code on a non-distributed system it works as expected, but when I run it on multiple processors (2 in my

Re: [petsc-users] DMPlex : Assemble global stiffness matrix problem

2016-04-05 Thread Morten Nobel-Jørgensen
outputs the matrix to a file called ‘Kmat.m'. Kind regards, Morten Fra: Matthew Knepley <knep...@gmail.com<mailto:knep...@gmail.com>> Dato: Tuesday 8 March 2016 at 19:47 Til: Morten Nobel-Jørgensen <m...@mek.dtu.dk<mailto:m...@mek.dtu.dk>> Cc: "petsc-users@mcs.anl.gov&l

[petsc-users] DMPlex : Assemble global stiffness matrix problem

2016-03-07 Thread Morten Nobel-Jørgensen
I have some problems using DMPlex on unstructured grids in 3D. After I have created the DMPlex and assigned dofs (3 dofs on each node), I run into some problems when assembling the global stiffness matrix. I have created a small example in the attached cc file. My problems are: * It seems

Re: [petsc-users] DMPlex: Ghost points after DMRefine

2015-12-01 Thread Morten Nobel-Jørgensen
I found a solution to my problem by using the global section instead. I still don’t quite understand what my problem ISLocalToGlobalMapping was. > Yes, that is the solution I use. Thanks – good to hear that I’m on the right track :) > I think I ignore nonlocal indices in the l2g mapping

Re: [petsc-users] DMPlex: Ghost points after DMRefine

2015-12-01 Thread Morten Nobel-Jørgensen
From: Morten Nobel-Jørgensen <m...@mek.dtu.dk<mailto:m...@mek.dtu.dk>> Date: Monday 30 November 2015 at 17:24 To: Matthew Knepley <knep...@gmail.com<mailto:knep...@gmail.com>> Cc: "petsc-users@mcs.anl.gov<mailto:petsc-users@mcs.anl.gov>" <petsc-users@mcs.anl

Re: [petsc-users] DMPlex: Ghost points after DMRefine

2015-11-30 Thread Morten Nobel-Jørgensen
DOUT); } From: Matthew Knepley <knep...@gmail.com<mailto:knep...@gmail.com>> Date: Monday 30 November 2015 at 14:08 To: Morten Nobel-Jørgensen <m...@mek.dtu.dk<mailto:m...@mek.dtu.dk>> Cc: "petsc-users@mcs.anl.gov<mailto:petsc-users@mcs.anl.gov>" <petsc-

[petsc-users] DMPlex: Ghost points after DMRefine

2015-11-30 Thread Morten Nobel-Jørgensen
I have a very simple unstructured mesh composed of two triangles (four vertices) with one shared edge using a DMPlex: /|\ / | \ \ | / \|/ After distributing this mesh to two processes, each process owns a triangle. However one process owns tree vertices, while the last vertex is owned by the

[petsc-users] Problem Iterating DMPlex

2015-11-17 Thread Morten Nobel-Jørgensen
After distributing a DMPlex it seems like my cells are appearing twice (or rather multiple cells maps onto the same vertices). I’m assuming the way I’m iterating the DMPlex is wrong. Essentially I iterate the DMPlex the following way after distribution (see code snippet below – or attached