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

1 f u n c t i o n r e t = N VWrmsNorm( 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 / l e n g t h ( x ) ) ;
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 ln = l e n g t h ( x ) ;
28 gn = 0 ;
29 MPI Allreduce ( ln , gn , ’SUM’ ,comm ) ;
30
31 r e t = s q r t ( gnrm/gn ) ;
32
33 end
56