User:Tom Allen/TensorExperimentsB11
Jump to navigation
Jump to search
...
...
...
... 6.25/1.75/8/2 ...
load 'plot'
vSgDd2M=.[([(0{])@(([mcn xRef)smx"2 1[(ydx((0|:[)smx])"2 1 B2kCEr1+SBkCEr1)yRef)hkxI),:@(S03xH{])
vSgDd2 =.gXB(-@(([([(0{])@ch2k xRef@hkxI),:@(S03xH{]))+/"1@((<0 2)|:])@((0|:[)smx])(]*/])@(S02d1L{]))+vSgDd2M)]
Initds2=.(gXB([(0{])@mcv xRef@hkxI),:@Initx)+/^:2@:*(]*/])@Initdx
gXcw =: 1.5 0 _ _ _"_
gXpAll =: 1 1 1 9 _1"_
gXq2All=: 1 1 1 9 _"_
gXq3All=: 0 0 1 1 0"_
gXq4All=: 1 1 1 1 _"_
gXq5All=: 2 0 1 1 _2"_
gXq6All=: 1 1 1 0 _"_
gXq7All=: 1 1 3 1 _2"_
Initx =: 6.25p1 , 0.5p1 , 1p1 , 1.75p1"_
Initdx=: (0 , 0 , _0.01 , 0 "_)*%@(2"_)
vS00xL =:Initx
vS01dxL =:0,0,0,0:
vS02d1L =:Initd1
vS06ds =:Initds
vS07count =:(0,8000"_)*2"_
pd vSgDplot pdDA=:vSgD@vS''
pd 'pdf'
pd 'show'
... geodesic ...
NB. ... execute (ijx) ...
NB. ... using pdDA from the plot above ...
$pdDA
16001 8 4
NB. ... x0 ...
A=:cx0"1@:((S00xL{])"2) pdDA
(6.25p1-<./) A
0.0674185
(>./-6.25p1"_) A
0.183833
NB. ... x1 ...
*./ (0=]) ((2^_48)qtsz]) (0.50p1-]) cx1"1@:((S00xL{])"2) pdDA
1
NB. ... x2 decreases continuously ...
*./ (0<]) (2-/\]) cx2"1@:((S00xL{])"2) pdDA
1
NB. ... x3 ...
*./ (0=]) ((2^_46)qtsz]) (1.75p1-]) cx3"1@:((S00xL{])"2) pdDA
1
load 'plot'
NB. ... using pdDA from the plot above ...
pd (-@cx2"1;cx0"1)@:((S00xL{])"2) (i.@(16000"_){]) pdDA
pd 'pdf'
pd 'show'
... x0 ...
...
...
...
... 6.25/1.75/8/3 ...
load 'plot'
vSgDd2M=.[([(0{])@(([mcn xRef)smx"2 1[(ydx((0|:[)smx])"2 1 B2kCEr1+SBkCEr1)yRef)hkxI),:@(S03xH{])
vSgDd2 =.gXB(-@(([([(0{])@ch2k xRef@hkxI),:@(S03xH{]))+/"1@((<0 2)|:])@((0|:[)smx])(]*/])@(S02d1L{]))+vSgDd2M)]
Initds2=.(gXB([(0{])@mcv xRef@hkxI),:@Initx)+/^:2@:*(]*/])@Initdx
gXcw =: 1.5 0 _ _ _"_
gXpAll =: 1 1 1 9 _1"_
gXq2All=: 1 1 1 9 _"_
gXq3All=: 0 0 1 1 0"_
gXq4All=: 1 1 1 1 _"_
gXq5All=: 2 0 1 1 _2"_
gXq6All=: 1 1 1 0 _"_
gXq7All=: 1 1 3 1 _2"_
Initx =: 6.25p1 , 0.5p1 , 1p1 , 1.75p1"_
Initdx=: (0 , 0 , _0.01 , 0 "_)*%@(3"_)
vS00xL =:Initx
vS01dxL =:0,0,0,0:
vS02d1L =:Initd1
vS06ds =:Initds
vS07count =:(0,8000"_)*3"_
pd vSgDplot pdDA=:vSgD@vS''
pd 'pdf'
pd 'show'
... geodesic ...
NB. ... execute (ijx) ...
NB. ... using pdDA from the plot above ...
$pdDA
24001 8 4
NB. ... x0 ...
A=:cx0"1@:((S00xL{])"2) pdDA
(6.25p1-<./) A
0.0392048
(>./-6.25p1"_) A
0.154201
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=]) ((2^_47)qtsz]) (1.75p1-]) cx3"1@:((S00xL{])"2) pdDA
1
load 'plot'
NB. ... using pdDA from the plot above ...
pd (-@cx2"1;cx0"1)@:((S00xL{])"2) (i.@(24000"_){]) pdDA
pd 'pdf'
pd 'show'
... x0 ...
...
...
...
... 6.25/1.75/8/6 ...
load 'plot'
vSgDd2M=.[([(0{])@(([mcn xRef)smx"2 1[(ydx((0|:[)smx])"2 1 B2kCEr1+SBkCEr1)yRef)hkxI),:@(S03xH{])
vSgDd2 =.gXB(-@(([([(0{])@ch2k xRef@hkxI),:@(S03xH{]))+/"1@((<0 2)|:])@((0|:[)smx])(]*/])@(S02d1L{]))+vSgDd2M)]
Initds2=.(gXB([(0{])@mcv xRef@hkxI),:@Initx)+/^:2@:*(]*/])@Initdx
gXcw =: 1.5 0 _ _ _"_
gXpAll =: 1 1 1 9 _1"_
gXq2All=: 1 1 1 9 _"_
gXq3All=: 0 0 1 1 0"_
gXq4All=: 1 1 1 1 _"_
gXq5All=: 2 0 1 1 _2"_
gXq6All=: 1 1 1 0 _"_
gXq7All=: 1 1 3 1 _2"_
Initx =: 6.25p1 , 0.5p1 , 1p1 , 1.75p1"_
Initdx=: (0 , 0 , _0.01 , 0 "_)*%@(6"_)
vS00xL =:Initx
vS01dxL =:0,0,0,0:
vS02d1L =:Initd1
vS06ds =:Initds
vS07count =:(0,8000"_)*6"_
pd vSgDplot pdDA=:vSgD@vS''
pd 'pdf'
pd 'show'
... geodesic ...
NB. ... execute (ijx) ...
NB. ... using pdDA from the plot above ...
$pdDA
48001 8 4
NB. ... x0 ...
A=:cx0"1@:((S00xL{])"2) pdDA
(6.25p1-<./) A
0.0170535
(>./-6.25p1"_) A
0.131537
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=]) ((2^_48)qtsz]) (1.75p1-]) cx3"1@:((S00xL{])"2) pdDA
1
load 'plot'
NB. ... using pdDA from the plot above ...
pd (-@cx2"1;cx0"1)@:((S00xL{])"2) (i.@(48000"_){]) pdDA
pd 'pdf'
pd 'show'
... x0 ...
...
...
...
... 6.25/1.75/8/8 ...
load 'plot'
vSgDd2M=.[([(0{])@(([mcn xRef)smx"2 1[(ydx((0|:[)smx])"2 1 B2kCEr1+SBkCEr1)yRef)hkxI),:@(S03xH{])
vSgDd2 =.gXB(-@(([([(0{])@ch2k xRef@hkxI),:@(S03xH{]))+/"1@((<0 2)|:])@((0|:[)smx])(]*/])@(S02d1L{]))+vSgDd2M)]
Initds2=.(gXB([(0{])@mcv xRef@hkxI),:@Initx)+/^:2@:*(]*/])@Initdx
gXcw =: 1.5 0 _ _ _"_
gXpAll =: 1 1 1 9 _1"_
gXq2All=: 1 1 1 9 _"_
gXq3All=: 0 0 1 1 0"_
gXq4All=: 1 1 1 1 _"_
gXq5All=: 2 0 1 1 _2"_
gXq6All=: 1 1 1 0 _"_
gXq7All=: 1 1 3 1 _2"_
Initx =: 6.25p1 , 0.5p1 , 1p1 , 1.75p1"_
Initdx=: (0 , 0 , _0.01 , 0 "_)*%@(8"_)
vS00xL =:Initx
vS01dxL =:0,0,0,0:
vS02d1L =:Initd1
vS06ds =:Initds
vS07count =:(0,8000"_)*8"_
pd vSgDplot pdDA=:vSgD@vS''
pd 'pdf'
pd 'show'
... geodesic ...
NB. ... execute (ijx) ...
NB. ... using pdDA from the plot above ...
$pdDA
64001 8 4
NB. ... x0 ...
A=:cx0"1@:((S00xL{])"2) pdDA
(6.25p1-<./) A
0.0123629
(>./-6.25p1"_) A
0.12676
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=]) ((2^_48)qtsz]) (1.75p1-]) cx3"1@:((S00xL{])"2) pdDA
1
load 'plot'
NB. ... using pdDA from the plot above ...
pd (-@cx2"1;cx0"1)@:((S00xL{])"2) (i.@(64000"_){]) pdDA
pd 'pdf'
pd 'show'
... x0 ...
Download MoinMoin source: File:TensorExperimentsB11.ijs