intercepting residues as array

Dominik Szczerba dominik at itis.ethz.ch
Sat Dec 12 07:48:00 CST 2009


Dear Barry and Jed,

Many thanks for the useful hints, the functions provide the 
functionality I needed.

Dominik

Barry Smith wrote:
>    See the manual page for KSPSetResidualHistory()
> 
>     Let us know if it does not provide the functionality you need.
> 
>    Barry
> 
> On Dec 11, 2009, at 1:36 PM, Dominik Szczerba wrote:
> 
>> Would it be easy to intercept the rnorms as computed in e.g. fgmres  
>> into an array? E.g. have
>>
>>  0 KSP Residual norm 5.955965225039e-04
>>  1 KSP Residual norm 8.064336982160e-05
>>  2 KSP Residual norm 2.550695121313e-05
>>  3 KSP Residual norm 7.872959324551e-06
>>  4 KSP Residual norm 2.084257872716e-06
>>
>> as {5.955965225039e-04, 8.064336982160e-05, ...} for independent  
>> logging etc.
>>
>> I know how I can do it the long way (e.g. writing my own convergence  
>> test where I have an access to rnorm), but I was wondering if there  
>> was a short way.
>>
>> with regards,
>> Dominik
> 



More information about the petsc-users mailing list