File tree Expand file tree Collapse file tree
Expand file tree Collapse file tree Original file line number Diff line number Diff line change @@ -319,25 +319,25 @@ def dotS(n,S):
319319 dS [4 ] = gdivl * (S [1 ]- 3 * S [3 ]) # d(dP/dt)/dt
320320 else :
321321 # Refer to comment block above.
322- #... compute Ainv
322+ # ... compute Ainv
323323 cosP = cos (S [3 ])
324324 Ainv = np .zeros ((2 ,2 ),dtype = np .float )
325325 Ainv [0 ,0 ] = - 1.0
326326 Ainv [0 ,1 ] = cosP + 1.0
327327 Ainv [1 ,0 ] = Ainv [0 ,1 ]
328328 Ainv [1 ,1 ] = - (2 * cosP + 3.0 )
329329 Ainv = Ainv / (cosP ** 2 - 2.0 )
330- #... compute G
330+ # ... compute G
331331 sinT = sin (S [1 ])
332332 sinP = sin (S [3 ])
333333 sinTpP = sin (S [1 ]+ S [3 ])
334334 G = np .zeros ((2 ,1 ),dtype = np .float )
335335 G [0 ,0 ] = gdivl * (sinTpP + 2 * sinT ) - sinP * S [4 ]* (S [4 ] + 2 * S [2 ])
336336 G [1 ,0 ] = gdivl * sinTpP + sinP * S [2 ]** 2
337337 G = - G
338- #... compute Y' = Ainv*G
338+ # ... compute Y' = Ainv*G
339339 AinvG = Ainv .dot (G )
340- #... load state vector values
340+ # ... load state vector values
341341 dS [1 ] = S [2 ]
342342 dS [2 ] = AinvG [0 ,0 ]
343343 dS [3 ] = S [4 ]
You can’t perform that action at this time.
0 commit comments