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

More documents

Recommendations

Info

161 162 i f NHL > 0 % a c c u r a t e l y account m u l t i p r o c e s s o r s 163 NCL = 0 ; % number o f CPUs in lamnodes 164 f o r i =1: s i z e (LNO, 1 ) % add the 2nd ”:”− separated 165 f i e l d s=s p l i t (LNO( i , : ) , ’ : ’ ) ; % f i e l d , ie , #CPUs 166 NCL = NCL + str2num ( f i e l d s ( 2 , : ) ) ; 167 end 168 i f NCL 0 % avoid lamhalt in weird c a s e s 191 disp ( ’ LAM Init : h a l t i n g LAM’ ) 192 system ( ’ lamhalt ’ ) ; % won ’ t get caught on t h i s 193 end 194 end 195 end 196 197 %−−−−−−−−−−−−−−−−−−−−−−−−−−−−−− 198 % LAMBOOT 199 %−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−− 200 % r e a s o n s <strong>to</strong> lamboot : % 201 % − not lambooted yet % s t a t ˜=0 202 % − lamhalted above ( or weird ) % NHL < NSLAVES+1 (0 i s
215 f p r i n t f ( f i d , ’%s \n ’ ,h ’ ) ; % w r i t e s l a v e s ’ hostnames 216 end 217 f c l o s e ( f i d ) ; 218 disp ( ’ LAM Init : booting LAM’ ) 219 220 s t a t = system ( ’ lamboot −s −v bhost ’ ) ; 221 222 i f s t a t % again , t h i s shouldn ’ t happen 223 f p r i n t f ( ’ pushing stubborn MATLAB ” system ” c a l l ( lamboot ) : ’ ) ; 224 while s t a t 225 f p r i n t f ( ’ . ’ ) ; s t a t = system ( ’ lamboot −s −v bhost ’ ) ; 226 end 227 f p r i n t f ( ’ \n ’ ) ; 228 end 229 230 system ( ’rm −f bhost ’ ) ; % don ’ t need bhost anymore 231 end % won ’ t wipe on e x i t / could lamhalt 232 233 %−−−−−−−−−−−−−−−−−−−−−−−−−−−−−− 234 % RPI CHECK 235 %−−−−−−−−−−−−−−−−−−−−−−−−−−−−−− 236 237 [ i n f I f l g I ] = M P I I n i t i a l i z e d ; % I n i t ? 238 [ infF f l g F ] = MPI Finalized ; % F i n a l i z e ? 239 240 i f i n f I | | infF 241 e r r o r ( ’ LAM Start : e r r o r c a l l i n g I n i t i a l i z e d / F i n a l i z e d ? ’ ) 242 end 243 244 i f f l g I && ˜ f l g F % Perfect , ready <strong>to</strong> s t a r t 245 e l s e % something we could f i x ? 246 i f f l g I | | f l g F % MPI used , w i l l break 247 c l e a r MPI ∗ % u n l e s s we c l e a r MPITB 248 disp ( ’ LAM Start : MPI a l r e a d y used− c l e a r i n g ’ ) % must s t a r t over 249 end 250 251 MPI Init ; 252 end 253 254 %−−−−−−−−−−−−−−−−−−−−−−−−−−−−−− 255 % NSLAVES CHECK 256 %−−−−−−−−−−−−−−−−−−−−−−−−−−−−−− 257 258 [ i n f o a t t r f l a g ] = MPI Attr get (MPI COMM WORLD, MPI UNIVERSE SIZE ) ; 259 i f i n f o | ˜ f l a g 260 e r r o r ( ’ LAM Init : a t t r i b u t e MPI UNIVERSE SIZE does not e x i s t ? ’ ) 261 end 262 i f attr
Page 1 and 2:
UCRL-SM-212121 sundialsTB, a Matlab
Page 3:
sundialsTB, a matlab Interface to s
Page 6 and 7:
1 Introduction sundials [2], SUite
Page 8 and 9:
2 matlab Interface to cvodes The ma
Page 10 and 11:
The relative and absolute tolerance
Page 12 and 13:
(see PrecModule), it specifies the
Page 14 and 15:
to (centered or forward) finite dif
Page 16 and 17:
If quadratures were computed (see C
Page 18 and 19: o RootInfo - strucutre with rootfin
Page 20 and 21: Jacobian approximation Fields in LS
Page 22 and 23: o cntr [ true | false ] If true, CV
Page 24 and 25: 123 f o r j = 1 :N 124 y ( j , i )
Page 26 and 27: 231 box on ; 232 g r i d on ; 233 x
Page 28 and 29: 339 s e t ( hc ( 4 ) , ’XData ’
Page 30 and 31: 447 end 448 f p r i n t f ( ’−
Page 32 and 33: Purpose CVDenseJacFn - type for use
Page 34 and 35: The function GLOCFUN must be define
Page 36 and 37: evaluated on the backward phase. Se
Page 38 and 39: See also CVodeSetOptions NOTE: ODES
Page 40 and 41: M = I - gamma*J. Here J is the syst
Page 42 and 43: Adjoint Problem The function PSOLFU
Page 44 and 45: 3.1 Interface functions Purpose KIN
Page 46 and 47: Jacobian-vector product. PrecModule
Page 48 and 49: then KINSol attempts to solve the s
Page 50 and 51: 3.2 Function types Purpose KINDense
Page 52 and 53: Purpose KINJacTimesVecFn KINJacTime
Page 54 and 55: The function PSOLFUN must be define
Page 56 and 57: 4.1 nvector functions Purpose N_VDo
Page 58 and 59: 13 14 i f nargin == 1 15 16 r e t =
Page 60 and 61: N_VWL2Norm returns the weighted Euc
Page 62 and 63: 4.2 Parallel utilities Purpose mpir
Page 64 and 65: 19 f p r i n t f ( ’ waiting f o
Page 66 and 67: 53 e r r o r ( ’ LAM Init : : can
Page 70 and 71: 269 %RPI STR F u l l LAM SSI RPI s
show all

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

You also want an ePaper? Increase the reach of your titles

Delete template?

Save as template?