Groups | Search | Server Info | Keyboard shortcuts | Login | Register [http] [https] [nntp] [nntps]

Groups > comp.soft-sys.math.mathematica > #16429 > unrolled thread

compensated summation in norm

Back to article view | Back to comp.soft-sys.math.mathematica

Contents

  compensated summation in norm Alex Krasnov <akrasnov@cory.eecs.berkeley.edu> - 2014-01-18 10:01 +0000

#16429 — compensated summation in norm

The documentation states that Norm[v, p] is equivalent to 
Total[Abs[v]^p]^(1/p). Total has a documented option Method -> 
Automatic|"CompensatedSummation". Norm has an undocumented option Method. 
Is this option supposed to correspond to the option of Total? If so, then 
it appears to have no effect. Furthermore, Norm is not consistent with the 
use of Total, having much larger numerical error:

In:	Norm[ConstantArray[0.1, 10^6], 1, Method -> Automatic]-10^5
Out:	2.2354652173817158*^-7

In:	Norm[ConstantArray[0.1, 10^6], 1, Method -> "CompensatedSummation"]-10^5
Out:	2.2354652173817158*^-7

In:	Total[ConstantArray[0.1, 10^6], Method -> Automatic]-10^5
Out:	8.731149137020111*^-10

In:	Total[ConstantArray[0.1, 10^6], Method -> "CompensatedSummation"]-10^5
Out:	0.

This behavior should be made consistent, and the documentation should be 
corrected.

Alex

[toc] | [standalone]

Back to top | Article view | comp.soft-sys.math.mathematica

csiph-web