On Tue, Nov 18, 2008 at 1:22 PM, Rafael Santos Coelho <rafaelsantoscoelho at gmail.com> wrote: > Originally, my idea is to implement this preconditioner. As far as I`m > concerned, there`s no need to form and store any matrices, only the diagonal > elements of the jacobian like I said. This way, I believe I`d have to stick > with option #2. > > With regard to "The matrix it will receive will be the MFFD that is created > by -mat_mf.", what matrix is that actually? The one that is created by the > `MatCreateMFFD` routine?
Yes. > "Then you would have to extract the KSP and set your matrix-free PC there." > > Do the calls to PCSetType and PCShellSetApply plus my own implementation of > the `MatrixFreePreconditioner` function do that? SNESGetKSP KSPGetPC then your functions Matt > Best regards, > > Rafael -- What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead. -- Norbert Wiener
