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'
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'
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