prehistoric, via dongarra From wmc Mon Aug 29 07:33:09 EDT 1988 Kincaid said that he would be happy to have netlib distribute their iterative codes. From cs.utexas.edu!oppe Tue Jan 16 16:50:15 CST 1990 Received: from cactus.cs.utexas.edu by cs.utexas.edu (5.59/1.46) id AA20731; Tue, 16 Jan 90 16:50:23 CST Posted-Date: Tue, 16 Jan 90 16:50:15 CST Message-Id: <9001162250.AA01064@cactus.cs.utexas.edu> Received: by cactus.cs.utexas.edu (5.59/1.4-Client) id AA01064; Tue, 16 Jan 90 16:50:16 CST From: oppe@cs.utexas.edu (Tom Oppe) Date: Tue, 16 Jan 90 16:50:15 CST X-Mailer: Mail User's Shell (6.5.6 6/30/89) To: ehg@research.att.com Dear Dr. Grosse, Would you do two favors for me? (1) Would you replace the "index" file of the itpack directory on netlib with the following file? To be consistent with the new index, please rename three of the files in the itpack directory as follows: Old Name New Name -------- -------- src2cv.f srcv2d.f tst2cv.f tstv2d.f user2cv.tex userv2d.tex ================new index file for itpack=================================== The ITPACK directory contains the following files: dsrc2c1.f - d.p. ITPACK 2C source code, part 1 dsrc2c2.f - d.p. ITPACK 2C source code, part 2 dtst2c.f - testing routine for d.p. ITPACK 2C index - this file that you are reading info.tex - information about CNA iterative packages (formated with LaTeX) nspcg1.f - NSPCG source code, part 1 nspcg2.f - NSPCG source code, part 2 nspcg3.f - NSPCG source code, part 3 nspcg4.f - NSPCG source code, part 4 nspcg5.f - NSPCG source code, part 5 src2c1.f - s.p. ITPACK 2C source code, part 1 src2c2.f - s.p. ITPACK 2C source code, part 2 srcv2d.f - ITPACKV 2D source code (vectorized for Cray X-MP) tst2c.f - testing routine for s.p. ITPACK 2C tstv2d.f - testing routine for ITPACKV 2D tstnsp1.f - testing program 1 for NSPCG tstnsp2.f - testing program 2 for NSPCG tstnsp3.f - testing program 3 for NSPCG tstnsp4.f - testing program 4 for NSPCG tstnsp5.f - testing program 5 for NSPCG user2c.tex - ITPACK 2C User's Guide (formated with LaTeX) userv2d.tex - ITPACKV 2D User's Guide (formated with LaTeX) usernsp.tex - NSPCG User's Guide (formated with LaTeX) There are four packages in the ITPACK directory for solving large sparse linear systems by iterative methods: ITPACK 2C (single precision), ITPACK 2C (double precision), ITPACKV 2D (a vectorized version of ITPACK 2C for the Cray X-MP and similar vector computers), and NSPCG. ITPACK 2C and ITPACKV 2D are intended for symmetric and positive-definite matrix problems. NSPCG has preconditioners and polynomial accelerators for nonsymmetric matrix problems as well. Only single precision versions are available for ITPACKV 2D and NSPCG. Machine-dependent constants and the timing function may need to be modified by the user when installing the packages for a particular computer. They are located in routines DFAULT and TIMER in all four packages. Currently, the codes assume a 32-bit word length and the UNIX etime timing facility. Comments and suggestions are welcome and should be addressed to: Center for Numerical Analysis RLM Bldg. 13.150 University of Texas at Austin Austin, Texas 78713-8510 Ph: (512) 471-1242 email: oppe@cs.utexas.edu or kincaid@cs.utexas.edu ================end of new index file for itpack=============================== (2) The routine SUBQ25 in the file nspcg1.f of the itpack directory needs to be replaced with the following version. Apparently, the current version is missing a few lines. =============beginning of subq25 file========================================== subroutine subq25 (coef,jcoef,wksp,iwksp,nn,r,z) c c ... subq25 calls ibsl, for bic preconditioning. c c common / dscons / ndim, mdim, maxnz common / itcom6 / method, iscale, iperm, nstore, a ifact, kblsz, lvfill, ltrunc, ndeg, a ipropa, isymm, ifctv logical propa common / cblock / propa, ncolor, maxd, nc, ipt, maxnew, a jcnew, lbhb, iblock, ncmax common / intern / ndt, ndb, maxt, maxb, ivers, irwise common / cfactr / nfactr, nfacti, ifactr, ifacti, timfac common / point / iptscl, iwkpt1, iwkpt2, iwkpt3 integer jcoef(2), iwksp(1) dimension r(1), z(1), coef(1), wksp(1) c n = nn ipt2 = ifactr + n*iwksp(ifacti+2) if (lvfill .gt. 0) go to 10 nwdiag = iwksp(ifacti+2) - ltrunc if (propa) call ibsl a (n,ndim,n,kblsz,kblsz,lbhb,iwksp(ifacti), a wksp(ifactr),coef(ndim*nwdiag+1), a jcoef(nwdiag+1),r,z,ivers,wksp(iwkpt1)) if (.not. propa) call ibsl a (n,n,n,kblsz,kblsz,lbhb,iwksp(ifacti), a wksp(ifactr),wksp(ipt2), a jcoef(nwdiag+1),r,z,ivers,wksp(iwkpt1)) return 10 ipt1 = ifacti + 3*lbhb + iwksp(ifacti+2) call ibsl (n,n,n,kblsz,kblsz,lbhb,iwksp(ifacti), a wksp(ifactr),wksp(ipt2),iwksp(ipt1),r,z, a ivers,wksp(iwkpt1)) return end =================end of subq25 file===================================== I believe that these changes need to be made at both the ornl.gov and research@att.com locations of netlib. Thank you very much for any help you can provide. Thomas C. Oppe (oppe@cs.utexas.edu) .