diff options
Diffstat (limited to 'math/Revisions')
| -rw-r--r-- | math/Revisions | 406 |
1 files changed, 406 insertions, 0 deletions
diff --git a/math/Revisions b/math/Revisions new file mode 100644 index 00000000..44110034 --- /dev/null +++ b/math/Revisions @@ -0,0 +1,406 @@ +.help revisions Jun88 math +.nf + +gsurfit/gs_deval.gx +iminterp/mrider.x +iminterp/mrieval.x + Fixed some procedure calls being closed with a ']' insted of + a ')' (2/17/08, MJF) + +gsurfit/gs_chomat.gx + The test for singularity would fail with certain kinds of problems + because the test used EPSILON (should have been EPSILOND for the + double precision) but this is for distinguishing numbers small + numbers from 1 and not from each other. The test is now done + with a comparison against the smallest real or double difference. + The place where this was found to be a problem was with CCSETWCS. + (8/31/04, Valdes) + +======== +V2.12.2a +======== + +curfit/cvinit.gx + If one of the error checks caused an error return the cv pointer + would have been allocated (cv != NULL) but some of the pointer + fields could have garbage values since a malloc was used instead + of a calloc. A later call to cvfree could result in a segmentation + error. This was changed so that 1) a null cv pointer is returned + in the initial error checks cause an error return and 2) the + cv pointer is initially allocated with calloc so that no pointer + fields will be non-NULL until explicitly set. + (11/18/02, Valdes) + +======= +V2.12.1 +======= + +gsurfit/gsder.gx +gsurfit/gsderr.x +gsurfit/gsderd.x +gsurfit/gs_fder.gx +gsurfit/gs_fderr.x +gsurfit/gs_fderd.x +gsurfit/gs_deval.gx +gsurfit/gs_devalr.x +gsurfit/gs_devald.x + Rewrote the routines which compute the derivatives of the polynomial + surfaces. There were some renormalization problems in this code + for non-linear chebyshev and legendre polynomial surfaces. The new + code uses recursion relations for computation the derivatives which + could cause small errors if they are used to do flux conservation. + (27/10/99 LED) + +curfit/mkpkg +deboor/mkpkg +iminterp/mkpkg + Added some missing file dependencies to the above mkpkg files. + (20/9/99 LED) + +slalib + Upgraded slalib from Version 1.6.3 to Version 2.3.0. Seventeen + new routines were added including 4 dealing with FK5 to ICRS + coordinate conversions (17/6/99 LED). + +arbpix.x +arider.x +arieval.x +asider.x +asieval.x +asifit.x +asifree.x +asigeti.x +asigetr.x +asigrl.x +asiinit.x +asirestore.x +asisave.x +asisinit.x +asitype.x +asivector.x +ii_1dinteg.x +ii_bieval.x +ii_cubspl.f +ii_eval.x +ii_greval.x +ii_pc1deval.x +ii_pc2deval.x +ii_polterp.x +ii_sinctable.x +ii_spline.x +ii_spline2d.x +im1interpdef.h +im2interpdef.h +mkpkg +mrider.x +mrieval.x +msider.x +msieval.x +msifit.x +msifree.x +msigeti.x +msigetr.x +msigrid.x +msigrl.x +msiinit.x +msirestore.x +msisave.x +msisinit.x +msisqgrl.x +msitype.x +msivector.x + The 1D interpolation routines were replaced with new versions which + support look-up table since and drizzle interpolation. Minor + modifications were also made to the existing 1D sinc routines. The 2D + image interpolation routines were replaced with new versions modified + which support sinc, look-up table sinc and drizzle interpolation. + (12/28/98 Davis) + +iminterp.hd +iminterp.men +iminterp.spc +arbpix.hlp +arider.hlp +arieval.hlp +asider.hlp +asieval.hlp +asifit.hlp +asifree.hlp +asigeti.hlp +asigetr.hlp +asigrl.hlp +asiinit.hlp +asirestore.hlp +asisave.hlp +asisinit.hlp +asitype.hlp +asivector.hlp +iminterp.hlp +mrider.hlp +mrieval.hlp +msider.hlp +msieval.hlp +msifit.hlp +msifree.hlp +msigeti.hlp +msigetr.hlp +msigrid.hlp +msigrl.hlp +msiinit.hlp +msirestore.hlp +msisave.hlp +msisinit.hlp +msisqgrl.hlp +msitype.hlp +msivector.hlp + The 1D and 2D image interpolation routines help pages were updated + to reflect the changes in both packages. (12/28/98 Davis) + +math$nlfit/nlchomat.gx +math$nlfit/nlchomatr.x +math$nlfit/nlchomatd.x + Modified the singular matrix test to make a comparison against + EPSILON instead of 0.0 to avoid floating point problems. + + (8/1/98 Davis) + +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 |