Differences between revisions 18 and 19
 ⇤ ← Revision 18 as of 2013-08-19 04:12:28 → Size: 11756 Editor: TomAllen Comment: ← Revision 19 as of 2013-08-19 04:15:16 → ⇥ Size: 11757 Editor: TomAllen Comment: Deletions are marked like this. Additions are marked like this. Line 100: Line 100: f_{r^{b}}=(f_r)^{b10} f_{r^b}=(f_r)^{b_{10}}

NB. ... script experimentalZ.ijs ...

qfr2=:(cz0^2:)+(cz1^2:)+cz2^2:

NB. ... script experimentalZ.ijs ...

qfr2UUdz        =:(2*cz0),(2*cz1),(2*cz2),0:
qfr2UUdzdz      =:(2,0,0,0:),(0,2,0,0:),(0,0,2,0:),:0,0,0,0:
qfr2UUdzdzdz    =:(3$4:)$0:
qfr2UUdzdzdzdz  =:(4$4:)$0:
qfr2UUdzdzdzdzdz=:(5$4:)$0:

qfr2XX=:qfr2;qfr2UUdz;qfr2UUdzdz;qfr2UUdzdzdz;qfr2UUdzdzdzdz;qfr2UUdzdzdzdzdz

NB. ... script experimentalZ.ijs ...

qfr=:qfr2^0.5"_

NB. ... script experimentalZ.ijs ...

qfr2PWR=:0.5"_

qfrU0=:                                                                qfr2^qfr2PWR
qfrU1=:                                                        qfr2PWR*qfr2^qfr2PWR+_1:
qfrU2=:                                          (qfr2PWR+_1:)*qfr2PWR*qfr2^qfr2PWR+_2:
qfrU3=:                            (qfr2PWR+_2:)*(qfr2PWR+_1:)*qfr2PWR*qfr2^qfr2PWR+_3:
qfrU4=:              (qfr2PWR+_3:)*(qfr2PWR+_2:)*(qfr2PWR+_1:)*qfr2PWR*qfr2^qfr2PWR+_4:
qfrU5=:(qfr2PWR+_4:)*(qfr2PWR+_3:)*(qfr2PWR+_2:)*(qfr2PWR+_1:)*qfr2PWR*qfr2^qfr2PWR+_5:

qfrXX=:[cfFN5(qfrU0;qfrU1;qfrU2;qfrU3;qfrU4;qfrU5),:qfr2XX

NB. ... script experimentalZ.ijs ...

qfrb=:qfr^gXb10

NB. ... script experimentalZ.ijs ...

qfrbU0=:                                                      qfr^gXb10
qfrbU1=:                                                gXb10*qfr^gXb10+_1:
qfrbU2=:                                    (gXb10+_1:)*gXb10*qfr^gXb10+_2:
qfrbU3=:                        (gXb10+_2:)*(gXb10+_1:)*gXb10*qfr^gXb10+_3:
qfrbU4=:            (gXb10+_3:)*(gXb10+_2:)*(gXb10+_1:)*gXb10*qfr^gXb10+_4:
qfrbU5=:(gXb10+_4:)*(gXb10+_3:)*(gXb10+_2:)*(gXb10+_1:)*gXb10*qfr^gXb10+_5:

qfrbXX=:[cfFN5(qfrbU0;qfrbU1;qfrbU2;qfrbU3;qfrbU4;qfrbU5),:qfrXX

NB. ... script experimentalZ.ijs ...

qfbr=:gXb11*qfr

NB. ... script experimentalZ.ijs ...

qfbrXX=:[cfMP5 gXb11XX,:qfrXX

NB. ... script experimentalZ.ijs ...

qfsin=:sin@qfbr

NB. ... script experimentalZ.ijs ...

qfsinXX=:[cfFN5(sin@qfbr;cos@qfbr;-@sin@qfbr;-@cos@qfbr;sin@qfbr;cos@qfbr),:qfbrXX

NB. ... script experimentalZ.ijs ...

qfsinP=:qfsin^gXb12

NB. ... script experimentalZ.ijs ...

qfsinPU0=:                                                      qfsin^gXb12
qfsinPU1=:                                                gXb12*qfsin^gXb12+_1:
qfsinPU2=:                                    (gXb12+_1:)*gXb12*qfsin^gXb12+_2:
qfsinPU3=:                        (gXb12+_2:)*(gXb12+_1:)*gXb12*qfsin^gXb12+_3:
qfsinPU4=:            (gXb12+_3:)*(gXb12+_2:)*(gXb12+_1:)*gXb12*qfsin^gXb12+_4:
qfsinPU5=:(gXb12+_4:)*(gXb12+_3:)*(gXb12+_2:)*(gXb12+_1:)*gXb12*qfsin^gXb12+_5:

qfsinPXX=:[cfFN5(qfsinPU0;qfsinPU1;qfsinPU2;qfsinPU3;qfsinPU4;qfsinPU5),:qfsinXX

NB. ... script experimentalZ.ijs ...

qfbt=:gXb13*cz3

NB. ... script experimentalZ.ijs ...

qfbtXX=:[cfMP5 gXb13XX,:cz3XX

NB. ... script experimentalZ.ijs ...

qfsinbt=:sin@qfbt

NB. ... script experimentalZ.ijs ...

qfsinbtXX=:[cfFN5(sin@qfbt;cos@qfbt;-@sin@qfbt;-@cos@qfbt;sin@qfbt;cos@qfbt),:qfbtXX

NB. ... script experimentalZ.ijs ...

qfcosbt=:cos@qfbt

NB. ... script experimentalZ.ijs ...

qfcosbtXX=:[cfFN5(cos@qfbt;-@sin@qfbt;-@cos@qfbt;sin@qfbt;cos@qfbt;-@sin@qfbt),:qfbtXX

NB. ... script experimentalZ.ijs ...

qffz0=:gXb14*cz0*qfrb*qfsinP
qffz1=:gXb14*cz1*qfrb*qfsinP
qffz2=:gXb14*cz2*qfrb*qfsinP

qffz0XX=:[cfMP5 gXb14XX,:[cfMP5 cz0XX,:[cfMP5 qfrbXX,:qfsinPXX
qffz1XX=:[cfMP5 gXb14XX,:[cfMP5 cz1XX,:[cfMP5 qfrbXX,:qfsinPXX
qffz2XX=:[cfMP5 gXb14XX,:[cfMP5 cz2XX,:[cfMP5 qfrbXX,:qfsinPXX

OzC0=:qffz0*qfsinbt
OzC1=:qffz1*qfsinbt
OzC2=:qffz2*qfsinbt

OzC0XX  =:[cfMP4 qffz0XX,:qfsinbtXX
OzC1XX  =:[cfMP4 qffz1XX,:qfsinbtXX
OzC2XX  =:[cfMP4 qffz2XX,:qfsinbtXX

gXb13RXX=:%@gXb13;((1$4:)$0:);((2$4:)$0:);((3$4:)$0:);((4$4:)$0:);(5$4:)$0:
OzC3t0XX=:[cfMP4((0{])&.>@}.@qffz0XX,<@(_"_)),:[cfMP4 gXb13RXX,:qfcosbtXX
OzC3t1XX=:[cfMP4((1{])&.>@}.@qffz1XX,<@(_"_)),:[cfMP4 gXb13RXX,:qfcosbtXX
OzC3t2XX=:[cfMP4((2{])&.>@}.@qffz2XX,<@(_"_)),:[cfMP4 gXb13RXX,:qfcosbtXX
OzC3XX  =:OzC3t2XX+&.>OzC3t1XX+&.>OzC3t0XX

OzCXXS=:<@:>"1@(0|:])@(OzC0XX,OzC1XX,OzC2XX,:OzC3XX)

load 'plot'

gXcw  =:  _  _  _  _  _  _"_
gXbAll=: _1  1  7  1  1  _"_

zCpts=:_ 0 2p1 300,_ 0 0 0,_ 0 0 0,:_ 0.00p1 0.00p1 0"_

pd (gXZ(cz0"_ 1;OzC0"_ 1)]) pTtotal zCpts''

zCpts=:_ 0 2p1 300,_ 0 0 0,_ 0 0 0,:_ 0.07p1 0.07p1 0"_

pd (gXZ(cz0"_ 1;OzC0"_ 1)]) pTtotal zCpts''

zCpts=:_ 0 2p1 300,_ 0 0 0,_ 0 0 0,:_ 0.15p1 0.15p1 0"_

pd (gXZ(cz0"_ 1;OzC0"_ 1)]) pTtotal zCpts''

zCpts=:_ 0 2p1 300,_ 0 0 0,_ 0 0 0,:_ 0.25p1 0.25p1 0"_

pd (gXZ(cz0"_ 1;OzC0"_ 1)]) pTtotal zCpts''

zCpts=:_ 0 2p1 300,_ 0 0 0,_ 0 0 0,:_ 0.50p1 0.50p1 0"_

pd (gXZ(cz0"_ 1;OzC0"_ 1)]) pTtotal zCpts''

zCpts=:_ 0 2p1 300,_ 0 0 0,_ 0 0 0,:_ 0.75p1 0.75p1 0"_

pd (gXZ(cz0"_ 1;OzC0"_ 1)]) pTtotal zCpts''

zCpts=:_ 0 2p1 300,_ 0 0 0,_ 0 0 0,:_ 1.00p1 1.00p1 0"_

pd (gXZ(cz0"_ 1;OzC0"_ 1)]) pTtotal zCpts''

zCpts=:_ 0 2p1 300,_ 0 0 0,_ 0 0 0,:_ 1.07p1 1.07p1 0"_

pd (gXZ(cz0"_ 1;OzC0"_ 1)]) pTtotal zCpts''

zCpts=:_ 0 2p1 300,_ 0 0 0,_ 0 0 0,:_ 1.15p1 1.15p1 0"_

pd (gXZ(cz0"_ 1;OzC0"_ 1)]) pTtotal zCpts''

zCpts=:_ 0 2p1 300,_ 0 0 0,_ 0 0 0,:_ 1.25p1 1.25p1 0"_

pd (gXZ(cz0"_ 1;OzC0"_ 1)]) pTtotal zCpts''

zCpts=:_ 0 2p1 300,_ 0 0 0,_ 0 0 0,:_ 1.50p1 1.50p1 0"_

pd (gXZ(cz0"_ 1;OzC0"_ 1)]) pTtotal zCpts''

zCpts=:_ 0 2p1 300,_ 0 0 0,_ 0 0 0,:_ 1.75p1 1.75p1 0"_

pd (gXZ(cz0"_ 1;OzC0"_ 1)]) pTtotal zCpts''

pd 'pdf'
pd 'reset'
 ... a standing wave ...

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

gXcw  =:  _    _    _  _    _    _"_
gXbAll=: _0.6  3.5  3  1.3  1.2  _"_

zCpts=:_ 0.1p1 0.5p1 27,_ 0.1p1 0.5p1 27,_ 0.1p1 0.5p1 27,:_ 0.1p1 0.9p1 27"_

NB. ... derivatives ...

(gXZ(cfXX1@OzCXXS"_ 1(((2^_44);2^_14)qteq[;])"2(0|:[:(gXZ cfXX0@OzCXXS])D.1])"_ 1)]) ((0;50"_)pTRandom zCpts)''
1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1

(gXZ(cfXX2@OzCXXS"_ 1(((2^_44);2^_10)qteq[;])"3(0|:[:(gXZ cfXX1@OzCXXS])D.1])"_ 1)]) ((0;50"_)pTRandom zCpts)''
1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1

(gXZ(cfXX3@OzCXXS"_ 1(((2^_44);2^_06)qteq[;])"4(0|:[:(gXZ cfXX2@OzCXXS])D.1])"_ 1)]) ((0;50"_)pTRandom zCpts)''
1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1

(gXZ(cfXX4@OzCXXS"_ 1(((2^_44);2^_06)qteq[;])"5(0|:[:(gXZ cfXX3@OzCXXS])D.1])"_ 1)]) ((0;50"_)pTRandom zCpts)''
1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1

NB. ... divergence is 0                          ...
NB. ... (compare Derek F. Lawden equation 28.18) ...

(0=]) ((2^_44)qtsz]) +/@((<0 1)|:])"2 (gXZ cfXX1@OzCXXS"_ 1]) ((0;50"_)pTRandom zCpts)''
1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1

NB. ... script experimentalZ.ijs ...

OzCD4        =:>@((<4 0){])
OzCD4dz      =:>@((<4 1){])
OzCD4dzdz    =:>@((<4 2){])
OzCD4dzdzdz  =:>@((<4 3){])
OzCD4dzdzdzdz=:>@((<4 4){])

OzCXX=:<@:>"1@(0|:])@:(OzCXXS"_ 1)

hkxL=:[(],[OzCXX zpts)hkxK
hkyL=:[(],[OzCXX zpts)hkyK

TomAllen/TensorExperimentsZ02 (last edited 2013-08-19 04:15:16 by TomAllen)