User:Tom Allen/TensorExperimentsR14
Jump to navigation
Jump to search
NB. ... script experimentalR.ijs ... NB. ... these verbs are valid only in Cartesian coordinates ... rfEEMcn0=:-@FEMcvcv30 rfEEMcn1=:-@FEMcvcv31 rfEEMcn2=:-@FEMcvcv32 rfEEMcn=:(0|:])@(rfEEMcn0,rfEEMcn1,:rfEEMcn2)
NB. ... script experimentalR.ijs ... NB. ... these verbs are valid only in Cartesian coordinates ... rfPEMcn0=:(rfEEMcn1*HEMcn2)+-@(rfEEMcn2*HEMcn1) rfPEMcn1=:(rfEEMcn2*HEMcn0)+-@(rfEEMcn0*HEMcn2) rfPEMcn2=:(rfEEMcn0*HEMcn1)+-@(rfEEMcn1*HEMcn0) rfPEMcn=:(0|:])@(rfPEMcn0,rfPEMcn1,:rfPEMcn2)
NB. ... script experimentalR.ijs ... NB. ... these verbs are valid only in Cartesian coordinates ... rfJEMcn =:+/"1@((<1 2)|:])"3@(-@FEMcvcvdn) rfJEMcndn=:+/"1@((<1 2)|:])"4@(-@FEMcvcvdndn)
NB. ... execute (ijx) ... gXRAll=:,:@(1 1 2 8 1 1 1 _1 0 0"_) xCpts=:_ 0.7 1.3 26,_ 0.4p1 0.6p1 26,_ 0 2p1 26,:_ 0 2p1 26"_ NB. ... for Maxwell's equations ... *./^:3"3 (0=]) ((2^_44)qtsz]) (gXR([(]+0|:])"2@:((0 1|:])"3)@(+/"1@((<0 1)|:])"5@R)xRef@hkxR03)]) ((0;50"_)pTRandom xCpts)'' 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 *./^:3"3 (0=]) ((2^_44)qtsz]) (]+(0|:])+0 1|:])"3 (gXR([+/"1@((<0 1)|:])"5@B2kdc yRef@hkxR03)]) ((0;50"_)pTRandom xCpts)'' 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
load 'plot' gXRAll=:,:@(0.05 1 2 8 1 1 1 _0.1 0 0"_) xCpts=:_ 0.4 2.0 100,_ 0.5p1 0.5p1 0,_ 0.25p1 0.25p1 0,:(_ 0.00p1 0.00p1 0"_)*10"_ pd (gXR(cx0"_ 1;[(0{])"1@(}:"1@}:"2@xdy smx"2 1 rfEEMcn)yRef@hkxR03)]) pTtotal xCpts'' xCpts=:_ 0.4 2.0 100,_ 0.5p1 0.5p1 0,_ 0.25p1 0.25p1 0,:(_ 0.25p1 0.25p1 0"_)*10"_ pd (gXR(cx0"_ 1;[(0{])"1@(}:"1@}:"2@xdy smx"2 1 rfEEMcn)yRef@hkxR03)]) pTtotal xCpts'' xCpts=:_ 0.4 2.0 100,_ 0.5p1 0.5p1 0,_ 0.25p1 0.25p1 0,:(_ 0.40p1 0.40p1 0"_)*10"_ pd (gXR(cx0"_ 1;[(0{])"1@(}:"1@}:"2@xdy smx"2 1 rfEEMcn)yRef@hkxR03)]) pTtotal xCpts'' xCpts=:_ 0.4 2.0 100,_ 0.5p1 0.5p1 0,_ 0.25p1 0.25p1 0,:(_ 0.50p1 0.50p1 0"_)*10"_ pd (gXR(cx0"_ 1;[(0{])"1@(}:"1@}:"2@xdy smx"2 1 rfEEMcn)yRef@hkxR03)]) pTtotal xCpts'' xCpts=:_ 0.4 2.0 100,_ 0.5p1 0.5p1 0,_ 0.25p1 0.25p1 0,:(_ 0.60p1 0.60p1 0"_)*10"_ pd (gXR(cx0"_ 1;[(0{])"1@(}:"1@}:"2@xdy smx"2 1 rfEEMcn)yRef@hkxR03)]) pTtotal xCpts'' xCpts=:_ 0.4 2.0 100,_ 0.5p1 0.5p1 0,_ 0.25p1 0.25p1 0,:(_ 0.75p1 0.75p1 0"_)*10"_ pd (gXR(cx0"_ 1;[(0{])"1@(}:"1@}:"2@xdy smx"2 1 rfEEMcn)yRef@hkxR03)]) pTtotal xCpts'' xCpts=:_ 0.4 2.0 100,_ 0.5p1 0.5p1 0,_ 0.25p1 0.25p1 0,:(_ 1.00p1 1.00p1 0"_)*10"_ pd (gXR(cx0"_ 1;[(0{])"1@(}:"1@}:"2@xdy smx"2 1 rfEEMcn)yRef@hkxR03)]) pTtotal xCpts'' pd (((0.8723"_)@]"0);]"0) pTsteps _ _0.00004 0.00004 160 pd (((1.0912"_)@]"0);]"0) pTsteps _ _0.00004 0.00004 160 pd (((1.3202"_)@]"0);]"0) pTsteps _ _0.00004 0.00004 160 pd 'pdf' pd 'show'
load 'plot' gXRAll=:,:@(0.05 1 2 8 1 1 1 _0.1 0 0"_) xCpts=:_ 0.4 2.0 100,_ 0.5p1 0.5p1 0,_ 0.25p1 0.25p1 0,:(_ 0.00p1 0.00p1 0"_)*10"_ pd (gXR(cx0"_ 1;[(3{])"1@(xdy smx"2 1 rfJEMcn)yRef@hkxR03)]) pTtotal xCpts'' xCpts=:_ 0.4 2.0 100,_ 0.5p1 0.5p1 0,_ 0.25p1 0.25p1 0,:(_ 0.25p1 0.25p1 0"_)*10"_ pd (gXR(cx0"_ 1;[(3{])"1@(xdy smx"2 1 rfJEMcn)yRef@hkxR03)]) pTtotal xCpts'' xCpts=:_ 0.4 2.0 100,_ 0.5p1 0.5p1 0,_ 0.25p1 0.25p1 0,:(_ 0.40p1 0.40p1 0"_)*10"_ pd (gXR(cx0"_ 1;[(3{])"1@(xdy smx"2 1 rfJEMcn)yRef@hkxR03)]) pTtotal xCpts'' xCpts=:_ 0.4 2.0 100,_ 0.5p1 0.5p1 0,_ 0.25p1 0.25p1 0,:(_ 0.50p1 0.50p1 0"_)*10"_ pd (gXR(cx0"_ 1;[(3{])"1@(xdy smx"2 1 rfJEMcn)yRef@hkxR03)]) pTtotal xCpts'' xCpts=:_ 0.4 2.0 100,_ 0.5p1 0.5p1 0,_ 0.25p1 0.25p1 0,:(_ 0.60p1 0.60p1 0"_)*10"_ pd (gXR(cx0"_ 1;[(3{])"1@(xdy smx"2 1 rfJEMcn)yRef@hkxR03)]) pTtotal xCpts'' xCpts=:_ 0.4 2.0 100,_ 0.5p1 0.5p1 0,_ 0.25p1 0.25p1 0,:(_ 0.75p1 0.75p1 0"_)*10"_ pd (gXR(cx0"_ 1;[(3{])"1@(xdy smx"2 1 rfJEMcn)yRef@hkxR03)]) pTtotal xCpts'' xCpts=:_ 0.4 2.0 100,_ 0.5p1 0.5p1 0,_ 0.25p1 0.25p1 0,:(_ 1.00p1 1.00p1 0"_)*10"_ pd (gXR(cx0"_ 1;[(3{])"1@(xdy smx"2 1 rfJEMcn)yRef@hkxR03)]) pTtotal xCpts'' pd (((0.8723"_)@]"0);]"0) pTsteps _ _0.0002 0.0002 160 pd (((1.0912"_)@]"0);]"0) pTsteps _ _0.0002 0.0002 160 pd (((1.3202"_)@]"0);]"0) pTsteps _ _0.0002 0.0002 160 pd 'pdf' pd 'show'
load 'plot' gXRAll=:,:@(0.05 1 2 8 1 1 1 _0.1 0 0"_) xCpts=:_ 1.0912 1.0912 0,_ 0.5p1 0.5p1 0,_ 0.25p1 2.25p1 100,:_ 0p1 0p1 0"_ pd (gXR(cx2"_ 1;[(0{])"1@(}:"1@}:"2@xdy smx"2 1 rfEEMcn)yRef@hkxR03)]) pTtotal xCpts'' xCpts=:_ 1.0912 1.0912 0,_ 0.5p1 0.5p1 0,_ 0.25p1 2.25p1 100,:_ 0p1 0p1 0"_ pd (gXR(cx2"_ 1;[(1{])"1@(}:"1@}:"2@xdy smx"2 1 rfEEMcn)yRef@hkxR03)]) pTtotal xCpts'' xCpts=:_ 1.0912 1.0912 0,_ 0.5p1 0.5p1 0,_ 0.25p1 2.25p1 100,:_ 0p1 0p1 0"_ pd (gXR(cx2"_ 1;[(2{])"1@(}:"1@}:"2@xdy smx"2 1 rfEEMcn)yRef@hkxR03)]) pTtotal xCpts'' pd 'pdf' pd 'show'
load 'plot' gXRAll=:,:@(0.05 1 2 8 1 1 1 _0.1 0 0"_) xCpts=:_ 1.0912 1.0912 0,_ 0.5p1 0.5p1 0,_ 0.25p1 2.25p1 100,:_ 0p1 0p1 0"_ pd (gXR(cx2"_ 1;[(0{])"1@(}:"1@}:"2@xdy smx"2 1 HEMcn)yRef@hkxR03)]) pTtotal xCpts'' xCpts=:_ 1.0912 1.0912 0,_ 0.5p1 0.5p1 0,_ 0.25p1 2.25p1 100,:_ 0p1 0p1 0"_ pd (gXR(cx2"_ 1;[(1{])"1@(}:"1@}:"2@xdy smx"2 1 HEMcn)yRef@hkxR03)]) pTtotal xCpts'' xCpts=:_ 1.0912 1.0912 0,_ 0.5p1 0.5p1 0,_ 0.25p1 2.25p1 100,:_ 0p1 0p1 0"_ pd (gXR(cx2"_ 1;[(2{])"1@(}:"1@}:"2@xdy smx"2 1 HEMcn)yRef@hkxR03)]) pTtotal xCpts'' pd 'pdf' pd 'show'
load 'plot' gXRAll=:,:@(0.05 1 2 8 1 1 1 _0.1 0 0"_) xCpts=:_ 1.0912 1.0912 0,_ 0.5p1 0.5p1 0,_ 0.25p1 2.25p1 100,:_ 0p1 0p1 0"_ pd (gXR(cx2"_ 1;[(0{])"1@(}:"1@}:"2@xdy smx"2 1 rfPEMcn)yRef@hkxR03)]) pTtotal xCpts'' xCpts=:_ 1.0912 1.0912 0,_ 0.5p1 0.5p1 0,_ 0.25p1 2.25p1 100,:_ 0p1 0p1 0"_ pd (gXR(cx2"_ 1;[(1{])"1@(}:"1@}:"2@xdy smx"2 1 rfPEMcn)yRef@hkxR03)]) pTtotal xCpts'' xCpts=:_ 1.0912 1.0912 0,_ 0.5p1 0.5p1 0,_ 0.25p1 2.25p1 500,:_ 0p1 0p1 0"_ pd (gXR(cx2"_ 1;[(2{])"1@(}:"1@}:"2@xdy smx"2 1 rfPEMcn)yRef@hkxR03)]) pTtotal xCpts'' pd 'pdf' pd 'show'
load 'plot' gXRAll=:,:@(0.05 1 2 8 1 1 1 _0.1 0 0"_) xCpts=:_ 1.0912 1.0912 0,_ 0.5p1 0.5p1 0,_ 0.25p1 2.25p1 100,:_ 0p1 0p1 0"_ pd (gXR(cx2"_ 1;[(3{])"1@(xdy smx"2 1 rfJEMcn)yRef@hkxR03)]) pTtotal xCpts'' xCpts=:_ 1.0912 1.0912 0,_ 0.5p1 0.5p1 0,_ 0.25p1 2.25p1 100,:_ 2p1 2p1 0"_ pd (gXR(cx2"_ 1;[(3{])"1@(xdy smx"2 1 rfJEMcn)yRef@hkxR03)]) pTtotal xCpts'' xCpts=:_ 1.0912 1.0912 0,_ 0.5p1 0.5p1 0,_ 0.25p1 2.25p1 100,:_ 4p1 4p1 0"_ pd (gXR(cx2"_ 1;[(3{])"1@(xdy smx"2 1 rfJEMcn)yRef@hkxR03)]) pTtotal xCpts'' pd 'pdf' pd 'show'
Download: File:LTER1400C.txt