User:Tom Allen/TensorExperimentsZ14

From J Wiki
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'


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


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


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


Zgraphic1404.jpg
... x2 and x3 ...



Download MoinMoin source: File:TensorExperimentsZ14.ijs