[petsc-dev] Semi-smooth VI solver and matrix-free SNES

2014-08-01 Thread Dmitry Karpeyev
Currently the semi-smooth VI solver (VINEWTONSSLS) doesn't play well with the matrix-free SNES (aka JFNK). This is partly because the calculation of the merit function gradient requires the action of the transpose of the Jacobian, which isn't implemented for MATMFFD. Presumably, the merit

Re: [petsc-dev] Semi-smooth VI solver and matrix-free SNES

2014-08-01 Thread Dmitry Karpeyev
Todd, What's the semi-smooth version in TAO that respect the bounds at the iterates? Maybe that would be a good place to attempt a unification with SNESVI. The current implementation could benefit from refactoring to enable more composability: could we use a VIRS method inside VISS and allow for

Re: [petsc-dev] Semi-smooth VI solver and matrix-free SNES

2014-08-01 Thread Jason Sarich
The ssfls and asfls solvers respect feasibility, ssils and asils do not (the 'f' stands for feasibility, 'i' for infeasibility) ssfls - semismooth, feasible, with line search asfls - active set semismooth, feasible, with line search On Fri, Aug 1, 2014 at 10:45 AM, Karpeyev, Dmitry