User:Tom Allen/TensorExperimentsB10

From J Wiki
Jump to: navigation, search


...

...

...

... 6.25/1.75/8/4 ...


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     "_)*%@(4"_)

vS00xL    =:Initx
vS01dxL   =:0,0,0,0:
vS02d1L   =:Initd1
vS06ds    =:Initds
vS07count =:(0,8000"_)*4"_
pd vSgDplot pdDA=:vSgD@vS''

pd 'pdf'
pd 'show'


Egraphic100g.jpg
... geodesic ...


NB. ... execute (ijx) ...

NB. ... using pdDA from the plot above ...

   $pdDA
32001 8 4

NB. ... x0 ...

   A=:cx0"1@:((S00xL{])"2) pdDA

   (6.25p1-<./) A
0.0274124

   (>./-6.25p1"_) A
0.142113

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.@(20000"_){]) pdDA

pd 'pdf'
pd 'show'


Egraphic1000.jpg
... x0 ...


...

...

...

... 6.375/1.75/8/4 ...


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.375p1 , 0.5p1 ,  1p1  , 1.75p1"_
Initdx=:  (0       , 0     , _0.01 , 0     "_)*%@(4"_)

vS00xL    =:Initx
vS01dxL   =:0,0,0,0:
vS02d1L   =:Initd1
vS06ds    =:Initds
vS07count =:(0,8000"_)*4"_
pd vSgDplot pdDA=:vSgD@vS''

pd 'pdf'
pd 'show'


Egraphic104g.jpg
... geodesic ...


NB. ... execute (ijx) ...

NB. ... using pdDA from the plot above ...

   $pdDA
32001 8 4

NB. ... x0 ...

   A=:cx0"1@:((S00xL{])"2) pdDA

   (6.25p1-<./) A
0.414795

   (>./-6.25p1"_) A
0.583729

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.@(27600"_){]) pdDA

pd 'pdf'
pd 'show'


Egraphic1040.jpg
... x0 ...


...

...

...

... 6.25/2.00/8/4 ...


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  , 2.00p1"_
Initdx=:  (0      , 0     , _0.01 , 0     "_)*%@(4"_)

vS00xL    =:Initx
vS01dxL   =:0,0,0,0:
vS02d1L   =:Initd1
vS06ds    =:Initds
vS07count =:(0,8000"_)*4"_
pd vSgDplot pdDA=:vSgD@vS''

pd 'pdf'
pd 'show'


Egraphic108g.jpg
... geodesic ...


NB. ... execute (ijx) ...

NB. ... using pdDA from the plot above ...

   $pdDA
32001 8 4

NB. ... x0 ...

   A=:cx0"1@:((S00xL{])"2) pdDA

   (6.25p1-<./) A
0.0108606

   (>./-6.25p1"_) A
0.0789094

NB. ... x1 ...

   *./ (0=]) (0.50p1-]) cx1"1@:((S00xL{])"2) pdDA
1

NB. ... x2 decreases continuously ...

   *./ (0<]) (2-/\]) cx2"1@:((S00xL{])"2) pdDA
1


load 'plot'

NB. ... using pdDA from the plot above ...

pd (-@cx2"1;cx0"1)@:((S00xL{])"2) (i.@(18000"_){]) pdDA

pd 'pdf'
pd 'show'


Egraphic1080.jpg
... x0 ...


load 'plot'

NB. ... using pdDA from the plot above ...

pd (-@cx2"1;cx3"1)@:((S00xL{])"2) (i.@(32000"_){]) pdDA

pd 'pdf'
pd 'show'


Egraphic1083.jpg
... x3 ...



Download MoinMoin source: File:TensorExperimentsB10.ijs