Re: [petsc-users] PetscViewer with 64bit

2023-02-14 Thread Dave May
On Tue 14. Feb 2023 at 21:27, Jed Brown wrote: > Dave May writes: > > > On Tue 14. Feb 2023 at 17:17, Jed Brown wrote: > > > >> Can you share a reproducer? I think I recall the format requiring > certain > >> things to be Int32. > > > > > > By default, the byte offset used with the appended

Re: [petsc-users] PetscViewer with 64bit

2023-02-14 Thread Jed Brown
Dave May writes: > On Tue 14. Feb 2023 at 17:17, Jed Brown wrote: > >> Can you share a reproducer? I think I recall the format requiring certain >> things to be Int32. > > > By default, the byte offset used with the appended data format is UInt32. I > believe that’s where the sizeof(int) is

Re: [petsc-users] PetscViewer with 64bit

2023-02-14 Thread Dave May
On Tue 14. Feb 2023 at 21:03, Dave May wrote: > > > On Tue 14. Feb 2023 at 17:17, Jed Brown wrote: > >> Can you share a reproducer? I think I recall the format requiring certain >> things to be Int32. > > > By default, the byte offset used with the appended data format is UInt32. > I believe

Re: [petsc-users] PetscViewer with 64bit

2023-02-14 Thread Dave May
On Tue 14. Feb 2023 at 17:17, Jed Brown wrote: > Can you share a reproducer? I think I recall the format requiring certain > things to be Int32. By default, the byte offset used with the appended data format is UInt32. I believe that’s where the sizeof(int) is coming from. This default is

Re: [petsc-users] PetscViewer with 64bit

2023-02-14 Thread Jed Brown
Can you share a reproducer? I think I recall the format requiring certain things to be Int32. Mike Michell writes: > Thanks for the note. > I understood that PETSc calculates the offsets for me through "boffset" > variable in plexvtu.c file. Please correct me if it is wrong. > > If plexvtu.c

Re: [petsc-users] PetscViewer with 64bit

2023-02-14 Thread Mike Michell
Thanks for the note. I understood that PETSc calculates the offsets for me through "boffset" variable in plexvtu.c file. Please correct me if it is wrong. If plexvtu.c has a bug, it could be around "write file header" part in which the boffset is also computed. Is this correct? I am not using

Re: [petsc-users] PetscViewer with 64bit

2023-02-14 Thread Matthew Knepley
On Tue, Feb 14, 2023 at 11:45 AM Mike Michell wrote: > I was trying to modify the header flags from "Int32" to "Int64", but the > problem was not resolved. Could I get any additional comments? > The calculated offsets are not correct I think. Matt > Thanks, > Mike > > >> Thanks for the

Re: [petsc-users] MatFDColoringSetUp with Periodic BC

2023-02-14 Thread Barry Smith
I have created an MR that documents this and moved the error checking to a more appropriate location https://gitlab.com/petsc/petsc/-/merge_requests/6070 > On Feb 14, 2023, at 4:33 AM, Pierre Bernigaud > wrote: > > Dear all, > > I hope this email finds you well. > We are currently

Re: [petsc-users] Help with fieldsplit performance

2023-02-14 Thread Edoardo alinovi
Hi Matt, So I have done some research these days and I have found out that I might try to assemble the SIMPLE for Schur approximation (myS = A11 - A10 inv(DIAGFORM(A00)) A01). Reading papers around, I come up with a doubt, which I believe to be a very silly one but worth asking... Is the way

Re: [petsc-users] PetscViewer with 64bit

2023-02-14 Thread Mike Michell
I was trying to modify the header flags from "Int32" to "Int64", but the problem was not resolved. Could I get any additional comments? Thanks, Mike > Thanks for the comments. > To be precise on the question, the entire part of the header of the .vtu > file is attached: > > > > > >

Re: [petsc-users] Problem setting Fieldsplit fields

2023-02-14 Thread Matthew Knepley
On Tue, Feb 14, 2023 at 8:36 AM Nicolas Barnafi wrote: > Hello Matt, > > After some discussions elsewhere (thanks @LMitchell!), we found out that > the problem is that the fields are setup with PCSetIS, without an attached > DM, which does not support this kind of nesting fields. > > I would

Re: [petsc-users] Segregated solvers in PETSc

2023-02-14 Thread Matthew Knepley
On Tue, Feb 14, 2023 at 9:25 AM Miguel Angel Salazar de Troya < miguel.sala...@corintis.com> wrote: > Hello, > > I am solving the Navier-Stokes equation and an advection-diffusion > equation to model the temperature. They are fully coupled because the > viscosity is temperature dependent. I plan

[petsc-users] Segregated solvers in PETSc

2023-02-14 Thread Miguel Angel Salazar de Troya
Hello, I am solving the Navier-Stokes equation and an advection-diffusion equation to model the temperature. They are fully coupled because the viscosity is temperature dependent. I plan to solve the fully-coupled problem with a segregated approach: I first solve the Navier-Stokes equation for a

Re: [petsc-users] Problem setting Fieldsplit fields

2023-02-14 Thread Nicolas Barnafi via petsc-users
Hello Matt, After some discussions elsewhere (thanks @LMitchell!), we found out that the problem is that the fields are setup with PCSetIS, without an attached DM, which does not support this kind of nesting fields. I would like to add this feature, meaning that during the setup of the

[petsc-users] MatFDColoringSetUp with Periodic BC

2023-02-14 Thread Pierre Bernigaud
Dear all, I hope this email finds you well. We are currently working on a solver which is employing DMDA with SNES. The jacobian is computed via FDColoring, ie: call DMDACreate1D(PETSC_COMM_WORLD, DM_BOUNDARY_GHOSTED, NY, NC, NGc, PETSC_NULL_INTEGER, dmF, ierr) ! - Some steps ... call