ruby-dcl 1.6.2

data/demo/math2/rnmlib/rnml01.rb

@@ -0,0 +1,21 @@
+#
+# rnml01.rb
+#
+# $Id: rnml01.rb,v 1.1.1.1 2011-02-23 07:21:28 horinout Exp $
+#
+
+require "narray"
+require "numru/dcl"
+
+include NumRu
+include Math
+
+
+N = 12
+
+rx = [0, 10, 20, 30, 40, 30, 20, 10, 0, -10, -20, -10].collect{|i| i.to_f}
+
+printf " LIST OF RX : %s\n", "%5.0f"*N % rx.to_a
+ry = DCL::vrrnm(rx, N, 1, 1, 3)
+print " AFTER CALLING VRRNM(RX,RY,12,1,1,3)\n"
+printf " LIST OF RY : %s\n", "%5.0f"*N % ry.to_a

data/demo/math2/vstlib/vstl01.rb

@@ -0,0 +1,29 @@
+#
+# vstl01.rb
+#
+# $Id: vstl01.rb,v 1.1.1.1 2011-02-23 07:21:28 horinout Exp $
+#
+
+require "narray"
+require "numru/dcl"
+
+include NumRu
+include Math
+
+
+IX = 3
+NN = 4
+
+x = NArray[1, 3,-1, 2,-4, 0, 2, 1, 3, 2, 1, 0]
+x.reshape!(IX, NN)
+
+print " *** TEST FOR VS1INT / VS1DIN / VS1OUT\n"
+wz1, nw1 = DCL::vs1int(IX)
+for j in 0..NN-1
+  printf " INPUT = %s\n", sprintf("%6.1f"*IX, x[0,j], x[1,j], x[2,j])
+  wz1, nw1 = DCL::vs1din(wz1, nw1, x[true,j])
+end
+wz1 = DCL::vs1out(wz1, nw1)
+printf " AVE.  = %s\n", "%6.1f"*IX % wz1[0..IX-1].to_a
+printf " VAR.  = %s\n", "%6.1f"*IX % wz1[IX..2*IX-1].to_a
+

data/demo/math2/vstlib/vstl02.rb

@@ -0,0 +1,38 @@
+#
+# vstl02.rb
+#
+# $Id: vstl02.rb,v 1.1.1.1 2011-02-23 07:21:28 horinout Exp $
+#
+
+require "narray"
+require "numru/dcl"
+
+include NumRu
+include Math
+
+
+IX = 3
+IY = 3
+NN = 4
+
+x = NArray[1, 3,-1, 2,-4, 0, 2, 1, 3, 2, 1, 0]
+y = NArray[0, 1, 0, 1,-2, 1, 3, 0, 2, 1,-1, 0]
+x.reshape!(IX, NN)
+y.reshape!(IY, NN)
+
+fmt1 = " INPUT(%s) = " + "%6.1f"*3 + "\n"
+fmt2 = "%6.1f"*3 + ("\n" + " "*12 + "%6.1f"*3)*2
+
+print " *** TEST FOR VS2INT / VS2DIN / VS2OUT\n"
+wz2, nw2 = DCL::vs2int(IX, IY)
+for j in 0..NN-1
+  printf fmt1, "X", x[0,j], x[1,j], x[2,j]
+  printf fmt1, "Y", y[0,j], y[1,j], y[2,j]
+  wz2, nw2 = DCL::vs2din(wz2, nw2, x[true,j], y[true,j])
+end
+wz2 = DCL::vs2out(wz2, nw2)
+printf " AVE. (X) = %s\n", fmt2 % wz2.to_a.flatten[0,IX*IY]
+printf " AVE. (Y) = %s\n", fmt2 % wz2.to_a.flatten[IX*IY,IX*IY]
+printf " VAR. (X) = %s\n", fmt2 % wz2.to_a.flatten[IX*IY*2,IX*IY]
+printf " VAR. (Y) = %s\n", fmt2 % wz2.to_a.flatten[IX*IY*3,IX*IY]
+printf " COV(X,Y) = %s\n", fmt2 % wz2.to_a.flatten[IX*IY*4,IX*IY]

data/demo/rakuraku/color/color1.rb

@@ -0,0 +1,61 @@
+#
+# color1.rb
+#   $Id: color1.rb,v 1.1.1.1 2011-02-23 07:21:26 horinout Exp $
+#
+
+
+require "narray"
+require "numru/dcl"
+include NumRu
+include Math
+
+
+np = 4
+nt = 10
+dd = 0.07
+ds = 0.003
+dx = dd+ds*2
+dy = 1.0/nt
+
+#-- data ---
+xbox = NArray[0.0, 1.0, 1.0, 0.0]
+ybox = NArray[0.0, 0.0, 1.0, 1.0]
+
+#-- graph ---
+iws = (ARGV[0] || (puts ' WORKSTATION ID (I)  ? ;'; DCL::sgpwsn; gets)).to_i
+DCL::gliset('MAXMSG', 300)
+
+DCL::sgopn(-(iws).abs)
+
+DCL::sglset('LFULL', true)
+DCL::slrat(0.85, 1.0)
+DCL::slmgn(0.0, 0.0, 0.05, 0.1)
+DCL::slsttl('TEST OF COLORMAP', 'T', 0.0, -0.5, 0.025, 1)
+
+DCL::sgfrm
+DCL::sgswnd(0.0, 1.0, 0.0, 1.0)
+DCL::sgstrn(1)
+
+DCL::sgstxs(0.01)
+DCL::sgstxc(1)
+for j in 0..9
+  for i in 0..9
+    x1 = dx*i     + ds + (1-dx*nt)/2
+    y1 = dy*(9-j) + ds
+    x2 = x1 + dd
+    y2 = y1 + dd
+    itpat = (i+j*10)*1000 + 999
+
+    DCL::sgsvpt(x1, x2, y1, y2)
+    DCL::sgstrf
+
+    DCL::slpvpr(1)
+    DCL::sgstnp(itpat)
+    DCL::sgtnu(xbox, ybox)
+    chr = format("%5d", itpat)
+    DCL::sgtxv(x2, y2+(dy-dd)/2, chr)
+  end
+end
+
+DCL::sgcls
+

data/demo/rakuraku/color/color2.rb

@@ -0,0 +1,57 @@
+#
+# color2.rb
+#   $Id: color2.rb,v 1.1.1.1 2011-02-23 07:21:26 horinout Exp $
+#
+
+
+require "narray"
+require "numru/dcl"
+include NumRu
+include Math
+
+
+nmax = 50
+x = NArray.sfloat(nmax+1)
+y = NArray.sfloat(nmax+1)
+
+#-- data ---
+x.indgen(1950.0, 1.0)
+r    = 3.7
+y[0] = 0.5
+for n in 0..nmax-1
+  y[n+1] = r*y[n]*(1.0-y[n])
+end
+yave = DCL::rave(y, nmax+1, 1)
+y = -4 * (y - yave)
+
+#-- graph ---
+iws = (ARGV[0] || (puts ' WORKSTATION ID (I)  ? ;'; DCL::sgpwsn; gets)).to_i
+DCL::gropn iws
+
+DCL::grfrm
+
+DCL::grswnd(x[0], x[nmax], -1.5, 2.0)
+DCL::grsvpt(0.2, 0.9, 0.2, 0.9)
+DCL::grstrn(1)
+DCL::grstrf
+
+DCL::sgiset('IFONT', 2)
+DCL::uziset('INDEXT1', 402)
+DCL::uziset('INDEXT2', 404)
+DCL::uziset('INDEXL1', 303)
+DCL::uziset('INDEXL2', 123)
+DCL::uxaxdv('B', 2.0, 10.0)
+DCL::uxaxdv('T', 2.0, 10.0)
+DCL::uxsttl('B', 'YEAR', 0.0)
+
+DCL::uyaxdv('L', 0.1, 0.5)
+DCL::uyaxdv('R', 0.1, 0.5)
+cttl = DCL::csgi(131) + 'T [K]'
+DCL::uysttl('L', cttl, 0.0)
+
+DCL::uxmttl('T', 'INTERANNUAL VARIATION', 0.0)
+
+DCL::sgspli(24)
+DCL::sgplu(x, y)
+DCL::grcls
+

data/demo/rakuraku/color/color3.rb

@@ -0,0 +1,81 @@
+#
+# color3.rb
+#   $Id: color3.rb,v 1.1.1.1 2011-02-23 07:21:26 horinout Exp $
+#
+
+
+require "narray"
+require "numru/dcl"
+include NumRu
+include Math
+
+
+nx = 21
+ny = 21
+xmin = -10
+xmax = 10
+ymin = -10
+ymax = 10
+dx1 = 1
+dx2 = 5
+dy1 = 1
+dy2 = 5
+kmax = 56
+pmin = 0
+pmax = 1
+p = NArray.sfloat(nx, ny)
+pi = NArray.sfloat(2, kmax+1)
+
+#-- data ---
+for j in 0..ny-1
+  for i in 0..nx-1
+    x = xmin + (xmax-xmin)*i/(nx-1).to_f
+    y = ymin + (ymax-ymin)*j/(ny-1).to_f
+    p[i,j] = exp(-x**2/64 -y**2/25)
+  end
+end
+
+#-- graph ---
+iws = (ARGV[0] || (puts ' WORKSTATION ID (I)  ? ;'; DCL::sgpwsn; gets)).to_i
+DCL::gropn iws
+
+DCL::grfrm
+
+DCL::grswnd(xmin, xmax, ymin, ymax)
+DCL::grsvpt(0.2, 0.8, 0.2, 0.8)
+DCL::grstrn(1)
+DCL::grstrf
+
+dp = (pmax-pmin)/kmax.to_f
+for k in 1..kmax
+  tlev1 = (k-1)*dp
+  tlev2 = tlev1 + dp
+  ipat  = (29+k)*1000 + 999
+  DCL::uestlv(tlev1, tlev2, ipat)
+end
+
+DCL::uetone(p)
+DCL::usdaxs
+
+#-- �ȡ���С� ----
+DCL::grswnd(0.0, 1.0, pmin, pmax)
+DCL::grsvpt(0.85, 0.9, 0.2, 0.8)
+DCL::grstrn(1)
+DCL::grstrf
+
+for k in 0..kmax
+  pi[0,k] = pmin + k*dp
+  pi[1,k] = pmin + k*dp
+end
+
+DCL::uetone(pi)
+
+DCL::slpvpr(3)
+DCL::uzlset('LABELYR', true)
+
+DCL::uzfact(0.8)
+DCL::uysfmt('(F4.1)')
+DCL::uyaxdv('R', 0.1, 0.2)
+
+DCL::grcls
+

data/demo/rakuraku/kihon/kihon1.rb

@@ -0,0 +1,39 @@
+#
+# kihon1.rb
+#   $Id: kihon1.rb,v 1.1.1.1 2011-02-23 07:21:26 horinout Exp $
+#
+
+require "narray"
+require "numru/dcl"
+
+include NumRu
+include NMath
+
+nmax = 50
+
+x  = NArray.sfloat(nmax+1)
+y1 = NArray.sfloat(nmax+1)
+y2 = NArray.sfloat(nmax+1)
+y3 = NArray.sfloat(nmax+1)
+
+#-- data ----
+x = x.indgen * 1.0/nmax
+y1 = x**3
+y2 = x**2
+y3 = sqrt(x)
+
+#-- graph ----
+iws = (ARGV[0] || (puts ' WORKSTATION ID (I)  ? ;'; DCL::sgpwsn; gets)).to_i
+DCL::sgopn iws
+DCL::sgfrm
+
+DCL::slpvpr(1)
+
+DCL::sgplv(x, y1)
+DCL::sgpmv(x, y2)
+DCL::sgtxv(0.5, 0.5, 'SGTXV')
+DCL::sglset('LSOFTF', true)
+DCL::sgtnv(x, y3)
+
+DCL::sgcls
+

data/demo/rakuraku/kihon/kihon2.rb

@@ -0,0 +1,56 @@
+#
+# kihon2.rb
+#   $Id: kihon2.rb,v 1.1.1.1 2011-02-23 07:21:26 horinout Exp $
+#
+
+require "narray"
+require "numru/dcl"
+
+include NumRu
+include Math
+
+nmax = 40
+imax = 4
+dt = 4.0*PI/nmax
+
+x = NArray.sfloat(nmax+1)
+y = NArray.sfloat(nmax+1, imax)
+
+#-- data ----
+x = x.indgen * 1.0/nmax
+for n in 0..nmax
+  for i in 0..imax-1
+    y[n,i] = 0.2*sin(n*dt) + 0.8 - 0.2*i
+  end
+end
+
+#-- graph ----
+iws = (ARGV[0] || (puts ' WORKSTATION ID (I)  ? ;'; DCL::sgpwsn; gets)).to_i
+DCL::sgopn iws
+
+#-- line index : frame 1 --
+DCL::sgfrm
+DCL::slpvpr(1) 
+
+DCL::sgplv(x, y[true,0])
+
+for i in 1..imax-1
+  DCL::sgspli(2*(i+1))
+  DCL::sgplv(x, y[true,i])
+end
+
+DCL::sgspli(1)
+
+#-- line type : frame 2 --
+DCL::sgfrm
+DCL::slpvpr(1)
+
+DCL::sgplv(x, y[true,0])
+
+for i in 1..imax-1
+  DCL::sgsplt(i+1)
+  DCL::sgplv(x, y[true,i])
+end
+
+DCL::sgcls
+

data/demo/rakuraku/kihon/kihon3.rb

@@ -0,0 +1,57 @@
+#
+# kihon3.rb
+#   $Id: kihon3.rb,v 1.1.1.1 2011-02-23 07:21:26 horinout Exp $
+#
+
+require "narray"
+require "numru/dcl"
+
+include NumRu
+include Math
+
+nmax = 40
+imax = 7
+dt = 4.0*PI/nmax
+
+x = NArray.sfloat(nmax+1)
+y = NArray.sfloat(nmax+1, imax)
+
+#-- data ----
+x = x.indgen * 1.0/nmax
+for n in 0..nmax
+  for i in 0..imax-1
+    y[n,i] = 0.2*sin(n*dt) + 0.8 - 0.1*i
+  end
+end
+
+#-- graph ----
+iws = (ARGV[0] || (puts ' WORKSTATION ID (I)  ? ;'; DCL::sgpwsn; gets)).to_i
+DCL::sgopn iws
+
+#-- marker type : frame 1 --
+DCL::sgfrm
+DCL::slpvpr(1) 
+
+DCL::sgpmv(x, y[true,0])
+
+for i in 1..imax-1
+  itype = i + 1
+  DCL::sgspmt(itype)
+  DCL::sgpmv(x, y[true,i])
+end
+
+#-- marker size : frame 2 --
+DCL::sgfrm
+DCL::slpvpr(1) 
+
+DCL::sgspmt(10)
+DCL::sgpmv(x, y[true,0])
+
+for i in 1..imax-1
+  rsize = 0.005*i
+  DCL::sgspms(rsize)
+  DCL::sgpmv(x, y[true,i])
+end
+
+DCL::sgcls
+