User:Tom Allen/TensorExperimentsB06

From J Wiki
Jump to navigation Jump to search


GTEB0600C.jpg


load 'plot'

vSgDd2M=.[([(0{])@(mcn smx"2 1 B2kCEr1+SBkCEr1)xRef@hkxH),:@(S03xH{])
vSgDd2 =.gXB(-@(([([(0{])@ch2k xRef@hkxH),:@(S03xH{]))+/"1@((<0 2)|:])@((0|:[)smx])(]*/])@(S02d1L{]))+-@vSgDd2M)]

Initds2=.(gXB([(0{])@mcv xRef@hkxH),:@Initx)+/^:2@:*(]*/])@Initdx

gXcw   =:  1.5  0  _  _              _"_
gXpAll =:  0    0  0  0              0"_
gXq2All=:  1    1  1  0              _"_
gXq3All=:  1    0  0  0              0"_
gXq4All=:  1    1  1  0              _"_
gXq5All=:  1    1  3  1             _3"_
gXq6All=:  1    1  1  020.62867905   _"_
gXq7All=:  2    0  1  1.125         _1"_

Initx =:  7 , 0.5p1 ,  1p1  , 0             "_
Initdx=:  0 , 0     , _0.01 , 0.047935260915"_

vS00xL    =:Initx
vS01dxL   =:0,0,0,0:
vS02d1L   =:Initd1
vS06ds    =:Initds
vS07count =:0,628"_
pd vSgDplot pdDA=:vSgD@vS''

pd 'pdf'
pd 'show'


Egraphic0600.jpg
020.62867905*(1/((c^2)*0.5*r))


...

...

...


load 'plot'

gXcw   =:  1.2  0    _    _     _"_
gXpAll =:  0    0    0    0     0"_
gXq2All=:  0.1  0.1  0.1  0     _"_
gXq3All=:  1    0    0    0     0"_
gXq4All=:  0.1  0.1  0.1  0     _"_
gXq5All=:  1    1    3    1    _3"_
gXq6All=:  0.1  0.1  0.1  3     _"_
gXq7All=:  1    0.5  3    0.5  _2"_

yCpts=:_ 14 23 300,_ 0 0 0,_ 0 0 0,:_ 0 0 0"_

pd (gXB(cy0"_ 1;fbF"_ 1)]) pTtotal yCpts''

pd 'pdf'
pd 'show'


Egraphic0601.jpg
... time-independent component of OyB3 ...


GTEB0601C.jpg


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

gXcw   =:  1.2  0    _    _     _"_
gXpAll =:  0    0    0    0     0"_
gXq2All=:  0.1  0.1  0.1  0     _"_
gXq3All=:  1    0    0    0     0"_
gXq4All=:  0.1  0.1  0.1  0     _"_
gXq5All=:  1    1    3    1    _3"_
gXq6All=:  0.1  0.1  0.1  3     _"_
gXq7All=:  1    0.5  3    0.5  _2"_

xCpts=:_ 14 23 26,_ 0 1p1 26,_ 0 2p1 26,:_ 0 2 26"_

NB. ... "divergence" does not include x3 ...

   ((2^_44)qtsz]) (gXB([(mcn smx"2 1 B2kCEr1+SBkCEr1)xRef@hkxH)]) ((0.1;10"_)pTRandom xCpts)''
 _0.00636245 0 0 0
  3.49883e_5 0 0 0
 _0.00584808 0 0 0
 _0.00584808 0 0 0
 0.000886009 0 0 0
  0.00129016 0 0 0
 _0.00480622 0 0 0
 _0.00532779 0 0 0
_0.000861438 0 0 0
_0.000861438 0 0 0

NB. ... no EM ...

   ((2^_44)qtsz]) (gXB([(0 1|:])"3@(+/"1@((<0 1)|:])"5@R)xRef@hkxH)]) ((0.1;1"_)pTRandom xCpts)''
0 0 0 0
0 0 0 0
0 0 0 0
0 0 0 0

0 0 0 0
0 0 0 0
0 0 0 0
0 0 0 0

0 0 0 0
0 0 0 0
0 0 0 0
0 0 0 0

0 0 0 0
0 0 0 0
0 0 0 0
0 0 0 0


...

...

...


load 'plot'

gXcw   =:  1.2   0     _     _   _"_
gXpAll =:  0     0     0     0   0"_
gXq2All=:  0.06  0.06  0.06  3   _"_
gXq3All=:  1     0     0     0   0"_
gXq4All=:  0.06  0.06  0.06  1   _"_
gXq5All=:  1     1     3     1  _3"_
gXq6All=:  0.06  0.06  0.06  0   _"_
gXq7All=:  1     1     3     1  _3"_

yCpts=:_ 0.1 60 300,_ 0 0 0,_ 0 0 0,:_ 0 0 0"_

pd (gXB(cy0"_ 1;OyBcK0"_ 1)]) pTtotal yCpts''

pd 'pdf'
pd 'show'


Egraphic0602.jpg
... steady state (no wave) ...


load 'plot'

gXcw   =:  1.2   0     _     _   _"_
gXpAll =:  0     0     0     0   0"_
gXq2All=:  0.06  0.06  0.06  3   _"_
gXq3All=:  1     0     0     0   0"_
gXq4All=:  0.06  0.06  0.06  1   _"_
gXq5All=:  1     1     3     1  _3"_
gXq6All=:  0.06  0.06  0.06  0   _"_
gXq7All=:  1     1     3     1  _3"_

xCpts=:_ 0.1 35 500,_ 0.25p1 0.25p1 0,_ 0.25p1 0.25p1 0,:_ 1.0 1.0 0"_

pd (gXB(cx0"_ 1;[(0{])"1@(mcn smx"2 1 B2kCEr1+SBkCEr1)xRef@hkxH)]) pTtotal xCpts''

xCpts=:_ 0.1 35 500,_ 0.25p1 0.25p1 0,_ 0.25p1 0.25p1 0,:_ 2.0 2.0 0"_

pd (gXB(cx0"_ 1;[(0{])"1@(mcn smx"2 1 B2kCEr1+SBkCEr1)xRef@hkxH)]) pTtotal xCpts''

xCpts=:_ 0.1 35 500,_ 0.25p1 0.25p1 0,_ 0.25p1 0.25p1 0,:_ 3.0 3.0 0"_

pd (gXB(cx0"_ 1;[(0{])"1@(mcn smx"2 1 B2kCEr1+SBkCEr1)xRef@hkxH)]) pTtotal xCpts''

xCpts=:_ 0.1 35 500,_ 0.25p1 0.25p1 0,_ 0.25p1 0.25p1 0,:_ 4.0 4.0 0"_

pd (gXB(cx0"_ 1;[(0{])"1@(mcn smx"2 1 B2kCEr1+SBkCEr1)xRef@hkxH)]) pTtotal xCpts''

xCpts=:_ 0.1 35 500,_ 0.25p1 0.25p1 0,_ 0.25p1 0.25p1 0,:_ 5.0 5.0 0"_

pd (gXB(cx0"_ 1;[(0{])"1@(mcn smx"2 1 B2kCEr1+SBkCEr1)xRef@hkxH)]) pTtotal xCpts''

pd 'pdf'
pd 'show'


Egraphic0603.jpg
... "divergence" x0 (early) ...


load 'plot'

gXcw   =:  1.2   0     _     _   _"_
gXpAll =:  0     0     0     0   0"_
gXq2All=:  0.06  0.06  0.06  3   _"_
gXq3All=:  1     0     0     0   0"_
gXq4All=:  0.06  0.06  0.06  1   _"_
gXq5All=:  1     1     3     1  _3"_
gXq6All=:  0.06  0.06  0.06  0   _"_
gXq7All=:  1     1     3     1  _3"_

xCpts=:_ 0.1 35 500,_ 0.25p1 0.25p1 0,_ 0.25p1 0.25p1 0,:_ 0.1 0.1 0"_

pd (gXB(cx0"_ 1;[(3{])"1@(mcn smx"2 1 B2kCEr1+SBkCEr1)xRef@hkxH)]) pTtotal xCpts''

xCpts=:_ 0.1 35 500,_ 0.25p1 0.25p1 0,_ 0.25p1 0.25p1 0,:_ 0.5 0.5 0"_

pd (gXB(cx0"_ 1;[(3{])"1@(mcn smx"2 1 B2kCEr1+SBkCEr1)xRef@hkxH)]) pTtotal xCpts''

xCpts=:_ 0.1 35 500,_ 0.25p1 0.25p1 0,_ 0.25p1 0.25p1 0,:_ 1.0 1.0 0"_

pd (gXB(cx0"_ 1;[(3{])"1@(mcn smx"2 1 B2kCEr1+SBkCEr1)xRef@hkxH)]) pTtotal xCpts''

xCpts=:_ 0.1 35 500,_ 0.25p1 0.25p1 0,_ 0.25p1 0.25p1 0,:_ 1.5 1.5 0"_

pd (gXB(cx0"_ 1;[(3{])"1@(mcn smx"2 1 B2kCEr1+SBkCEr1)xRef@hkxH)]) pTtotal xCpts''

xCpts=:_ 0.1 35 500,_ 0.25p1 0.25p1 0,_ 0.25p1 0.25p1 0,:_ 2.0 2.0 0"_

pd (gXB(cx0"_ 1;[(3{])"1@(mcn smx"2 1 B2kCEr1+SBkCEr1)xRef@hkxH)]) pTtotal xCpts''

pd 'pdf'
pd 'show'


Egraphic0604.jpg
... "divergence" x3 (early) ...


load 'plot'

gXcw   =:  1.2   0     _     _   _"_
gXpAll =:  0     0     0     0   0"_
gXq2All=:  0.06  0.06  0.06  3   _"_
gXq3All=:  1     0     0     0   0"_
gXq4All=:  0.06  0.06  0.06  1   _"_
gXq5All=:  1     1     3     1  _3"_
gXq6All=:  0.06  0.06  0.06  0   _"_
gXq7All=:  1     1     3     1  _3"_

xCpts=:_ 0.1 30 500,_ 0.25p1 0.25p1 0,_ 0.25p1 0.25p1 0,:_ 20 20 0"_

pd (gXB(cx0"_ 1;[(B2kCEr0+-@LamCEr0)xRef@hkxH)]) pTtotal xCpts''

xCpts=:_ 0.1 30 500,_ 0.25p1 0.25p1 0,_ 0.25p1 0.25p1 0,:_ 30 30 0"_

pd (gXB(cx0"_ 1;[(B2kCEr0+-@LamCEr0)xRef@hkxH)]) pTtotal xCpts''

xCpts=:_ 0.1 30 500,_ 0.25p1 0.25p1 0,_ 0.25p1 0.25p1 0,:_ 40 40 0"_

pd (gXB(cx0"_ 1;[(B2kCEr0+-@LamCEr0)xRef@hkxH)]) pTtotal xCpts''

xCpts=:_ 0.1 30 500,_ 0.25p1 0.25p1 0,_ 0.25p1 0.25p1 0,:_ 50 50 0"_

pd (gXB(cx0"_ 1;[(B2kCEr0+-@LamCEr0)xRef@hkxH)]) pTtotal xCpts''

xCpts=:_ 0.1 30 500,_ 0.25p1 0.25p1 0,_ 0.25p1 0.25p1 0,:_ 60 60 0"_

pd (gXB(cx0"_ 1;[(B2kCEr0+-@LamCEr0)xRef@hkxH)]) pTtotal xCpts''

pd 'pdf'
pd 'show'


Egraphic0605.jpg
... (B2kCEr0+-@LamCEr0) at times 20 30 40 50 60 (later) ...


load 'plot'

gXcw   =:  1.2   0     _     _   _"_
gXpAll =:  0     0     0     0   0"_
gXq2All=:  0.06  0.06  0.06  3   _"_
gXq3All=:  1     0     0     0   0"_
gXq4All=:  0.06  0.06  0.06  1   _"_
gXq5All=:  1     1     3     1  _3"_
gXq6All=:  0.06  0.06  0.06  0   _"_
gXq7All=:  1     1     3     1  _3"_

xCpts=:_ 0.1 30 500,_ 0.25p1 0.25p1 0,_ 0.25p1 0.25p1 0,:_ 20 20 0"_

pd (gXB(cx0"_ 1;[(SBkCEr0+-@XiCEr0)xRef@hkxH)]) pTtotal xCpts''

xCpts=:_ 0.1 30 500,_ 0.25p1 0.25p1 0,_ 0.25p1 0.25p1 0,:_ 30 30 0"_

pd (gXB(cx0"_ 1;[(SBkCEr0+-@XiCEr0)xRef@hkxH)]) pTtotal xCpts''

xCpts=:_ 0.1 30 500,_ 0.25p1 0.25p1 0,_ 0.25p1 0.25p1 0,:_ 40 40 0"_

pd (gXB(cx0"_ 1;[(SBkCEr0+-@XiCEr0)xRef@hkxH)]) pTtotal xCpts''

xCpts=:_ 0.1 30 500,_ 0.25p1 0.25p1 0,_ 0.25p1 0.25p1 0,:_ 50 50 0"_

pd (gXB(cx0"_ 1;[(SBkCEr0+-@XiCEr0)xRef@hkxH)]) pTtotal xCpts''

xCpts=:_ 0.1 30 500,_ 0.25p1 0.25p1 0,_ 0.25p1 0.25p1 0,:_ 60 60 0"_

pd (gXB(cx0"_ 1;[(SBkCEr0+-@XiCEr0)xRef@hkxH)]) pTtotal xCpts''

pd 'pdf'
pd 'show'


Egraphic0606.jpg
... (SBkCEr0+-@XiCEr0) at times 20 30 40 50 60 (later) ...


load 'plot'

gXcw   =:  1.2   0     _     _   _"_
gXpAll =:  0     0     0     0   0"_
gXq2All=:  0.06  0.06  0.06  3   _"_
gXq3All=:  1     0     0     0   0"_
gXq4All=:  0.06  0.06  0.06  1   _"_
gXq5All=:  1     1     3     1  _3"_
gXq6All=:  0.06  0.06  0.06  0   _"_
gXq7All=:  1     1     3     1  _3"_

xCpts=:_ 0.1 30 500,_ 0.25p1 0.25p1 0,_ 0.25p1 0.25p1 0,:_ 20 20 0"_

pd (gXB(cx0"_ 1;[LamCEr0 xRef@hkxH)]) pTtotal xCpts''

xCpts=:_ 0.1 30 500,_ 0.25p1 0.25p1 0,_ 0.25p1 0.25p1 0,:_ 30 30 0"_

pd (gXB(cx0"_ 1;[LamCEr0 xRef@hkxH)]) pTtotal xCpts''

xCpts=:_ 0.1 30 500,_ 0.25p1 0.25p1 0,_ 0.25p1 0.25p1 0,:_ 40 40 0"_

pd (gXB(cx0"_ 1;[LamCEr0 xRef@hkxH)]) pTtotal xCpts''

xCpts=:_ 0.1 30 500,_ 0.25p1 0.25p1 0,_ 0.25p1 0.25p1 0,:_ 50 50 0"_

pd (gXB(cx0"_ 1;[LamCEr0 xRef@hkxH)]) pTtotal xCpts''

xCpts=:_ 0.1 30 500,_ 0.25p1 0.25p1 0,_ 0.25p1 0.25p1 0,:_ 60 60 0"_

pd (gXB(cx0"_ 1;[LamCEr0 xRef@hkxH)]) pTtotal xCpts''

NB. ...

xCpts=:_ 0.1 30 500,_ 0.25p1 0.25p1 0,_ 0.25p1 0.25p1 0,:_ 20 20 0"_

pd (gXB(cx0"_ 1;[XiCEr0  xRef@hkxH)]) pTtotal xCpts''

xCpts=:_ 0.1 30 500,_ 0.25p1 0.25p1 0,_ 0.25p1 0.25p1 0,:_ 30 30 0"_

pd (gXB(cx0"_ 1;[XiCEr0  xRef@hkxH)]) pTtotal xCpts''

xCpts=:_ 0.1 30 500,_ 0.25p1 0.25p1 0,_ 0.25p1 0.25p1 0,:_ 40 40 0"_

pd (gXB(cx0"_ 1;[XiCEr0  xRef@hkxH)]) pTtotal xCpts''

xCpts=:_ 0.1 30 500,_ 0.25p1 0.25p1 0,_ 0.25p1 0.25p1 0,:_ 50 50 0"_

pd (gXB(cx0"_ 1;[XiCEr0  xRef@hkxH)]) pTtotal xCpts''

xCpts=:_ 0.1 30 500,_ 0.25p1 0.25p1 0,_ 0.25p1 0.25p1 0,:_ 60 60 0"_

pd (gXB(cx0"_ 1;[XiCEr0  xRef@hkxH)]) pTtotal xCpts''

pd 'pdf'
pd 'show'


Egraphic0607.jpg
... LamCEr0 and XiCEr0 at times 20 30 40 50 60 (later) ...


load 'plot'

gXcw   =:  1.2   0     _     _   _"_
gXpAll =:  0     0     0     0   0"_
gXq2All=:  0.06  0.06  0.06  3   _"_
gXq3All=:  1     0     0     0   0"_
gXq4All=:  0.06  0.06  0.06  1   _"_
gXq5All=:  1     1     3     1  _3"_
gXq6All=:  0.06  0.06  0.06  0   _"_
gXq7All=:  1     1     3     1  _3"_

xCpts=:_ 0.1 30 500,_ 0.25p1 0.25p1 0,_ 0.25p1 0.25p1 0,:_ 20 20 0"_

pd (gXB(cx0"_ 1;[(B2kEMr0+SBkEMr0)xRef@hkxH)]) pTtotal xCpts''

xCpts=:_ 0.1 30 500,_ 0.25p1 0.25p1 0,_ 0.25p1 0.25p1 0,:_ 30 30 0"_

pd (gXB(cx0"_ 1;[(B2kEMr0+SBkEMr0)xRef@hkxH)]) pTtotal xCpts''

xCpts=:_ 0.1 30 500,_ 0.25p1 0.25p1 0,_ 0.25p1 0.25p1 0,:_ 40 40 0"_

pd (gXB(cx0"_ 1;[(B2kEMr0+SBkEMr0)xRef@hkxH)]) pTtotal xCpts''

xCpts=:_ 0.1 30 500,_ 0.25p1 0.25p1 0,_ 0.25p1 0.25p1 0,:_ 50 50 0"_

pd (gXB(cx0"_ 1;[(B2kEMr0+SBkEMr0)xRef@hkxH)]) pTtotal xCpts''

xCpts=:_ 0.1 30 500,_ 0.25p1 0.25p1 0,_ 0.25p1 0.25p1 0,:_ 60 60 0"_

pd (gXB(cx0"_ 1;[(B2kEMr0+SBkEMr0)xRef@hkxH)]) pTtotal xCpts''

pd 'pdf'
pd 'show'


Egraphic0608.jpg
... EM at times 20 30 40 50 60 (later) ...



Download: File:LTEB0600C.txt
Download: File:LTEB0601C.txt




Download MoinMoin source: File:TensorExperimentsB06.ijs