Fortran Inverse Error Function
type shall erfinv c++ be REAL. Return value:The return value is of type REAL, of the same kind as X and lies in the range -1 \leq erf (x) \leq 1 . Example: program test_erf real(8) :: x = 0.17_8 x = erf(x) end program test_erf Specific names: Name Argument Return type Standard DERF(X) REAL(8) X REAL(8) GNU extension
here for a quick overview of the site Help Center Detailed answers to any questions you might have Meta Discuss the workings and policies of this site About Us Learn more about Stack Overflow the company Business Learn more about hiring developers or posting ads with us Stack Overflow Questions Jobs Documentation Tags Users Badges Ask Question x Dismiss Join the Stack Overflow Community Stack Overflow is a community of 4.7 million programmers, just like you, helping each other. Join them; it only takes a minute: Sign up Need code for Inverse Error Function up https://gcc.gnu.org/onlinedocs/gfortran/ERF.html vote 5 down vote favorite 4 Does anyone know where I could find code for the "Inverse Error Function?" Freepascal/Delphi would be preferable but C/C++ would be fine too. The TMath/DMath library did not have it :( delphi math pascal share|improve this question asked May 11 '11 at 23:46 Mike Furlender 1,52422053 do you mean inverse erf? –soandos May 11 '11 at 23:47 http://stackoverflow.com/questions/5971830/need-code-for-inverse-error-function Do you refer to erfinv? –David Heffernan May 12 '11 at 0:02 @soandosDavid Yup thats what I'm talking about –Mike Furlender May 12 '11 at 0:07 1 If you can find fortran code convert it to c with f2c -a . If you can find c code great. Compile the c with bcc32 and link with $L that's how I always do it! –David Heffernan May 12 '11 at 0:10 This is a really nifty piece of math. If you find an implementation, I'll make sure it gets added to JEDI Math. It's planned for JEDI Math library in the future, but it looks like it's not in there yet! –Warren P May 12 '11 at 1:49 | show 1 more comment 7 Answers 7 active oldest votes up vote 4 down vote accepted Here's an implementation of erfinv(). Note that for it to work well, you also need a good implementation of erf(). function erfinv(const y: Double): Double; //rational approx coefficients const a: array [0..3] of Double = ( 0.886226899, -1.645349621, 0.914624893, -0.140543331); b: array [0..3] of Double = (-2.118377725, 1.442710462, -0.329097515, 0.012229801); c: array [0..3] of Double = (-1.970840454, -
missing. Despite inverf can be expressed directly with series, we can not use this serie evalution because it is error function slow and lack of precision. Therefore, I found another implementation based on rational function approximation. It is quicker and more precise than the previous implementation from the old inverse error function correlator. Original implementation in Fortan Code (look for MERFI subroutine) --Blanchet 17:50, 17 March 2008 (CET) Retrieved from "https://www.iram.fr/computer/wiki/widex/index.php/Inverse_Error_Function" Views Article Discussion View source History Personal tools Log in / create account Navigation Main Page Community portal Current events Recent changes Random page Help Search Toolbox What links here Related changes Upload file Special pages Printable version Permanent link This page was last modified 16:50, 17 March 2008. This page has been accessed 11,172 times. Privacy policy About Widex Wiki Disclaimers
be down. Please try the request again. Your cache administrator is webmaster. Generated Sun, 16 Oct 2016 00:49:18 GMT by s_ac15 (squid/3.5.20)