help-glpk
[Top][All Lists]
Advanced

[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

Re: [Help-glpk] Working with larger numbers


From: Andrew Makhorin
Subject: Re: [Help-glpk] Working with larger numbers
Date: Wed, 30 Jul 2008 13:48:47 +0400

To fix the defect causing round-off errors please insert the following
fragment in the routine set_d_eps (file glplpx10.c, lines 63-71):

63    int s, n, j;
64    double f, p, q, eps = 1e-9;
65    mpq_t temp;
66    xassert(-DBL_MAX <= val && val <= +DBL_MAX);
+     if (val == floor(val))
+     {  /* if val is integral, do not approximate */
+        mpq_set_d(x, val);
+        goto done;
+     }
67    if (val > 0.0)
68       s = +1;
69    else if (val < 0.0)
70       s = -1;
71    else

and then recompile the package as usual.

Note that these changes will appear in a next release of the package.





reply via email to

[Prev in Thread] Current Thread [Next in Thread]