User:Tom Allen/TensorExperimentsZ13
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'
... 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'
... 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