[petsc-dev] Semismooth Newton

2012-08-02 Thread Todd Munson
> How does one precondition a system like that? Presumably, the user may at > best be able to provide a preconditioner for JF. > > If JF responds well to multigrid, then it may be possible to use AMG-KKT > to coarsen the constraint Hessian, but how does MG perform for the sum > JF + Hc? There i

[petsc-dev] Semismooth Newton

2012-08-02 Thread Dmitry Karpeev
On Wed, Aug 1, 2012 at 5:01 PM, Todd Munson wrote: > > On Aug 1, 2012, at 4:13 PM, Jed Brown wrote: > > > On Wed, Aug 1, 2012 at 11:10 AM, Todd Munson > wrote: > > > > Yes, I am all in favor of actually specifying a VI externally and then > > internally formulating the appropriate augmented syst

[petsc-dev] Semismooth Newton

2012-08-01 Thread Todd Munson
On Aug 1, 2012, at 4:13 PM, Jed Brown wrote: > On Wed, Aug 1, 2012 at 11:10 AM, Todd Munson wrote: > > Yes, I am all in favor of actually specifying a VI externally and then > internally formulating the appropriate augmented system. This is what > I do in other contexts (i.e. PATH for GAMS/AMP

[petsc-dev] Semismooth Newton

2012-08-01 Thread Jed Brown
On Wed, Aug 1, 2012 at 11:10 AM, Todd Munson wrote: > > Yes, I am all in favor of actually specifying a VI externally and then > internally formulating the appropriate augmented system. This is what > I do in other contexts (i.e. PATH for GAMS/AMPL). Some things you > probably want to do with s

[petsc-dev] Semismooth Newton

2012-08-01 Thread Todd Munson
Yes, I am all in favor of actually specifying a VI externally and then internally formulating the appropriate augmented system. This is what I do in other contexts (i.e. PATH for GAMS/AMPL). Some things you probably want to do with such an interface is differentiate between linear and nonlinea

[petsc-dev] Semismooth Newton

2012-08-01 Thread Jed Brown
On Tue, Jul 31, 2012 at 5:38 AM, Todd Munson wrote: > Knowing the constraint structure and forming the augmented system is fine > and can be > helpful information to know when determining the active set. I would > concentrate on > linear constraints and inequalities though; they are nicer to dea

[petsc-dev] Semismooth Newton

2012-07-31 Thread Todd Munson
On Jul 31, 2012, at 2:36 AM, Jed Brown wrote: > On Mon, Jul 30, 2012 at 5:40 PM, Todd Munson wrote: > > 1. change SNESVI to support general complementarity constraints because > > only doing box constraints on state variables is lame > > Okay; and how do you want to do this? You need to write

[petsc-dev] Semismooth Newton

2012-07-31 Thread Jed Brown
On Mon, Jul 30, 2012 at 5:40 PM, Todd Munson wrote: > > 1. change SNESVI to support general complementarity constraints because > only doing box constraints on state variables is lame > > Okay; and how do you want to do this? You need to write out the > corresponding KT system > and then you jus

[petsc-dev] Semismooth Newton

2012-07-30 Thread Todd Munson
> 1. change SNESVI to support general complementarity constraints because only > doing box constraints on state variables is lame Okay; and how do you want to do this? You need to write out the corresponding KT system and then you just have a box constrained problem. One can make it easier to

[petsc-dev] Semismooth Newton

2012-07-30 Thread Jed Brown
I think we should: 1. change SNESVI to support general complementarity constraints because only doing box constraints on state variables is lame 2. have a mode to solve the system in the full space, instead of eliminating the semi-smooth variables (which produces the bad conditioning) 3. eventua