User:Tom Allen/TensorExperimentsB11

From J Wiki
Jump to: navigation, 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'


Egraphic1182g.jpg
... 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'


Egraphic11820.jpg
... 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'


Egraphic1183g.jpg
... 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'


Egraphic11830.jpg
... 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'


Egraphic1186g.jpg
... 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'


Egraphic11860.jpg
... 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'


Egraphic1188g.jpg
... 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'


Egraphic11880.jpg
... x0 ...



Download MoinMoin source: File:TensorExperimentsB11.ijs