User:Tom Allen/TensorExperimentsZ32
Jump to navigation
Jump to search
NB. ... script experimentalZ.ijs ... hkyGmetric=:<@:>"1@(0|:])@:(((1 0 0 0,0 1 0 0,0 0 1 0,:0 0 0 _1"_);((3$4:)$0:);((4$4:)$0:);((5$4:)$0:);((6$4:)$0:);1$_"_)"_ 1) hkyG=:(6$a:"_),(6$a:"_),hkyGmetric,(6$a:"_),<@:>"1@(0|:])@:(OyAXXS"_ 1),:6$a:"_
gXcw =. 1 0 _ _ _ "_ gXpAll=. 1 1 1 2 _0.6"_ Inity =. _5.000 , 0.0000 , 0 , 0 "_ Initdy=. ( 0.000 , 0.1600 , 0 , 0.375"_)*%@(64"_) Initds=.%:@|@((gXT([(0{])@mcv yRef@hkyG),:@Inity)+/^:2@:*(]*/])@Initdy) vS00xL =.Inity vS01dxL =.0,0,0,0: vS02d1L =.Initdy%Initds vS06ds =.Initds vS07count=.(0,123"_)*64"_ vSgDd2=.gXT([([(0{])@E1cn yRef@hkyG),:@(S03xH{]))] pdDU=:vSgD@vS'' Inity =. _5.000 , 0.0000 , 0 , 0 "_ Initdy=. ( 0.000 , 0.1600 , 0 , 0.375"_)*%@(128"_) Initds=.%:@|@((gXT([(0{])@mcv yRef@hkyG),:@Inity)+/^:2@:*(]*/])@Initdy) vS00xL =.Inity vS01dxL =.0,0,0,0: vS02d1L =.Initdy%Initds vS06ds =.Initds vS07count=.(0,123"_)*128"_ vSgDd2=.gXT([([(0{])@E1cn yRef@hkyG),:@(S03xH{]))] pdDV=:vSgD@vS'' NB. ... overwrite ... nCcn3 =: cn3"1@:((S00xL{])"2) pdDV nCOcn =:(gXT([Ocn yRef@hkyG)])@:((S00xL{])"2) pdDV nCOcndc =:(gXT([Ocndc yRef@hkyG)])@:((S00xL{])"2) pdDV nCE1cn =:(gXT([E1cn yRef@hkyG)])@:((S00xL{])"2) pdDV nCEEMcn =:(gXT([EEMcn yRef@hkyG)])@:((S00xL{])"2) pdDV nCHEMcn =:(gXT([HEMcn yRef@hkyG)])@:((S00xL{])"2) pdDV nCPEMcn =:(gXT([PEMcn yRef@hkyG)])@:((S00xL{])"2) pdDV nCJEMcn =:(gXT([JEMcn yRef@hkyG)])@:((S00xL{])"2) pdDV nCJEMcndn=:(gXT([JEMcndn yRef@hkyG)])@:((S00xL{])"2) pdDV nCX1cncv =:(gXT([X1cncv yRef@hkyG)])@:((S00xL{])"2) pdDV pdDU arrayPut 'pdDU' arrayClose 'pdDU' pdDV arrayPut 'pdDV' arrayClose 'pdDV' nCcn3 arrayPut 'nCcn3' arrayClose 'nCcn3' nCOcn arrayPut 'nCOcn' arrayClose 'nCOcn' nCOcndc arrayPut 'nCOcndc' arrayClose 'nCOcndc' nCE1cn arrayPut 'nCE1cn' arrayClose 'nCE1cn' nCEEMcn arrayPut 'nCEEMcn' arrayClose 'nCEEMcn' nCHEMcn arrayPut 'nCHEMcn' arrayClose 'nCHEMcn' nCPEMcn arrayPut 'nCPEMcn' arrayClose 'nCPEMcn' nCJEMcn arrayPut 'nCJEMcn' arrayClose 'nCJEMcn' nCJEMcndn arrayPut 'nCJEMcndn' arrayClose 'nCJEMcndn' nCX1cncv arrayPut 'nCX1cncv' arrayClose 'nCX1cncv'
NB. ... establish arrays if necessary ... arrayOpen 'pdDU' pdDU=: 7873 8 4 arrayGet 'pdDU' NB. ... shape required ... arrayClose 'pdDU' arrayOpen 'pdDV' pdDV=: 15745 8 4 arrayGet 'pdDV' arrayClose 'pdDV' arrayOpen 'nCcn3' nCcn3=: 15745 arrayGet 'nCcn3' arrayClose 'nCcn3' arrayOpen 'nCOcn' nCOcn=: 15745 4 arrayGet 'nCOcn' arrayClose 'nCOcn' arrayOpen 'nCOcndc' nCOcndc=: 15745 4 4 arrayGet 'nCOcndc' arrayClose 'nCOcndc' arrayOpen 'nCE1cn' nCE1cn=: 15745 4 arrayGet 'nCE1cn' arrayClose 'nCE1cn' arrayOpen 'nCEEMcn' nCEEMcn=: 15745 3 arrayGet 'nCEEMcn' arrayClose 'nCEEMcn' arrayOpen 'nCHEMcn' nCHEMcn=: 15745 3 arrayGet 'nCHEMcn' arrayClose 'nCHEMcn' arrayOpen 'nCPEMcn' nCPEMcn=: 15745 3 arrayGet 'nCPEMcn' arrayClose 'nCPEMcn' arrayOpen 'nCJEMcn' nCJEMcn=: 15745 4 arrayGet 'nCJEMcn' arrayClose 'nCJEMcn' arrayOpen 'nCJEMcndn' nCJEMcndn=:15745 4 4 arrayGet 'nCJEMcndn' arrayClose 'nCJEMcndn' arrayOpen 'nCX1cncv' nCX1cncv=: 15745 4 4 arrayGet 'nCX1cncv' arrayClose 'nCX1cncv'
load 'plot' pd <"1@(0|:])@:((cy0,cy1)@(S00xL{])"2) pdDU pd <"1@(0|:])@:((cy0,cy1)@(S00xL{])"2) pdDV pd 'pdf' pd 'show'
...
... valid forms ...
...
... (+/-)energy flow with (-/+)energy density ...
...
NB. ... execute (ijx) ... gXcw =: 1 0 _ _ _ "_ gXpAll=: 1 1 1 2 _0.6"_ xCpts=:_ 1 3 26,_ 0 1p1 26,_ 0 2p1 26,:_ 0 20 26"_ NB. ... reversed ... NB. ... (+/-)'energy flow' with (-/+)'energy density' ... (gXT(([(3{])"1@Ocn xRef@hkxF)-:"0-@((gXp3*gXp4^_1:)*(sin@(cx0+gXp4*cx3)*cx0^_1:)+-@(cos@(cx0+gXp4*cx3)*cx0^_2:))"_ 1)]) ((0.1;50"_)pTRandom xCpts)'' 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. ... execute (ijx) ... gXcw =: 1 0 _ _ _ "_ gXpAll=: 1 1 1 2 _0.6"_ xCpts=:_ 1 3 26,_ 0 1p1 26,_ 0 2p1 26,:_ 0 20 26"_ (gXT(([(3{])"1@Ocn yRef@hkxF)-:"0[(3{])"1@Ocn xRef@hkxF)]) ((0.1;50"_)pTRandom xCpts)'' 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
Download: File:LTEZ3200C.txt
Download: File:LTEZ3201D.txt
Download: File:LTEZ3202C.txt
Download: File:LTEZ3204C.txt
Download MoinMoin source: File:TensorExperimentsZ32.ijs