User:Tom Allen/TensorExperimentsZ11
Jump to navigation
Jump to search
...
... J1cn ...
...
...
... one sign ...
...
load 'plot'
vSgDd2M=.[([(0{])@([(xdz smx"2 1 J1cn)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 =: 2p1 , 0.5p1 , 1p1 , 0.5p1"_
Initdx=: 0 , 0 , _0.01 , 0 "_
vS00xL =:Initx
vS01dxL =:0,0,0,0:
vS02d1L =:Initd1
vS06ds =:Initds
vS07count=:0,628"_
pd vSgDplot pdDA=:vSgD@vS''
Initx =: 3p1 , 0.5p1 , 1p1 , 0.5p1"_
Initdx=: 0 , 0 , _0.01 , 0 "_
vS00xL =:Initx
vS01dxL =:0,0,0,0:
vS02d1L =:Initd1
vS06ds =:Initds
vS07count=:0,628"_
pd vSgDplot pdDB=:vSgD@vS''
Initx =: 4p1 , 0.5p1 , 1p1 , 0.5p1"_
Initdx=: 0 , 0 , _0.01 , 0 "_
vS00xL =:Initx
vS01dxL =:0,0,0,0:
vS02d1L =:Initd1
vS06ds =:Initds
vS07count=:0,628"_
pd vSgDplot pdDC=:vSgD@vS''
pd 'aspect 1'
pd 'pdf'
pd 'show'
...
NB. ... execute (ijx) ...
NB. ... for pdDA ...
$pdDA
629 8 4
NB. ... x0 ...
*./ (0=]) (2p1-]) cx0"1@:((S00xL{])"2) pdDA
1
NB. ... x1 ...
*./ (0=]) (0.50p1-]) cx1"1@:((S00xL{])"2) pdDA
1
NB. ... x2 decreases continuously ...
*./ (0<]) (2-/\]) cx2"1@:((S00xL{])"2) pdDA
1
NB. ... x3 ...
*./ (0=]) (0.50p1-]) cx3"1@:((S00xL{])"2) pdDA
1
NB. ... for pdDB ...
$pdDB
629 8 4
NB. ... etc ...
...
... the other sign ...
...
load 'plot'
vSgDd2M=.[([(0{])@([(xdz smx"2 1 J1cn)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 =: 2p1 , 0.5p1 , 1p1 , 0"_
Initdx=: 0 , 0 , _0.01 , 0"_
vS00xL =:Initx
vS01dxL =:0,0,0,0:
vS02d1L =:Initd1
vS06ds =:Initds
vS07count=:0,628"_
pd vSgDplot pdDD=:vSgD@vS''
Initx =: 3p1 , 0.5p1 , 1p1 , 0"_
Initdx=: 0 , 0 , _0.01 , 0"_
vS00xL =:Initx
vS01dxL =:0,0,0,0:
vS02d1L =:Initd1
vS06ds =:Initds
vS07count=:0,628"_
pd vSgDplot pdDE=:vSgD@vS''
Initx =: 4p1 , 0.5p1 , 1p1 , 0"_
Initdx=: 0 , 0 , _0.01 , 0"_
vS00xL =:Initx
vS01dxL =:0,0,0,0:
vS02d1L =:Initd1
vS06ds =:Initds
vS07count=:0,628"_
pd vSgDplot pdDF=:vSgD@vS''
pd 'aspect 1'
pd 'pdf'
pd 'show'
...
NB. ... execute (ijx) ...
NB. ... for pdDD ...
$pdDD
629 8 4
NB. ... x0 ...
*./ (0=]) (2p1-]) cx0"1@:((S00xL{])"2) pdDD
1
NB. ... x1 ...
*./ (0=]) (0.50p1-]) cx1"1@:((S00xL{])"2) pdDD
1
NB. ... x2 decreases continuously ...
*./ (0<]) (2-/\]) cx2"1@:((S00xL{])"2) pdDD
1
NB. ... x3 ...
*./ (0=]) cx3"1@:((S00xL{])"2) pdDD
1
NB. ... for pdDE ...
$pdDE
629 8 4
NB. ... etc ...
...
... stability (orientation and position) ...
...
load 'plot'
vSgDd2M=.[([(0{])@([(xdz smx"2 1 J1cn)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 =: 1.92p1 , 0.5p1 , 1p1 , 0.5p1"_
Initdx=: (0 , 0 , _0.01 , 0 "_)*%@(10"_)
vS00xL =:Initx
vS01dxL =:0,0,0,0:
vS02d1L =:Initd1
vS06ds =:Initds
vS07count=:(0,3140"_)*10"_
pd vSgDplot pdDA=:vSgD@vS''
pd 'aspect 1'
pd 'pdf'
pd 'show'
...
NB. ... execute (ijx) ...
$pdDA
31401 8 4
NB. ... x0 ...
A=: cx0"1@:((S00xL{])"2) pdDA
(2p1-<./) A
0.261011
(>./-2p1"_) A
0.129868
NB. ... x1 ...
*./ (0=]) (0.50p1-]) cx1"1@:((S00xL{])"2) pdDA
1
NB. ... x2 decreases continuously ...
*./ (0<]) (2-/\]) cx2"1@:((S00xL{])"2) pdDA
1
NB. ... x3 ...
*./ (0=]) (0.50p1-]) cx3"1@:((S00xL{])"2) pdDA
1
Download MoinMoin source: File:TensorExperimentsZ11.ijs