User:Tom Allen/TensorExperimentsZ32

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'

... geodesic ...

...

... 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

---