I'm checking two different procedures that will eventually be used with regular floating point arithmetic. I need to estimate numerical error arising from floating point arithmetic.

Is there a way to force computation on floating point numbers? Model problem could be error estimation for:

Erf[x] - Erf[x+h]

for small vaues of $h$ and large values of $x$.

  • 4
    $\begingroup$ How about if you try setting $MinPrecision = $MachinePrecision; $MaxPrecision = $MachinePrecision; before the computation. Will this do what you want? if you give a MWE and what output you expect, I could try it. You can put the above inside a Block also. $\endgroup$
    – Nasser
    Apr 28 '19 at 10:19
  • $\begingroup$ @MichaelE2 that's the whole point here. Sometimes it's more accurate to find this difference directly and sometimes it's better to approximate by series in $h$. I'm trying to investigate when each method is better, so I need both arbitrary precision arithmetic for correct result and floating to check numerical errors. $\endgroup$
    – Radost
    Apr 28 '19 at 15:40

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service, privacy policy and cookie policy

Browse other questions tagged or ask your own question.