Re: [petsc-users] SNESSetFunctionDomainError

2014-08-31 Thread Barry Smith
We don’t currently have the concept of “matrix-vector product failed due to funky thing like domain error” therefor end the current linear solve and also the nonlinear solve but without an error condition. I’m open to any suggestions, but not comfortable with reorganizing PETSc errors

Re: [petsc-users] SNESSetFunctionDomainError

2014-08-31 Thread Jed Brown
Derek Gaston fried...@gmail.com writes: As a workaround I have tried returning various diverged statuses in my linear and nonlinear convergence checks... but it still has a pesky problem where there are a couple more residual evaluations at the end of the nonlinear solve. I think this is

Re: [petsc-users] SNESSetFunctionDomainError

2014-08-31 Thread Dmitry Karpeyev
Handling this at the KSP level (I actually think the Mat level is more appropriate, since the operator, not the solver, knows its domain), we have to fix MatMFFD, I would think. Yes, that thing has the concept of a nonlinear function, but it doesn't have any clue about or access to a SNES. I did,

Re: [petsc-users] SNESSetFunctionDomainError

2014-08-31 Thread Jed Brown
Dmitry Karpeyev karp...@mcs.anl.gov writes: Handling this at the KSP level (I actually think the Mat level is more appropriate, since the operator, not the solver, knows its domain), We are dynamically discovering the domain, but I don't think it's appropriate for Mat to refuse to evaluate