[comp.lang.fortran] Routines for output of matrices

bernhold@red8.qtp.ufl.edu (David E. Bernholdt) (04/25/91)

Some time ago, while writing a code which has occasion to output many
matrices in many different formats. I was tired of working with codes
where the format, line width, etc. are hardwired all over the program.
I decided to write some fairly general routines to handle the task.

Input to these routines includes matrix to be output (of course), the
unit to write on, a title to put above the matrix, the format
descriptor to be used for the elements of the array, and the width of
the output.  The routine parses the format descriptor, computes how
many values will fit on each line, and writes out the matrix in chunks
as wide as will fit in the allowed line width.

The naming and the call sequence for these routines is based on the
BLAS conventions, and each routine is 100-200 lines long (50%
comments, though).  The set is not complete -- I've only written the
ones I needed, but it does handle single and double precision matrices
and a few other things (no complex yet).

If anyone is interested, I can mail them.  I'd be interested to hear
comments on the routines and ideas for better ways to do things.
-- 
David Bernholdt			bernhold@qtp.ufl.edu
Quantum Theory Project		bernhold@ufpine.bitnet
University of Florida
Gainesville, FL  32611		904/392 6365