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

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

compensated summation in norm

Show all headers | View raw

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

Back to comp.soft-sys.math.mathematica | Previous | Next | Find similar | Unroll thread

Thread

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

csiph-web