[petsc-dev] Adding VecErrorNormWRMS()

Abhyankar, Shrirang G. abhyshr at anl.gov
Fri May 1 07:51:38 CDT 2015

> On Apr 30, 2015, at 2:34 PM, "Jed Brown" <jed at jedbrown.org> wrote:
> 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?  
Thanks Jed, I didn't think about that. It certainly seems cleaner to just set the entries of vatol to infinity rather than creating index sets. 

> That's why we have vector weighted norms
> after all.
>> * Implement TSErrorNormWRMS() in terms of VecErrorNormWRMS().

Please note that there is also a variation that uses infinity norm for error calculation.

>> 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.


More information about the petsc-dev mailing list