# User:Tom Allen/TensorExperimentsZ23

Jump to navigation Jump to search

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

gXm00=:(<0 0){[
gXm01=:(<0 1){[
gXm02=:(<0 2){[
gXm03=:(<0 3){[
gXm10=:(<1 0){[
gXm11=:(<1 1){[
gXm12=:(<1 2){[
gXm13=:(<1 3){[

gXM=:gXm0All,:gXm1All

gXm00XX=:gXm00;((1\$4:)\$0:);((2\$4:)\$0:);((3\$4:)\$0:);((4\$4:)\$0:);(5\$4:)\$0:
gXm01XX=:gXm01;((1\$4:)\$0:);((2\$4:)\$0:);((3\$4:)\$0:);((4\$4:)\$0:);(5\$4:)\$0:
gXm02XX=:gXm02;((1\$4:)\$0:);((2\$4:)\$0:);((3\$4:)\$0:);((4\$4:)\$0:);(5\$4:)\$0:
gXm03XX=:gXm03;((1\$4:)\$0:);((2\$4:)\$0:);((3\$4:)\$0:);((4\$4:)\$0:);(5\$4:)\$0:
gXm10XX=:gXm10;((1\$4:)\$0:);((2\$4:)\$0:);((3\$4:)\$0:);((4\$4:)\$0:);(5\$4:)\$0:
gXm11XX=:gXm11;((1\$4:)\$0:);((2\$4:)\$0:);((3\$4:)\$0:);((4\$4:)\$0:);(5\$4:)\$0:
gXm12XX=:gXm12;((1\$4:)\$0:);((2\$4:)\$0:);((3\$4:)\$0:);((4\$4:)\$0:);(5\$4:)\$0:
gXm13XX=:gXm13;((1\$4:)\$0:);((2\$4:)\$0:);((3\$4:)\$0:);((4\$4:)\$0:);(5\$4:)\$0:

cn0=:0{]
cn1=:1{]
cn2=:2{]
cn3=:3{]

cn0XX=:cn0;1 0 0 0;((2\$4:)\$0:);((3\$4:)\$0:);((4\$4:)\$0:);(5\$4:)\$0:
cn1XX=:cn1;0 1 0 0;((2\$4:)\$0:);((3\$4:)\$0:);((4\$4:)\$0:);(5\$4:)\$0:
cn2XX=:cn2;0 0 1 0;((2\$4:)\$0:);((3\$4:)\$0:);((4\$4:)\$0:);(5\$4:)\$0:
cn3XX=:cn3;0 0 0 1;((2\$4:)\$0:);((3\$4:)\$0:);((4\$4:)\$0:);(5\$4:)\$0:
```

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

fez=:^@((gXm10*cn0)+gXm11*cn3)+^@-@((gXm10*cn0)+gXm11*cn3)

fezn0=:gXm10*^@((gXm10*cn0)+gXm11*cn3)-^@-@((gXm10*cn0)+gXm11*cn3)
fezn3=:gXm11*^@((gXm10*cn0)+gXm11*cn3)-^@-@((gXm10*cn0)+gXm11*cn3)
fezdn=:fezn0,0,0,fezn3

fezn0n0=:(gXm10^2:)*fez
fezn0n3=:gXm10*gXm11*fez
fezn3n0=:gXm11*gXm10*fez
fezn3n3=:(gXm11^2:)*fez
fezdndn=:(fezn0n0,0,0,fezn0n3),0 0 0 0,0 0 0 0,:fezn3n0,0,0,fezn3n3
```

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

gXm0All=: _    _    _    _  "_
gXm1All=: 1.4  1.5  1.6  1.7"_

nCpts=:_ _4 4 161,_ 0 0 0,_ 0 0 0,:_ _4 4 161"_

(gXM(fezdn"_ 1  (((2^_44);2^_18)qteq[;])"1(0|:[:(gXM fez  ])D.1])"_ 1)]) ((0;50"_)pTRandom nCpts)''
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

(gXM(fezdndn"_ 1(((2^_44);2^_21)qteq[;])"2(0|:[:(gXM fezdn])D.1])"_ 1)]) ((0;50"_)pTRandom nCpts)''
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 ...

fex=:gXm13*fez^_1:

fexn0=:-@gXm13*(fez^_2:)*fezn0
fexn3=:-@gXm13*(fez^_2:)*fezn3
fexdn=:fexn0,0,0,fexn3

fexn0n0=:(-@gXm13*(fez^_2:)*fezn0n0)+2*gXm13*(fez^_3:)*fezn0^2:
fexn0n3=:(-@gXm13*(fez^_2:)*fezn0n3)+2*gXm13*(fez^_3:)*fezn3*fezn0
fexn3n0=:(-@gXm13*(fez^_2:)*fezn3n0)+2*gXm13*(fez^_3:)*fezn0*fezn3
fexn3n3=:(-@gXm13*(fez^_2:)*fezn3n3)+2*gXm13*(fez^_3:)*fezn3^2:
fexdndn=:(fexn0n0,0,0,fexn0n3),0 0 0 0,0 0 0 0,:fexn3n0,0,0,fexn3n3
```

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

gXm0All=: _    _    _    _  "_
gXm1All=: 1.4  1.5  1.6  1.7"_

nCpts=:_ _4 4 161,_ 0 0 0,_ 0 0 0,:_ _4 4 161"_

(gXM(fexdn"_ 1  (((2^_44);2^_18)qteq[;])"1(0|:[:(gXM fex  ])D.1])"_ 1)]) ((0;50"_)pTRandom nCpts)''
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

(gXM(fexdndn"_ 1(((2^_44);2^_18)qteq[;])"2(0|:[:(gXM fexdn])D.1])"_ 1)]) ((0;50"_)pTRandom nCpts)''
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 ...

OnD0=:(-@gXm12*fex*sin@(gXm12*cn3))+fexn3*cos@(gXm12*cn3)
OnD3=:-@fexn0*cos@(gXm12*cn3)
OnD =:OnD0,0,0,OnD3

OnD0n0=:(-@gXm12*fexn0*sin@(gXm12*cn3))+fexn3n0*cos@(gXm12*cn3)
OnD0n3=:(-@(gXm12^2:)*fex*cos@(gXm12*cn3))+(_2*gXm12*fexn3*sin@(gXm12*cn3))+fexn3n3*cos@(gXm12*cn3)
OnD3n0=:-@fexn0n0*cos@(gXm12*cn3)
OnD3n3=:(gXm12*fexn0*sin@(gXm12*cn3))+-@fexn0n3*cos@(gXm12*cn3)
OnDdn =:(OnD0n0,0,0,OnD0n3),0 0 0 0,0 0 0 0,:OnD3n0,0,0,OnD3n3
```

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

gXm0All=: _    _    _    _  "_
gXm1All=: 1.4  1.5  1.6  1.7"_

nCpts=:_ _4 4 161,_ 0 0 0,_ 0 0 0,:_ _4 4 161"_

(gXM(OnDdn"_ 1(((2^_44);2^_15)qteq[;])"2(0|:[:(gXM OnD])D.1])"_ 1)]) ((0;50"_)pTRandom nCpts)''
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 ...

(0=]) ((2^_44)qtsz]) +/@((<0 1)|:])"2 (gXM OnDdn"_ 1]) ((0;50"_)pTRandom nCpts)''
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
```

```load 'plot'

gXm0All=: _  _  _  _"_
gXm1All=: 1  0  1  1"_

nCpts=:_ _2p1 2p1 161,_ 0 0 0,_ 0 0 0,:_ 0.00p1 0.00p1 0"_

pd (gXM(cn0"_ 1;OnD0"_ 1)]) pTtotal nCpts''

nCpts=:_ _2p1 2p1 161,_ 0 0 0,_ 0 0 0,:_ 0.10p1 0.10p1 0"_

pd (gXM(cn0"_ 1;OnD0"_ 1)]) pTtotal nCpts''

nCpts=:_ _2p1 2p1 161,_ 0 0 0,_ 0 0 0,:_ 0.20p1 0.20p1 0"_

pd (gXM(cn0"_ 1;OnD0"_ 1)]) pTtotal nCpts''

nCpts=:_ _2p1 2p1 161,_ 0 0 0,_ 0 0 0,:_ 0.30p1 0.30p1 0"_

pd (gXM(cn0"_ 1;OnD0"_ 1)]) pTtotal nCpts''

nCpts=:_ _2p1 2p1 161,_ 0 0 0,_ 0 0 0,:_ 0.40p1 0.40p1 0"_

pd (gXM(cn0"_ 1;OnD0"_ 1)]) pTtotal nCpts''

nCpts=:_ _2p1 2p1 161,_ 0 0 0,_ 0 0 0,:_ 0.50p1 0.50p1 0"_

pd (gXM(cn0"_ 1;OnD0"_ 1)]) pTtotal nCpts''

nCpts=:_ _2p1 2p1 161,_ 0 0 0,_ 0 0 0,:_ 0.60p1 0.60p1 0"_

pd (gXM(cn0"_ 1;OnD0"_ 1)]) pTtotal nCpts''

nCpts=:_ _2p1 2p1 161,_ 0 0 0,_ 0 0 0,:_ 0.70p1 0.70p1 0"_

pd (gXM(cn0"_ 1;OnD0"_ 1)]) pTtotal nCpts''

nCpts=:_ _2p1 2p1 161,_ 0 0 0,_ 0 0 0,:_ 0.80p1 0.80p1 0"_

pd (gXM(cn0"_ 1;OnD0"_ 1)]) pTtotal nCpts''

nCpts=:_ _2p1 2p1 161,_ 0 0 0,_ 0 0 0,:_ 0.90p1 0.90p1 0"_

pd (gXM(cn0"_ 1;OnD0"_ 1)]) pTtotal nCpts''

nCpts=:_ _2p1 2p1 161,_ 0 0 0,_ 0 0 0,:_ 1.00p1 1.00p1 0"_

pd (gXM(cn0"_ 1;OnD0"_ 1)]) pTtotal nCpts''

nCpts=:_ _2p1 2p1 161,_ 0 0 0,_ 0 0 0,:_ 1.10p1 1.10p1 0"_

pd (gXM(cn0"_ 1;OnD0"_ 1)]) pTtotal nCpts''

nCpts=:_ _2p1 2p1 161,_ 0 0 0,_ 0 0 0,:_ 1.20p1 1.20p1 0"_

pd (gXM(cn0"_ 1;OnD0"_ 1)]) pTtotal nCpts''

nCpts=:_ _2p1 2p1 161,_ 0 0 0,_ 0 0 0,:_ 1.30p1 1.30p1 0"_

pd (gXM(cn0"_ 1;OnD0"_ 1)]) pTtotal nCpts''

nCpts=:_ _2p1 2p1 161,_ 0 0 0,_ 0 0 0,:_ 1.40p1 1.40p1 0"_

pd (gXM(cn0"_ 1;OnD0"_ 1)]) pTtotal nCpts''

nCpts=:_ _2p1 2p1 161,_ 0 0 0,_ 0 0 0,:_ 1.50p1 1.50p1 0"_

pd (gXM(cn0"_ 1;OnD0"_ 1)]) pTtotal nCpts''

nCpts=:_ _2p1 2p1 161,_ 0 0 0,_ 0 0 0,:_ 1.60p1 1.60p1 0"_

pd (gXM(cn0"_ 1;OnD0"_ 1)]) pTtotal nCpts''

nCpts=:_ _2p1 2p1 161,_ 0 0 0,_ 0 0 0,:_ 1.70p1 1.70p1 0"_

pd (gXM(cn0"_ 1;OnD0"_ 1)]) pTtotal nCpts''

nCpts=:_ _2p1 2p1 161,_ 0 0 0,_ 0 0 0,:_ 1.80p1 1.80p1 0"_

pd (gXM(cn0"_ 1;OnD0"_ 1)]) pTtotal nCpts''

nCpts=:_ _2p1 2p1 161,_ 0 0 0,_ 0 0 0,:_ 1.90p1 1.90p1 0"_

pd (gXM(cn0"_ 1;OnD0"_ 1)]) pTtotal nCpts''

nCpts=:_ _2p1 2p1 161,_ 0 0 0,_ 0 0 0,:_ 2.00p1 2.00p1 0"_

pd (gXM(cn0"_ 1;OnD0"_ 1)]) pTtotal nCpts''

pd 'pdf'
pd 'show'
```

 ... energy flow (string) ...

```load 'plot'

gXm0All=: _  _  _  _"_
gXm1All=: 1  0  1  1"_

nCpts=:_ _2p1 2p1 161,_ 0 0 0,_ 0 0 0,:_ 0.00p1 0.00p1 0"_

pd (gXM(cn0"_ 1;OnD3"_ 1)]) pTtotal nCpts''

nCpts=:_ _2p1 2p1 161,_ 0 0 0,_ 0 0 0,:_ 0.10p1 0.10p1 0"_

pd (gXM(cn0"_ 1;OnD3"_ 1)]) pTtotal nCpts''

nCpts=:_ _2p1 2p1 161,_ 0 0 0,_ 0 0 0,:_ 0.20p1 0.20p1 0"_

pd (gXM(cn0"_ 1;OnD3"_ 1)]) pTtotal nCpts''

nCpts=:_ _2p1 2p1 161,_ 0 0 0,_ 0 0 0,:_ 0.30p1 0.30p1 0"_

pd (gXM(cn0"_ 1;OnD3"_ 1)]) pTtotal nCpts''

nCpts=:_ _2p1 2p1 161,_ 0 0 0,_ 0 0 0,:_ 0.40p1 0.40p1 0"_

pd (gXM(cn0"_ 1;OnD3"_ 1)]) pTtotal nCpts''

nCpts=:_ _2p1 2p1 161,_ 0 0 0,_ 0 0 0,:_ 0.50p1 0.50p1 0"_

pd (gXM(cn0"_ 1;OnD3"_ 1)]) pTtotal nCpts''

nCpts=:_ _2p1 2p1 161,_ 0 0 0,_ 0 0 0,:_ 0.60p1 0.60p1 0"_

pd (gXM(cn0"_ 1;OnD3"_ 1)]) pTtotal nCpts''

nCpts=:_ _2p1 2p1 161,_ 0 0 0,_ 0 0 0,:_ 0.70p1 0.70p1 0"_

pd (gXM(cn0"_ 1;OnD3"_ 1)]) pTtotal nCpts''

nCpts=:_ _2p1 2p1 161,_ 0 0 0,_ 0 0 0,:_ 0.80p1 0.80p1 0"_

pd (gXM(cn0"_ 1;OnD3"_ 1)]) pTtotal nCpts''

nCpts=:_ _2p1 2p1 161,_ 0 0 0,_ 0 0 0,:_ 0.90p1 0.90p1 0"_

pd (gXM(cn0"_ 1;OnD3"_ 1)]) pTtotal nCpts''

nCpts=:_ _2p1 2p1 161,_ 0 0 0,_ 0 0 0,:_ 1.00p1 1.00p1 0"_

pd (gXM(cn0"_ 1;OnD3"_ 1)]) pTtotal nCpts''

nCpts=:_ _2p1 2p1 161,_ 0 0 0,_ 0 0 0,:_ 1.10p1 1.10p1 0"_

pd (gXM(cn0"_ 1;OnD3"_ 1)]) pTtotal nCpts''

nCpts=:_ _2p1 2p1 161,_ 0 0 0,_ 0 0 0,:_ 1.20p1 1.20p1 0"_

pd (gXM(cn0"_ 1;OnD3"_ 1)]) pTtotal nCpts''

nCpts=:_ _2p1 2p1 161,_ 0 0 0,_ 0 0 0,:_ 1.30p1 1.30p1 0"_

pd (gXM(cn0"_ 1;OnD3"_ 1)]) pTtotal nCpts''

nCpts=:_ _2p1 2p1 161,_ 0 0 0,_ 0 0 0,:_ 1.40p1 1.40p1 0"_

pd (gXM(cn0"_ 1;OnD3"_ 1)]) pTtotal nCpts''

nCpts=:_ _2p1 2p1 161,_ 0 0 0,_ 0 0 0,:_ 1.50p1 1.50p1 0"_

pd (gXM(cn0"_ 1;OnD3"_ 1)]) pTtotal nCpts''

nCpts=:_ _2p1 2p1 161,_ 0 0 0,_ 0 0 0,:_ 1.60p1 1.60p1 0"_

pd (gXM(cn0"_ 1;OnD3"_ 1)]) pTtotal nCpts''

nCpts=:_ _2p1 2p1 161,_ 0 0 0,_ 0 0 0,:_ 1.70p1 1.70p1 0"_

pd (gXM(cn0"_ 1;OnD3"_ 1)]) pTtotal nCpts''

nCpts=:_ _2p1 2p1 161,_ 0 0 0,_ 0 0 0,:_ 1.80p1 1.80p1 0"_

pd (gXM(cn0"_ 1;OnD3"_ 1)]) pTtotal nCpts''

nCpts=:_ _2p1 2p1 161,_ 0 0 0,_ 0 0 0,:_ 1.90p1 1.90p1 0"_

pd (gXM(cn0"_ 1;OnD3"_ 1)]) pTtotal nCpts''

nCpts=:_ _2p1 2p1 161,_ 0 0 0,_ 0 0 0,:_ 2.00p1 2.00p1 0"_

pd (gXM(cn0"_ 1;OnD3"_ 1)]) pTtotal nCpts''

pd 'pdf'
pd 'show'
```

 ... energy density ...

```load 'plot'

gXm0All=: _   _  _  _"_
gXm1All=: 1  _1  1  1"_

nCpts=:_ _2p1 3p1 161,_ 0 0 0,_ 0 0 0,:_ 0.00p1 0.00p1 0"_

pd (gXM(cn0"_ 1;OnD0"_ 1)]) pTtotal nCpts''                  NB. ... blue ...

nCpts=:_ _2p1 3p1 161,_ 0 0 0,_ 0 0 0,:_ 0.15p1 0.15p1 0"_

pd (gXM(cn0"_ 1;OnD0"_ 1)]) pTtotal nCpts''                  NB. ... red ...

nCpts=:_ _2p1 3p1 161,_ 0 0 0,_ 0 0 0,:_ 0.50p1 0.50p1 0"_

pd (gXM(cn0"_ 1;OnD0"_ 1)]) pTtotal nCpts''                  NB. ... green ...

nCpts=:_ _2p1 3p1 161,_ 0 0 0,_ 0 0 0,:_ 0.85p1 0.85p1 0"_

pd (gXM(cn0"_ 1;OnD0"_ 1)]) pTtotal nCpts''                  NB. ... purple ...

nCpts=:_ _2p1 3p1 161,_ 0 0 0,_ 0 0 0,:_ 1.00p1 1.00p1 0"_

pd (gXM(cn0"_ 1;OnD0"_ 1)]) pTtotal nCpts''                  NB. ... pink ...

pd 'pdf'
pd 'show'
```

```load 'plot'

gXm0All=: _   _  _  _"_
gXm1All=: 1  _1  1  1"_

nCpts=:_ _1p1 4p1 161,_ 0 0 0,_ 0 0 0,:_ 1.00p1 1.00p1 0"_

pd (gXM(cn0"_ 1;OnD0"_ 1)]) pTtotal nCpts''                  NB. ... blue ...

nCpts=:_ _1p1 4p1 161,_ 0 0 0,_ 0 0 0,:_ 1.15p1 1.15p1 0"_

pd (gXM(cn0"_ 1;OnD0"_ 1)]) pTtotal nCpts''                  NB. ... red ...

nCpts=:_ _1p1 4p1 161,_ 0 0 0,_ 0 0 0,:_ 1.50p1 1.50p1 0"_

pd (gXM(cn0"_ 1;OnD0"_ 1)]) pTtotal nCpts''                  NB. ... green ...

nCpts=:_ _1p1 4p1 161,_ 0 0 0,_ 0 0 0,:_ 1.85p1 1.85p1 0"_

pd (gXM(cn0"_ 1;OnD0"_ 1)]) pTtotal nCpts''                  NB. ... purple ...

nCpts=:_ _1p1 4p1 161,_ 0 0 0,_ 0 0 0,:_ 2.00p1 2.00p1 0"_

pd (gXM(cn0"_ 1;OnD0"_ 1)]) pTtotal nCpts''                  NB. ... pink ...

pd 'pdf'
pd 'show'
```

 ... energy flow ...

 ... energy flow ...

```load 'plot'

gXm0All=: _   _  _  _"_
gXm1All=: 1  _1  1  1"_

nCpts=:_ _2p1 3p1 161,_ 0 0 0,_ 0 0 0,:_ 0.00p1 0.00p1 0"_

pd (gXM(cn0"_ 1;OnD3"_ 1)]) pTtotal nCpts''                  NB. ... blue ...

nCpts=:_ _2p1 3p1 161,_ 0 0 0,_ 0 0 0,:_ 0.25p1 0.25p1 0"_

pd (gXM(cn0"_ 1;OnD3"_ 1)]) pTtotal nCpts''                  NB. ... red ...

nCpts=:_ _2p1 3p1 161,_ 0 0 0,_ 0 0 0,:_ 0.40p1 0.40p1 0"_

pd (gXM(cn0"_ 1;OnD3"_ 1)]) pTtotal nCpts''                  NB. ... green ...

nCpts=:_ _2p1 3p1 161,_ 0 0 0,_ 0 0 0,:_ 0.60p1 0.60p1 0"_

pd (gXM(cn0"_ 1;OnD3"_ 1)]) pTtotal nCpts''                  NB. ... purple ...

nCpts=:_ _2p1 3p1 161,_ 0 0 0,_ 0 0 0,:_ 0.75p1 0.75p1 0"_

pd (gXM(cn0"_ 1;OnD3"_ 1)]) pTtotal nCpts''                  NB. ... pink ...

nCpts=:_ _2p1 3p1 161,_ 0 0 0,_ 0 0 0,:_ 1.00p1 1.00p1 0"_

pd (gXM(cn0"_ 1;OnD3"_ 1)]) pTtotal nCpts''                  NB. ... dark green ...

pd 'pdf'
pd 'show'
```

```load 'plot'

gXm0All=: _   _  _  _"_
gXm1All=: 1  _1  1  1"_

nCpts=:_ _1p1 4p1 161,_ 0 0 0,_ 0 0 0,:_ 1.00p1 1.00p1 0"_

pd (gXM(cn0"_ 1;OnD3"_ 1)]) pTtotal nCpts''                  NB. ... blue ...

nCpts=:_ _1p1 4p1 161,_ 0 0 0,_ 0 0 0,:_ 1.25p1 1.25p1 0"_

pd (gXM(cn0"_ 1;OnD3"_ 1)]) pTtotal nCpts''                  NB. ... red ...

nCpts=:_ _1p1 4p1 161,_ 0 0 0,_ 0 0 0,:_ 1.40p1 1.40p1 0"_

pd (gXM(cn0"_ 1;OnD3"_ 1)]) pTtotal nCpts''                  NB. ... green ...

nCpts=:_ _1p1 4p1 161,_ 0 0 0,_ 0 0 0,:_ 1.60p1 1.60p1 0"_

pd (gXM(cn0"_ 1;OnD3"_ 1)]) pTtotal nCpts''                  NB. ... purple ...

nCpts=:_ _1p1 4p1 161,_ 0 0 0,_ 0 0 0,:_ 1.75p1 1.75p1 0"_

pd (gXM(cn0"_ 1;OnD3"_ 1)]) pTtotal nCpts''                  NB. ... pink ...

nCpts=:_ _1p1 4p1 161,_ 0 0 0,_ 0 0 0,:_ 2.00p1 2.00p1 0"_

pd (gXM(cn0"_ 1;OnD3"_ 1)]) pTtotal nCpts''                  NB. ... dark green ...

pd 'pdf'
pd 'show'
```

 ... energy density ...

 ... energy density ...

...

... also ...

...

...

... energy density adjusted for velocity ...

... (compare with above) ...

...

```load 'plot'

gXm0All=: _   _  _  _"_
gXm1All=: 1  _1  1  1"_

nCpts=:_ _2p1 3p1 161,_ 0 0 0,_ 0 0 0,:_ 0.00p1 0.00p1 0"_

pd (gXM((cn0+gXm11*cn3)"_ 1;OnD3"_ 1)]) pTtotal nCpts''

nCpts=:_ _2p1 3p1 161,_ 0 0 0,_ 0 0 0,:_ 0.10p1 0.10p1 0"_

pd (gXM((cn0+gXm11*cn3)"_ 1;OnD3"_ 1)]) pTtotal nCpts''

nCpts=:_ _2p1 3p1 161,_ 0 0 0,_ 0 0 0,:_ 0.20p1 0.20p1 0"_

pd (gXM((cn0+gXm11*cn3)"_ 1;OnD3"_ 1)]) pTtotal nCpts''

nCpts=:_ _2p1 3p1 161,_ 0 0 0,_ 0 0 0,:_ 0.30p1 0.30p1 0"_

pd (gXM((cn0+gXm11*cn3)"_ 1;OnD3"_ 1)]) pTtotal nCpts''

nCpts=:_ _2p1 3p1 161,_ 0 0 0,_ 0 0 0,:_ 0.40p1 0.40p1 0"_

pd (gXM((cn0+gXm11*cn3)"_ 1;OnD3"_ 1)]) pTtotal nCpts''

nCpts=:_ _2p1 3p1 161,_ 0 0 0,_ 0 0 0,:_ 0.50p1 0.50p1 0"_

pd (gXM((cn0+gXm11*cn3)"_ 1;OnD3"_ 1)]) pTtotal nCpts''

nCpts=:_ _2p1 3p1 161,_ 0 0 0,_ 0 0 0,:_ 0.60p1 0.60p1 0"_

pd (gXM((cn0+gXm11*cn3)"_ 1;OnD3"_ 1)]) pTtotal nCpts''

nCpts=:_ _2p1 3p1 161,_ 0 0 0,_ 0 0 0,:_ 0.70p1 0.70p1 0"_

pd (gXM((cn0+gXm11*cn3)"_ 1;OnD3"_ 1)]) pTtotal nCpts''

nCpts=:_ _2p1 3p1 161,_ 0 0 0,_ 0 0 0,:_ 0.80p1 0.80p1 0"_

pd (gXM((cn0+gXm11*cn3)"_ 1;OnD3"_ 1)]) pTtotal nCpts''

nCpts=:_ _2p1 3p1 161,_ 0 0 0,_ 0 0 0,:_ 0.90p1 0.90p1 0"_

pd (gXM((cn0+gXm11*cn3)"_ 1;OnD3"_ 1)]) pTtotal nCpts''

nCpts=:_ _2p1 3p1 161,_ 0 0 0,_ 0 0 0,:_ 1.00p1 1.00p1 0"_

pd (gXM((cn0+gXm11*cn3)"_ 1;OnD3"_ 1)]) pTtotal nCpts''

nCpts=:_ _1p1 4p1 161,_ 0 0 0,_ 0 0 0,:_ 1.10p1 1.10p1 0"_

pd (gXM((cn0+gXm11*cn3)"_ 1;OnD3"_ 1)]) pTtotal nCpts''

nCpts=:_ _1p1 4p1 161,_ 0 0 0,_ 0 0 0,:_ 1.20p1 1.20p1 0"_

pd (gXM((cn0+gXm11*cn3)"_ 1;OnD3"_ 1)]) pTtotal nCpts''

nCpts=:_ _1p1 4p1 161,_ 0 0 0,_ 0 0 0,:_ 1.30p1 1.30p1 0"_

pd (gXM((cn0+gXm11*cn3)"_ 1;OnD3"_ 1)]) pTtotal nCpts''

nCpts=:_ _1p1 4p1 161,_ 0 0 0,_ 0 0 0,:_ 1.40p1 1.40p1 0"_

pd (gXM((cn0+gXm11*cn3)"_ 1;OnD3"_ 1)]) pTtotal nCpts''

nCpts=:_ _1p1 4p1 161,_ 0 0 0,_ 0 0 0,:_ 1.50p1 1.50p1 0"_

pd (gXM((cn0+gXm11*cn3)"_ 1;OnD3"_ 1)]) pTtotal nCpts''

nCpts=:_ _1p1 4p1 161,_ 0 0 0,_ 0 0 0,:_ 1.60p1 1.60p1 0"_

pd (gXM((cn0+gXm11*cn3)"_ 1;OnD3"_ 1)]) pTtotal nCpts''

nCpts=:_ _1p1 4p1 161,_ 0 0 0,_ 0 0 0,:_ 1.70p1 1.70p1 0"_

pd (gXM((cn0+gXm11*cn3)"_ 1;OnD3"_ 1)]) pTtotal nCpts''

nCpts=:_ _1p1 4p1 161,_ 0 0 0,_ 0 0 0,:_ 1.80p1 1.80p1 0"_

pd (gXM((cn0+gXm11*cn3)"_ 1;OnD3"_ 1)]) pTtotal nCpts''

nCpts=:_ _1p1 4p1 161,_ 0 0 0,_ 0 0 0,:_ 1.90p1 1.90p1 0"_

pd (gXM((cn0+gXm11*cn3)"_ 1;OnD3"_ 1)]) pTtotal nCpts''

nCpts=:_ _1p1 4p1 161,_ 0 0 0,_ 0 0 0,:_ 2.00p1 2.00p1 0"_

pd (gXM((cn0+gXm11*cn3)"_ 1;OnD3"_ 1)]) pTtotal nCpts''

pd 'pdf'
pd 'show'
```

 ... energy density (adjusted for velocity) ...

Download: File:LTEZ2300C.txt
Download: File:LTEZ2301C.txt
Download: File:LTEZ2302C.txt

Download MoinMoin source: File:TensorExperimentsZ23.ijs