dir    templates
for     Base Directory

file    BiCG.c
for     /* BiCG.f -- translated by f2c (version of 20 August 1993  13:15:44).
,          You must link the resulting object file with the libraries:
,       	-lf2c -lm   (in that order)
,       */
,       
,       #include "f2c.h"
,       
,       /* Table of constant values */
,       
,       static integer c__1 = 1;
size    9748

file    BiCGSTAB.c
for     /* BiCGSTAB.f -- translated by f2c (version of 20 August 1993  13:15:44).
,          You must link the resulting object file with the libraries:
,       	-lf2c -lm   (in that order)
,       */
,       
,       #include "f2c.h"
,       
,       /* Table of constant values */
,       
,       static integer c__1 = 1;
size    9662

file    CG.c
for     /* CG.f -- translated by f2c (version of 20 August 1993  13:15:44).
,          You must link the resulting object file with the libraries:
,       	-lf2c -lm   (in that order)
,       */
,       
,       #include "f2c.h"
,       
,       /* Table of constant values */
,       
,       static integer c__1 = 1;
size    7165

file    CGS.c
for     /* CGS.f -- translated by f2c (version of 20 August 1993  13:15:44).
,          You must link the resulting object file with the libraries:
,       	-lf2c -lm   (in that order)
,       */
,       
,       #include "f2c.h"
,       
,       /* Table of constant values */
,       
,       static integer c__1 = 1;
size    9573

file    Cheby.c
for     /* Cheby.f -- translated by f2c (version of 20 August 1993  13:15:44).
,          You must link the resulting object file with the libraries:
,       	-lf2c -lm   (in that order)
,       */
,       
,       #include "f2c.h"
,       
,       /* Table of constant values */
,       
,       static integer c__1 = 1;
size    7162

file    Notice
for     /****************************************************************
,       Copyright 1990, 1991, 1992, 1993 by AT&T Bell Laboratories and Bellcore.
,       
,       Permission to use, copy, modify, and distribute this software
,       and its documentation for any purpose and without fee is hereby
,       granted, provided that the above copyright notice appear in all
,       copies and that both that the copyright notice and this
,       permission notice and warranty disclaimer appear in supporting
,       documentation, and that the names of AT&T Bell Laboratories or
,       Bellcore or any of their entities not be used in advertising or
size    1195

file    GMRES.c
for     /* GMRES.f -- translated by f2c (version of 20 August 1993  13:15:44).
,          You must link the resulting object file with the libraries:
,       	-lf2c -lm   (in that order)
,       */
,       
,       #include "f2c.h"
,       
,       /* Table of constant values */
,       
,       static integer c__1 = 1;
size    12284

file    Jacobi.c
for     /* Jacobi.f -- translated by f2c (version of 20 August 1993  13:15:44).
,          You must link the resulting object file with the libraries:
,       	-lf2c -lm   (in that order)
,       */
,       
,       #include "f2c.h"
,       
,       /* Table of constant values */
,       
,       static real c_b2 = (float)1.;
size    6070

file    Makefile
for     FC = f77
,       
,       FFLAGS= -g -u
,       
,       libs = slapack.o sblas.o
,       
,       tester_o  = Tester.o MatGen.o sspdchk.o snsychk.o \
,                   Utils.o \
,                   MatVec.o PSolve.o \
,       			CG.o CGREVCOM.o \
size    560

file    MatGen.c
for     /* MatGen.f -- translated by f2c (version of 20 August 1993  13:15:44).
,          You must link the resulting object file with the libraries:
,       	-lf2c -lm   (in that order)
,       */
,       
,       #include "f2c.h"
,       
,       /* Table of constant values */
,       
,       static doublereal c_b14 = 2.7182817459106;
size    21747

file    MatVec.c
for     /* MatVec.f -- translated by f2c (version of 20 August 1993  13:15:44).
,          You must link the resulting object file with the libraries:
,       	-lf2c -lm   (in that order)
,       */
,       
,       #include "f2c.h"
,       
,       /* Common Block Declarations */
,       
,       struct {
size    1767

file    PSolve.c
for     /* PSolve.f -- translated by f2c (version of 20 August 1993  13:15:44).
,          You must link the resulting object file with the libraries:
,       	-lf2c -lm   (in that order)
,       */
,       
,       #include "f2c.h"
,       
,       /* Common Block Declarations */
,       
,       struct {
size    9866

file    QMR.c
for     /* QMR.f -- translated by f2c (version of 20 August 1993  13:15:44).
,          You must link the resulting object file with the libraries:
,       	-lf2c -lm   (in that order)
,       */
,       
,       #include "f2c.h"
,       
,       /* Table of constant values */
,       
,       static integer c__1 = 1;
size    14420

file    SOR.c
for     /* SOR.f -- translated by f2c (version of 20 August 1993  13:15:44).
,          You must link the resulting object file with the libraries:
,       	-lf2c -lm   (in that order)
,       */
,       
,       #include "f2c.h"
,       
,       /* Table of constant values */
,       
,       static integer c__1 = 1;
size    7609

file    Tester.c
for     /* Tester.f -- translated by f2c (version of 20 August 1993  13:15:44).
,          You must link the resulting object file with the libraries:
,       	-lf2c -lm   (in that order)
,       */
,       
,       #include "f2c.h"
,       
,       /* Common Block Declarations */
,       
,       struct {
size    50225

file    Utils.c
for     /* Utils.f -- translated by f2c (version of 20 August 1993  13:15:44).
,          You must link the resulting object file with the libraries:
,       	-lf2c -lm   (in that order)
,       */
,       
,       #include "f2c.h"
,       
,       /* Common Block Declarations */
,       
,       struct {
size    10835

file    f2c.h
for     /* f2c.h  --  Standard Fortran to C header file */
,       
,       /**  barf  [ba:rf]  2.  "He suggested using FORTRAN, and everybody barfed."
,       
,       	- From The Shogakukan DICTIONARY OF NEW ENGLISH (Second edition) */
,       
,       #ifndef F2C_INCLUDE
,       #define F2C_INCLUDE
,       
,       typedef long int integer;
size    4206

file    sblas.c
for     /* sblas.f -- translated by f2c (version of 20 August 1993  13:15:44).
,          You must link the resulting object file with the libraries:
,       	-lf2c -lm   (in that order)
,       */
,       
,       #include "f2c.h"
,       
,       /* Table of constant values */
,       
,       static real c_b100 = (float)1.;
size    89107

file    slapack.c
for     /* slapack.f -- translated by f2c (version of 20 August 1993  13:15:44).
,          You must link the resulting object file with the libraries:
,       	-lf2c -lm   (in that order)
,       */
,       
,       #include "f2c.h"
,       
,       /* Table of constant values */
,       
,       static integer c__1 = 1;
size    198998

file    snsychk.c
for     /* snsychk.f -- translated by f2c (version of 20 August 1993  13:15:44).
,          You must link the resulting object file with the libraries:
,       	-lf2c -lm   (in that order)
,       */
,       
,       #include "f2c.h"
,       
,       /* Common Block Declarations */
,       
,       struct {
size    12582

file    sspdchk.c
for     /* sspdchk.f -- translated by f2c (version of 20 August 1993  13:15:44).
,          You must link the resulting object file with the libraries:
,       	-lf2c -lm   (in that order)
,       */
,       
,       #include "f2c.h"
,       
,       /* Common Block Declarations */
,       
,       struct {
size    14797

file    test.data
for     1.0E-7                             CONVERGENCE TOLERANCE
,       1.0E+1                             SCALED RESIDUAL TOLERANCE
,       CG     T PUT F FOR NO TEST.        ALGORITHMS TO BE TESTED
,       CHEBY  T PUT F FOR NO TEST.
,       SOR    T PUT F FOR NO TEST.
,       BICG   T PUT F FOR NO TEST.
,       CGS    T PUT F FOR NO TEST.
,       BICGS  T PUT F FOR NO TEST.
,       GMRES  T PUT F FOR NO TEST.
,       QMR    T PUT F FOR NO TEST.
size    1564

file    README
for     =======================
,       SFTEMPLATES README FILE
,       =======================
,       
,       
,       This directory contains SINGLE PRECISION code for the algorithms
,       discussed in the Templates book, implemented using the Reverse
,       Communicaiton mechanism.
,       
,       Unpack the shar file as:
size    1714

file    Version.c
for     static char junk[] = "\n@(#)LIBF77 VERSION 2.01 12 March 1993\n";
,       
,       /*
,       2.00	11 June 1980.  File version.c added to library.
,       2.01	31 May 1988.  s_paus() flushes stderr; names of hl_* fixed
,       	[ d]erf[c ] added
,       	 8 Aug. 1989: #ifdefs for f2c -i2 added to s_cat.c
,       	29 Nov. 1989: s_cmp returns long (for f2c)
,       	30 Nov. 1989: arg types from f2c.h
,       	12 Dec. 1989: s_rnge allows long names
size    1247

file    backspace.c
for     #include "f2c.h"
,       #include "fio.h"
,       #ifdef KR_headers
,       integer f_back(a) alist *a;
,       #else
,       integer f_back(alist *a)
,       #endif
,       {	unit *b;
,       	int i, n, ndec;
,       	long x;
size    1341

file    close.c
for     #include "f2c.h"
,       #include "fio.h"
,       #ifdef KR_headers
,       integer f_clos(a) cllist *a;
,       #else
,       #undef abs
,       #undef min
,       #undef max
,       #include "stdlib.h"
,       #ifdef NON_UNIX_STDIO
size    1362

file    dfe.c
for     #include "f2c.h"
,       #include "fio.h"
,       #include "fmt.h"
,       
,       y_rsk(Void)
,       {
,       	if(f__curunit->uend || f__curunit->url <= f__recpos
,       		|| f__curunit->url == 1) return 0;
,       	do {
,       		getc(f__cf);
size    2903

file    dolio.c
for     #include "f2c.h"
,       
,       #ifdef __cplusplus
,       extern "C" {
,       #endif
,       #ifdef KR_headers
,       extern int (*f__lioproc)();
,       
,       integer do_lio(type,number,ptr,len) ftnint *number,*type; char *ptr; ftnlen len;
,       #else
size    404

file    due.c
for     #include "f2c.h"
,       #include "fio.h"
,       
,       #ifdef KR_headers
,       c_due(a) cilist *a;
,       #else
,       c_due(cilist *a)
,       #endif
,       {
,       	if(!f__init) f_init();
size    1430

file    endfile.c
for     #include "f2c.h"
,       #include "fio.h"
,       #include "sys/types.h"
,       #include "rawio.h"
,       
,       #ifdef KR_headers
,       extern char *strcpy();
,       #else
,       #undef abs
,       #undef min
size    3370

file    err.c
for     #ifndef NON_UNIX_STDIO
,       #include "sys/types.h"
,       #include "sys/stat.h"
,       #endif
,       #include "f2c.h"
,       #include "fio.h"
,       #include "fmt.h"	/* for struct syl */
,       #include "rawio.h"	/* for fcntl.h, fdopen */
,       #ifdef NON_UNIX_STDIO
,       #ifdef KR_headers
size    6262

file    f2ch.add
for     /* If you are using a C++ compiler, append the following to f2c.h
,          for compiling libF77 and libI77. */
,       
,       #ifdef __cplusplus
,       extern "C" {
,       extern int abort_(void);
,       extern double c_abs(complex *);
,       extern void c_cos(complex *, complex *);
,       extern void c_div(complex *, complex *, complex *);
,       extern void c_exp(complex *, complex *);
size    6056

file    fio.h
for     #include "stdio.h"
,       #include "errno.h"
,       #ifndef NULL
,       /* ANSI C */
,       #include "stddef.h"
,       #endif
,       
,       #ifndef SEEK_SET
,       #define SEEK_SET 0
,       #define SEEK_CUR 1
size    2309

file    fmt.c
for     #include "f2c.h"
,       #include "fio.h"
,       #include "fmt.h"
,       #define skip(s) while(*s==' ') s++
,       #ifdef interdata
,       #define SYLMX 300
,       #endif
,       #ifdef pdp11
,       #define SYLMX 300
,       #endif
size    7869

file    fmt.h
for     struct f__syl
,       {	int op,p1,p2,p3;
,       };
,       #define RET1 1
,       #define REVERT 2
,       #define GOTO 3
,       #define X 4
,       #define SLASH 5
,       #define STACK 6
,       #define I 7
size    1649

file    fmtlib.c
for     /*	@(#)fmtlib.c	1.2	*/
,       #define MAXINTLENGTH 23
,       #ifdef KR_headers
,       char *f__icvt(value,ndigit,sign, base) long value; int *ndigit,*sign;
,        register int base;
,       #else
,       char *f__icvt(long value, int *ndigit, int *sign, int base)
,       #endif
,       {	static char buf[MAXINTLENGTH+1];
,       	register int i;
size    582

file    fp.h
for     #define FMAX 40
,       #define EXPMAXDIGS 8
,       #define EXPMAX 99999999
,       /* FMAX = max number of nonzero digits passed to atof() */
,       /* EXPMAX = 10^EXPMAXDIGS - 1 = largest allowed exponent absolute value */
,       
,       #include "local.h"
,       
,       /* MAXFRACDIGS and MAXINTDIGS are for wrt_F -- bounds (not necessarily
,          tight) on the maximum number of digits to the right and left of
size    613

file    iio.c
for     #include "f2c.h"
,       #include "fio.h"
,       #include "fmt.h"
,       extern char *f__icptr;
,       char *f__icend;
,       extern icilist *f__svic;
,       int f__icnum;
,       extern int f__hiwater;
,       z_getc(Void)
,       {
size    2287

file    ilnw.c
for     #include "f2c.h"
,       #include "fio.h"
,       #include "lio.h"
,       extern char *f__icptr;
,       extern char *f__icend;
,       extern icilist *f__svic;
,       extern int f__icnum;
,       #ifdef KR_headers
,       extern int z_putc();
,       #else
size    1049

file    inquire.c
for     #include "f2c.h"
,       #include "fio.h"
,       #ifdef KR_headers
,       integer f_inqu(a) inlist *a;
,       #else
,       #ifdef MSDOS
,       #undef abs
,       #undef min
,       #undef max
,       #include "string.h"
size    2536

file    libI77.xsum
for     Notice	1211689a	1195
,       README	140275a4	6531
,       Version.c	e5d880dd	7661
,       backspace.c	14a32c39	1341
,       close.c	1467fa93	1362
,       dfe.c	3c6b216	2903
,       dolio.c	17595b24	404
,       due.c	1bbe319b	1430
,       endfile.c	f266a8f6	3370
,       err.c	f411968b	6262
size    868

file    lio.h
for     /*	copy of ftypes from the compiler */
,       /* variable types
,        * numeric assumptions:
,        *	int < reals < complexes
,        *	TYDREAL-TYREAL = TYDCOMPLEX-TYCOMPLEX
,        */
,       
,       /* 0-10 retain their old (pre LOGICAL*1, etc.) */
,       /* values to allow mixing old and new objects. */
,       
size    1550

file    local.h
for     
size    1

file    lread.c
for     #include "f2c.h"
,       #include "fio.h"
,       #include "fmt.h"
,       #include "lio.h"
,       #include "ctype.h"
,       #include "fp.h"
,       
,       extern char *f__fmtbuf;
,       #ifdef KR_headers
,       extern double atof();
size    10998

file    lwrite.c
for     #include "f2c.h"
,       #include "fio.h"
,       #include "fmt.h"
,       #include "lio.h"
,       ftnint L_len;
,       
,       #ifdef KR_headers
,       t_putc(c)
,       #else
,       t_putc(int c)
size    3943

file    makefile
for     .SUFFIXES: .c .o
,       CC = gcc
,       SHELL = /bin/sh
,       CFLAGS = -O
,       
,       # compile, then strip unnecessary symbols
,       .c.o:
,       	$(CC) -c -DSkip_f2c_Undefs $(CFLAGS) $*.c
,       	ld -r -x -o $*.xxx $*.o
,       	mv $*.xxx $*.o
size    3097

file    open.c
for     #include "sys/types.h"
,       #include "sys/stat.h"
,       #include "f2c.h"
,       #include "fio.h"
,       #include "string.h"
,       #include "rawio.h"
,       
,       #ifdef KR_headers
,       extern char *malloc(), *mktemp();
,       extern integer f_clos();
size    4356

file    rawio.h
for     #ifdef KR_headers
,       extern FILE *fdopen();
,       #else
,       #ifdef MSDOS
,       #include "io.h"
,       #define close _close
,       #define creat _creat
,       #define open _open
,       #define read _read
,       #define write _write
size    807

file    rdfmt.c
for     #include "f2c.h"
,       #include "fio.h"
,       #include "fmt.h"
,       #include "fp.h"
,       
,       extern int f__cursor;
,       #ifdef KR_headers
,       extern double atof();
,       #else
,       #undef abs
size    8356

file    rewind.c
for     #include "f2c.h"
,       #include "fio.h"
,       #ifdef KR_headers
,       integer f_rew(a) alist *a;
,       #else
,       integer f_rew(alist *a)
,       #endif
,       {
,       	unit *b;
,       	if(a->aunit>=MXUNIT || a->aunit<0)
size    408

file    rsfe.c
for     /* read sequential formatted external */
,       #include "f2c.h"
,       #include "fio.h"
,       #include "fmt.h"
,       
,       xrd_SL(Void)
,       {	int ch;
,       	if(!f__curunit->uend)
,       		while((ch=getc(f__cf))!='\n' && ch!=EOF);
,       	f__cursor=f__recpos=0;
size    1299

file    rsli.c
for     #include "f2c.h"
,       #include "fio.h"
,       #include "lio.h"
,       #include "fmt.h" /* for f__doend */
,       
,       extern flag f__lquit;
,       extern int f__lcount;
,       extern char *f__icptr;
,       extern char *f__icend;
,       extern icilist *f__svic;
size    1688

file    rsne.c
for     #include "f2c.h"
,       #include "fio.h"
,       #include "lio.h"
,       
,       #define MAX_NL_CACHE 3	/* maximum number of namelist hash tables to cache */
,       #define MAXDIM 20	/* maximum number of subscripts */
,       
,        struct dimen {
,       	ftnlen extent;
,       	ftnlen curval;
size    10597

file    sfe.c
for     /* sequential formatted external common routines*/
,       #include "f2c.h"
,       #include "fio.h"
,       
,       extern char *f__fmtbuf;
,       
,       integer e_rsfe(Void)
,       {	int n;
,       	n=en_fio();
,       	if (f__cf == stdout)
size    638

file    sue.c
for     #include "f2c.h"
,       #include "fio.h"
,       extern uiolen f__reclen;
,       long f__recloc;
,       
,       #ifdef KR_headers
,       c_sue(a) cilist *a;
,       #else
,       c_sue(cilist *a)
,       #endif
size    1740

file    typesize.c
for     #include "f2c.h"
,       
,       ftnlen f__typesize[] = { 0, 0, sizeof(shortint), sizeof(integer),
,       			sizeof(real), sizeof(doublereal),
,       			sizeof(complex), sizeof(doublecomplex),
,       			sizeof(logical), sizeof(char) };
size    200

file    uio.c
for     #include "f2c.h"
,       #include "fio.h"
,       uiolen f__reclen;
,       
,       #ifdef KR_headers
,       do_us(number,ptr,len) ftnint *number; char *ptr; ftnlen len;
,       #else
,       do_us(ftnint *number, char *ptr, ftnlen len)
,       #endif
,       {
size    1265

file    util.c
for     #ifndef NON_UNIX_STDIO
,       #include "sys/types.h"
,       #include "sys/stat.h"
,       #endif
,       #include "f2c.h"
,       #include "fio.h"
,       
,        VOID
,       #ifdef KR_headers
,       g_char(a,alen,b) char *a,*b; ftnlen alen;
size    1354

file    wref.c
for     #include "f2c.h"
,       #include "fio.h"
,       #include "fmt.h"
,       #include "fp.h"
,       #ifndef VAX
,       #include "ctype.h"
,       #endif
,       
,       #ifndef KR_headers
,       #undef abs
size    4028

file    wrtfmt.c
for     #include "f2c.h"
,       #include "fio.h"
,       #include "fmt.h"
,       extern int f__cursor;
,       #ifdef KR_headers
,       extern char *f__icvt();
,       #else
,       extern char *f__icvt(long, int*, int*, int);
,       #endif
,       int f__hiwater;
size    8023

file    wsfe.c
for     /*write sequential formatted external*/
,       #include "f2c.h"
,       #include "fio.h"
,       #include "fmt.h"
,       extern int f__hiwater;
,       
,       #ifdef KR_headers
,       x_putc(c)
,       #else
,       x_putc(int c)
size    1581

file    wsle.c
for     #include "f2c.h"
,       #include "fio.h"
,       #include "fmt.h"
,       #include "lio.h"
,       
,       #ifdef KR_headers
,       integer s_wsle(a) cilist *a;
,       #else
,       integer s_wsle(cilist *a)
,       #endif
size    611

file    wsne.c
for     #include "f2c.h"
,       #include "fio.h"
,       #include "lio.h"
,       
,        integer
,       #ifdef KR_headers
,       s_wsne(a) cilist *a;
,       #else
,       s_wsne(cilist *a)
,       #endif
size    438

file    xwsne.c
for     #include "f2c.h"
,       #include "fio.h"
,       #include "lio.h"
,       #include "fmt.h"
,       
,        static VOID
,       nl_donewrec(Void)
,       {
,       	(*f__donewrec)();
,       	PUT(' ');
size    1080

file    abort_.c
for     #include "stdio.h"
,       #include "f2c.h"
,       
,       #ifdef KR_headers
,       extern VOID sig_die();
,       
,       int abort_()
,       #else
,       extern void sig_die(char*,int);
,       
size    239

file    c_abs.c
for     #include "f2c.h"
,       
,       #ifdef KR_headers
,       extern double f__cabs();
,       
,       double c_abs(z) complex *z;
,       #else
,       extern double f__cabs(double, double);
,       
,       double c_abs(complex *z)
size    205

file    c_cos.c
for     #include "f2c.h"
,       
,       #ifdef KR_headers
,       extern double sin(), cos(), sinh(), cosh();
,       
,       VOID c_cos(r, z) complex *r, *z;
,       #else
,       #undef abs
,       #include "math.h"
,       
size    260

file    c_div.c
for     #include "f2c.h"
,       
,       #ifdef KR_headers
,       extern VOID sig_die();
,       VOID c_div(c, a, b)
,       complex *a, *b, *c;
,       #else
,       extern void sig_die(char*,int);
,       void c_div(complex *c, complex *a, complex *b)
,       #endif
size    665

file    c_exp.c
for     #include "f2c.h"
,       
,       #ifdef KR_headers
,       extern double exp(), cos(), sin();
,       
,        VOID c_exp(r, z) complex *r, *z;
,       #else
,       #undef abs
,       #include "math.h"
,       
size    270

file    c_log.c
for     #include "f2c.h"
,       
,       #ifdef KR_headers
,       extern double log(), f__cabs(), atan2();
,       VOID c_log(r, z) complex *r, *z;
,       #else
,       #undef abs
,       #include "math.h"
,       extern double f__cabs(double, double);
,       
size    292

file    c_sin.c
for     #include "f2c.h"
,       
,       #ifdef KR_headers
,       extern double sin(), cos(), sinh(), cosh();
,       
,       VOID c_sin(r, z) complex *r, *z;
,       #else
,       #undef abs
,       #include "math.h"
,       
size    258

file    c_sqrt.c
for     #include "f2c.h"
,       
,       #ifdef KR_headers
,       extern double sqrt(), f__cabs();
,       
,       VOID c_sqrt(r, z) complex *r, *z;
,       #else
,       #undef abs
,       #include "math.h"
,       extern double f__cabs(double, double);
size    505

file    cabs.c
for     #ifdef KR_headers
,       extern double sqrt();
,       double f__cabs(real, imag) double real, imag;
,       #else
,       #undef abs
,       #include "math.h"
,       double f__cabs(double real, double imag)
,       #endif
,       {
,       double temp;
size    427

file    d_abs.c
for     #include "f2c.h"
,       
,       #ifdef KR_headers
,       double d_abs(x) doublereal *x;
,       #else
,       double d_abs(doublereal *x)
,       #endif
,       {
,       if(*x >= 0)
,       	return(*x);
size    151

file    d_acos.c
for     #include "f2c.h"
,       
,       #ifdef KR_headers
,       double acos();
,       double d_acos(x) doublereal *x;
,       #else
,       #undef abs
,       #include "math.h"
,       double d_acos(doublereal *x)
,       #endif
size    178

file    d_asin.c
for     #include "f2c.h"
,       
,       #ifdef KR_headers
,       double asin();
,       double d_asin(x) doublereal *x;
,       #else
,       #undef abs
,       #include "math.h"
,       double d_asin(doublereal *x)
,       #endif
size    178

file    d_atan.c
for     #include "f2c.h"
,       
,       #ifdef KR_headers
,       double atan();
,       double d_atan(x) doublereal *x;
,       #else
,       #undef abs
,       #include "math.h"
,       double d_atan(doublereal *x)
,       #endif
size    178

file    d_atn2.c
for     #include "f2c.h"
,       
,       #ifdef KR_headers
,       double atan2();
,       double d_atn2(x,y) doublereal *x, *y;
,       #else
,       #undef abs
,       #include "math.h"
,       double d_atn2(doublereal *x, doublereal *y)
,       #endif
size    204

file    d_cnjg.c
for     #include "f2c.h"
,       
,        VOID
,       #ifdef KR_headers
,       d_cnjg(r, z) doublecomplex *r, *z;
,       #else
,       d_cnjg(doublecomplex *r, doublecomplex *z)
,       #endif
,       {
,       r->r = z->r;
size    165

file    d_cos.c
for     #include "f2c.h"
,       
,       #ifdef KR_headers
,       double cos();
,       double d_cos(x) doublereal *x;
,       #else
,       #undef abs
,       #include "math.h"
,       double d_cos(doublereal *x)
,       #endif
size    174

file    d_cosh.c
for     #include "f2c.h"
,       
,       #ifdef KR_headers
,       double cosh();
,       double d_cosh(x) doublereal *x;
,       #else
,       #undef abs
,       #include "math.h"
,       double d_cosh(doublereal *x)
,       #endif
size    178

file    d_dim.c
for     #include "f2c.h"
,       
,       #ifdef KR_headers
,       double d_dim(a,b) doublereal *a, *b;
,       #else
,       double d_dim(doublereal *a, doublereal *b)
,       #endif
,       {
,       return( *a > *b ? *a - *b : 0);
,       }
size    165

file    d_exp.c
for     #include "f2c.h"
,       
,       #ifdef KR_headers
,       double exp();
,       double d_exp(x) doublereal *x;
,       #else
,       #undef abs
,       #include "math.h"
,       double d_exp(doublereal *x)
,       #endif
size    174

file    d_imag.c
for     #include "f2c.h"
,       
,       #ifdef KR_headers
,       double d_imag(z) doublecomplex *z;
,       #else
,       double d_imag(doublecomplex *z)
,       #endif
,       {
,       return(z->i);
,       }
size    134

file    d_int.c
for     #include "f2c.h"
,       
,       #ifdef KR_headers
,       double floor();
,       double d_int(x) doublereal *x;
,       #else
,       #undef abs
,       #include "math.h"
,       double d_int(doublereal *x)
,       #endif
size    202

file    d_lg10.c
for     #include "f2c.h"
,       
,       #define log10e 0.43429448190325182765
,       
,       #ifdef KR_headers
,       double log();
,       double d_lg10(x) doublereal *x;
,       #else
,       #undef abs
,       #include "math.h"
size    224

file    d_log.c
for     #include "f2c.h"
,       
,       #ifdef KR_headers
,       double log();
,       double d_log(x) doublereal *x;
,       #else
,       #undef abs
,       #include "math.h"
,       double d_log(doublereal *x)
,       #endif
size    174

file    d_mod.c
for     #include "f2c.h"
,       
,       #ifdef KR_headers
,       #ifdef IEEE_drem
,       double drem();
,       #else
,       double floor();
,       #endif
,       double d_mod(x,y) doublereal *x, *y;
,       #else
size    621

file    d_nint.c
for     #include "f2c.h"
,       
,       #ifdef KR_headers
,       double floor();
,       double d_nint(x) doublereal *x;
,       #else
,       #undef abs
,       #include "math.h"
,       double d_nint(doublereal *x)
,       #endif
size    214

file    d_prod.c
for     #include "f2c.h"
,       
,       #ifdef KR_headers
,       double d_prod(x,y) real *x, *y;
,       #else
,       double d_prod(real *x, real *y)
,       #endif
,       {
,       return( (*x) * (*y) );
,       }
size    140

file    d_sign.c
for     #include "f2c.h"
,       
,       #ifdef KR_headers
,       double d_sign(a,b) doublereal *a, *b;
,       #else
,       double d_sign(doublereal *a, doublereal *b)
,       #endif
,       {
,       double x;
,       x = (*a >= 0 ? *a : - *a);
size    199

file    d_sin.c
for     #include "f2c.h"
,       
,       #ifdef KR_headers
,       double sin();
,       double d_sin(x) doublereal *x;
,       #else
,       #undef abs
,       #include "math.h"
,       double d_sin(doublereal *x)
,       #endif
size    174

file    d_sinh.c
for     #include "f2c.h"
,       
,       #ifdef KR_headers
,       double sinh();
,       double d_sinh(x) doublereal *x;
,       #else
,       #undef abs
,       #include "math.h"
,       double d_sinh(doublereal *x)
,       #endif
size    178

file    d_sqrt.c
for     #include "f2c.h"
,       
,       #ifdef KR_headers
,       double sqrt();
,       double d_sqrt(x) doublereal *x;
,       #else
,       #undef abs
,       #include "math.h"
,       double d_sqrt(doublereal *x)
,       #endif
size    178

file    d_tan.c
for     #include "f2c.h"
,       
,       #ifdef KR_headers
,       double tan();
,       double d_tan(x) doublereal *x;
,       #else
,       #undef abs
,       #include "math.h"
,       double d_tan(doublereal *x)
,       #endif
size    174

file    d_tanh.c
for     #include "f2c.h"
,       
,       #ifdef KR_headers
,       double tanh();
,       double d_tanh(x) doublereal *x;
,       #else
,       #undef abs
,       #include "math.h"
,       double d_tanh(doublereal *x)
,       #endif
size    178

file    derf_.c
for     #include "f2c.h"
,       
,       #ifdef KR_headers
,       double erf();
,       double derf_(x) doublereal *x;
,       #else
,       extern double erf(double);
,       double derf_(doublereal *x)
,       #endif
,       {
size    172

file    derfc_.c
for     #include "f2c.h"
,       
,       #ifdef KR_headers
,       extern double erfc();
,       
,       double derfc_(x) doublereal *x;
,       #else
,       extern double erfc(double);
,       
,       double derfc_(doublereal *x)
size    186

file    ef1asc_.c
for     /* EFL support routine to copy string b to string a */
,       
,       #include "f2c.h"
,       
,       
,       #define M	( (long) (sizeof(long) - 1) )
,       #define EVEN(x)	( ( (x)+ M) & (~M) )
,       
,       #ifdef KR_headers
,       extern VOID s_copy();
size    453

file    ef1cmc_.c
for     /* EFL support routine to compare two character strings */
,       
,       #include "f2c.h"
,       
,       #ifdef KR_headers
,       extern integer s_cmp();
,       integer ef1cmc_(a, la, b, lb) ftnint *a, *b; ftnlen *la, *lb;
,       #else
,       extern integer s_cmp(char*,char*,ftnlen,ftnlen);
,       integer ef1cmc_(ftnint *a, ftnlen *la, ftnint *b, ftnlen *lb)
size    360

file    erf_.c
for     #include "f2c.h"
,       
,       #ifdef KR_headers
,       double erf();
,       double erf_(x) real *x;
,       #else
,       extern double erf(double);
,       double erf_(real *x)
,       #endif
,       {
size    158

file    erfc_.c
for     #include "f2c.h"
,       
,       #ifdef KR_headers
,       double erfc();
,       double erfc_(x) real *x;
,       #else
,       extern double erfc(double);
,       double erfc_(real *x)
,       #endif
,       {
size    163

file    getarg_.c
for     #include "f2c.h"
,       
,       /*
,        * subroutine getarg(k, c)
,        * returns the kth unix command argument in fortran character
,        * variable argument c
,       */
,       
,       #ifdef KR_headers
,       VOID getarg_(n, s, ls) ftnint *n; register char *s; ftnlen ls;
size    495

file    getenv_.c
for     #include "f2c.h"
,       
,       /*
,        * getenv - f77 subroutine to return environment variables
,        *
,        * called by:
,        *	call getenv (ENV_NAME, char_var)
,        * where:
,        *	ENV_NAME is the name of an environment variable
,        *	char_var is a character variable which will receive
size    975

file    h_abs.c
for     #include "f2c.h"
,       
,       #ifdef KR_headers
,       shortint h_abs(x) shortint *x;
,       #else
,       shortint h_abs(shortint *x)
,       #endif
,       {
,       if(*x >= 0)
,       	return(*x);
size    151

file    h_dim.c
for     #include "f2c.h"
,       
,       #ifdef KR_headers
,       shortint h_dim(a,b) shortint *a, *b;
,       #else
,       shortint h_dim(shortint *a, shortint *b)
,       #endif
,       {
,       return( *a > *b ? *a - *b : 0);
,       }
size    163

file    h_dnnt.c
for     #include "f2c.h"
,       
,       #ifdef KR_headers
,       double floor();
,       shortint h_dnnt(x) doublereal *x;
,       #else
,       #undef abs
,       #include "math.h"
,       shortint h_dnnt(doublereal *x)
,       #endif
size    218

file    h_indx.c
for     #include "f2c.h"
,       
,       #ifdef KR_headers
,       shortint h_indx(a, b, la, lb) char *a, *b; ftnlen la, lb;
,       #else
,       shortint h_indx(char *a, char *b, ftnlen la, ftnlen lb)
,       #endif
,       {
,       ftnlen i, n;
,       char *s, *t, *bend;
size    375

file    h_len.c
for     #include "f2c.h"
,       
,       #ifdef KR_headers
,       shortint h_len(s, n) char *s; ftnlen n;
,       #else
,       shortint h_len(char *s, ftnlen n)
,       #endif
,       {
,       return(n);
,       }
size    138

file    h_mod.c
for     #include "f2c.h"
,       
,       #ifdef KR_headers
,       shortint h_mod(a,b) short *a, *b;
,       #else
,       shortint h_mod(short *a, short *b)
,       #endif
,       {
,       return( *a % *b);
,       }
size    140

file    h_nint.c
for     #include "f2c.h"
,       
,       #ifdef KR_headers
,       double floor();
,       shortint h_nint(x) real *x;
,       #else
,       #undef abs
,       #include "math.h"
,       shortint h_nint(real *x)
,       #endif
size    206

file    h_sign.c
for     #include "f2c.h"
,       
,       #ifdef KR_headers
,       shortint h_sign(a,b) shortint *a, *b;
,       #else
,       shortint h_sign(shortint *a, shortint *b)
,       #endif
,       {
,       shortint x;
,       x = (*a >= 0 ? *a : - *a);
size    199

file    hl_ge.c
for     #include "f2c.h"
,       
,       #ifdef KR_headers
,       extern integer s_cmp();
,       shortlogical hl_ge(a,b,la,lb) char *a, *b; ftnlen la, lb;
,       #else
,       extern integer s_cmp(char *, char *, ftnlen, ftnlen);
,       shortlogical hl_ge(char *a, char *b, ftnlen la, ftnlen lb)
,       #endif
,       {
size    279

file    hl_gt.c
for     #include "f2c.h"
,       
,       #ifdef KR_headers
,       extern integer s_cmp();
,       shortlogical hl_gt(a,b,la,lb) char *a, *b; ftnlen la, lb;
,       #else
,       extern integer s_cmp(char *, char *, ftnlen, ftnlen);
,       shortlogical hl_gt(char *a, char *b, ftnlen la, ftnlen lb)
,       #endif
,       {
size    278

file    hl_le.c
for     #include "f2c.h"
,       
,       #ifdef KR_headers
,       extern integer s_cmp();
,       shortlogical hl_le(a,b,la,lb) char *a, *b; ftnlen la, lb;
,       #else
,       extern integer s_cmp(char *, char *, ftnlen, ftnlen);
,       shortlogical hl_le(char *a, char *b, ftnlen la, ftnlen lb)
,       #endif
,       {
size    279

file    hl_lt.c
for     #include "f2c.h"
,       
,       #ifdef KR_headers
,       extern integer s_cmp();
,       shortlogical hl_lt(a,b,la,lb) char *a, *b; ftnlen la, lb;
,       #else
,       extern integer s_cmp(char *, char *, ftnlen, ftnlen);
,       shortlogical hl_lt(char *a, char *b, ftnlen la, ftnlen lb)
,       #endif
,       {
size    278

file    i_abs.c
for     #include "f2c.h"
,       
,       #ifdef KR_headers
,       integer i_abs(x) integer *x;
,       #else
,       integer i_abs(integer *x)
,       #endif
,       {
,       if(*x >= 0)
,       	return(*x);
size    147

file    i_dim.c
for     #include "f2c.h"
,       
,       #ifdef KR_headers
,       integer i_dim(a,b) integer *a, *b;
,       #else
,       integer i_dim(integer *a, integer *b)
,       #endif
,       {
,       return( *a > *b ? *a - *b : 0);
,       }
size    158

file    i_dnnt.c
for     #include "f2c.h"
,       
,       #ifdef KR_headers
,       double floor();
,       integer i_dnnt(x) doublereal *x;
,       #else
,       #undef abs
,       #include "math.h"
,       integer i_dnnt(doublereal *x)
,       #endif
size    216

file    i_indx.c
for     #include "f2c.h"
,       
,       #ifdef KR_headers
,       integer i_indx(a, b, la, lb) char *a, *b; ftnlen la, lb;
,       #else
,       integer i_indx(char *a, char *b, ftnlen la, ftnlen lb)
,       #endif
,       {
,       ftnlen i, n;
,       char *s, *t, *bend;
size    363

file    i_len.c
for     #include "f2c.h"
,       
,       #ifdef KR_headers
,       integer i_len(s, n) char *s; ftnlen n;
,       #else
,       integer i_len(char *s, ftnlen n)
,       #endif
,       {
,       return(n);
,       }
size    136

file    i_mod.c
for     #include "f2c.h"
,       
,       #ifdef KR_headers
,       integer i_mod(a,b) integer *a, *b;
,       #else
,       integer i_mod(integer *a, integer *b)
,       #endif
,       {
,       return( *a % *b);
,       }
size    144

file    i_nint.c
for     #include "f2c.h"
,       
,       #ifdef KR_headers
,       double floor();
,       integer i_nint(x) real *x;
,       #else
,       #undef abs
,       #include "math.h"
,       integer i_nint(real *x)
,       #endif
size    204

file    i_sign.c
for     #include "f2c.h"
,       
,       #ifdef KR_headers
,       integer i_sign(a,b) integer *a, *b;
,       #else
,       integer i_sign(integer *a, integer *b)
,       #endif
,       {
,       integer x;
,       x = (*a >= 0 ? *a : - *a);
size    193

file    iargc_.c
for     #include "f2c.h"
,       
,       #ifdef KR_headers
,       ftnint iargc_()
,       #else
,       ftnint iargc_(void)
,       #endif
,       {
,       extern int xargc;
,       return ( xargc - 1 );
size    129

file    l_ge.c
for     #include "f2c.h"
,       
,       #ifdef KR_headers
,       extern integer s_cmp();
,       logical l_ge(a,b,la,lb) char *a, *b; ftnlen la, lb;
,       #else
,       extern integer s_cmp(char *, char *, ftnlen, ftnlen);
,       logical l_ge(char *a, char *b, ftnlen la, ftnlen lb)
,       #endif
,       {
size    267

file    l_gt.c
for     #include "f2c.h"
,       
,       #ifdef KR_headers
,       extern integer s_cmp();
,       logical l_gt(a,b,la,lb) char *a, *b; ftnlen la, lb;
,       #else
,       extern integer s_cmp(char *, char *, ftnlen, ftnlen);
,       logical l_gt(char *a, char *b, ftnlen la, ftnlen lb)
,       #endif
,       {
size    266

file    l_le.c
for     #include "f2c.h"
,       
,       #ifdef KR_headers
,       extern integer s_cmp();
,       logical l_le(a,b,la,lb) char *a, *b; ftnlen la, lb;
,       #else
,       extern integer s_cmp(char *, char *, ftnlen, ftnlen);
,       logical l_le(char *a, char *b, ftnlen la, ftnlen lb)
,       #endif
,       {
size    267

file    l_lt.c
for     #include "f2c.h"
,       
,       #ifdef KR_headers
,       extern integer s_cmp();
,       logical l_lt(a,b,la,lb) char *a, *b; ftnlen la, lb;
,       #else
,       extern integer s_cmp(char *, char *, ftnlen, ftnlen);
,       logical l_lt(char *a, char *b, ftnlen la, ftnlen lb)
,       #endif
,       {
size    266

file    libF77.xsum
for     Notice	1211689a	1195
,       README	1b1001bf	3633
,       Version.c	fe780238	1247
,       abort_.c	eaf90dc0	239
,       c_abs.c	ecce7a47	205
,       c_cos.c	f2338a46	260
,       c_div.c	f780c50e	665
,       c_exp.c	e1b005d5	270
,       c_log.c	4050533	292
,       c_sin.c	f19855c9	258
size    2518

file    main.c
for     /* STARTUP PROCEDURE FOR UNIX FORTRAN PROGRAMS */
,       
,       #include "stdio.h"
,       #include "signal.h"
,       
,       #ifndef SIGIOT
,       #define SIGIOT SIGABRT
,       #endif
,       
,       #ifndef KR_headers
size    1969

file    pow_ci.c
for     #include "f2c.h"
,       
,       #ifdef KR_headers
,       VOID pow_ci(p, a, b) 	/* p = a**b  */
,        complex *p, *a; integer *b;
,       #else
,       extern void pow_zi(doublecomplex*, doublecomplex*, integer*);
,       void pow_ci(complex *p, complex *a, integer *b) 	/* p = a**b  */
,       #endif
,       {
size    345

file    pow_dd.c
for     #include "f2c.h"
,       
,       #ifdef KR_headers
,       double pow();
,       double pow_dd(ap, bp) doublereal *ap, *bp;
,       #else
,       #undef abs
,       #include "math.h"
,       double pow_dd(doublereal *ap, doublereal *bp)
,       #endif
size    209

file    pow_di.c
for     #include "f2c.h"
,       
,       #ifdef KR_headers
,       double pow_di(ap, bp) doublereal *ap; integer *bp;
,       #else
,       double pow_di(doublereal *ap, integer *bp)
,       #endif
,       {
,       double pow, x;
,       integer n;
size    360

file    pow_hh.c
for     #include "f2c.h"
,       
,       #ifdef KR_headers
,       shortint pow_hh(ap, bp) shortint *ap, *bp;
,       #else
,       shortint pow_hh(shortint *ap, shortint *bp)
,       #endif
,       {
,       	shortint pow, x, n;
,       
size    401

file    pow_ii.c
for     #include "f2c.h"
,       
,       #ifdef KR_headers
,       integer pow_ii(ap, bp) integer *ap, *bp;
,       #else
,       integer pow_ii(integer *ap, integer *bp)
,       #endif
,       {
,       	integer pow, x, n;
,       
size    395

file    pow_qq.c
for     #include "f2c.h"
,       
,       #ifdef KR_headers
,       longint pow_qq(ap, bp) longint *ap, *bp;
,       #else
,       longint pow_qq(longint *ap, longint *bp)
,       #endif
,       {
,       	longint pow, x, n;
,       
size    395

file    pow_ri.c
for     #include "f2c.h"
,       
,       #ifdef KR_headers
,       double pow_ri(ap, bp) real *ap; integer *bp;
,       #else
,       double pow_ri(real *ap, integer *bp)
,       #endif
,       {
,       double pow, x;
,       integer n;
size    347

file    pow_zi.c
for     #include "f2c.h"
,       
,       #ifdef KR_headers
,       VOID pow_zi(p, a, b) 	/* p = a**b  */
,        doublecomplex *p, *a; integer *b;
,       #else
,       extern void z_div(doublecomplex*, doublecomplex*, doublecomplex*);
,       void pow_zi(doublecomplex *p, doublecomplex *a, integer *b) 	/* p = a**b  */
,       #endif
,       {
size    694

file    pow_zz.c
for     #include "f2c.h"
,       
,       #ifdef KR_headers
,       double log(), exp(), cos(), sin(), atan2(), f__cabs();
,       VOID pow_zz(r,a,b) doublecomplex *r, *a, *b;
,       #else
,       #undef abs
,       #include "math.h"
,       extern double f__cabs(double,double);
,       void pow_zz(doublecomplex *r, doublecomplex *a, doublecomplex *b)
size    482

file    r_abs.c
for     #include "f2c.h"
,       
,       #ifdef KR_headers
,       double r_abs(x) real *x;
,       #else
,       double r_abs(real *x)
,       #endif
,       {
,       if(*x >= 0)
,       	return(*x);
size    139

file    r_acos.c
for     #include "f2c.h"
,       
,       #ifdef KR_headers
,       double acos();
,       double r_acos(x) real *x;
,       #else
,       #undef abs
,       #include "math.h"
,       double r_acos(real *x)
,       #endif
size    166

file    r_asin.c
for     #include "f2c.h"
,       
,       #ifdef KR_headers
,       double asin();
,       double r_asin(x) real *x;
,       #else
,       #undef abs
,       #include "math.h"
,       double r_asin(real *x)
,       #endif
size    166

file    r_atan.c
for     #include "f2c.h"
,       
,       #ifdef KR_headers
,       double atan();
,       double r_atan(x) real *x;
,       #else
,       #undef abs
,       #include "math.h"
,       double r_atan(real *x)
,       #endif
size    166

file    r_atn2.c
for     #include "f2c.h"
,       
,       #ifdef KR_headers
,       double atan2();
,       double r_atn2(x,y) real *x, *y;
,       #else
,       #undef abs
,       #include "math.h"
,       double r_atn2(real *x, real *y)
,       #endif
size    186

file    r_cnjg.c
for     #include "f2c.h"
,       
,       #ifdef KR_headers
,       VOID r_cnjg(r, z) complex *r, *z;
,       #else
,       VOID r_cnjg(complex *r, complex *z)
,       #endif
,       {
,       r->r = z->r;
,       r->i = - z->i;
size    151

file    r_cos.c
for     #include "f2c.h"
,       
,       #ifdef KR_headers
,       double cos();
,       double r_cos(x) real *x;
,       #else
,       #undef abs
,       #include "math.h"
,       double r_cos(real *x)
,       #endif
size    162

file    r_cosh.c
for     #include "f2c.h"
,       
,       #ifdef KR_headers
,       double cosh();
,       double r_cosh(x) real *x;
,       #else
,       #undef abs
,       #include "math.h"
,       double r_cosh(real *x)
,       #endif
size    166

file    r_dim.c
for     #include "f2c.h"
,       
,       #ifdef KR_headers
,       double r_dim(a,b) real *a, *b;
,       #else
,       double r_dim(real *a, real *b)
,       #endif
,       {
,       return( *a > *b ? *a - *b : 0);
,       }
size    147

file    r_exp.c
for     #include "f2c.h"
,       
,       #ifdef KR_headers
,       double exp();
,       double r_exp(x) real *x;
,       #else
,       #undef abs
,       #include "math.h"
,       double r_exp(real *x)
,       #endif
size    162

file    r_imag.c
for     #include "f2c.h"
,       
,       #ifdef KR_headers
,       double r_imag(z) complex *z;
,       #else
,       double r_imag(complex *z)
,       #endif
,       {
,       return(z->i);
,       }
size    122

file    r_int.c
for     #include "f2c.h"
,       
,       #ifdef KR_headers
,       double floor();
,       double r_int(x) real *x;
,       #else
,       #undef abs
,       #include "math.h"
,       double r_int(real *x)
,       #endif
size    190

file    r_lg10.c
for     #include "f2c.h"
,       
,       #define log10e 0.43429448190325182765
,       
,       #ifdef KR_headers
,       double log();
,       double r_lg10(x) real *x;
,       #else
,       #undef abs
,       #include "math.h"
size    212

file    r_log.c
for     #include "f2c.h"
,       
,       #ifdef KR_headers
,       double log();
,       double r_log(x) real *x;
,       #else
,       #undef abs
,       #include "math.h"
,       double r_log(real *x)
,       #endif
size    162

file    r_mod.c
for     #include "f2c.h"
,       
,       #ifdef KR_headers
,       #ifdef IEEE_drem
,       double drem();
,       #else
,       double floor();
,       #endif
,       double r_mod(x,y) real *x, *y;
,       #else
size    611

file    r_nint.c
for     #include "f2c.h"
,       
,       #ifdef KR_headers
,       double floor();
,       double r_nint(x) real *x;
,       #else
,       #undef abs
,       #include "math.h"
,       double r_nint(real *x)
,       #endif
size    202

file    r_sign.c
for     #include "f2c.h"
,       
,       #ifdef KR_headers
,       double r_sign(a,b) real *a, *b;
,       #else
,       double r_sign(real *a, real *b)
,       #endif
,       {
,       double x;
,       x = (*a >= 0 ? *a : - *a);
size    181

file    r_sin.c
for     #include "f2c.h"
,       
,       #ifdef KR_headers
,       double sin();
,       double r_sin(x) real *x;
,       #else
,       #undef abs
,       #include "math.h"
,       double r_sin(real *x)
,       #endif
size    162

file    r_sinh.c
for     #include "f2c.h"
,       
,       #ifdef KR_headers
,       double sinh();
,       double r_sinh(x) real *x;
,       #else
,       #undef abs
,       #include "math.h"
,       double r_sinh(real *x)
,       #endif
size    166

file    r_sqrt.c
for     #include "f2c.h"
,       
,       #ifdef KR_headers
,       double sqrt();
,       double r_sqrt(x) real *x;
,       #else
,       #undef abs
,       #include "math.h"
,       double r_sqrt(real *x)
,       #endif
size    166

file    r_tan.c
for     #include "f2c.h"
,       
,       #ifdef KR_headers
,       double tan();
,       double r_tan(x) real *x;
,       #else
,       #undef abs
,       #include "math.h"
,       double r_tan(real *x)
,       #endif
size    162

file    r_tanh.c
for     #include "f2c.h"
,       
,       #ifdef KR_headers
,       double tanh();
,       double r_tanh(x) real *x;
,       #else
,       #undef abs
,       #include "math.h"
,       double r_tanh(real *x)
,       #endif
size    166

file    s_cat.c
for     #include "f2c.h"
,       
,       #ifdef KR_headers
,       VOID s_cat(lp, rpp, rnp, np, ll) char *lp, *rpp[]; ftnlen rnp[], *np, ll;
,       #else
,       VOID s_cat(char *lp, char *rpp[], ftnlen rnp[], ftnlen *np, ftnlen ll)
,       #endif
,       {
,       ftnlen i, n, nc;
,       char *f__rp;
size    413

file    s_cmp.c
for     #include "f2c.h"
,       
,       /* compare two strings */
,       
,       #ifdef KR_headers
,       integer s_cmp(a, b, la, lb) register char *a, *b; ftnlen la, lb;
,       #else
,       integer s_cmp(register char *a, register char *b, ftnlen la, ftnlen lb)
,       #endif
,       {
size    609

file    s_copy.c
for     #include "f2c.h"
,       
,       /* assign strings:  a = b */
,       
,       #ifdef KR_headers
,       VOID s_copy(a, b, la, lb) register char *a, *b; ftnlen la, lb;
,       #else
,       void s_copy(register char *a, register char *b, ftnlen la, ftnlen lb)
,       #endif
,       {
size    397

file    z_abs.c
for     #include "f2c.h"
,       
,       #ifdef KR_headers
,       double f__cabs();
,       double z_abs(z) doublecomplex *z;
,       #else
,       double f__cabs(double, double);
,       double z_abs(doublecomplex *z)
,       #endif
,       {
size    201

file    s_paus.c
for     #include "stdio.h"
,       #include "f2c.h"
,       #define PAUSESIG 15
,       
,       #ifdef KR_headers
,       #define Void /* void */
,       #define Int /* int */
,       #else
,       #define Void void
,       #define Int int
size    1150

file    s_rnge.c
for     #include "stdio.h"
,       #include "f2c.h"
,       
,       /* called when a subscript is out of range */
,       
,       #ifdef KR_headers
,       extern VOID sig_die();
,       integer s_rnge(varn, offset, procn, line) char *varn, *procn; ftnint offset, line;
,       #else
,       extern VOID sig_die(char*,int);
size    680

file    s_stop.c
for     #include "stdio.h"
,       #include "f2c.h"
,       
,       #ifdef KR_headers
,       extern void f_exit();
,       VOID s_stop(s, n) char *s; ftnlen n;
,       #else
,       #undef abs
,       #include "stdlib.h"
,       #ifdef __cplusplus
size    465

file    sig_die.c
for     #include "stdio.h"
,       #include "signal.h"
,       
,       #ifndef SIGIOT
,       #define SIGIOT SIGABRT
,       #endif
,       
,       #ifdef KR_headers
,       void sig_die(s, kill) register char *s; int kill;
,       #else
size    546

file    signal_.c
for     #include "f2c.h"
,       
,       #ifdef KR_headers
,       typedef int (*sig_type)();
,       extern sig_type signal();
,       
,       ftnint signal_(sigp, proc) integer *sigp; sig_type proc;
,       #else
,       #include "signal.h"
,       typedef void (*sig_type)(int);
size    327

file    system_.c
for     /* f77 interface to system routine */
,       
,       #include "f2c.h"
,       
,       #ifdef KR_headers
,       system_(s, n) register char *s; ftnlen n;
,       #else
,       #undef abs
,       #include "stdlib.h"
,       system_(register char *s, ftnlen n)
size    374

file    z_cos.c
for     #include "f2c.h"
,       
,       #ifdef KR_headers
,       double sin(), cos(), sinh(), cosh();
,       VOID z_cos(r, z) doublecomplex *r, *z;
,       #else
,       #undef abs
,       #include "math.h"
,       void z_cos(doublecomplex *r, doublecomplex *z)
,       #endif
size    269

file    z_div.c
for     #include "f2c.h"
,       
,       #ifdef KR_headers
,       extern void sig_die();
,       VOID z_div(c, a, b) doublecomplex *a, *b, *c;
,       #else
,       extern void sig_die(char*, int);
,       void z_div(doublecomplex *c, doublecomplex *a, doublecomplex *b)
,       #endif
,       {
size    675

file    z_exp.c
for     #include "f2c.h"
,       
,       #ifdef KR_headers
,       double exp(), cos(), sin();
,       VOID z_exp(r, z) doublecomplex *r, *z;
,       #else
,       #undef abs
,       #include "math.h"
,       void z_exp(doublecomplex *r, doublecomplex *z)
,       #endif
size    278

file    z_log.c
for     #include "f2c.h"
,       
,       #ifdef KR_headers
,       double log(), f__cabs(), atan2();
,       VOID z_log(r, z) doublecomplex *r, *z;
,       #else
,       #undef abs
,       #include "math.h"
,       extern double f__cabs(double, double);
,       void z_log(doublecomplex *r, doublecomplex *z)
size    305

file    z_sin.c
for     #include "f2c.h"
,       
,       #ifdef KR_headers
,       double sin(), cos(), sinh(), cosh();
,       VOID z_sin(r, z) doublecomplex *r, *z;
,       #else
,       #undef abs
,       #include "math.h"
,       void z_sin(doublecomplex *r, doublecomplex *z)
,       #endif
size    267

file    z_sqrt.c
for     #include "f2c.h"
,       
,       #ifdef KR_headers
,       double sqrt(), f__cabs();
,       VOID z_sqrt(r, z) doublecomplex *r, *z;
,       #else
,       #undef abs
,       #include "math.h"
,       extern double f__cabs(double, double);
,       void z_sqrt(doublecomplex *r, doublecomplex *z)
size    492

file    BiCG.f
for     *
,             SUBROUTINE BICG( N, B, X, WORK, LDW, ITER, RESID, MATVEC, 
,            $                 MATVECTRANS, PSOLVE, PSOLVETRANS, INFO )
,       *
,       *
,       *  -- Iterative template routine --
,       *     Univ. of Tennessee and Oak Ridge National Laboratory
,       *     October 1, 1993
,       *     Details of this algorithm are described in "Templates for the 
,       *     Solution of Linear Systems: Building Blocks for Iterative 
size    8109

file    BiCGREVCOM.f
for     *
,             SUBROUTINE BICGREVCOM( N, B, X, WORK, LDW, ITER, RESID, INFO,
,            $                       NDX1, NDX2, SCLR1, SCLR2, IJOB)
,       *
,       *
,       *  -- Iterative template routine --
,       *     Univ. of Tennessee and Oak Ridge National Laboratory
,       *     October 1, 1993
,       *     Details of this algorithm are described in "Templates for the 
,       *     Solution of Linear Systems: Building Blocks for Iterative 
size    10373

file    BiCGSTAB.f
for     *
,             SUBROUTINE BICGSTAB( N, B, X, WORK, LDW, ITER, RESID, MATVEC,
,            $                     PSOLVE, INFO )
,       *
,       *  -- Iterative template routine --
,       *     Univ. of Tennessee and Oak Ridge National Laboratory
,       *     October 1, 1993
,       *     Details of this algorithm are described in "Templates for the 
,       *     Solution of Linear Systems: Building Blocks for Iterative 
,       *     Methods", Barrett, Berry, Chan, Demmel, Donato, Dongarra, 
size    6987

file    BiCGSTABREVCOM.f
for     *
,             SUBROUTINE BICGSTABREVCOM(N, B, X, WORK, LDW, ITER, RESID, INFO,
,            $                    NDX1, NDX2, SCLR1, SCLR2, IJOB)
,       *
,       *  -- Iterative template routine --
,       *     Univ. of Tennessee and Oak Ridge National Laboratory
,       *     October 1, 1993
,       *     Details of this algorithm are described in "Templates for the 
,       *     Solution of Linear Systems: Building Blocks for Iterative 
,       *     Methods", Barrett, Berry, Chan, Demmel, Donato, Dongarra, 
size    11714

file    CG.f
for     *
,             SUBROUTINE CG( N, B, X, WORK, LDW, ITER, RESID, MATVEC, 
,            $               PSOLVE, INFO )
,       *
,       *  -- Iterative template routine --
,       *     Univ. of Tennessee and Oak Ridge National Laboratory
,       *     October 1, 1993
,       *     Details of this algorithm are described in "Templates for the 
,       *     Solution of Linear Systems: Building Blocks for Iterative 
,       *     Methods", Barrett, Berry, Chan, Demmel, Donato, Dongarra, 
size    5861

file    CGREVCOM.f
for     *
,             SUBROUTINE CGREVCOM( N, B, X, WORK, LDW, ITER, RESID, INFO,
,            $                     NDX1, NDX2, SCLR1, SCLR2, IJOB)
,       *
,       *  -- Iterative template routine --
,       *     Univ. of Tennessee and Oak Ridge National Laboratory
,       *     October 1, 1993
,       *     Details of this algorithm are described in "Templates for the 
,       *     Solution of Linear Systems: Building Blocks for Iterative 
,       *     Methods", Barrett, Berry, Chan, Demmel, Donato, Dongarra, 
size    8143

file    CGS.f
for     *
,             SUBROUTINE CGS( N, B, X, WORK, LDW, ITER, RESID, MATVEC, PSOLVE,
,            $                INFO )
,       *
,       *  -- Iterative template routine --
,       *     Univ. of Tennessee and Oak Ridge National Laboratory
,       *     October 1, 1993
,       *     Details of this algorithm are described in "Templates for the 
,       *     Solution of Linear Systems: Building Blocks for Iterative 
,       *     Methods", Barrett, Berry, Chan, Demmel, Donato, Dongarra, 
size    6841

file    CGSREVCOM.f
for     *
,             SUBROUTINE CGSREVCOM(N, B, X, WORK, LDW, ITER, RESID, INFO,
,            $                     NDX1, NDX2, SCLR1, SCLR2, IJOB)
,       *
,       *  -- Iterative template routine --
,       *     Univ. of Tennessee and Oak Ridge National Laboratory
,       *     October 1, 1993
,       *     Details of this algorithm are described in "Templates for the 
,       *     Solution of Linear Systems: Building Blocks for Iterative 
,       *     Methods", Barrett, Berry, Chan, Demmel, Donato, Dongarra, 
size    11578

file    Cheby.f
for     *
,             SUBROUTINE CHEBY( N, B, X, WORK, LDW, ITER, RESID, MATVEC, 
,            $                  PSOLVE, INFO )
,       *
,       *  -- Iterative template routine --
,       *     Univ. of Tennessee and Oak Ridge National Laboratory
,       *     October 1, 1993
,       *     Details of this algorithm are described in "Templates for the 
,       *     Solution of Linear Systems: Building Blocks for Iterative 
,       *     Methods", Barrett, Berry, Chan, Demmel, Donato, Dongarra, 
size    6448

file    ChebyREVCOM.f
for     *
,             SUBROUTINE CHEBYREVCOM( N, B, X, WORK, LDW, ITER, RESID, INFO,
,            $                        NDX1, NDX2, SCLR1, SCLR2, IJOB)
,       *
,       *  -- Iterative template routine --
,       *     Univ. of Tennessee and Oak Ridge National Laboratory
,       *     October 1, 1993
,       *     Details of this algorithm are described in "Templates for the 
,       *     Solution of Linear Systems: Building Blocks for Iterative 
,       *     Methods", Barrett, Berry, Chan, Demmel, Donato, Dongarra, 
size    8147

file    GMRES.f
for     *
,             SUBROUTINE GMRES( N, B, X, RESTRT, WORK, LDW, WORK2, LDW2, 
,            $                  ITER, RESID, MATVEC, PSOLVE, INFO )
,       *
,       *  -- Iterative template routine --
,       *     Univ. of Tennessee and Oak Ridge National Laboratory
,       *     October 1, 1993
,       *     Details of this algorithm are described in "Templates for the
,       *     Solution of Linear Systems: Building Blocks for Iterative
,       *     Methods", Barrett, Berry, Chan, Demmel, Donato, Dongarra,
size    6401

file    GMRESREVCOM.f
for     *
,              SUBROUTINE GMRESREVCOM(N, B, X, RESTRT, WORK, LDW, WORK2,
,            $                  LDW2, ITER, RESID, INFO, NDX1, NDX2, SCLR1, 
,            $                  SCLR2, IJOB)
,       *
,       *  -- Iterative template routine --
,       *     Univ. of Tennessee and Oak Ridge National Laboratory
,       *     October 1, 1993
,       *     Details of this algorithm are described in "Templates for the
,       *     Solution of Linear Systems: Building Blocks for Iterative
size    15783

file    Jacobi.f
for     *
,             SUBROUTINE JACOBI( N, B, X, WORK, LDW, ITER, RESID, MATVEC, 
,            $                   INFO )
,       *
,       *  -- Iterative template routine --
,       *     Univ. of Tennessee and Oak Ridge National Laboratory
,       *     October 1, 1993
,       *     Details of this algorithm are described in "Templates for the 
,       *     Solution of Linear Systems: Building Blocks for Iterative 
,       *     Methods", Barrett, Berry, Chan, Demmel, Donato, Dongarra, 
size    5620

file    QMR.f
for     *
,             SUBROUTINE QMR( N, B, X, WORK, LDW, ITER, RESID, MATVEC, 
,            $                MATVECTRANS, PSOLVEQ, PSOLVETRANSQ, INFO )
,       *
,       *
,       *  -- Iterative template routine --
,       *     Univ. of Tennessee and Oak Ridge National Laboratory
,       *     October 1, 1993
,       *     Details of this algorithm are described in "Templates for the 
,       *     Solution of Linear Systems: Building Blocks for Iterative 
size    9617

file    JacobiREVCOM.f
for     *
,             SUBROUTINE JACOBIREVCOM(N, B, X, WORK, LDW, ITER, RESID,
,            $                     INFO, NDX1, NDX2, SCLR1, SCLR2, IJOB)
,       *
,       *  -- Iterative template routine --
,       *     Univ. of Tennessee and Oak Ridge National Laboratory
,       *     October 1, 1993
,       *     Details of this algorithm are described in "Templates for the 
,       *     Solution of Linear Systems: Building Blocks for Iterative 
,       *     Methods", Barrett, Berry, Chan, Demmel, Donato, Dongarra, 
size    7227

file    MatGen.f
for     *
,       *     This file contains routines for forming matrices that result from
,       *     a 5- or 7-point discretization of elliptic PDEs with Dirichlet
,       *     boundary conditions, and a consistent mass matrix "Wathen".
,       *
,       *        GEN57PT and GETSTEN are from SPARSEKIT. They actually form the
,       *        row compressed matrix.
,       *
,       *        COEFF provides the functions for computing the coefficients 
,       *        of the PDE.
size    18899

file    MatVec.f
for     *
,             SUBROUTINE MATVEC( ALPHA, X, BETA, Y )
,       *
,       *     This MatVec routine assumes the matrix is in dense format,
,       *     and uses the BLAS SGEMV.
,       *
,       *     .. Scalar Arguments ..
,             REAL               ALPHA, BETA
,       *     ..
,       *     .. Array Arguments ..
size    1601

file    PSolve.f
for     *
,       *  This file contains the preconditioner solve routines:
,       *
,       *  PSOLVE and PSOLVETRANS call the appropriate solver:
,       *
,       *     PSOLVENONE and PSOLVENONETRANS for using no preconditioning.
,       *
,       *     PSOLVEJAC and PSOLVEJACTRANS for Jacobi preconditioning.
,       *
,       *     Also included are the solvers for QMR which require left and right
size    7328

file    QMRREVCOM.f
for     *
,             SUBROUTINE QMRREVCOM(N, B, X, WORK, LDW, ITER, RESID, INFO,
,            $                     NDX1, NDX2, SCLR1, SCLR2, IJOB)
,       *
,       *
,       *  -- Iterative template routine --
,       *     Univ. of Tennessee and Oak Ridge National Laboratory
,       *     October 1, 1993
,       *     Details of this algorithm are described in "Templates for the 
,       *     Solution of Linear Systems: Building Blocks for Iterative 
size    15379

file    SOR.f
for     *
,             SUBROUTINE SOR( N, B, X, WORK, LDW, ITER, RESID, MATVEC, 
,            $                BACKSOLVE, INFO )
,       *
,       *  -- Iterative template routine --
,       *     Univ. of Tennessee and Oak Ridge National Laboratory
,       *     October 1, 1993
,       *     Details of this algorithm are described in "Templates for the 
,       *     Solution of Linear Systems: Building Blocks for Iterative 
,       *     Methods", Barrett, Berry, Chan, Demmel, Donato, Dongarra, 
size    6483

file    SORREVCOM.f
for     *
,             SUBROUTINE SORREVCOM(N, B, X, WORK, LDW, ITER, RESID, INFO,
,            $                     NDX1, NDX2, SCLR1, SCLR2, IJOB)
,       *
,       *  -- Iterative template routine --
,       *     Univ. of Tennessee and Oak Ridge National Laboratory
,       *     October 1, 1993
,       *     Details of this algorithm are described in "Templates for the 
,       *     Solution of Linear Systems: Building Blocks for Iterative 
,       *     Methods", Barrett, Berry, Chan, Demmel, Donato, Dongarra, 
size    8609

file    STOPTEST2.f
for     *
,             SUBROUTINE STOPTEST2( N, R, B, BNRM2, RESID, TOL, INFO )
,       *
,       *     .. Scalar Arguments ..
,             INTEGER            N, INFO
,             REAL               RESID, TOL, BNRM2
,       *     ..
,       *     .. Array Arguments ..
,             REAL               R( * ), B( * )
,       *     ..
size    1840

file    Tester.f
for     *
,             PROGRAM TEMPLATESTESTER
,       *
,       *  Test program for the SINGLE PRECISION iterative templates.
,       *
,       *  The program must be driven by a short data file. The first 18 records
,       *  of the file are read using list-directed input, the last 16 records
,       *  are read using the format ( A6, L2 ). An annotated example of a data
,       *  file s a follows:
,       *
size    25480

file    Utils.f
for     *
,       *  This file contains routines used by Jacobi, SOR, and Chebyshev:
,       *
,       *  Jacobi/SOR:
,       *
,       *     MATSPLIT  calls specific matrix splitting routine
,       *     JACSPLIT
,       *     SORSPLIT
,       *     BACKSOLVE
,       *
size    7171

file    sblas.f
for     *
,       *  This file contains the single precision fortran BLAS routines needed 
,       *  by the Templates. Therefore, if you can link to a BLAS library, this
,       *  file is not needed.
,       *
,             subroutine saxpy(n,sa,sx,incx,sy,incy)
,       c
,       c     constant times a vector plus a vector.
,       c     uses unrolled loop for increments equal to one.
,       c     jack dongarra, linpack, 3/11/78.
size    102471

file    slapack.f
for     *
,       *  The routines in this file are from LAPACK. SSYEV is used by the
,       *  Chebyshev method for computing the extremal eigenvalues of the
,       *  matrix. All the other routines in this file support this routine.
,       *  If you can link to the LAPACK library, this file is not needed.
,       *
,             SUBROUTINE SSYEV( JOBZ, UPLO, N, A, LDA, W, WORK, LWORK, INFO )
,       *
,       *  -- LAPACK driver routine (version 1.1) --
,       *     Univ. of Tennessee, Univ. of California Berkeley, NAG Ltd.,
size    184522

file    snsychk.f
for     *
,             SUBROUTINE SNSYCHK( X, LDX, B, X0, WORK, LDW, PFORM, MATVEC,
,            $                    MATVECTRANS, PSOLVE, PSOLVETRANS, PSOLVEQ,
,            $                    PSOLVETRANSQ, BACKSOLVE, TOL, SCALEDTOL,
,            $                    LTEST, NSYRES, NUMTESTS, NUMSUSP, CRITERR )
,       *
,       *     .. Scalar Arguments ..
,             INTEGER            LDW, LDX, NUMTESTS, NUMSUSP, CRITERR
,             REAL               TOL, SCALEDTOL
,             LOGICAL            NSYRES
size    8716

file    sspdchk.f
for     *
,             SUBROUTINE SSPDCHK( X, LDX, B, X0, WORK, LDW, PFORM, MATVEC, 
,            $                    MATVECTRANS, PSOLVE, PSOLVETRANS, PSOLVEQ,
,            $                    PSOLVETRANSQ, BACKSOLVE, TOL, SCALEDTOL,
,            $                    LTEST, SPDRES, NUMTESTS, NUMSUSP, CRITERR )
,       *
,       *     .. Scalar Arguments ..
,             INTEGER            LDW, LDX, NUMTESTS, NUMSUSP, CRITERR
,             REAL               TOL, SCALEDTOL
,             LOGICAL            SPDRES
size    10728

file    sutils.f
for     *  This file contains various routines needed by the tester for generating
,       *  the test matrices. They are all copied from the LAPACK test suite,
,       *  so if you can link to this tester, this file is not needed.
,       *
,             REAL FUNCTION SLARAN( ISEED )
,       *
,       *  -- LAPACK auxiliary routine (version 1.1) --
,       *     Univ. of Tennessee, Univ. of California Berkeley, NAG Ltd.,
,       *     Courant Institute, Argonne National Lab, and Rice University
,       *     February 29, 1992
size    2849

file    test.results
for      
,        DETAILS OF ITERATIVE TEMPLATES TEST:
,        
,           Univ. of Tennessee and Oak Ridge National Laboratory
,           October 1, 1993
,           Details of these algorithms are described in "Templates
,           for the Solution of Linear Systems: Building Blocks for
,           Iterative Methods", Barrett, Berry, Chan, Demmel, Donato,
,           Dongarra, Eijkhout, Pozo, Romine, and van der Vorst,
,           SIAM Publications, 1993.
size    7206