Re: [petsc-users] Periodic BC in petsc DMPlex / firedrake

2018-10-25 Thread Maximilian Hartig
>> Ahh, thanks. I was missing the option " -dm_plex_gmsh_periodic “. But using >> this option I now generate a segmentation fault error when calling VecView() >> on the solution vector with vtk and hdf5 viewers. Any suggestions? >> > Small example? VTK is deprecated. HDF5 should work,

Re: [petsc-users] Periodic BC in petsc DMPlex / firedrake

2018-10-25 Thread Stephen Wornom
Did you forget to attach "The nice picture attached" Stephen > From: "Stefano Zampini" > To: "petsc-maint" > Cc: "imilian hartig" , "PETSc users list" > > Sent: Thursday, October 25, 2018 4:47:43 PM > Subject: Re: [p

Re: [petsc-users] Periodic BC in petsc DMPlex / firedrake

2018-10-25 Thread Stefano Zampini
Opened the PR https://bitbucket.org/petsc/petsc/pull-requests/1203/fix-dump-vtk-field-with-periodic-meshes/diff Il giorno gio 25 ott 2018 alle ore 17:44 Matthew Knepley ha scritto: > Good catch Stefano. > > Matt > > On Thu, Oct 25, 2018 at 9:36 AM Stefano Zampini > wrote: > >> Maybe this is

Re: [petsc-users] Periodic BC in petsc DMPlex / firedrake

2018-10-25 Thread Matthew Knepley
Good catch Stefano. Matt On Thu, Oct 25, 2018 at 9:36 AM Stefano Zampini wrote: > Maybe this is a fix > > diff --git a/src/dm/impls/plex/plexvtu.c b/src/dm/impls/plex/plexvtu.c > index acdea12c2f..1a8bbada6a 100644 > --- a/src/dm/impls/plex/plexvtu.c > +++ b/src/dm/impls/plex/plexvtu.c > @@

Re: [petsc-users] Periodic BC in petsc DMPlex / firedrake

2018-10-25 Thread Stefano Zampini
Matt, you can reproduce it via $ valgrind ./ex12 -quiet -run_type test -interpolate 1 -bc_type dirichlet -petscspace_degree 1 -vec_view vtk:test.vtu:vtk_vtu -f ${PETSC_DIR}/share/petsc/datafiles/meshes/square_periodic.msh -dm_plex_gmsh_periodic Long time ago I added support for viewing meshes

Re: [petsc-users] Periodic BC in petsc DMPlex / firedrake

2018-10-24 Thread Matthew Knepley
On Wed, Oct 24, 2018 at 11:36 AM Maximilian Hartig wrote: > > > On 24. Oct 2018, at 12:49, Matthew Knepley wrote: > > On Wed, Oct 24, 2018 at 6:29 AM Lawrence Mitchell wrote: > >> Hi Max, >> >> (I'm cc'ing in the petsc-users mailing list which may have more advice, >> if you are using PETSc

Re: [petsc-users] Periodic BC in petsc DMPlex / firedrake

2018-10-24 Thread Maximilian Hartig
> On 24. Oct 2018, at 12:49, Matthew Knepley wrote: > > On Wed, Oct 24, 2018 at 6:29 AM Lawrence Mitchell > wrote: > Hi Max, > > (I'm cc'ing in the petsc-users mailing list which may have more advice, if > you are using PETSc you should definitely subscribe! > > > On

Re: [petsc-users] Periodic BC in petsc DMPlex / firedrake

2018-10-24 Thread Stefano Zampini
To read periodic meshes from GMSH, you need to use the option -dm_plex_gmsh_periodic and DMPlexCreateFromFile See src/dm/impls/plex/examples/tests/ex1.c. An example runs $ ./ex1 -filename ${PETSC_DIR}/share/petsc/datafiles/meshes/cube_periodic_bin.msh -dm_plex_gmsh_periodic -dm_view

Re: [petsc-users] Periodic BC in petsc DMPlex / firedrake

2018-10-24 Thread Matthew Knepley
On Wed, Oct 24, 2018 at 6:29 AM Lawrence Mitchell wrote: > Hi Max, > > (I'm cc'ing in the petsc-users mailing list which may have more advice, if > you are using PETSc you should definitely subscribe! > > > On 24 Oct 2018, at 09:27, Maximilian Hartig > wrote: > > > > Hello Lawrence, > > > >

Re: [petsc-users] Periodic BC in petsc DMPlex / firedrake

2018-10-24 Thread Lawrence Mitchell
Hi Max, (I'm cc'ing in the petsc-users mailing list which may have more advice, if you are using PETSc you should definitely subscribe! > On 24 Oct 2018, at 09:27, Maximilian Hartig wrote: > > Hello Lawrence, > > sorry to message you out of the blue. My name is Max and I found your post on