home *** CD-ROM | disk | FTP | other *** search
/ ProfitPress Mega CDROM2 …eeware (MSDOS)(1992)(Eng) / ProfitPress-MegaCDROM2.B6I / MISC / MATH / ULTRA101.ZIP / DEMO.FOR < prev    next >
Encoding:
Text File  |  1992-03-25  |  2.0 KB  |  95 lines
  1. c  ULTRA DEMO
  2. c
  3. c  This calls the various subroutines of ultra as a test
  4.  
  5.     external uni,vni,duni,dvni
  6.     external i32bit,i31bit,i16bit,i15bit,i8bit,i7bit,i1bit
  7.  
  8.     real*4     uni,     vni
  9.     real*8     duni,    dvni,    d(4)
  10.     integer*4  i32bit,  i31bit,  k(32)
  11.     integer*2  i16bit,  i15bit
  12.     integer*2  i8bit,   i7bit,   i1bit
  13.  
  14.     print 1
  15.     1    format (' A Sample Printout:')
  16.  
  17.     call rinit(1234567,7654321)
  18.     do 2 i=1,10000
  19.     j=min(2,i)
  20.     2    d(j)=dvni()
  21.         print 12,(d(j),j=1,2)
  22.    12    format (' Signed Double: ',2f32.15)
  23.  
  24.     call rinit(1234567,7654321)
  25.     do 3 i=1,10000
  26.     j=min(2,i)
  27.     3    d(j)=duni()
  28.     print 13,(d(j),j=1,2)
  29.    13    format (' Unsigned Dble: ',2f32.15)
  30.  
  31.     call rinit(1234567,7654321)
  32.     do 4 i=1,10000
  33.     j=min(4,i)
  34.     4    d(j)=vni()
  35.     print 14,(d(j),j=1,4)
  36.    14    format (' Signed Single: ',4f16.7)
  37.  
  38.     call rinit(1234567,7654321)
  39.     do 5 i=1,10000
  40.     j=min(4,i)
  41.     5    d(j)=uni()
  42.     print 15,(d(j),j=1,4)
  43.    15    format (' Uns. Single:   ',4f16.7)
  44.  
  45.     call rinit(1234567,7654321)
  46.     do 6 i=1,10000
  47.     j=min(4,i)
  48.     6    k(j)=i32bit()
  49.     print 16,(k(j),j=1,4)
  50.    16    format (' Signed 4 byte: ',4i16)
  51.  
  52.     call rinit(1234567,7654321)
  53.     do 7 i=1,10000
  54.     j=min(4,i)
  55.     7    k(j)=i31bit()
  56.     print 17,(k(j),j=1,4)
  57.    17    format (' Uns. 4 byte:   ',4i16)
  58.  
  59.     call rinit(1234567,7654321)
  60.     do 8 i=1,10000
  61.     j=min(8,i)
  62.     8    k(j)=i16bit()
  63.     print 18,(k(j),j=1,8)
  64.    18    format (' Signed 2 byte: ',8i8)
  65.  
  66.     call rinit(1234567,7654321)
  67.     do 9 i=1,10000
  68.     j=min(8,i)
  69.     9    k(j)=i15bit()
  70.     print 19,(k(j),j=1,8)
  71.    19    format (' Uns. 2 byte:   ',8i8)
  72.  
  73.     call rinit(1234567,7654321)
  74.     do 10 i=1,10000
  75.     j=min(16,i)
  76.    10    k(j)=i8bit()
  77.     print 20,(k(j),j=1,16)
  78.    20    format (' Signed byte:   ',16i4)
  79.  
  80.     call rinit(1234567,7654321)
  81.     do 11 i=1,10000
  82.     j=min(16,i)
  83.    11    k(j)=i7bit()
  84.     print 21,(k(j),j=1,16)
  85.    21    format (' Uns. byte:     ',16i4)
  86.  
  87.     call rinit(1234567,7654321)
  88.     do 22 i=1,10000
  89.     j=min(32,i)
  90.    22    k(j)=i1bit()
  91.     print 23,(k(j),j=1,32)
  92.    23    format (' Random bits:   ',32i2)
  93.  
  94.         end
  95.