User:Tom Allen/TensorExperimentsZ20

gXsAll=. _1  1  1  1  1  0"_

Initu =.  _2.05p1 , 0       , 0 , 0.5p1"_
Initdu=. ( 0      , 0.01331 , 0 , 0.1  "_)*%@(1"_)

Initds=.%:@|@((gXS([(0{])@mcv uRef@hkuR),:@Initu)+/^:2@:*(]*/])@Initdu)

vS00xL   =.Initu
vS01dxL  =.0,0,0,0:
vS02d1L  =.Initdu%Initds
vS06ds   =.Initds
vS07count=.(0,3724"_)*1"_

vSgDd2=.gXS([-@([(0{])@E1cn uRef@hkuR),:@(S03xH{]))]

pdDA=:vSgD@vS''

NB. ... pdDA to disk file ...

pdDA arrayPut   'pdDA'
     arrayClose 'pdDA'


Initu =.  _2.20p1 , 0       , 0 , 0.5p1"_
Initdu=. ( 0      , 0.01375 , 0 , 0.1  "_)*%@(1"_)

Initds=.%:@|@((gXS([(0{])@mcv uRef@hkuR),:@Initu)+/^:2@:*(]*/])@Initdu)

vS00xL   =.Initu
vS01dxL  =.0,0,0,0:
vS02d1L  =.Initdu%Initds
vS06ds   =.Initds
vS07count=.(0,3362"_)*1"_

vSgDd2=.gXS([-@([(0{])@E1cn uRef@hkuR),:@(S03xH{]))]

pdDB=:vSgD@vS''

NB. ... pdDB to disk file ...

pdDB arrayPut   'pdDB'
     arrayClose 'pdDB'

load 'plot'

pd ((2.05p1*cos@]"0);2.05p1*sin@]"0) pTsteps _ 0p1 2p1 160

pd <"1@(0|:])@:((cu0,cu1)@(S00xL{])"2) pdDA

pd ((2.20p1*cos@]"0);2.20p1*sin@]"0) pTsteps _ 0p1 2p1 160

pd <"1@(0|:])@:((cu0,cu1)@(S00xL{])"2) pdDB

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

... waves ...

---

Download MoinMoin source: File:TensorExperimentsZ20.ijs

---