[petsc-dev] vector inner products
Munson, Todd
tmunson at mcs.anl.gov
Thu Apr 12 12:21:04 CDT 2018
There is a bit of code in TAO that allows the user to change the norm to
a matrix norm. This was introduced to get some mesh independent
behavior in one example (tao/examples/tutorials/ex3.c). That
norm, however, does not propagate down into the KSP methods
and is only used for testing convergence of the nonlinear
problem.
A few questions then: Is similar functionality needed in SNES? Are
TAO and SNES even the right place for this functionality? Should
it belong to the Vector class so that you can change the inner
products and have all the KSP methods (hopefully) work
correctly?
Note: that this discussion brings us to the brink of supporting an
optimize-then-discretize approach. I am not convinced we should
go down that rabbit hole.
Thanks, Todd.
More information about the petsc-dev
mailing list