User:Tom Allen/TensorExperimentsZ13

From J Wiki
Jump to navigation Jump to search


...

... decreasing interval size ...

...

... pdDA pdDB pdDC pdDD pdDE pdDF ...

...


load 'plot'

pd ((2.05p1*cos@cx0"1);2.05p1*sin@cx0"1) pTtotal (_ 0 2p1 160,_ 0 0 0,_ 0 0 0,:_ 0 0 0"_)''

vSgDd2M=.[([(0{])@([(xdz smx"2 1 E1cn)zRef)hkxM),:@(S03xH{])
vSgDd2 =.gXZ(-@(([([(0{])@ch2k xRef@hkxM),:@(S03xH{]))+/"1@((<0 2)|:])@((0|:[)smx])(]*/])@(S02d1L{]))+-@vSgDd2M)]

Initds2=.(gXZ([(0{])@mcv xRef@hkxM),:@Initx)+/^:2@:*(]*/])@Initdx

gXcw  =:  1  0  _  _  _  _"_
gXbAll=: _1  1  1  1  1  0"_

Initx =:  2.05p1 , 0.5p1 ,  1p1  , 0.5p1 "_
Initdx=: (0      , 0     , _0.01 , 0.4809"_)*%@(1"_)

vS00xL   =:Initx
vS01dxL  =:0,0,0,0:
vS02d1L  =:Initd1
vS06ds   =:Initds
vS07count=:(0,18648"_)*1"_
pd vSgDplot pdDA=:vSgD@vS''


vSgDd2M=.[([(0{])@([(xdz smx"2 1 E1cn)zRef)hkxM),:@(S03xH{])
vSgDd2 =.gXZ(-@(([([(0{])@ch2k xRef@hkxM),:@(S03xH{]))+/"1@((<0 2)|:])@((0|:[)smx])(]*/])@(S02d1L{]))+-@vSgDd2M)]

Initds2=.(gXZ([(0{])@mcv xRef@hkxM),:@Initx)+/^:2@:*(]*/])@Initdx

gXcw  =:  1  0  _  _  _  _"_
gXbAll=: _1  1  1  1  1  0"_

Initx =:  2.05p1 , 0.5p1 ,  1p1  , 0.5p1 "_
Initdx=: (0      , 0     , _0.01 , 0.4809"_)*%@(2"_)

vS00xL   =:Initx
vS01dxL  =:0,0,0,0:
vS02d1L  =:Initd1
vS06ds   =:Initds
vS07count=:(0,18648"_)*2"_
pd vSgDplot pdDB=:vSgD@vS''


vSgDd2M=.[([(0{])@([(xdz smx"2 1 E1cn)zRef)hkxM),:@(S03xH{])
vSgDd2 =.gXZ(-@(([([(0{])@ch2k xRef@hkxM),:@(S03xH{]))+/"1@((<0 2)|:])@((0|:[)smx])(]*/])@(S02d1L{]))+-@vSgDd2M)]

Initds2=.(gXZ([(0{])@mcv xRef@hkxM),:@Initx)+/^:2@:*(]*/])@Initdx

gXcw  =:  1  0  _  _  _  _"_
gXbAll=: _1  1  1  1  1  0"_

Initx =:  2.05p1 , 0.5p1 ,  1p1  , 0.5p1 "_
Initdx=: (0      , 0     , _0.01 , 0.4809"_)*%@(4"_)

vS00xL   =:Initx
vS01dxL  =:0,0,0,0:
vS02d1L  =:Initd1
vS06ds   =:Initds
vS07count=:(0,18648"_)*4"_
pd vSgDplot pdDC=:vSgD@vS''


vSgDd2M=.[([(0{])@([(xdz smx"2 1 E1cn)zRef)hkxM),:@(S03xH{])
vSgDd2 =.gXZ(-@(([([(0{])@ch2k xRef@hkxM),:@(S03xH{]))+/"1@((<0 2)|:])@((0|:[)smx])(]*/])@(S02d1L{]))+-@vSgDd2M)]

Initds2=.(gXZ([(0{])@mcv xRef@hkxM),:@Initx)+/^:2@:*(]*/])@Initdx

gXcw  =:  1  0  _  _  _  _"_
gXbAll=: _1  1  1  1  1  0"_

Initx =:  2.05p1 , 0.5p1 ,  1p1  , 0.5p1 "_
Initdx=: (0      , 0     , _0.01 , 0.4809"_)*%@(8"_)

vS00xL   =:Initx
vS01dxL  =:0,0,0,0:
vS02d1L  =:Initd1
vS06ds   =:Initds
vS07count=:(0,18648"_)*8"_
pd vSgDplot pdDD=:vSgD@vS''


vSgDd2M=.[([(0{])@([(xdz smx"2 1 E1cn)zRef)hkxM),:@(S03xH{])
vSgDd2 =.gXZ(-@(([([(0{])@ch2k xRef@hkxM),:@(S03xH{]))+/"1@((<0 2)|:])@((0|:[)smx])(]*/])@(S02d1L{]))+-@vSgDd2M)]

Initds2=.(gXZ([(0{])@mcv xRef@hkxM),:@Initx)+/^:2@:*(]*/])@Initdx

gXcw  =:  1  0  _  _  _  _"_
gXbAll=: _1  1  1  1  1  0"_

Initx =:  2.05p1 , 0.5p1 ,  1p1  , 0.5p1 "_
Initdx=: (0      , 0     , _0.01 , 0.4809"_)*%@(16"_)

vS00xL   =:Initx
vS01dxL  =:0,0,0,0:
vS02d1L  =:Initd1
vS06ds   =:Initds
vS07count=:(0,18648"_)*16"_
pd vSgDplot pdDE=:vSgD@vS''


vSgDd2M=.[([(0{])@([(xdz smx"2 1 E1cn)zRef)hkxM),:@(S03xH{])
vSgDd2 =.gXZ(-@(([([(0{])@ch2k xRef@hkxM),:@(S03xH{]))+/"1@((<0 2)|:])@((0|:[)smx])(]*/])@(S02d1L{]))+-@vSgDd2M)]

Initds2=.(gXZ([(0{])@mcv xRef@hkxM),:@Initx)+/^:2@:*(]*/])@Initdx

gXcw  =:  1  0  _  _  _  _"_
gXbAll=: _1  1  1  1  1  0"_

Initx =:  2.05p1 , 0.5p1 ,  1p1  , 0.5p1 "_
Initdx=: (0      , 0     , _0.01 , 0.4809"_)*%@(32"_)

vS00xL   =:Initx
vS01dxL  =:0,0,0,0:
vS02d1L  =:Initd1
vS06ds   =:Initds
vS07count=:(0,18648"_)*32"_
pd vSgDplot pdDF=:vSgD@vS''

pd 'aspect 1'
pd 'pdf'
pd 'show'


pdDA arrayPut   'pdDA'
     arrayClose 'pdDA'
pdDB arrayPut   'pdDB'
     arrayClose 'pdDB'
pdDC arrayPut   'pdDC'
     arrayClose 'pdDC'
pdDD arrayPut   'pdDD'
     arrayClose 'pdDD'
pdDE arrayPut   'pdDE'
     arrayClose 'pdDE'
pdDF arrayPut   'pdDF'
     arrayClose 'pdDF'


NB. ... establish arrays if necessary ...

                      arrayOpen  'pdDA'
pdDA=:     18649  8 4 arrayGet   'pdDA'   NB. ... shape required ...
                      arrayClose 'pdDA'

                      arrayOpen  'pdDB'
pdDB=:     37297  8 4 arrayGet   'pdDB'   NB. ... shape required ...
                      arrayClose 'pdDB'

                      arrayOpen  'pdDC'
pdDC=:     74593  8 4 arrayGet   'pdDC'   NB. ... shape required ...
                      arrayClose 'pdDC'

                      arrayOpen  'pdDD'
pdDD=:     149185 8 4 arrayGet   'pdDD'   NB. ... shape required ...
                      arrayClose 'pdDD'

                      arrayOpen  'pdDE'
pdDE=:     298369 8 4 arrayGet   'pdDE'   NB. ... shape required ...
                      arrayClose 'pdDE'

                      arrayOpen  'pdDF'
pdDF=:     596737 8 4 arrayGet   'pdDF'   NB. ... shape required ...
                      arrayClose 'pdDF'


load 'plot'

pd ((2.05p1*cos@cx0"1);2.05p1*sin@cx0"1) pTtotal (_ 0 2p1 160,_ 0 0 0,_ 0 0 0,:_ 0 0 0"_)''

pd vSgDplot pdDF

pd 'aspect 1'
pd 'pdf'
pd 'show'


Zgraphic1300.jpg
... plot shows pdDF only ...


load 'plot'

pd (cx3"1;-@cx2"1)@:((S00xL{])"2) pdDA
pd (cx3"1;-@cx2"1)@:((S00xL{])"2) pdDB
pd (cx3"1;-@cx2"1)@:((S00xL{])"2) pdDC
pd (cx3"1;-@cx2"1)@:((S00xL{])"2) pdDD
pd (cx3"1;-@cx2"1)@:((S00xL{])"2) pdDE
pd (cx3"1;-@cx2"1)@:((S00xL{])"2) pdDF

pd 'pdf'
pd 'show'


Zgraphic1301.jpg
... x3 and x2 ...


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

   $pdDF
596737 8 4

NB. ... x1 ...

   *./ (0=]) (0.50p1-])   cx1"1@:((S00xL{])"2) pdDF
1

NB. ... x2 decreases continuously ...

   *./ (0<]) (2-/\])      cx2"1@:((S00xL{])"2) pdDF
1

NB. ... x3 ...

NB. ... increases continuously ...

   *./ (0>]) (2-/\])      cx3"1@:((S00xL{])"2) pdDF
1

NB. ... min/max interval ...

   (<./@],>./@])@(2-~/\]) cx3"1@:((S00xL{])"2) pdDF
0.0147052 0.0154749

NB. ... last value ...

   {:                     cx3"1@:((S00xL{])"2) pdDF
9010.36



Download MoinMoin source: File:TensorExperimentsZ13.ijs