More problems with SetPrecision[] and/or $MinPrecision,...

*To*: mathgroup at smc.vnet.net*Subject*: [mg56954] More problems with SetPrecision[] and/or $MinPrecision,...*From*: terryisnow at yahoo.com*Date*: Wed, 11 May 2005 05:24:09 -0400 (EDT)*Sender*: owner-wri-mathgroup at wolfram.com

Dear Mathematica Guru's, Lately I seem to be finding a lot of SetPrecision[]/$MinPrecision bugs. This appears to be another one: In[1]:= $MinPrecision = 20; savedData1 = 669151541.9328941107875; savedData2 = 0.99960472897267246018765; savedData3 = Log[savedData1^savedData2/savedData1^savedData2]; Print[Apply[Plus, savedData3 - Log[savedData1^savedData2/savedData1^savedData2]] ]; x = 0.99960472897267245645; y = 669151541.932894110712; Do[ z = SetPrecision[y,24]^x; ,{50000}]; Print[Apply[Plus, savedData3 - Log[savedData1^savedData2/savedData1^savedData2]] ]; The output looks like this on my Mathematica 3.0 system: 0 1.9586707534418188626 x 10^-30 Obviously the question is why is the second Print statement not also returning a 0, after all the code that runs between the two Print statements neither changes any of the variables referenced by the two Print statements, nor even references any of the same variables. The problem maybe reproducible without the $MinPrecision but with this specific test case it is needed (there might however be other constants that can reproduce the problem without $MinPrecision). Can someone also run this on a more recent version of Mathematica to see if it reproduces? All help is very much appreciated! Terry

**Follow-Ups**:**Re: More problems with SetPrecision[] and/or $MinPrecision,...***From:*Chris Chiasson <chris.chiasson@gmail.com>

**Re: More problems with SetPrecision[] and/or $MinPrecision,...***From:*Chris Chiasson <chris.chiasson@gmail.com>