- Cuchulainn
**Posts:**60865**Joined:****Location:**Amsterdam-
**Contact:**

QuoteOriginally posted by: outrunYes,.. that what's need to be done.I'm using this "casio" site a lot (the one from the 80s digital watches with a calculator in them!), it has all sorts of high precision calculators:http://keisan.casio.com/exec/system/1180573449erf(0.01) = 0.011283415555849616915909523548069210949100954830116erf(0.1) = 0.1124629160182848922032750717439683832216962991597erf(0.2) = 0.22270258921047845414013900680014381638826903843023erf(0.8) = 0.74210096470766048616711058650294587731768957991471erf(2.0) = 0.99532226501895273416206925636725292861089179704006erf(4.0) = 0.99999998458274209971998114784032651311595142785475Why is accuracy important in this case?

http://www.datasimfinancial.com

http://www.datasim.nl

Approach your problem from the right end and begin with the answers. Then one day, perhaps you will find the final question..

R. van Gulik

http://www.datasim.nl

Approach your problem from the right end and begin with the answers. Then one day, perhaps you will find the final question..

R. van Gulik

- Cuchulainn
**Posts:**60865**Joined:****Location:**Amsterdam-
**Contact:**

QuoteOriginally posted by: outrunQuoteOriginally posted by: CuchulainnQuoteOriginally posted by: outrunYes,.. that what's need to be done.I'm using this "casio" site a lot (the one from the 80s digital watches with a calculator in them!), it has all sorts of high precision calculators:http://keisan.casio.com/exec/system/1180573449erf(0.01) = 0.011283415555849616915909523548069210949100954830116erf(0.1) = 0.1124629160182848922032750717439683832216962991597erf(0.2) = 0.22270258921047845414013900680014381638826903843023erf(0.8) = 0.74210096470766048616711058650294587731768957991471erf(2.0) = 0.99532226501895273416206925636725292861089179704006erf(4.0) = 0.99999998458274209971998114784032651311595142785475Why is accuracy important in this case?It's for some high precision hobby project that I want to work on at some point in the future when I have more time. The idea is to sell it and retire of course.I see. Have you tried "Dragons' Den"?

http://www.datasimfinancial.com

http://www.datasim.nl

Approach your problem from the right end and begin with the answers. Then one day, perhaps you will find the final question..

R. van Gulik

http://www.datasim.nl

Approach your problem from the right end and begin with the answers. Then one day, perhaps you will find the final question..

R. van Gulik

- Cuchulainn
**Posts:**60865**Joined:****Location:**Amsterdam-
**Contact:**

QuoteOriginally posted by: CuchulainnQuoteI also want to know this: * is erf slow compared to boost?* does is give valid (accurate, 80 bit) results for "long double"?* does boost multiprecision overload the std::erf for the multiprecision type? Or do I need to use boost:erf?Should be easy to determine. Boost mailing. A statistical regression framework is needed, does such a framework exist? My feeling is NO.So, press a button, let code run for 2 hours and look at report "std erf 20% slower than boost erf" etc.upper estimate possible. As well as accuracy estimate.These tests take about [8,24] manhours?

http://www.datasimfinancial.com

http://www.datasim.nl

Approach your problem from the right end and begin with the answers. Then one day, perhaps you will find the final question..

R. van Gulik

http://www.datasim.nl

Approach your problem from the right end and begin with the answers. Then one day, perhaps you will find the final question..

R. van Gulik

- Cuchulainn
**Posts:**60865**Joined:****Location:**Amsterdam-
**Contact:**

In 2006, BoostQuoteRewrite erf/erfc as rational approximations (valid to 128-bit precision). ?

http://www.datasim.nl

Approach your problem from the right end and begin with the answers. Then one day, perhaps you will find the final question..

R. van Gulik

- Cuchulainn
**Posts:**60865**Joined:****Location:**Amsterdam-
**Contact:**

AFAIK, Boost multiprecision does not support erf(x) (but it does support exp(x):)) but if you want high precision you could try1) Convert erf(x) to an initial value problem du/dt = .. exp(-x^2), u(-infinity) = 02) Solve IVP using Boost odeint that supports multiprecision. Have you tried any of these libraries

Last edited by Cuchulainn on January 21st, 2016, 11:00 pm, edited 1 time in total.

http://www.datasim.nl

Approach your problem from the right end and begin with the answers. Then one day, perhaps you will find the final question..

R. van Gulik

- Cuchulainn
**Posts:**60865**Joined:****Location:**Amsterdam-
**Contact:**

QuoteOriginally posted by: outrunThanks. Ik think eg 2) here (pdf) is good. Looks simple.I notices that boost multiprecision has a set of eval_cos() eval_atan2() routines in the detail namespace, it would be nice to add erf to that list.Looks doable.Maybe accelerate the series summation with Shanks transformation?

Last edited by Cuchulainn on January 21st, 2016, 11:00 pm, edited 1 time in total.

http://www.datasim.nl

Approach your problem from the right end and begin with the answers. Then one day, perhaps you will find the final question..

R. van Gulik

- Cuchulainn
**Posts:**60865**Joined:****Location:**Amsterdam-
**Contact:**

QuoteOriginally posted by: outrunQuoteOriginally posted by: CuchulainnQuoteOriginally posted by: outrunThanks. Ik think eg 2) here (pdf) is good. Looks simple.I notices that boost multiprecision has a set of eval_cos() eval_atan2() routines in the detail namespace, it would be nice to add erf to that list.Looks doable.Maybe accelerate the series summation with Shanks transformation?That would be a nice innovation!I have a feeling that something like Kahan summation is also needed as well?

http://www.datasim.nl

Approach your problem from the right end and begin with the answers. Then one day, perhaps you will find the final question..

R. van Gulik

- Cuchulainn
**Posts:**60865**Joined:****Location:**Amsterdam-
**Contact:**

QuoteOriginally posted by: outrunQuoteOriginally posted by: CuchulainnQuoteOriginally posted by: outrunQuoteOriginally posted by: CuchulainnQuoteOriginally posted by: outrunThanks. Ik think eg 2) here (pdf) is good. Looks simple.I notices that boost multiprecision has a set of eval_cos() eval_atan2() routines in the detail namespace, it would be nice to add erf to that list.Looks doable.Maybe accelerate the series summation with Shanks transformation?That would be a nice innovation!I have a feeling that something like Kahan summation is also needed as well?Full precision (up to the number of digits of the float) is actually the main topic of that paper.Clear. But can round-off error creep into the computations and hence compromise the accuracy?

Last edited by Cuchulainn on January 21st, 2016, 11:00 pm, edited 1 time in total.

http://www.datasim.nl

Approach your problem from the right end and begin with the answers. Then one day, perhaps you will find the final question..

R. van Gulik

GZIP: On