The first result below is correct; the second is close but has fewer significant digits. The strings have 211 characters. I can't figure out what my precision specification has to do with the output.
It's the mapm m_apm_pow() function, specifically the second ("decimal_places") argument. If it's set to fewer characters than the result, you'll get rounding errors. (TCC is already adding 20 digits on its own, which is why you didn't see it until you got a > 30 digit result.)
Maybe it's a bug, or a lack of documentation, (or me) but when computing base**exponent, (MAPM::pow() or m_apm_pow()) and base>1, "decimal_places" actually specifies the number of significant digits (which certainly includes what's to the left of the decimal point) and not the number of places after the decimal point.
If you add the number below (instead of 20) to the desired number of decimal places (when base>1) you should be safe. If you wish, I'll translate it into the non-C++ version (i.e., using m_apm_func() instead of MAPM::func()).
Code:
// this is before the decimal point; add the desired number of decimal places
INT extra (MAPM base, MAPM exponent)
{
MAPM u, mm1 = MM_One;
INT x;
u = ++((base.log10() * exponent).ceil());
CHAR *szdp = (CHAR*) AllocMem(u.significant_digits()+1);
u.toIntegerString(szdp);
x = atoi(szdp);
FreeMem(szdp);
return x;
}
I reproduced (and fixed, as above) the problem with 2**103. But I don't need the fix for .5**(-103), which is the same number. @EVAL, however, has the same problem with both.
FWIW, extra() is a pretty tight estimate of what's needed. Below, extra is 46 when the number has 45 digits before the decimal point. I think there are special cases when extra()-1 won't be enough.
Code:
v:\> echo %@pow[1.01 10310.99]^r^n%_ruler
extra = 46
361111948879059149982803812474402661922976059.62768572890231244276
....+....1....+....2....+....3....+....4....+....5....+....6....+....7....+....8....+....9....+....0....+....1...
FWIW, this (below) is more than 20 times faster than extra() and should work in reasonable cases (intergalactic distances and the national debt being possible exceptions).
Code:
double b = ... , // base
e = ... , // exponent
extra = 1 + ceil(e * abs(log10(b)));
This site uses cookies to help personalise content, tailor your experience and to keep you logged in if you register.
By continuing to use this site, you are consenting to our use of cookies.