User:Tom Allen/TensorExperimentsZ09
Jump to navigation
Jump to search
...
... E1cn (geodesics) ...
...
...
... one sign ...
...
load 'plot' vSgDd2M=.[([(0{])@([(xdz smx"2 1 E1cn)zRef)hkxM),:@(S03xH{]) vSgDd2 =.gXZ(-@(([([(0{])@ch2k xRef@hkxM),:@(S03xH{]))+/"1@((<0 2)|:])@((0|:[)smx])(]*/])@(S02d1L{]))+vSgDd2M)] Initds2=.(gXZ([(0{])@mcv xRef@hkxM),:@Initx)+/^:2@:*(]*/])@Initdx gXcw =: 1 0 _ _ _ _"_ gXbAll=: _1 1 1 1 1 0"_ Initx =: 2p1 , 0.5p1 , 1p1 , 0.5p1"_ Initdx=: 0 , 0 , _0.01 , 0 "_ vS00xL =:Initx vS01dxL =:0,0,0,0: vS02d1L =:Initd1 vS06ds =:Initds vS07count=:0,628"_ pd vSgDplot pdDA=:vSgD@vS'' Initx =: 3p1 , 0.5p1 , 1p1 , 0.5p1"_ Initdx=: 0 , 0 , _0.01 , 0 "_ vS00xL =:Initx vS01dxL =:0,0,0,0: vS02d1L =:Initd1 vS06ds =:Initds vS07count=:0,628"_ pd vSgDplot pdDB=:vSgD@vS'' Initx =: 4p1 , 0.5p1 , 1p1 , 0.5p1"_ Initdx=: 0 , 0 , _0.01 , 0 "_ vS00xL =:Initx vS01dxL =:0,0,0,0: vS02d1L =:Initd1 vS06ds =:Initds vS07count=:0,628"_ pd vSgDplot pdDC=:vSgD@vS'' pd 'aspect 1' pd 'pdf' pd 'show'
NB. ... execute (ijx) ... NB. ... for pdDA ... $pdDA 629 8 4 NB. ... x0 ... *./ (0=]) (2p1-]) cx0"1@:((S00xL{])"2) pdDA 1 NB. ... x1 ... *./ (0=]) (0.50p1-]) cx1"1@:((S00xL{])"2) pdDA 1 NB. ... x2 decreases continuously ... *./ (0<]) (2-/\]) cx2"1@:((S00xL{])"2) pdDA 1 NB. ... x3 ... *./ (0=]) (0.50p1-]) cx3"1@:((S00xL{])"2) pdDA 1 NB. ... for pdDB ... $pdDB 629 8 4 NB. ... etc ...
Download MoinMoin source: File:TensorExperimentsZ09.ijs