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'
... waves ...
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'
... x2 and x3 ...
...
... 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'
... waves ...
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'
... x2 and x3 ...
Download MoinMoin source: File:TensorExperimentsZ14.ijs