Re: [petsc-dev] obscure changes in TSGetStages_Theta

2021-01-24 Thread Matthew Knepley
On Sun, Jan 24, 2021 at 11:48 AM Zhang, Hong via petsc-dev < petsc-dev@mcs.anl.gov> wrote: > Some TS methods such as TSRK do have an array of vectors like this to > store the stage values. But not all TS methods have it. I am fine adding > the scratch for TSTheta and any other method missing it. A

Re: [petsc-dev] Memory problem with OpenMP and Fieldsplit sub solvers

2021-01-24 Thread Mark Adams
Thanks Pierre, that works. On Sun, Jan 24, 2021 at 11:42 AM Pierre Jolivet wrote: > > > On 24 Jan 2021, at 4:54 PM, Mark Adams wrote: > > Hi Sherry, I have this running with OMP, with cuSparse solves (PETSc CPU > factorizations) > > Building SuperLU_dist w/o _OPENMP was not easy for me. > > > E

Re: [petsc-dev] obscure changes in TSGetStages_Theta

2021-01-24 Thread Zhang, Hong via petsc-dev
Some TS methods such as TSRK do have an array of vectors like this to store the stage values. But not all TS methods have it. I am fine adding the scratch for TSTheta and any other method missing it. A little drawback is that it is used only by TSGetStages and the TSStep implementation does not

Re: [petsc-dev] Memory problem with OpenMP and Fieldsplit sub solvers

2021-01-24 Thread Pierre Jolivet
> On 24 Jan 2021, at 4:54 PM, Mark Adams wrote: > > Hi Sherry, I have this running with OMP, with cuSparse solves (PETSc CPU > factorizations) > > Building SuperLU_dist w/o _OPENMP was not easy for me. Expanding on Sherry’s answer on Thu Jan 21, it should be as easy as adding to your config

Re: [petsc-dev] Memory problem with OpenMP and Fieldsplit sub solvers

2021-01-24 Thread Mark Adams
Hi Sherry, I have this running with OMP, with cuSparse solves (PETSc CPU factorizations) Building SuperLU_dist w/o _OPENMP was not easy for me. We need to get a better way to do this. (Satish or Barry?) SuperLU works with one thread and two subdomains. With two threads I see this (appended). So t

Re: [petsc-dev] obscure changes in TSGetStages_Theta

2021-01-24 Thread Lisandro Dalcin
On Sun, 24 Jan 2021 at 10:08, Stefano Zampini wrote: > You actually do not need a RestoreStages if you use a scratch > > Vec stages[2]; > > in TS_theta, and pass it back with TSGetStages_Theta. I understand PETSc > philosophy is that every Get should have a Restore method, but is this > really ne