ruby-dcl 1.6.2

data/demo/grph2/udegpk/u2df09e.rb

@@ -0,0 +1,88 @@
+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 = 5
+pmin = 0
+pmax = 1
+u = NArray.sfloat(nx, ny)
+v = NArray.sfloat(nx, ny)
+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
+    u[i,j] =  x
+    v[i,j] = -y
+    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::sglset('LSOFTF', true)
+DCL::grfrm
+
+DCL::grswnd(xmin, xmax, ymin, ymax)
+DCL::grsvpt(0.2, 0.8, 0.2, 0.8)
+DCL::grstrn(1)
+DCL::grstrf
+
+DCL::uxaxdv('B', dx1, dx2)
+DCL::uxaxdv('T', dx1, dx2)
+DCL::uxsttl('B', 'X-axis', 0.0)
+DCL::uxmttl('T', 'FIGURE TITLE', 0.0)
+
+DCL::uyaxdv('L', dy1, dy2)
+DCL::uyaxdv('R', dy1, dy2)
+DCL::uysttl('L', 'Y-axis', 0.0)
+
+dp = (pmax-pmin)/kmax.to_f
+for k in 1..kmax
+  tlev1 = (k-1)*dp
+  tlev2 = tlev1 + dp
+  ipat  = 600 + k - 1
+  DCL::uestlv(tlev1, tlev2, ipat)
+end
+DCL::uetone(p)
+
+#-- �ȡ���С� ----
+DCL::grswnd(0.0, 1.0, pmin, pmax)
+DCL::grsvpt(0.85, 0.9, 0.45, 0.75)
+DCL::grstrn(1)
+DCL::grstrf
+
+for k in 1..kmax+1
+  pi[0,k-1] = pmin + (k-1)*dp
+  pi[1,k-1] = pmin + (k-1)*dp
+end
+
+
+DCL::udlset('LMSG', false)
+DCL::udcntr(pi)
+DCL::uetone(pi)
+
+DCL::slpvpr(3)
+DCL::uzlset('LABELYR', true)
+DCL::uzfact(0.8)
+DCL::uysfmt('(F4.1)')
+DCL::uyaxdv('R', dp, dp)
+
+DCL::grcls

data/demo/grph2/udegpk/u2df09f.rb

@@ -0,0 +1,89 @@
+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 = 5
+pmin = 0
+pmax = 1
+u = NArray.sfloat(nx, ny)
+v = NArray.sfloat(nx, ny)
+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
+    u[i,j] =  x
+    v[i,j] = -y
+    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::sgscmn(5)
+DCL::gropn iws
+
+DCL::sglset('LSOFTF', false)
+DCL::grfrm
+
+DCL::grswnd(xmin, xmax, ymin, ymax)
+DCL::grsvpt(0.2, 0.8, 0.2, 0.8)
+DCL::grstrn(1)
+DCL::grstrf
+
+DCL::uxaxdv('B', dx1, dx2)
+DCL::uxaxdv('T', dx1, dx2)
+DCL::uxsttl('B', 'X-axis', 0.0)
+DCL::uxmttl('T', 'FIGURE TITLE', 0.0)
+
+DCL::uyaxdv('L', dy1, dy2)
+DCL::uyaxdv('R', dy1, dy2)
+DCL::uysttl('L', 'Y-axis', 0.0)
+
+dp = (pmax-pmin)/kmax.to_f
+for k in 1..kmax
+  tlev1 = (k-1)*dp
+  tlev2 = tlev1 + dp
+  ipat  = 999 + (99 - k*15)*1000
+  DCL::uestlv(tlev1, tlev2, ipat)
+end
+DCL::uetonf(p)
+
+#-- �ȡ���С� ----
+DCL::grswnd(0.0, 1.0, pmin, pmax)
+DCL::grsvpt(0.85, 0.9, 0.45, 0.75)
+DCL::grstrn(1)
+DCL::grstrf
+
+for k in 1..kmax+1
+  pi[0,k-1] = pmin + (k-1)*dp
+  pi[1,k-1] = pmin + (k-1)*dp
+end
+
+
+DCL::udlset('LMSG', false)
+DCL::udcntr(pi)
+DCL::uetonf(pi)
+
+DCL::slpvpr(3)
+DCL::uzlset('LABELYR', true)
+DCL::uzfact(0.8)
+DCL::uysfmt('(F4.1)')
+DCL::uyaxdv('R', dp, dp)
+
+DCL::grcls

data/demo/grph2/ulpack/ulpk01.rb

@@ -0,0 +1,111 @@
+#
+# ulpk01.rb
+#
+# $Id: ulpk01.rb,v 1.1.1.1 2011-02-23 07:21:27 horinout Exp $
+#
+
+require "narray"
+require "numru/dcl"
+
+include NumRu
+include Math
+
+
+nbl = 2
+
+#-- data ---
+bl = NArray[1.0, 3.0]
+
+#-- graph ---
+iws = (ARGV[0] || (puts ' WORKSTATION ID (I)  ? ;'; DCL::sgpwsn; gets)).to_i
+DCL::sgopn(-(iws.abs))
+
+DCL::sglset('LFULL', true)
+DCL::uzfact(0.7)
+DCL::slrat(0.75, 1.0)
+DCL::sldiv('T', 1, 7)
+
+DCL::sgfrm
+DCL::uzinit
+DCL::sgswnd(0.15, 200.0, 0.1, 1e5)
+DCL::sgsvpt(0.1, 0.9, 0.1, 0.11)
+DCL::sgstrn(4)
+DCL::sgstrf
+DCL::uliset('IXTYPE', 1)
+DCL::ulxlog('B', 1, 9)
+ctl = 'ULXLOG (IXTYPE=1,NLBL=1,NTCK=9)'
+DCL::uxsttl('B', ctl, 0.0)
+
+DCL::sglset('LCNTL', true)
+
+DCL::sgfrm
+DCL::uzinit
+DCL::sgswnd(0.15, 200.0, 0.1, 1e5)
+DCL::sgsvpt(0.1, 0.9, 0.1, 0.11)
+DCL::sgstrn(4)
+DCL::sgstrf
+DCL::uliset('IXTYPE', 1)
+DCL::ulxlog('B', 3, 9)
+ctl = 'ULXLOG (IXTYPE=1,NLBL=3,NTCK=9)'
+DCL::uxsttl('B', ctl, 0.0)
+
+DCL::sgfrm
+DCL::uzinit
+DCL::sgswnd(0.15, 200.0, 0.1, 1e5)
+DCL::sgsvpt(0.1, 0.9, 0.1, 0.11)
+DCL::sgstrn(4)
+DCL::sgstrf
+DCL::uliset('IXTYPE', 2)
+DCL::ulxlog('B', 3, 9)
+ctl = 'ULXLOG (IXTYPE=2,NLBL=3,NTCK=9)'
+DCL::uxsttl('B', ctl, 0.0)
+
+DCL::sgfrm
+DCL::uzinit
+DCL::sgswnd(0.15, 200.0, 0.1, 1e5)
+DCL::sgsvpt(0.1, 0.9, 0.1, 0.11)
+DCL::sgstrn(4)
+DCL::sgstrf
+DCL::uliset('IXTYPE', 3)
+DCL::ulxlog('B', 3, 9)
+ctl = 'ULXLOG (IXTYPE=3,NLBL=3,NTCK=9)'
+DCL::uxsttl('B', ctl, 0.0)
+
+DCL::sgfrm
+DCL::uzinit
+DCL::sgswnd(0.15, 200.0, 0.1, 1e5)
+DCL::sgsvpt(0.1, 0.9, 0.1, 0.11)
+DCL::sgstrn(4)
+DCL::sgstrf
+DCL::uliset('IXTYPE', 4)
+DCL::ulxlog('B', 3, 9)
+ctl = 'ULXLOG (IXTYPE=4,NLBL=3,NTCK=9)'
+DCL::uxsttl('B', ctl, 0.0)
+
+DCL::sgfrm
+DCL::uzinit
+DCL::sgswnd(0.15, 200.0, 0.1, 1e5)
+DCL::sgsvpt(0.1, 0.9, 0.1, 0.11)
+DCL::sgstrn(4)
+DCL::sgstrf
+DCL::ulxsfm('(F5.1)')
+DCL::uliset('IXTYPE', 3)
+DCL::ulxlog('B', 3, 9)
+ctl = 'ULXLOG (IXTYPE=3,NLBL=3,NTCK=9)'
+DCL::uxsttl('B', ctl, 0.0)
+
+DCL::sgfrm
+DCL::uzinit
+DCL::sgswnd(0.15, 200.0, 0.1, 1e5)
+DCL::sgsvpt(0.1, 0.9, 0.1, 0.11)
+DCL::sgstrn(4)
+DCL::sgstrf
+DCL::uliset('IXTYPE', 1)
+DCL::uliset('IXCHR', 194)
+DCL::ulsxbl(bl)
+DCL::ulxlog('B', 4, 5)
+ctl = 'ULXLOG (IXTYPE=1,NLBL=4,NTCK=5)'
+DCL::uxsttl('B', ctl, 0.0)
+
+DCL::sgcls
+

data/demo/grph2/ulpack/ulpk01n.rb

@@ -0,0 +1,111 @@
+#
+# ulpk01.rb
+#
+# $Id: ulpk01n.rb,v 1.1.1.1 2011-02-23 07:21:27 horinout Exp $
+#
+
+require "narray"
+require "numru/dcl"
+
+include NumRu
+include Math
+
+
+nbl = 2
+
+#-- data ---
+bl = NArray[1.0, 3.0]
+
+#-- graph ---
+iws = (ARGV[0] || (puts ' WORKSTATION ID (I)  ? ;'; DCL::sgpwsn; gets)).to_i
+DCL::sgopn(-(iws.abs))
+
+DCL::sglset('LFULL', true)
+DCL::uzfact(0.6)
+DCL::slrat(0.75, 1.0)
+DCL::sldiv('T', 1, 7)
+
+DCL::sgfrm
+DCL::uzinit
+DCL::sgswnd(-0.15, -200.0, 0.1, 1e5)
+DCL::sgsvpt(0.1, 0.9, 0.1, 0.11)
+DCL::sgstrn(4)
+DCL::sgstrf
+DCL::uliset('IXTYPE', 1)
+DCL::ulxlog('B', 1, 9)
+ctl = 'ULXLOG (IXTYPE=1,NLBL=1,NTCK=9)'
+DCL::uxsttl('B', ctl, 0.0)
+
+DCL::sglset('LCNTL', true)
+
+DCL::sgfrm
+DCL::uzinit
+DCL::sgswnd(-0.15, -200.0, 0.1, 1e5)
+DCL::sgsvpt(0.1, 0.9, 0.1, 0.11)
+DCL::sgstrn(4)
+DCL::sgstrf
+DCL::uliset('IXTYPE', 1)
+DCL::ulxlog('B', 3, 9)
+ctl = 'ULXLOG (IXTYPE=1,NLBL=3,NTCK=9)'
+DCL::uxsttl('B', ctl, 0.0)
+
+DCL::sgfrm
+DCL::uzinit
+DCL::sgswnd(-0.15, -200.0, 0.1, 1e5)
+DCL::sgsvpt(0.1, 0.9, 0.1, 0.11)
+DCL::sgstrn(4)
+DCL::sgstrf
+DCL::uliset('IXTYPE', 2)
+DCL::ulxlog('B', 3, 9)
+ctl = 'ULXLOG (IXTYPE=2,NLBL=3,NTCK=9)'
+DCL::uxsttl('B', ctl, 0.0)
+
+DCL::sgfrm
+DCL::uzinit
+DCL::sgswnd(-0.15, -200.0, 0.1, 1e5)
+DCL::sgsvpt(0.1, 0.9, 0.1, 0.11)
+DCL::sgstrn(4)
+DCL::sgstrf
+DCL::uliset('IXTYPE', 3)
+DCL::ulxlog('B', 3, 9)
+ctl = 'ULXLOG (IXTYPE=3,NLBL=3,NTCK=9)'
+DCL::uxsttl('B', ctl, 0.0)
+
+DCL::sgfrm
+DCL::uzinit
+DCL::sgswnd(-0.15, -200.0, 0.1, 1e5)
+DCL::sgsvpt(0.1, 0.9, 0.1, 0.11)
+DCL::sgstrn(4)
+DCL::sgstrf
+DCL::uliset('IXTYPE', 4)
+DCL::ulxlog('B', 3, 9)
+ctl = 'ULXLOG (IXTYPE=4,NLBL=3,NTCK=9)'
+DCL::uxsttl('B', ctl, 0.0)
+
+DCL::sgfrm
+DCL::uzinit
+DCL::sgswnd(-0.15, -200.0, 0.1, 1e5)
+DCL::sgsvpt(0.1, 0.9, 0.1, 0.11)
+DCL::sgstrn(4)
+DCL::sgstrf
+DCL::ulxsfm('(F5.1)')
+DCL::uliset('IXTYPE', 3)
+DCL::ulxlog('B', 3, 9)
+ctl = 'ULXLOG (IXTYPE=3,NLBL=3,NTCK=9)'
+DCL::uxsttl('B', ctl, 0.0)
+
+DCL::sgfrm
+DCL::uzinit
+DCL::sgswnd(-0.15, -200.0, 0.1, 1e5)
+DCL::sgsvpt(0.1, 0.9, 0.1, 0.11)
+DCL::sgstrn(4)
+DCL::sgstrf
+DCL::uliset('IXTYPE', 1)
+DCL::uliset('IXCHR', 194)
+DCL::ulsxbl(bl)
+DCL::ulxlog('B', 4, 5)
+ctl = 'ULXLOG (IXTYPE=1,NLBL=4,NTCK=5)'
+DCL::uxsttl('B', ctl, 0.0)
+
+DCL::sgcls
+

data/demo/grph2/ulpack/ulpk02.rb

@@ -0,0 +1,111 @@
+#
+# ulpk02.rb
+#
+# $Id: ulpk02.rb,v 1.1.1.1 2011-02-23 07:21:27 horinout Exp $
+#
+
+require "narray"
+require "numru/dcl"
+
+include NumRu
+include Math
+
+
+nbl = 2
+
+#-- data ---
+bl = [1.0, 3.0]
+
+#-- graph ---
+iws = (ARGV[0] || (puts ' WORKSTATION ID (I)  ? ;'; DCL::sgpwsn; gets)).to_i
+DCL::sgopn(+iws.abs)
+
+DCL::sglset('LFULL', true)
+DCL::uzfact(0.7)
+DCL::slrat(1.0, 0.75)
+DCL::sldiv('Y', 7, 1)
+
+DCL::sgfrm
+DCL::uzinit
+DCL::sgswnd(0.2, 200.0, 0.15, 200.0)
+DCL::sgsvpt(0.13, 0.14, 0.1, 0.9)
+DCL::sgstrn(4)
+DCL::sgstrf
+DCL::uliset('IYTYPE', 1)
+DCL::ulylog('L', 1, 9)
+ctl = 'ULYLOG (IYTYPE=1,NLBL=1,NTCK=9)'
+DCL::uysttl('L', ctl, 0.0)
+
+DCL::sglset('LCNTL', true)
+
+DCL::sgfrm
+DCL::uzinit
+DCL::sgswnd(0.2, 200.0, 0.15, 200.0)
+DCL::sgsvpt(0.13, 0.14, 0.1, 0.9)
+DCL::sgstrn(4)
+DCL::sgstrf
+DCL::uliset('IYTYPE', 1)
+DCL::ulylog('L', 3, 9)
+ctl = 'ULYLOG (IYTYPE=1,NLBL=3,NTCK=9)'
+DCL::uysttl('L', ctl, 0.0)
+
+DCL::sgfrm
+DCL::uzinit
+DCL::sgswnd(0.2, 200.0, 0.15, 200.0)
+DCL::sgsvpt(0.13, 0.14, 0.1, 0.9)
+DCL::sgstrn(4)
+DCL::sgstrf
+DCL::uliset('IYTYPE', 2)
+DCL::ulylog('L', 3, 9)
+ctl = 'ULYLOG (IYTYPE=2,NLBL=3,NTCK=9)'
+DCL::uysttl('L', ctl, 0.0)
+
+DCL::sgfrm
+DCL::uzinit
+DCL::sgswnd(0.2, 200.0, 0.15, 200.0)
+DCL::sgsvpt(0.13, 0.14, 0.1, 0.9)
+DCL::sgstrn(4)
+DCL::sgstrf
+DCL::uliset('IYTYPE', 3)
+DCL::ulylog('L', 3, 9)
+ctl = 'ULYLOG (IYTYPE=3,NLBL=3,NTCK=9)'
+DCL::uysttl('L', ctl, 0.0)
+
+DCL::sgfrm
+DCL::uzinit
+DCL::sgswnd(0.2, 200.0, 0.15, 200.0)
+DCL::sgsvpt(0.13, 0.14, 0.1, 0.9)
+DCL::sgstrn(4)
+DCL::sgstrf
+DCL::uliset('IYTYPE', 4)
+DCL::ulylog('L', 3, 9)
+ctl = 'ULYLOG (IYTYPE=4,NLBL=3,NTCK=9)'
+DCL::uysttl('L', ctl, 0.0)
+
+DCL::sgfrm
+DCL::uzinit
+DCL::sgswnd(0.2, 200.0, 0.15, 200.0)
+DCL::sgsvpt(0.13, 0.14, 0.1, 0.9)
+DCL::sgstrn(4)
+DCL::sgstrf
+DCL::ulysfm('(F5.1)')
+DCL::uliset('IYTYPE', 3)
+DCL::ulylog('L', 3, 9)
+ctl = 'ULYLOG (IYTYPE=3,NLBL=3,NTCK=9)'
+DCL::uysttl('L', ctl, 0.0)
+
+DCL::sgfrm
+DCL::uzinit
+DCL::sgswnd(0.2, 200.0, 0.15, 200.0)
+DCL::sgsvpt(0.13, 0.14, 0.1, 0.9)
+DCL::sgstrn(4)
+DCL::sgstrf
+DCL::uliset('IYTYPE', 1)
+DCL::uliset('IYCHR', 194)
+DCL::ulsybl(bl)
+DCL::ulylog('L', 4, 5)
+ctl = 'ULYLOG (IYTYPE=1,NLBL=4,NTCK=5)'
+DCL::uysttl('L', ctl, 0.0)
+
+DCL::sgcls
+