User:Tom Allen/TensorExperimentsR24

...

... ADDENDUM ...

...

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

   gXRAll=:,:@(1 1 2 8 1 1 1 _1 0 0"_)
   
   xCpts=:_ 0.7 1.3 26,_ 0.4p1 0.6p1 26,_ 0 2p1 26,:_ 0 2p1 26"_

   A=:(gXR([((<0      ){])@(2{])"2@Ocndc                         xRef@hkxR03)])                ((0;50"_)pTRandom xCpts)''
   B=:(gXR  ((<0      ){])@(2{])"2@:((0|:[:(gXR([(0{])@Ocn       xRef@hkxR03),:@])D.1])"_ 1)]) ((0;50"_)pTRandom xCpts)''
   
   A (((2^_44);2^_17)qteq[;])"0 B
1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1

   A=:(gXR([((<0 0    ){])@(2{])"3@Ocndcdc                       xRef@hkxR03)])                ((0;50"_)pTRandom xCpts)''
   B=:(gXR  ((<0 0    ){])@(2{])"3@:((0|:[:(gXR([(0{])@Ocndc     xRef@hkxR03),:@])D.1])"_ 1)]) ((0;50"_)pTRandom xCpts)''
   
   A (((2^_44);2^_14)qteq[;])"0 B
1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1

   A=:(gXR([((<0 0 0  ){])@(2{])"4@Ocndcdcdc                     xRef@hkxR03)])                ((0;50"_)pTRandom xCpts)''
   B=:(gXR  ((<0 0 0  ){])@(2{])"4@:((0|:[:(gXR([(0{])@Ocndcdc   xRef@hkxR03),:@])D.1])"_ 1)]) ((0;50"_)pTRandom xCpts)''
   
   A (((2^_44);2^_14)qteq[;])"0 B
1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1

   A=:(gXR([((<0 0 0 0){])@(2{])"5@Ocndcdcdcdc                   xRef@hkxR03)])                ((0;50"_)pTRandom xCpts)''
   B=:(gXR  ((<0 0 0 0){])@(2{])"5@:((0|:[:(gXR([(0{])@Ocndcdcdc xRef@hkxR03),:@])D.1])"_ 1)]) ((0;50"_)pTRandom xCpts)''
   
   A (((2^_44);2^_16)qteq[;])"0 B
1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1

... with ...
...
... a2=2 ...
... a3=7 ...
...
... at the origin the value is zero ...
... but ...
... at the origin the derivative is not zero (therefore not continuous through the origin) ...

load 'plot'

gXRAll=:,:@(1 1 2 7 1 1 1 _1 0 0"_)

xCpts=:_ 0.005 0.014 200,_ 0.5p1 0.5p1 0,_ 0.25p1 0.25p1 0,:_ 0 0 0"_

pd (gXR(cx0"_ 1;[((<0 0 0 0){])@(2{])"5@Ocndcdcdcdc xRef@hkxR03)]) pTtotal xCpts''

pd (0 0.01;0 0.8485)

pd 'pdf'
pd 'show'

... with ...
...
... a2=2 ...
... a3=8 ...
...
... at the origin the value is zero ...
... and ...
... at the origin the derivative is zero (therefore continuous through the origin) ...

load 'plot'

gXRAll=:,:@(1 1 2 8 1 1 1 _1 0 0"_)

xCpts=:_ 0 0.02 200,_ 0.5p1 0.5p1 0,_ 0.25p1 0.25p1 0,:_ 0 0 0"_

pd (gXR(cx0"_ 1;[((<0 0 0 0){])@(2{])"5@Ocndcdcdcdc xRef@hkxR03)]) pTtotal xCpts''

pd 'pdf'
pd 'show'

... an invariant ...

load 'plot'

gXRAll=:,:@(1 1 2 8 1 1 1 _1 0 0"_)

xCpts=:_ 0.001 1.6 200,_ 0.5p1 0.5p1 0,_ 0.25p1 0.25p1 0,:_ 0 0 0"_

pd (gXR(cx0"_ 1;[B2kCEr0 xRef@hkxR03)]) pTtotal xCpts''

pd 'pdf'
pd 'show'

---

Download: File:LTER2400D.txt
Download: File:LTER2401C.txt
Download: File:LTER2402C.txt

---