User:Tom Allen/TensorExperimentsZ14
Jump to navigation
Jump to search
...
... E1cn (geodesics) ...
...
...
... one sign ...
...
load 'plot' pd ((2.05p1*cos@cx0"1);2.05p1*sin@cx0"1) pTtotal (_ 0 2p1 160,_ 0 0 0,_ 0 0 0,:_ 0 0 0"_)'' 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 =: 2.05p1 , 0.5p1 , 1p1 , 0.5p1 "_ Initdx=: (0 , 0 , _0.01 , 0.4809"_)*%@(12"_) vS00xL =:Initx vS01dxL =:0,0,0,0: vS02d1L =:Initd1 vS06ds =:Initds vS07count=:(0,777"_)*12"_ pd vSgDplot pdDA=:vSgD@vS'' pd ((2.20p1*cos@cx0"1);2.20p1*sin@cx0"1) pTtotal (_ 0 2p1 160,_ 0 0 0,_ 0 0 0,:_ 0 0 0"_)'' 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 =: 2.20p1 , 0.5p1 , 1p1 , 0.5p1 "_ Initdx=: (0 , 0 , _0.01 , 0.5022"_)*%@(12"_) vS00xL =:Initx vS01dxL =:0,0,0,0: vS02d1L =:Initd1 vS06ds =:Initds vS07count=:(0,669"_)*12"_ pd vSgDplot pdDB=:vSgD@vS'' pd 'aspect 1' pd 'pdf' pd 'show'
NB. ... execute (ijx) ... $pdDA 9325 8 4 $pdDB 8029 8 4 NB. ... x1 ... *./ (0=]) (0.50p1-]) cx1"1@:((S00xL{])"2) pdDA 1 *./ (0=]) (0.50p1-]) cx1"1@:((S00xL{])"2) pdDB 1 NB. ... x2 decreases continuously ... *./ (0<]) (2-/\]) cx2"1@:((S00xL{])"2) pdDA 1 *./ (0<]) (2-/\]) cx2"1@:((S00xL{])"2) pdDB 1 NB. ... x3 ... NB. ... increases continuously ... *./ (0>]) (2-/\]) cx3"1@:((S00xL{])"2) pdDA 1 *./ (0>]) (2-/\]) cx3"1@:((S00xL{])"2) pdDB 1 NB. ... min/max interval ... (<./@],>./@])@(2-~/\]) cx3"1@:((S00xL{])"2) pdDA 0.0392337 0.0412501 (<./@],>./@])@(2-~/\]) cx3"1@:((S00xL{])"2) pdDB 0.0405035 0.0423874 NB. ... last value ... {: cx3"1@:((S00xL{])"2) pdDA 376.955 {: cx3"1@:((S00xL{])"2) pdDB 334.555
...
... also ...
...
load 'plot' NB. ... x2 and x3 ... pd (-@cx2"1;cx3"1)@:((S00xL{])"2) pdDA NB. blue pd (-@cx2"1;cx3"1)@:((S00xL{])"2) pdDB NB. red pd 'pdf' pd 'show'
...
... the other sign ...
...
load 'plot' pd ((2.55p1*cos@cx0"1);2.55p1*sin@cx0"1) pTtotal (_ 0 2p1 160,_ 0 0 0,_ 0 0 0,:_ 0 0 0"_)'' 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 =: 2.55p1 , 0.5p1 , 1p1 , 0.5p1 "_ Initdx=: (0 , 0 , _0.01 , 0.5471"_)*%@(12"_) vS00xL =:Initx vS01dxL =:0,0,0,0: vS02d1L =:Initd1 vS06ds =:Initds vS07count=:(0,742"_)*12"_ pd vSgDplot pdDC=:vSgD@vS'' pd ((2.70p1*cos@cx0"1);2.70p1*sin@cx0"1) pTtotal (_ 0 2p1 160,_ 0 0 0,_ 0 0 0,:_ 0 0 0"_)'' 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 =: 2.70p1 , 0.5p1 , 1p1 , 0.5p1 "_ Initdx=: (0 , 0 , _0.01 , 0.5410"_)*%@(12"_) vS00xL =:Initx vS01dxL =:0,0,0,0: vS02d1L =:Initd1 vS06ds =:Initds vS07count=:(0,658"_)*12"_ pd vSgDplot pdDD=:vSgD@vS'' pd 'aspect 1' pd 'pdf' pd 'show'
NB. ... execute (ijx) ... $pdDC 8905 8 4 $pdDD 7897 8 4 NB. ... x1 ... *./ (0=]) (0.50p1-]) cx1"1@:((S00xL{])"2) pdDC 1 *./ (0=]) (0.50p1-]) cx1"1@:((S00xL{])"2) pdDD 1 NB. ... x2 decreases continuously ... *./ (0<]) (2-/\]) cx2"1@:((S00xL{])"2) pdDC 1 *./ (0<]) (2-/\]) cx2"1@:((S00xL{])"2) pdDD 1 NB. ... x3 ... NB. ... increases continuously ... *./ (0>]) (2-/\]) cx3"1@:((S00xL{])"2) pdDC 1 *./ (0>]) (2-/\]) cx3"1@:((S00xL{])"2) pdDD 1 NB. ... min/max interval ... (<./@],>./@])@(2-~/\]) cx3"1@:((S00xL{])"2) pdDC 0.0455584 0.0474994 (<./@],>./@])@(2-~/\]) cx3"1@:((S00xL{])"2) pdDD 0.0450811 0.0461152 NB. ... last value ... {: cx3"1@:((S00xL{])"2) pdDC 415.777 {: cx3"1@:((S00xL{])"2) pdDD 361.245
...
... also ...
...
load 'plot' NB. ... x2 and x3 ... pd (-@cx2"1;cx3"1)@:((S00xL{])"2) pdDC NB. blue pd (-@cx2"1;cx3"1)@:((S00xL{])"2) pdDD NB. red pd 'pdf' pd 'show'
Download MoinMoin source: File:TensorExperimentsZ14.ijs