.help revisions Jun88 math
.nf
math$gsurfit/
	Installed a completely new version of gsurfit which support "half"
	or "diagonal" cross terms. All the .gx and .x files were replaced.

	(4/30/97 Davis)

math$gsurfit/gsrestore.gx
math$gsurfit/gsrestorer.x
math$gsurfit/gsrestored.x
	Changed the type declaration of the xmin, xmax, ymin, ymax variables
	from real to PIXEL to avoid machine precision problems.

	(3/21/96 Davis)

math$nlfit/nlfit.gx
math$nlfit/nlfitr.x
math$nlfit/nlfitd.x
	Fixed a divide by zero error in the nlscatter routine which occurs
	if the fit has no degrees of freedom and the weighting type is
	WTS_SCATTER.

	(6/13/94 Davis)

math/gsurfit/gssub.gx
math/gsurfit/gssubd.x
	The gsurfit double precision routine gssubd was incorrectly calling
	the real precision routine gsgeti instead of the double precision
	routine dgsgeti due to an error in the gssub.gx file. This bug
	results in the wrong value of the parameter GSNSAVE being returned to
	the calling routine gssubd in the double precision case.

	(1/17/94 Davis)

math/iminterp/arbpix.x
math/interp/arbpix.x
        Changed all if (x == INDEF) constructs to IS_INDEF(x) constructs.

	(1/15/93 Davis)

math/surfit/iseval.x
	Two indices were not being initialized before being used in a pointer
	offset computation in the spline1 and spline3 single input point code
	producing incorrect results.

	(12/17/92 Davis)

math/nlfit/nlfit.gx
math/nlfit/nlfitr.x
math/nlfit/nlfitd.x
math/nlfit/nliter.gx
math/nlfit/nliterr.x
math/nlfit/nliterd.x
	Added a test for lambda<=0.0 in the nlfit loop and in the iteration
	loop. Precision problems could cause the code to go into an infinite
	loop in some circumstances. The code returns a NOT_DONE status
	message if the convergence criteria have not been met in this
	circumstance. Made some other minor changes to the convergence
	criteria.

	(7/20/92 Davis)

math/nlfit/nlacpts.gx
math/nlfit/nlerrors.gx
	Fixed a bug in the routine nlacpts where the variable sum was
	not being initialized properly for each iteration. This could
	cause a slow drift in repeated fits. Also fixed a related problem
	in the error computation. Note this was not a problem
	with the version of nlfit in apphot in versions 2.9 or in 
	testphot.

	(1/9/92 Davis)

math/nlfit/nlfit.gx
math/nlfit/nlfitr.x
math/nlfit/nlfitd.x
	Added some code to force the nlfit package to do at least MINITER
	iterations whether or not the solution actually converged; MINITER
	+ 1 if the weighting type is WTS_SCATTER.
	(1/3/91 Davis)

math/nlfit/nliter.gx
math/nlfit/nliterr.x
math/nlfit/nliterd.x
	Simplied the test for determining whether the lambda factor needs
	to be increased or decreased.
	(1/3/91 Davis)

math/nlfit/nlzero.gx
math/nlfit/nlfitdefr.h
math/nlfit/nlfitdefd.h
	Removed an extra argument from two aclr$t routines. Also these routines
	were being called with the pointer instead of the array as the 
	argument. (10/3/91 Davis)

lib/math/iminterp.h
math/iminterp/arbpix.x
math/iminterp/arider.x
math/iminterp/arieval.x
math/iminterp/asider.x
math/iminterp/asieval.x
math/iminterp/asifit.x
math/iminterp/asiffree.x
math/iminterp/asigrl.x
math/iminterp/asiinit.x
math/iminterp/asirestore.x
math/iminterp/asisave.x
math/iminterp/asivector.x
math/iminterp/ii_1dinteg.x
math/iminterp/ii_eval.x
math/iminterp/ii_pc1deval.x
math/iminterp/im1interpdef.h
math/iminterp/doc/arbpix.hlp
math/iminterp/doc/arider.hlp
math/iminterp/doc/arieval.hlp
math/iminterp/doc/asider.hlp
math/iminterp/doc/asieval.hlp
math/iminterp/doc/asifit.hlp
math/iminterp/doc/asigrl.hlp
math/iminterp/doc/asiinit.hlp
math/iminterp/doc/asivector.hlp
	The 1D image interpolation routines were replaced with new versions
	modified by Frank Valdes to support sinc interpolation. (7/26/91 Davis)

lib/math/curfit.h
lib/math/gsurfit.h
lib/math/iminterp.h
lib/math/interp.h
lib/math/surfit.h
	Added dictionary string definitions for the interpolator types and
	weights.  Use of these strings insures that strdic/clgwrd will
	return the correct type code.  Currently one has to assume the code
	definitions will not change or put a lookup table with a data
	statement containing the macro definitions.  The code would also
	have to be modified to add new types.  The dates were restored to
	their original values to avoid recompilation.  (5/9, Valdes)

math$curfit/cverrors.gx: Davis, May 6, 1991
math$curfit/cvpower.gx
math$curfit/cvrefit.gx
math$curfit/cvrestore.gx
    Finished cleanup of curfit package .gx files.

math$gsurfit/gs_f1deval.gx: Davis, May 6, 1991
    Changed the line "call amulk$t (,,,2.,,,)" to "call amulk$t (,,,2$f,,,)
    to remove a type dependence problem found by ftoc on the Mac.

math$curfit/cvpower.gx: Davis, May 6, 1991
    Changed the constant from INDEFR to INDEF in the amov$t call in cvpower.gx.
    This was causing a problem for the Mac compiler.

math$curfit/cv_b1eval.gx: Davis, April 23, 1991
math$curfit/cv_beval.gx
math$curfit/cv_feval.gx
math$curfit/cvaccum.gx
math$curfit/cvacpts.gx
math$curfit/cvchomat.gx
math$curfit/cvfree.gx
math$curfit/cvinit.gx
    Did some cleaning up in the following .gx files to make the code easier
    to read.

math$nlfit/: Davis, Jan 31, 1991
    Added the nlfit package to the math package. 

math$surfit/isreplace.x: Davis, September 18, 1990
    1. Changed the int calls in isreplace.x to nint calls. This is a totally
    safe way to do the conversion from floating point to integer
    quantities in the surfit package and removes any potential precision
    problems for task which must read the surfit structure back from a
    text database file.

math$gsurfit/gsrestore.gx,gsrestorer.x,gsrestored.x: Davis, September 18, 1990
    1. Changed the int calls in gsrestore.gx to nint calls. This is a totally
    safe way to do the conversion from floating point to integer
    quantities in the gsurfit package and removes any potential precision
    problems for task which must read the gsurfit structure back from a
    text database file.

math$curfit/cvrestore.gx,cvrestorer.x,cvrestored.x: Davis, September 18, 1990
    1. Changed the int calls in cvrestore.gx to nint calls. This is a totally
    safe way to do the conversion from floating point to integer
    quantities in the curfit package and removes any potential precision
    problems for task which must read the curfit structure back from a
    text database file.

math$gsurfit/gsder.gx,gsderd.x: Davis, August 1, 1990
    1. Fixed a typo in the routine which computes the double precision
    derivatives of a surface. Due to this type the salloc routine was being
    passed a pointer instead of size. When the value of this pointer
    got very large this could cause an out of memory allocation error.

math$iminterp/msigrl.x: Davis, Feb 8, 1989
    1. There was a bug in the shift and wrap routines in the 2D routine
    which integrates in polygonal apertures. If the vertices of the
    polygon  were listed in certain orders the routine would produce
    wrong results. The routine has been totally rewritten to use
    the polygon clipping code in polyphot.

math$interp/asigrl.x: Davis, Jan 27, 1988
    1. The routine ii_getpcoeff was duplicated in the files asigrl.x and
    ii_1dinteg.x. I removed the copy in asigrl.x and rebuilt the library.

math$interp/arider.x: Davis, February 10, 1988
    1. The routines iidr_poly3, iidr_poly5 and iidr_spine3 were declared
    as functions but called everywhere as subroutines. I removed the
    function declarations.

math$gsurfit/gsder.gx: Davis, October 16, 1987
    1. changed a real type conversion call to double in the gsder routine.
    This error was causing the derivative evaluations to fail on the sun.

math$gsurfit/gsder.gx: Davis, August 11, 1987
    1. A call in gsder to gscoeff was not being set to double precision.

math$gsurfit: Davis, January 22, 1986
    1. A new version of the gsurfit package was installed. A double precision
    gsurfit has  been added. The double precision calls are the same as
    the single precision calls but prefixed by a d. Some of the evaluation
    routines were vectorized and run significantly faster than the old version.
    A power series polynomial function was added to the package for convenience
    in evaluating the simple rotation transformation functions.

math$iminterp/arider.x: Davis, September 7, 1986
    1. The routines ii_pcpoly3, ii_pcpoly5 and ii_pcspline3 were declared
    as functions but called everywhere as subroutines. Removed the function
    declaration. I also moved these routines into a file of their own
    ii_pc1deval.x and changed the prefix to ia_ from ii_ to avoid conflict
    with the 2D equivalents.

math$interp/arider.x: Davis, September 7, 1986
    1. The routines ii_pcpoly3, ii_pcpoly5 and ii_pcspline3 were declared
    as functions but called everywhere as subroutines. Removed the function
    declaration.

From Davis Jan 8, 1986

1. Added documentation for surfit package.

----------------------------------------------------------------------------
.endhelp
