User:Tom Allen/TensorExperimentsR09
Jump to navigation
Jump to search
NB. ... script experimentalR.ijs ... rfhcvXXS=:(1 0 0 0,0 1 0 0,0 0 1 0,:0 0 0 _1"_);((3$4:)$0:);((4$4:)$0:);((5$4:)$0:);((6$4:)$0:);1$_"_ rfhcvXX=:<@:>"1@(0|:])@:(rfhcvXXS"_ 1) hkxR01=:[(],[rfhcvXX ypts)hkxR00 hkyR01=:[(],[rfhcvXX ypts)hkyR00
... tensor mcv ...
... (lower indices on tensors) ...
... (generate Christoffel symbols for Spherical coordinates) ...
NB. ... execute (ijx) ... gXRAll=:,:@(0 0 0 0 0 0 0 0 0 0"_) xCpts=:_ 1 3 26,_ 0.1p1 0.9p1 26,_ 0 2p1 26,:_ 0 20 26"_ (gXR([(hcvD4(((2^_44);2^_44)qteq[;])"2 xdy([((0|:[)smx])"2((0|:[)smx])"2)gcv)hkxR01)]) ((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
NB. ... execute (ijx) ...
gXRAll=:,:@(0 0 0 0 0 0 0 0 0 0"_)
xCpts=:_ 1 3 26,_ 0.1p1 0.9p1 26,_ 0 2p1 26,:_ 0 20 26"_
A=:(gXR([gcvdx hkxR01)]) ((0;50"_)pTRandom xCpts)''
B=:(gXR((0|:[:(gXR([(0{])@gcv hkxR01),:@])D.1])"_ 1)]) ((0;50"_)pTRandom xCpts)''
A (((2^_24);2^_20)qteq[;])"3 B
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
NB. ... execute (ijx) ...
gXRAll=:,:@(0 0 0 0 0 0 0 0 0 0"_)
xCpts=:_ 1 3 26,_ 0.1p1 0.9p1 26,_ 0 2p1 26,:_ 0 20 26"_
A=:(gXR([gcvdxdx hkxR01)]) ((0;50"_)pTRandom xCpts)''
B=:(gXR((0|:[:(gXR([(0{])@gcvdx hkxR01),:@])D.1])"_ 1)]) ((0;50"_)pTRandom xCpts)''
A (((2^_24);2^_16)qteq[;])"4 B
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
NB. ... execute (ijx) ...
gXRAll=:,:@(0 0 0 0 0 0 0 0 0 0"_)
xCpts=:_ 1 3 26,_ 0.1p1 0.9p1 26,_ 0 2p1 26,:_ 0 20 26"_
A=:(gXR([gcvdxdxdx hkxR01)]) ((0;50"_)pTRandom xCpts)''
B=:(gXR((0|:[:(gXR([(0{])@gcvdxdx hkxR01),:@])D.1])"_ 1)]) ((0;50"_)pTRandom xCpts)''
A (((2^_23);2^_16)qteq[;])"5 B
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
NB. ... execute (ijx) ...
gXRAll=:,:@(0 0 0 0 0 0 0 0 0 0"_)
xCpts=:_ 1 3 26,_ 0.1p1 0.9p1 26,_ 0 2p1 26,:_ 0 20 26"_
A=:(gXR([gcvdxdxdxdx hkxR01)]) ((0;50"_)pTRandom xCpts)''
B=:(gXR((0|:[:(gXR([(0{])@gcvdxdxdx hkxR01),:@])D.1])"_ 1)]) ((0;50"_)pTRandom xCpts)''
A (((2^_22);2^_16)qteq[;])"6 B
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
NB. ... script experimentalR.ijs ... rfgcvXX=:gcv;gcvdx;gcvdxdx;gcvdxdxdx;gcvdxdxdxdx;[(1$_"_)"_ 1 xpts hkxR02=:[(],rfgcvXX)hkxR01 hkyR02=:[(],rfgcvXX)hkyR01
... tensor mcn ...
... (raise indices on tensors) ...
... (generate Christoffel symbols for Spherical coordinates) ...
NB. ... execute (ijx) ...
gXRAll=:,:@(0 0 0 0 0 0 0 0 0 0"_)
xCpts=:_ 1 3 26,_ 0.1p1 0.9p1 26,_ 0 2p1 26,:_ 0 20 26"_
A=:(gXR([mcndc xRef@(hkxR02,2 6$a:"_))]) ((0;50"_)pTRandom xCpts)''
B=:(gXR((0|:[:(gXR([(0{])@mcn xRef@(hkxR02,2 6$a:"_)),:@])D.1])"_ 1)]) ((0;50"_)pTRandom xCpts)''
A (((2^_25);2^_19)qteq[;])"3 B
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
A=:(gXR([mcndcdc xRef@(hkxR02,2 6$a:"_))]) ((0;50"_)pTRandom xCpts)''
B=:(gXR((0|:[:(gXR([(0{])@mcndc xRef@(hkxR02,2 6$a:"_)),:@])D.1])"_ 1)]) ((0;50"_)pTRandom xCpts)''
A (((2^_22);2^_19)qteq[;])"4 B
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
A=:(gXR([mcndcdcdc xRef@(hkxR02,2 6$a:"_))]) ((0;50"_)pTRandom xCpts)''
B=:(gXR((0|:[:(gXR([(0{])@mcndcdc xRef@(hkxR02,2 6$a:"_)),:@])D.1])"_ 1)]) ((0;50"_)pTRandom xCpts)''
A (((2^_18);2^_18)qteq[;])"5 B
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
... Christoffel symbols of the first kind ...
NB. ... execute (ijx) ...
gXRAll=:,:@(0 0 0 0 0 0 0 0 0 0"_)
xCpts=:_ 1 3 26,_ 0.1p1 0.9p1 26,_ 0 2p1 26,:_ 0 20 26"_
A=:(gXR([ch1kdc xRef@(hkxR02,2 6$a:"_))]) ((0;50"_)pTRandom xCpts)''
B=:(gXR((0|:[:(gXR([(0{])@ch1k xRef@(hkxR02,2 6$a:"_)),:@])D.1])"_ 1)]) ((0;50"_)pTRandom xCpts)''
A (((2^_25);2^_16)qteq[;])"4 B
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
A=:(gXR([ch1kdcdc xRef@(hkxR02,2 6$a:"_))]) ((0;50"_)pTRandom xCpts)''
B=:(gXR((0|:[:(gXR([(0{])@ch1kdc xRef@(hkxR02,2 6$a:"_)),:@])D.1])"_ 1)]) ((0;50"_)pTRandom xCpts)''
A (((2^_24);2^_16)qteq[;])"5 B
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
A=:(gXR([ch1kdcdcdc xRef@(hkxR02,2 6$a:"_))]) ((0;50"_)pTRandom xCpts)''
B=:(gXR((0|:[:(gXR([(0{])@ch1kdcdc xRef@(hkxR02,2 6$a:"_)),:@])D.1])"_ 1)]) ((0;50"_)pTRandom xCpts)''
A (((2^_23);2^_16)qteq[;])"6 B
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
... Christoffel symbols of the second kind ...
NB. ... execute (ijx) ...
gXRAll=:,:@(0 0 0 0 0 0 0 0 0 0"_)
xCpts=:_ 1 3 26,_ 0.1p1 0.9p1 26,_ 0 2p1 26,:_ 0 20 26"_
A=:(gXR([ch2kdc xRef@(hkxR02,2 6$a:"_))]) ((0;50"_)pTRandom xCpts)''
B=:(gXR((0|:[:(gXR([(0{])@ch2k xRef@(hkxR02,2 6$a:"_)),:@])D.1])"_ 1)]) ((0;50"_)pTRandom xCpts)''
A (((2^_25);2^_16)qteq[;])"4 B
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
A=:(gXR([ch2kdcdc xRef@(hkxR02,2 6$a:"_))]) ((0;50"_)pTRandom xCpts)''
B=:(gXR((0|:[:(gXR([(0{])@ch2kdc xRef@(hkxR02,2 6$a:"_)),:@])D.1])"_ 1)]) ((0;50"_)pTRandom xCpts)''
A (((2^_23);2^_16)qteq[;])"5 B
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
A=:(gXR([ch2kdcdcdc xRef@(hkxR02,2 6$a:"_))]) ((0;50"_)pTRandom xCpts)''
B=:(gXR((0|:[:(gXR([(0{])@ch2kdcdc xRef@(hkxR02,2 6$a:"_)),:@])D.1])"_ 1)]) ((0;50"_)pTRandom xCpts)''
A (((2^_19);2^_16)qteq[;])"6 B
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
... tensor mcn ...
... (a method to calculate) ...
NB. ... script experimentalR.ijs ... rfhcn=:1 0 0 0,0 1 0 0,0 0 1 0,:0 0 0 _1"_
NB. ... script experimentalR.ijs ... rfmcn=:xdy([smx"2 smx"2)rfhcn"_
NB. ... execute (ijx) ... gXRAll=:,:@(0 0 0 0 0 0 0 0 0 0"_) xCpts=:_ 1 3 26,_ 0.1p1 0.9p1 26,_ 0 2p1 26,:_ 0 20 26"_ (gXR([(rfmcn(((2^_44);2^_44)qteq[;])"2 mcn)xRef@(hkxR02,2 6$a:"_))]) ((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
NB. ... script experimentalR.ijs ... rfmcndc0=:(1|:])"3@(xdy(( [)smx ])"2 3 ydx ((0|:[)smx 1|:])"2 3 xdydy((1|:[)smx ])"3 2 rfhcn"_) rfmcndc1=:(0|:])"3@(ydx((0|:[)smx 1|:])"2 3 xdydy((1|:[)smx ])"3 2 xdy (( [)smx ])"2 2 rfhcn"_) rfmcndc=:rfmcndc0+rfmcndc1
NB. ... execute (ijx) ... gXRAll=:,:@(0 0 0 0 0 0 0 0 0 0"_) xCpts=:_ 1 3 26,_ 0.1p1 0.9p1 26,_ 0 2p1 26,:_ 0 20 26"_ (gXR([(rfmcndc(((2^_44);2^_44)qteq[;])"3 mcndc)xRef@(hkxR02,2 6$a:"_))]) ((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
NB. ... script experimentalR.ijs ... rfmcndcdc00=:(1 2|:])"4@(xdy(( [)smx ])"2 4 ydx((0|:[)smx 2|:])"2 4 ydx((0|:[)smx 2|:])"2 4 xdydydy((1|:[)smx ])"4 2 rfhcn"_) rfmcndcdc01=:(1 2|:])"4@( xdy(( [)smx ])"2 4 ydxdx((0|:[)smx 1|:])"3 3 xdydy((1|:[)smx ])"3 2 rfhcn"_) rfmcndcdc02=:(2 0|:])"4@(ydx((0|:[)smx 1|:])"2 4 xdydy((1|:[)smx ])"3 3 ydx((0|:[)smx 1|:])"2 3 xdydy((1|:[)smx ])"3 2 rfhcn"_) rfmcndcdc10=:(0 2|:])"4@(ydx((0|:[)smx 1|:])"2 4 xdydy((1|:[)smx ])"3 3 ydx((0|:[)smx 1|:])"2 3 xdydy((1|:[)smx ])"3 2 rfhcn"_) rfmcndcdc11=:(0 1|:])"4@(ydx((0|:[)smx 2|:])"2 4 ydx((0|:[)smx 2|:])"2 4 xdydydy((1|:[)smx ])"4 2 xdy((1|:[)smx ])"2 2 rfhcn"_) rfmcndcdc12=:(0 1|:])"4@( ydxdx((0|:[)smx 1|:])"3 3 xdydy((1|:[)smx ])"3 2 xdy((1|:[)smx ])"2 2 rfhcn"_) rfmcndcdc0=:rfmcndcdc00+rfmcndcdc01+rfmcndcdc02 rfmcndcdc1=:rfmcndcdc10+rfmcndcdc11+rfmcndcdc12 rfmcndcdc=:rfmcndcdc0+rfmcndcdc1
NB. ... execute (ijx) ... gXRAll=:,:@(0 0 0 0 0 0 0 0 0 0"_) xCpts=:_ 1 3 26,_ 0.1p1 0.9p1 26,_ 0 2p1 26,:_ 0 20 26"_ (gXR([(rfmcndcdc(((2^_44);2^_44)qteq[;])"4 mcndcdc)xRef@(hkxR02,2 6$a:"_))]) ((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
NB. ... script experimentalR.ijs ... rfmcndcdcdc000=:(1 2 3|:])"5@(xdy((1|:[)smx ])"2 5 ydx((0|:[)smx 3|:])"2 5 ydx((0|:[)smx 3|:])"2 5 ydx((0|:[)smx 3|:])"2 5 xdydydydy((1|:[)smx ])"5 2 rfhcn"_) rfmcndcdcdc001=:(1 2 3|:])"5@( xdy((1|:[)smx ])"2 5 ydx((0|:[)smx 3|:])"2 5 ydxdx((0|:[)smx 2|:])"3 4 xdydydy((1|:[)smx ])"4 2 rfhcn"_) rfmcndcdcdc002=:(1 3 2|:])"5@( xdy((1|:[)smx ])"2 5 ydxdx((0|:[)smx 2|:])"3 4 ydx((0|:[)smx 2|:])"2 4 xdydydy((1|:[)smx ])"4 2 rfhcn"_) rfmcndcdcdc003=:(2 3 0|:])"5@(ydx((0|:[)smx 1|:])"2 5 xdydy((1|:[)smx ])"3 4 ydx((0|:[)smx 2|:])"2 4 ydx((0|:[)smx 2|:])"2 4 xdydydy((1|:[)smx ])"4 2 rfhcn"_) rfmcndcdcdc010=:(1 2 3|:])"5@( xdy((1|:[)smx ])"2 5 ydxdx((0|:[)smx 2|:])"3 4 ydx((0|:[)smx 2|:])"2 4 xdydydy((1|:[)smx ])"4 2 rfhcn"_) rfmcndcdcdc011=:(1 2 3|:])"5@( xdy((1|:[)smx ])"2 5 ydxdxdx((0|:[)smx 1|:])"4 3 xdydy((1|:[)smx ])"3 2 rfhcn"_) rfmcndcdcdc012=:(2 3 0|:])"5@( ydx((0|:[)smx 1|:])"2 5 xdydy((1|:[)smx ])"3 4 ydxdx((0|:[)smx 1|:])"3 3 xdydy((1|:[)smx ])"3 2 rfhcn"_) rfmcndcdcdc020=:(2 0 3|:])"5@(ydx((0|:[)smx 1|:])"2 5 xdydy((1|:[)smx ])"3 4 ydx((0|:[)smx 2|:])"2 4 ydx((0|:[)smx 2|:])"2 4 xdydydy((1|:[)smx ])"4 2 rfhcn"_) rfmcndcdcdc021=:(2 0 3|:])"5@( ydx((0|:[)smx 1|:])"2 5 xdydy((1|:[)smx ])"3 4 ydxdx((0|:[)smx 1|:])"3 3 xdydy((1|:[)smx ])"3 2 rfhcn"_) rfmcndcdcdc022=:(3 0 1|:])"5@(ydx((0|:[)smx 2|:])"2 5 ydx((0|:[)smx 2|:])"2 5 xdydydy((1|:[)smx ])"4 3 ydx((0|:[)smx 1|:])"2 3 xdydy((1|:[)smx ])"3 2 rfhcn"_) rfmcndcdcdc023=:(3 0 1|:])"5@( ydxdx((0|:[)smx 1|:])"3 4 xdydy((1|:[)smx ])"3 3 ydx((0|:[)smx 1|:])"2 3 xdydy((1|:[)smx ])"3 2 rfhcn"_) rfmcndcdcdc100=:(0 2 3|:])"5@(ydx((0|:[)smx 1|:])"2 5 xdydy((1|:[)smx ])"3 4 ydx((0|:[)smx 2|:])"2 4 ydx((0|:[)smx 2|:])"2 4 xdydydy((1|:[)smx ])"4 2 rfhcn"_) rfmcndcdcdc101=:(0 2 3|:])"5@( ydx((0|:[)smx 1|:])"2 5 xdydy((1|:[)smx ])"3 4 ydxdx((0|:[)smx 1|:])"3 3 xdydy((1|:[)smx ])"3 2 rfhcn"_) rfmcndcdcdc102=:(0 3 1|:])"5@(ydx((0|:[)smx 2|:])"2 5 ydx((0|:[)smx 2|:])"2 5 xdydydy((1|:[)smx ])"4 3 ydx((0|:[)smx 1|:])"2 3 xdydy((1|:[)smx ])"3 2 rfhcn"_) rfmcndcdcdc103=:(0 3 1|:])"5@( ydxdx((0|:[)smx 1|:])"3 4 xdydy((1|:[)smx ])"3 3 ydx((0|:[)smx 1|:])"2 3 xdydy((1|:[)smx ])"3 2 rfhcn"_) rfmcndcdcdc110=:(0 1 3|:])"5@(ydx((0|:[)smx 2|:])"2 5 ydx((0|:[)smx 2|:])"2 5 xdydydy((1|:[)smx ])"4 3 ydx((0|:[)smx 1|:])"2 3 xdydy((1|:[)smx ])"3 2 rfhcn"_) rfmcndcdcdc111=:(0 1 2|:])"5@(ydx((0|:[)smx 3|:])"2 5 ydx((0|:[)smx 3|:])"2 5 ydx((0|:[)smx 3|:])"2 5 xdydydydy((1|:[)smx ])"5 2 xdy((1|:[)smx ])"2 2 rfhcn"_) rfmcndcdcdc112=:(0 1 2|:])"5@( ydx((0|:[)smx 3|:])"2 5 ydxdx((0|:[)smx 2|:])"3 4 xdydydy((1|:[)smx ])"4 2 xdy((1|:[)smx ])"2 2 rfhcn"_) rfmcndcdcdc113=:(0 2 1|:])"5@( ydxdx((0|:[)smx 2|:])"3 4 ydx((0|:[)smx 2|:])"2 4 xdydydy((1|:[)smx ])"4 2 xdy((1|:[)smx ])"2 2 rfhcn"_) rfmcndcdcdc120=:(0 1 3|:])"5@( ydxdx((0|:[)smx 1|:])"3 4 xdydy((1|:[)smx ])"3 3 ydx((0|:[)smx 1|:])"2 3 xdydy((1|:[)smx ])"3 2 rfhcn"_) rfmcndcdcdc121=:(0 1 2|:])"5@( ydxdx((0|:[)smx 2|:])"3 4 ydx((0|:[)smx 2|:])"2 4 xdydydy((1|:[)smx ])"4 2 xdy((1|:[)smx ])"2 2 rfhcn"_) rfmcndcdcdc122=:(0 1 2|:])"5@( ydxdxdx((0|:[)smx 1|:])"4 3 xdydy((1|:[)smx ])"3 2 xdy((1|:[)smx ])"2 2 rfhcn"_) rfmcndcdcdc00=:rfmcndcdcdc000+rfmcndcdcdc001+rfmcndcdcdc002+rfmcndcdcdc003 rfmcndcdcdc01=:rfmcndcdcdc010+rfmcndcdcdc011+rfmcndcdcdc012 rfmcndcdcdc02=:rfmcndcdcdc020+rfmcndcdcdc021+rfmcndcdcdc022+rfmcndcdcdc023 rfmcndcdcdc10=:rfmcndcdcdc100+rfmcndcdcdc101+rfmcndcdcdc102+rfmcndcdcdc103 rfmcndcdcdc11=:rfmcndcdcdc110+rfmcndcdcdc111+rfmcndcdcdc112+rfmcndcdcdc113 rfmcndcdcdc12=:rfmcndcdcdc120+rfmcndcdcdc121+rfmcndcdcdc122 rfmcndcdcdc0=:rfmcndcdcdc00+rfmcndcdcdc01+rfmcndcdcdc02 rfmcndcdcdc1=:rfmcndcdcdc10+rfmcndcdcdc11+rfmcndcdcdc12 rfmcndcdcdc=:rfmcndcdcdc0+rfmcndcdcdc1
NB. ... execute (ijx) ...
gXRAll=:,:@(0 0 0 0 0 0 0 0 0 0"_)
xCpts=:_ 1 3 26,_ 0.1p1 0.9p1 26,_ 0 2p1 26,:_ 0 20 26"_
A=:(gXR([rfmcndcdcdc00 xRef@(hkxR02,2 6$a:"_))]) ((0;50"_)pTRandom xCpts)''
B=:(gXR((0|:[:(gXR([(0{])@rfmcndcdc00 xRef@(hkxR02,2 6$a:"_)),:@])D.1])"_ 1)]) ((0;50"_)pTRandom xCpts)''
A (((2^_21);2^_16)qteq[;])"5 B
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
A=:(gXR([rfmcndcdcdc01 xRef@(hkxR02,2 6$a:"_))]) ((0;50"_)pTRandom xCpts)''
B=:(gXR((0|:[:(gXR([(0{])@rfmcndcdc01 xRef@(hkxR02,2 6$a:"_)),:@])D.1])"_ 1)]) ((0;50"_)pTRandom xCpts)''
A (((2^_23);2^_16)qteq[;])"5 B
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
A=:(gXR([rfmcndcdcdc02 xRef@(hkxR02,2 6$a:"_))]) ((0;50"_)pTRandom xCpts)''
B=:(gXR((0|:[:(gXR([(0{])@rfmcndcdc02 xRef@(hkxR02,2 6$a:"_)),:@])D.1])"_ 1)]) ((0;50"_)pTRandom xCpts)''
A (((2^_22);2^_16)qteq[;])"5 B
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
A=:(gXR([rfmcndcdcdc10 xRef@(hkxR02,2 6$a:"_))]) ((0;50"_)pTRandom xCpts)''
B=:(gXR((0|:[:(gXR([(0{])@rfmcndcdc10 xRef@(hkxR02,2 6$a:"_)),:@])D.1])"_ 1)]) ((0;50"_)pTRandom xCpts)''
A (((2^_22);2^_16)qteq[;])"5 B
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
A=:(gXR([rfmcndcdcdc11 xRef@(hkxR02,2 6$a:"_))]) ((0;50"_)pTRandom xCpts)''
B=:(gXR((0|:[:(gXR([(0{])@rfmcndcdc11 xRef@(hkxR02,2 6$a:"_)),:@])D.1])"_ 1)]) ((0;50"_)pTRandom xCpts)''
A (((2^_21);2^_16)qteq[;])"5 B
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
A=:(gXR([rfmcndcdcdc12 xRef@(hkxR02,2 6$a:"_))]) ((0;50"_)pTRandom xCpts)''
B=:(gXR((0|:[:(gXR([(0{])@rfmcndcdc12 xRef@(hkxR02,2 6$a:"_)),:@])D.1])"_ 1)]) ((0;50"_)pTRandom xCpts)''
A (((2^_23);2^_16)qteq[;])"5 B
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
NB. ... execute (ijx) ... gXRAll=:,:@(0 0 0 0 0 0 0 0 0 0"_) xCpts=:_ 1 3 26,_ 0.1p1 0.9p1 26,_ 0 2p1 26,:_ 0 20 26"_ (gXR([(rfmcndcdcdc(((2^_41);2^_44)qteq[;])"5 mcndcdcdc)xRef@(hkxR02,2 6$a:"_))]) ((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
Download: File:LTER0900D.txt
Download: File:LTER0901D.txt
Download: File:LTER0902D.txt
Download: File:LTER0903H.txt