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
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
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
> 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
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
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