For questions send mail to: fkrogh@mathalacarte.com
A new option has been added to DIVA/SIVA to allow the user to indicate that DIVADB/SIVADB is being called, so that variables that might be printed prior to being initialized will get initialized in this case.
DRSPOLZ has been changed to use a complex type for the roots, and new capability has been added to M77CON to simplify doing this with our conversion tools.
XERBLA now prints the value of INFO.
In addition, the following drivers were changed: DRDBESJ0/DRSBESJ0, DRDDASL2/DRSDASL2, DRDDASL3/DRSDASL3, DRDDASL6/DRSDASL6, DRDJACG2/DRSJACG2, DRCWOFC, DRDPOLZ/DRSPOLZ/DRZPOLZ, DRDPOLZ2/DRSPOLZ2, DRDIVDB/DRSIVDB, and DRDIVX/DRSIVX.
DINTA/SINTA (2009-10-17) Set AACUM to ACUM when nodes are coalescing and a few other changes.
DPLOT/SPLOT (2009-10-18) Added "save ADJIN" in dplota. Added a facility to read sets of points from a file.
DXRK8/SXRK8 (2009-15-09) Changed starting stepsize, many changes in the way the stepsize is selected, and a new method of checking for stiffness. There have also been some changes in usage and some bug fixes. Also some fixes in DXRK8G/SXRK8G.
DINTO/SINTO (2009-07-15) Fixes in some of the messages output.
DMESS/SMESS (2009-07-23) Many many minor fixes in outputs.
MESS (2009-02-28) Fixes for NAG compiler and fixed "f" format for C code, definitions.
DPRTSV/SPRTSV (2008-11-26) Fixed format to be Fortran standard.
DDASLX/SDASLX (2008-11-13) Fixed comments on starting, corrected leading dimension for banded problem, allowed INFO(11) = 0, 1, 2 for new startup, and added argument of leading dimension to ddasf/sdasf.
CSORTP (2008-10-29) Moved the string.h line for the C code up.
XERBLA (2008-10-27) Changed to use the MATH77/mathc90 error processor.
DDASLS/SDASLS (2008-10-26) Big changes in how we get starting values in solving differential algebraic equations. (The earlier code has some serious problems.)
DDASTP/SDASTP (2008-10-24) Declared dnrm2, added argument of leading dimension to ddasf/sdasf.
DDASF/SDASF (2008-10-24) Added row dimension argument, shortened the first line to 72 characters.
DDAS1/SDAS1, DDASJ/SDASJ, DDASLV/SDASLV (2008-08-26) Added argument for leading dimension to ddasf/sdasf.
DZERO/SZERO (2008-03-01) Minor changes in diagnostic print.
DINTMA/SINTMA (2008-01-11) Added a new error message.
DINTDU/SINTDU (2007-03-28) Bug fixed.
DTGPD/STGPD (2007-02-28) Big changes to eliminate assigned go to's.
DRPRL3/SRPRL3 (2007-02-28) Fixed a rounding problem in the labels.
DAXPY/SAXPY, DCOPY/SCOPY, DDOT,SDOT, DROTM/SROTM, DROTMG/SROTMG, DSWAP/SSWAP (2006-06-07) Removed arithmetic if's.
DRDDASL2/DRSDASL2 (2009-10-19) Fixed so 0 value of d(4,4) does not lead to divide by 0. Code changed to use new starting procedure for y', row dimension added to evaluators, minor changes for C conversion.
DRDJACG2/DRSJACG2 (2008-11-01) Added total energy constraint to the third example, and added row dimension parameter to evaluators.
DRDDASL4/DRSDASL4 (2008-10-29) New starting computation for y', added row dimension to evaluator call, multiple changes for C conversion.
DRDDASL7/DRSDASL7 (2008-10-26) New starting computation for y', changes for conversion to C.
DRDDASL5/DRSDASL5 (2008-10-26) Added total energy as a constraint, new code for getting initial y', a few changes for conversion to C.
DRDDASL3/DRSDASL3, DRDDASL1/DRSDASL1 (2008-10-26) Added row dimension for evaluators, new starting computation of y', minor changes for conversion to C.
DRDXRK8/DRSXRK8 (2008-02-24) Changed dimensions for new usage.
DRDCKDER/DRSCKDER (2007-01-02) Put commas around ':' in formats.
Bugs in the interpolation codes (D/S)ILUP have been corrected. They were flagging extrapolation in the case when interpolation was being done at the end point of a table when using indexed look ups. This led to getting incorrect results in the multi-dimensional codes (D/S)ILUPM when using indexed look ups. Also the multi-dimensional codes were getting slightly better accuracy when error estimates were requested. They have been changed so that results don't depend on whether error estimates are requested.
Significant changes to the DAE solvers, (D/S)DASLX, have been made, especially in the code for projecting on to constraints.
Code has been added, (D/S)ILUPMD, to print out the structure of tables used in multi-dimensional interpolation.
A sparse matrix solver using Gaussian Elimination, (D/S)SPGE, has been added. In connection with this the message routines (MESS,(D?S)MESS) have had code added to support the output of sparse vectors and indirectly the output of sparse matrices.
There has also been a change to M77CON to make it more portable and to handle the type conversions for code that is in the comments for the conversion to C. And some bugs in GLOCH1 (used mainly for the conversion to C) were corrected.
Another major change has been in the ddasslx/sddaslx codes for differential algebraic systems. A new method of stepsize selection due to Gustaf Söderlind has been incorporated into the code. The code does a much better job of getting accuracy that is smoothly dependent on the error tolerance and is slightly more efficient when getting a given accuracy. These changes required some changes in the logic for checking convergence of the iterations for solving the corrector equation. An option of getting numerical derivatives has also been added.
Other miscellaneous changes include,
Similarly for the drivers, the following codes were changed (with a "drs" file for every "drd" file listed): drcgam, drcgefsc, drdasinh, drdblas3, drdcdchi, drdcdnml, drdcdpoi, drdcon, drdcpval, drdfrenl, drdgami, drdge2, drdgefsc, drdhfti, drdilup, drdilupm, drdintmf, drdintmr, drdiva, drdivdb, drdivx, drdmlc01, drdnlafb, drdnlagb, drdnlsfb, drdnlsfu, drdnlsgb, drdnlsgu, drdpoi, drdrane, drdrang, drdrangv, drdranr, drdtgfi, drdtgfi1, drdtgfi2, drdtgrec, drdxrk8, drmatp, drmatpr, drvecp, drvecpr, drzgam. All such changes were minor.