An important use case for condition number estimation / calculating one-norms is repeated estimation of matrices with the same dimension.
To make that work efficiently, it's necessary to not constantly reallocate when calling normest1. There should thus be an object that contains all allocated temporaries, and normest1 should be redefined in terms of that struct.
An important use case for condition number estimation / calculating one-norms is repeated estimation of matrices with the same dimension.
To make that work efficiently, it's necessary to not constantly reallocate when calling
normest1
. There should thus be an object that contains all allocated temporaries, andnormest1
should be redefined in terms of that struct.