Re: [petsc-users] low dimensional sub-problem

Great, thanks for letting us know. Barry > On Feb 25, 2017, at 6:24 PM, Gideon Simpson wrote: > > Yea, I implemented it as suggested by Matt and Barry with the BCast of the > element that sits only on process 0, and that seems to work fine. > > -gideon > >>

Re: [petsc-users] Multigrid with defect correction

> On Feb 25, 2017, at 2:21 PM, Matt Landreman wrote: > > Thanks everyone for the help. > > On item 1 (using Amat different from Pmat with geometric multigrid), I tried > Barry's suggestion but it did not seem to resolve the issue. For example, in > ksp ex25.c, I

Re: [petsc-users] low dimensional sub-problem

Yea, I implemented it as suggested by Matt and Barry with the BCast of the element that sits only on process 0, and that seems to work fine. -gideon > On Feb 25, 2017, at 5:57 PM, Barry Smith wrote: > > >> On Feb 25, 2017, at 11:00 AM, Matthew Knepley

Re: [petsc-users] Multigrid with defect correction

Matt Landreman writes: > On Jed's comment, the application I have in mind is indeed a > convection-dominated equation (a steady linear 3D convection-diffusion > equation with smoothly varying anisotropic coefficients and recirculating > convection). Gamg and

Re: [petsc-users] low dimensional sub-problem

> On Feb 25, 2017, at 11:00 AM, Matthew Knepley wrote: > > On Sat, Feb 25, 2017 at 10:12 AM, Gideon Simpson > wrote: > > -gideon > >> On Feb 25, 2017, at 10:44 AM, Jed Brown wrote: >> >> Gideon Simpson

Re: [petsc-users] Multigrid with defect correction

On Sat, Feb 25, 2017 at 2:21 PM, Matt Landreman wrote: > Thanks everyone for the help. > > On item 1 (using Amat different from Pmat with geometric multigrid), I > tried Barry's suggestion but it did not seem to resolve the issue. For > example, in ksp ex25.c, I tried

Re: [petsc-users] Multigrid with defect correction

Thanks everyone for the help. On item 1 (using Amat different from Pmat with geometric multigrid), I tried Barry's suggestion but it did not seem to resolve the issue. For example, in ksp ex25.c, I tried adding the following lines after line 112: if (J == jac) { ierr =

Re: [petsc-users] low dimensional sub-problem

On Sat, Feb 25, 2017 at 10:12 AM, Gideon Simpson wrote: > > -gideon > > On Feb 25, 2017, at 10:44 AM, Jed Brown wrote: > > Gideon Simpson writes: > > I’ve been continuing working on implementing a projection method problem,

Re: [petsc-users] low dimensional sub-problem

-gideon > On Feb 25, 2017, at 10:44 AM, Jed Brown wrote: > > Gideon Simpson writes: > >> I’ve been continuing working on implementing a projection method problem, >> which, loosely, looks like the following. Vec y contains the state vector >>

Re: [petsc-users] low dimensional sub-problem

Gideon Simpson writes: > I’ve been continuing working on implementing a projection method problem, > which, loosely, looks like the following. Vec y contains the state vector > for my system, y’ = f(y) which is solved with a TS, using, for now, rk4. > > I have added

[petsc-users] low dimensional sub-problem

I’ve been continuing working on implementing a projection method problem, which, loosely, looks like the following. Vec y contains the state vector for my system, y’ = f(y) which is solved with a TS, using, for now, rk4. I have added to this a TSPostStep which does root finding on the