TOMS Algorithms 582, (508, 509)
Gibbs-Poole-Stockmeyer (GPS) bandwidth reduction algorithm, and Gibbs-King profile reduction algorithm.
The author of the GPS code
advises that there is a bug in TOMS 582 and suggests the following correction:
"In GPSKCE, modify the DO 500 loop by replacing the lower limit, now 1, by zero."
Roger Young and Ian MacPhedran