User:Tom Allen/ta/AddendumX01P00

From J Wiki
Jump to: navigation, search


NB. ... script experimentalX01.ijs ...

load jpath '~temp\experimentalT00.ijs'
load jpath '~temp\experimentalT01.ijs'
load jpath '~temp\experimentalT02.ijs'
load jpath '~temp\experimentalT03.ijs'
load jpath '~temp\experimentalU00.ijs'
load jpath '~temp\experimentalU01.ijs'
load jpath '~temp\experimentalU02.ijs'
load jpath '~temp\experimentalU03.ijs'
load jpath '~temp\experimentalU04.ijs'
load jpath '~temp\experimentalV00.ijs'
load jpath '~temp\experimentalW00.ijs'
load jpath '~temp\experimentalW01.ijs'


... check so far ...


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

   zS0=:,:@(0 0 0 0 0 0 8 _8 _5"_)

   xCpts=:_ 0.7 1.3 26,_ 0.4p1 0.6p1 26,_ 0 2p1 26,:_ 0 2p1 26"_

   (zS0(zSbQ17-:>@(0{])@zSbQ17XX)"_ 1])                                                        ((2^_44)qtsz])@(((0{])+5*3{]),1 2 3{])@zfyc"1@((0;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

   (zS0(>@(1{])@zSbQ17XX"_ 1(((2^_44);2^_12)qteq[;])"1(0|:[:(zS0>@(0{])@zSbQ17XX])D.1])"_ 1)]) ((2^_44)qtsz])@(((0{])+5*3{]),1 2 3{])@zfyc"1@((0;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

   (zS0(>@(2{])@zSbQ17XX"_ 1(((2^_44);2^_10)qteq[;])"2(0|:[:(zS0>@(1{])@zSbQ17XX])D.1])"_ 1)]) ((2^_44)qtsz])@(((0{])+5*3{]),1 2 3{])@zfyc"1@((0;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

   (zS0(>@(3{])@zSbQ17XX"_ 1(((2^_44);2^_09)qteq[;])"3(0|:[:(zS0>@(2{])@zSbQ17XX])D.1])"_ 1)]) ((2^_44)qtsz])@(((0{])+5*3{]),1 2 3{])@zfyc"1@((0;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

   (zS0(>@(4{])@zSbQ17XX"_ 1(((2^_44);2^_07)qteq[;])"4(0|:[:(zS0>@(3{])@zSbQ17XX])D.1])"_ 1)]) ((2^_44)qtsz])@(((0{])+5*3{]),1 2 3{])@zfyc"1@((0;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) ...

NB. ... for the calculation of derivatives using 'D.1' ...

   zS0=:,:@(0 0 0 0 0 0 8 _8 0"_)

NB. ... for y1=0 ...
   
   xCpts=:_ 1 1 0,_ 0.5p1 0.5p1 0,_ 1p1 1p1 0,:_ 1 1 0"_
   
   ((2^_44)qtsz])                               zfyc"1@((0;1"_)pTRandom xCpts)''    NB. ... y coordinates ...
_1 0 0 1
   
   (zS0>@(0{])@zSbQ17XX"_ 1])                   zfyc"1@((0;1"_)pTRandom xCpts)''    NB. ... value at _1 0 0 1 ...
_1.13486
   
   ((2^_44)qtsz]) (zS0>@(1{])@zSbQ17XX"_ 1])    zfyc"1@((0;1"_)pTRandom xCpts)''    NB. ... first derivatives ...
0 _6.75087 0 _6.75087
   
   (zS0(0|:[:(zS0>@(0{])@zSbQ17XX])D.1])"_ 1])  zfyc"1@((0;1"_)pTRandom xCpts)''    NB. ... note derivative w.r.t. y1 calculated to be 0 ('D.1' error) ...
0 0 0 _6.75086
 
NB. ... for y0=0 ...
   
   xCpts=:_ 1 1 0,_ 0.5p1 0.5p1 0,_ 0.5p1 0.5p1 0,:_ 1 1 0"_
   
   ((2^_44)qtsz])                               zfyc"1@((0;1"_)pTRandom xCpts)''    NB. ... y coordinates ...
0 1 0 1

   (zS0>@(0{])@zSbQ17XX"_ 1])                   zfyc"1@((0;1"_)pTRandom xCpts)''    NB. ... value at 0 1 0 1 ...
_1.13486

   ((2^_44)qtsz]) (zS0>@(1{])@zSbQ17XX"_ 1])    zfyc"1@((0;1"_)pTRandom xCpts)''    NB. ... first derivatives ...
_6.75087 0 0 _6.75087

   (zS0(0|:[:(zS0>@(0{])@zSbQ17XX])D.1])"_ 1])  zfyc"1@((0;1"_)pTRandom xCpts)''    NB. ... note derivative w.r.t. y0 calculated to be 0 ('D.1' error) ...
0 0 0 _6.75086


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

NB. ... check zSbQ17 consistent with rfaV31 ...

   gX0=:,:@(0 0 0 0 0 0 8 _8 0 0"_)

   zS0=:,:@(0 0 0 0 0 0 8 _8 0"_)

   xCpts=:_ 0.7 1.3 26,_ 0.4p1 0.6p1 26,_ 0 2p1 26,:_ 0 2p1 26"_

   ((gX0>@(0{])@([rfTT4([zfxXXS zfyc),:rfaV31XX)"_ 1])(((2^_44);2^_44)qteq[;])"0 zS0(([>@(0{])@zSbQ17XX zfyc)"_ 1)]) ((0;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
   
   ((gX0>@(1{])@([rfTT4([zfxXXS zfyc),:rfaV31XX)"_ 1])(((2^_44);2^_43)qteq[;])"1 zS0(([>@(1{])@zSbQ17XX zfyc)"_ 1)]) ((0;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
   
   ((gX0>@(2{])@([rfTT4([zfxXXS zfyc),:rfaV31XX)"_ 1])(((2^_44);2^_41)qteq[;])"2 zS0(([>@(2{])@zSbQ17XX zfyc)"_ 1)]) ((0;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
   
   ((gX0>@(3{])@([rfTT4([zfxXXS zfyc),:rfaV31XX)"_ 1])(((2^_44);2^_42)qteq[;])"3 zS0(([>@(3{])@zSbQ17XX zfyc)"_ 1)]) ((0;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
   
   ((gX0>@(4{])@([rfTT4([zfxXXS zfyc),:rfaV31XX)"_ 1])(((2^_44);2^_41)qteq[;])"4 zS0(([>@(4{])@zSbQ17XX zfyc)"_ 1)]) ((0;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