[petsc-dev] Adding VecErrorNormWRMS()

[Jed Brown]

Lisandro Dalcin <dalcinl at gmail.com> writes:

> Jed, I would like to make the following changes:
>
> *Add PetscErrorCode VecErrorNormWRMS(Vec U,Vec Y,PetscReal atol,Vec
> vatol,PetscReal rtol,Vec vrtol,PetscReal *norm).
>
> Any alternative name to suggest?

Fine, but how would you deal with Shri's ts->is_diff?  And Shri, why did
you implement that by creating these index sets instead of setting those
entries of vatol to +Infinity?  That's why we have vector weighted norms
after all.

> * Implement TSErrorNormWRMS() in terms of VecErrorNormWRMS().
>
> PS: These changes are preliminary ones, next I want to suggest a
> mechanism to let TS subtypes to decide how to compute the LTE error
> (or approximation to it).

Cool; the current solution is more complicated and less flexible than
necessary, IMO.
-------------- next part --------------
A non-text attachment was scrubbed...
Name: signature.asc
Type: application/pgp-signature
Size: 818 bytes
Desc: not available
URL: <http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20150430/d5750793/attachment.sig>