User:Tom Allen/TensorExperimentsB06
Jump to navigation
Jump to search
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'
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'
... time-independent component of OyB3 ...
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'
... 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'
... "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'
... "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'
... (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'
... (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'
... 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'
... EM at times 20 30 40 50 60 (later) ...
Download: File:LTEB0600C.txt
Download: File:LTEB0601C.txt
Download MoinMoin source: File:TensorExperimentsB06.ijs