Why does N not upgrade precision? [duplicate]

Why does N not upgrade precision? [duplicate]



This question already has an answer here:


Precision[N[1.0, 20]]
Precision[N[1, 20]]


MachinePrecision
20.



It would be so much more intuitive and less error prone, if Precision[N[1.0, 20]] would be 20 and not MachinePrecision. Why do I have to explicitly use N[Rationalized[1.0], 20] to upgrade the precision?


Precision[N[1.0, 20]]


20


MachinePrecision


N[Rationalized[1.0], 20]



If it is about the warning messages during the calculations, Mathematica could return a warning already at the Precision[N[1.0, 20]] step.


Precision[N[1.0, 20]]



Edit
I am not attempting to use N for rounding. I also understand the difference between MachinePrecision and arbitrary precision, which is very well described in this answer.


N


MachinePrecision



I want to know, why does Mathematica not consider the description of the number existing in the underling binary representation as exact if I apply N, and why I need to wrap it with Rationalize.


N


Rationalize



Is it performance?



Is it some deep semantic meaning of N vs SetPrecision?


N


SetPrecision



Is set SetPrecision any different from N@Rationalize@?


SetPrecision


N@Rationalize@



This question has been asked before and already has an answer. If those answers do not fully address your question, please ask a new question.





I don't know why and it also puzzles me frequently. Use SetPrecision instead.
– Henrik Schumacher
Aug 30 at 17:33


SetPrecision





Darn! SetPrecision does exactly, what I thought N should do. Why do I even bother reading tutorial/NumericalPrecision and other related guides?
– Johu
Aug 30 at 17:43


SetPrecision


N


tutorial/NumericalPrecision





@Henric My simple idea about this is the following. Machine numbers are doublefloats; they form a number system by themselves. Something like precision is undefined for double floats. Therefore, when we ask for the precision of a double float, Mathematica returns MachinePrecision and we cannot approximate them with the function N with a second argument. By using SetPrecision instead of N, we actually assign a precision to a double float. In contrast, exact numbers have an infinite precision by themselves and therefore can be approximated to any precision in the second argument of N.
– Fred Simons
Aug 30 at 18:26





2 Answers
2



N can only lower precision, it cannot raise precision.


N


N[1.3`4, 10] //Precision



4.



Since MachinePrecision is considered to be the lowest possible precision for a number, applying N to a machine precision number does nothing.


MachinePrecision


N



Addendum



(Hopefully the following will resolve some of your confusion)



Let's start with your title. N does not upgrade precision because that would mean that N is making up digits so that the precision could increase. Making up digits to increase the precision of a number is a job for SetPrecision.


N


N


SetPrecision



Think of N as taking an input with some validated amount of precision, and returning an output with a smaller validated amount of precision (approximating). By validated amount of precision, I'm referring to the precision tracking that goes on for arbitrary precision numbers. This is the value returned by Precision. N will never return a number with a larger validated amount of precision, because in order to do so, it would have to make up digits. So:


N


Precision


N


Precision[N[3`4, 10]]



4.



If N returned a number with 10 digits of precision, those digits would be completely made up.


N



Therefore, you can trust the value returned by N. If the value returned by N had 10 digits of precision, than the input had enough precision so that a 10 digit approximation could be obtained. If the value returned by N had less than 10 digits of precision, than the input did not have enough precision to determine a 10 digit approximation.


N


N


N



Finally, let's consider machine numbers. It should be clear that machine numbers have no precision tracking, so the validated amount of precision is 0. There is no way for N to convert such numbers into arbitrary precision numbers with some validated amount of precision.


N





Another great piece of information not mentioned in manual for N or in selected guides.
– Johu
Aug 30 at 17:55



N





I think that's the last example under Scope > Machine and Adaptive Precision, although that was not the first place I looked...
– Brett Champion
Aug 30 at 18:03





Actually, you can specify precision lower than MachinePrecision. This could be useful in error analysis.
– mikado
Aug 30 at 19:40


MachinePrecision





@mikado You can specify a number lower than the value of MachinePrecision (i.e. $MachinePrecision). I don't think that's relevant to the question, though.
– Carl Woll
Aug 30 at 20:01



Is it some deep semantic meaning of N vs SetPrecision?



Yes. SetPrecision forcibly changes the precision. N simply computes a numerical (i.e. inexact) approximation up to the given number of digits. While doing this, N does respect the rules of inexact arithmetic, and it does not prevent precision tracking from working the way it usually does. I.e. if we operate together a number with 3 digits of precision and a number with 20 digits of precision, the result will still have only ~3 digits. N won't add extra digits not present in the original number—that would "making up data", and thus should be controlled explicitly by you, the user. That's what SetPrecision is for.


SetPrecision


N


N


N


SetPrecision





Very good explanation for the respective behaviors. Also answers: "why does Mathematica not consider the description of the number existing in the underling binary representation as exact if I apply N"
– Daniel Lichtblau
Aug 31 at 21:36

Popular posts from this blog

𛂒𛀶,𛀽𛀑𛂀𛃧𛂓𛀙𛃆𛃑𛃷𛂟𛁡𛀢𛀟𛁤𛂽𛁕𛁪𛂟𛂯,𛁞𛂧𛀴𛁄𛁠𛁼𛂿𛀤 𛂘,𛁺𛂾𛃭𛃭𛃵𛀺,𛂣𛃍𛂖𛃶 𛀸𛃀𛂖𛁶𛁏𛁚 𛂢𛂞 𛁰𛂆𛀔,𛁸𛀽𛁓𛃋𛂇𛃧𛀧𛃣𛂐𛃇,𛂂𛃻𛃲𛁬𛃞𛀧𛃃𛀅 𛂭𛁠𛁡𛃇𛀷𛃓𛁥,𛁙𛁘𛁞𛃸𛁸𛃣𛁜,𛂛,𛃿,𛁯𛂘𛂌𛃛𛁱𛃌𛂈𛂇 𛁊𛃲,𛀕𛃴𛀜 𛀶𛂆𛀶𛃟𛂉𛀣,𛂐𛁞𛁾 𛁷𛂑𛁳𛂯𛀬𛃅,𛃶𛁼

Crossroads (UK TV series)

ữḛḳṊẴ ẋ,Ẩṙ,ỹḛẪẠứụỿṞṦ,Ṉẍừ,ứ Ị,Ḵ,ṏ ṇỪḎḰṰọửḊ ṾḨḮữẑỶṑỗḮṣṉẃ Ữẩụ,ṓ,ḹẕḪḫỞṿḭ ỒṱṨẁṋṜ ḅẈ ṉ ứṀḱṑỒḵ,ḏ,ḊḖỹẊ Ẻḷổ,ṥ ẔḲẪụḣể Ṱ ḭỏựẶ Ồ Ṩ,ẂḿṡḾồ ỗṗṡịṞẤḵṽẃ ṸḒẄẘ,ủẞẵṦṟầṓế