Re: [petsc-users] Question about PetscViewerHDF5SetTimestep usage

2018-12-14 Thread Sajid Ali via petsc-users
Hi Matthew, Thanks for the info. Could you give a specific commit that I should rebuild PETSc with ? I’m currently building with commit=83593104eb216a023c00ab8a3f9897d55e531412'. On Fri, Dec 14, 2018 at 12:56 PM Matthew Knepley wrote: > I believe this is now fixed in 'master', It was a bug in

[petsc-users] Combining petsc with libocca

2018-12-14 Thread Samuel Miller via petsc-users
I’m in the process of planning a basic fluids code based on the lattice Boltzmann method (LBM) and would like some input from the petsc user base. All the computation in LBM is local, so parallelization is straightforward and can be simply done in a threaded application; a simple google search

[petsc-users] Question about PetscViewerHDF5SetTimestep usage

2018-12-14 Thread Sajid Ali via petsc-users
Hi, I have a question about PetscViewerHDF5SetTimestep usage. I have tried the following minimal example which tries to write a vector twice to an hdf5 file and it fails with the error message shown below. Example : https://pastebin.com/cAeZsUgA Error : writing vector in hdf5 to vector.dat

Re: [petsc-users] Solution "jumps" after setting timestep

2018-12-14 Thread Zhang, Hong via petsc-users
> On Dec 14, 2018, at 10:28 AM, moritz.h...@isea.rwth-aachen.de wrote: > > Hello, > thank you this seems to be a reasonable explanation. > I will try this for my "steady" phases. > I have also phases where my Neumann boundary conditions change with time (in > my real application, here I uses

Re: [petsc-users] Solution "jumps" after setting timestep

2018-12-14 Thread Moritz.Huck--- via petsc-users
Hello, thank you this seems to be a reasonable explanation. I will try this for my "steady" phases. I have also phases where my Neumann boundary conditions change with time (in my real application, here I uses constant conditions), in this time I also need the results at exact time points. Can

Re: [petsc-users] Solution "jumps" after setting timestep

2018-12-14 Thread Zhang, Hong via petsc-users
> On Dec 14, 2018, at 5:16 AM, moritz.h...@isea.rwth-aachen.de wrote: > > Correction. > Not exactly more like this > While t<2 > t=TSGetTime > h=TSGetTimeStep > TSSetMaxTime(min(t+h,2) > TSSolve > t=TSGetTime > TSSetMaxTime(t+1e-3) > TSSolve > While t<10 >

Re: [petsc-users] Solution "jumps" after setting timestep

2018-12-14 Thread Moritz.Huck--- via petsc-users
Correction. Not exactly more like this While t<2 t=TSGetTime h=TSGetTimeStep TSSetMaxTime(min(t+h,2) TSSolve t=TSGetTime TSSetMaxTime(t+1e-3) TSSolve While t<10 t=TSGetTime h=TSGetTimeStep TSSetMaxTime(min(t+h,10) TSSolve @Hong No

Re: [petsc-users] Solution "jumps" after setting timestep

2018-12-14 Thread Moritz.Huck--- via petsc-users
In this example yes. -Ursprüngliche Nachricht- Von: Emil Constantinescu Gesendet: Donnerstag, 13. Dezember 2018 20:21 An: Zhang, Hong ; Huck, Moritz Cc: petsc-users Betreff: Re: [petsc-users] Solution "jumps" after setting timestep Okay, so you start your integration with