 Technical Briefs
ROHAN help Files
IMSL - Compiling and Executing the Program
IMPORTANT: Before compiling any program using IMSL functions or subroutines
you must issue the command
source /opt/imsl/ipt/bin/iptsetup.csh
Add this line to your .login file, so that it is run automatically every time you log in.
Below is a record of a ROHAN terminal session during which the program
is compiled and executed. To start recording a session on rohan, type
script; to end recording, type exit. Optionally, you can specify a
filename for the session record in the script command.
Note that the program duvsta.f specifies that results are to be printed
to the terminal (actually, to the file a.out which can then be viewed as
below). Alternatively, you could write the results to an output file;
the DUVSTA subroutine stores results in the two-dimensional STAT array.
Script started on Mon Oct 02 14:38:19 1995
rohan% f77 -limsl -lsocket duvsta.f
duvsta.f:
MAIN sumstats:
rohan% a.out
Univariate Statistics from UVSTA
Variable Mean Variance Std. Dev. Skewness Kurtosis
1 72.3873 1471.6868 38.3626 0.9082 0.3221
Variable Minimum Maximum Range Coef. Var. Count
1 17.8800 173.4000 155.5200 0.5300 22.0000
Variable Lower CLM Upper CLM Lower CLV Upper CLV
1 58.3134 86.4611 945.9712 2666.2590
rohan% exit
rohan%
script done on Mon Oct 02 14:54:29 1995
Writing the Main Program
Here is the listing of a Fortran program that invokes IMSL to compute
summary statistics for a set of data:
program sumstats
C -----------------------------------------------------------------
C This program computes basic univariate statistics for a data set
C contained in the input file DATA.DAT. The IMSL subroutine DUVSTA
C is called to perform the calculations in double precision.
C
C Written by John Q. Public -- February 1990.
C Most recent modification -- September 1995.
C -----------------------------------------------------------------
parameter(ldstat=15, ldx=22, nvar=1)
implicit double precision (a-h,o-r,t-z)
double precision x(ldx,nvar),stat(ldstat,nvar)
C -----------------------------------------------------------------
C Open the input file and read the data.
C -----------------------------------------------------------------
open(unit=13,file='data.dat',status='old')
do i=1,ldx
do j=1,nvar
read(13,*,end=99) x(i,j)
end do
end do
99 close(unit=13,dispose='keep')
C -----------------------------------------------------------------
C Specify values for the input arguments to DUVSTA, then issue
C the subroutine call to perform the computations and print
C the results out to the screen.
C -----------------------------------------------------------------
ido=0
nrow=ldx
ifrq=0
iwt=0
mopt=1
conprm=90.0d0
conprv=90.0d0
iprint=1
call duvsta(ido,nrow,nvar,x,ldx,ifrq,iwt,mopt,
+ conprm,conprv,iprint,stat,ldstat,nrmiss)
end
|
|
ROHAN Help Index