sundialsTB, a Matlab Interface to SUNDIALS - Site Index Page ...

N_VWL2Norm returns the weighted Euclidean L2 norm of x
with weight vector w:
sqrt [(sum (i = 0 to N-1) (x[i]*w[i])^2)]
Usage: RET = N_VWL2Norm ( X, W [, COMM] )
If COMM is not present, N_VWL2Norm returns the weighted L2
norm of the local portion of X. Otherwise, it returns the
global weighted L2 norm..
Source Code
1 f u n c t i o n r e t = N VWL2Norm( x ,w,comm)
11
12 % Radu Serban
13 % Copyright ( c ) 2005 , The Regents o f the U n i v e r s i t y o f C a l i f o r n i a .
14 % $Revision$Date$
15
16 i f nargin == 2
17
18 r e t = dot ( x . ˆ 2 ,w. ˆ 2 ) ;
19 r e t = s q r t ( r e t ) ;
20
21 e l s e
22
23 lnrm = dot ( x . ˆ 2 ,w. ˆ 2 ) ;
24 gnrm = 0 . 0 ;
25 MPI Allreduce ( lnrm , gnrm , ’SUM’ ,comm ) ;
26
27 r e t = s q r t ( gnrm ) ;
28
29 end
Purpose
N_VWrmsNorm
N_VWrmsNorm returns the weighted root mean square norm of x
Synopsis
function ret = N_VWrmsNorm(x,w,comm)
Description
N_VWrmsNorm returns the weighted root mean square norm of x
with weight vector w:
sqrt [(sum (i = 0 to N-1) (x[i]*w[i])^2)/N]
Usage: RET = N_VWrmsNorm ( X, W [, COMM] )
If COMM is not present, N_VWrmsNorm returns the WRMS norm
of the local portion of X. Otherwise, it returns the global
WRMS norm..
Source Code
55