C++ source code for compiled plugins (MEX files) to compute various error functions for complex arguments:
** Faddeeva_erf(z)  the error function
** Faddeeva_erfc(z) = 1  erf(z)  complementary error function
** Faddeeva_erfi(z) = i erf(iz)  imaginary error function
** Faddeeva_erfcx(z) = exp(z^2) erfc(z)  scaled complementary error function
** Faddeeva_w(z) = exp(z^2) erfc(iz)  Faddeeva function
** Faddeeva_Dawson(z) = 0.5 sqrt(pi) exp(z^2) erfi(z)  Dawson function
From e.g. the Faddeeva function, one can also obtain the Voigt functions and other related functions.
Assuming you have a C++ compiler (and have configured it in MATLAB with mex setup), compile by running the included Faddeeva_build.m script in MATLAB:
Faddeeva_build
All of the functions have usage of the form:
w = Faddeeva_w(z)
or optionally Faddeeva_w(z, relerr), where relerr is a desired relative error (default: machine precision). z may be an array or matrix of complex or real numbers.
This code may also be downloaded from
http://abinitio.mit.edu/Faddeeva
along with documentation and other versions. As described in the source code, this implementation uses a combination of algorithms for the Faddeeva function: a continuedfraction expansion for large z [similar to G. P. M. Poppe and C. M. J. Wijers, "More efficient computation of the complex error function," ACM Trans. Math. Soft. 16 (1), pp. 38–46 (1990)], and a completely different algorithm for smaller z [Mofreh R. Zaghloul and Ahmed N. Ali, "Algorithm 916: Computing the Faddeyeva and Voigt Functions," ACM Trans. Math. Soft. 38 (2), 15 (2011).]. Given the Faddeeva function, we can then compute the other error functions, although we must switch to Taylor expansions and use other tricks in certain regions of the complex plane to avoid cancellation errors or other floatingpoint problems.
1.5.0.0  portability fixes, slight accuracy improvements 

1.4.0.0  Now includes separate plugins for all of the error functions. 

1.3.0.0  note how to compute erfi using Faddeeva function 

1.1.0.0  Improve accuracy in Re[w(z)] taken by itself. 
Inspired: Voigt model fit
Create scripts with code, output, and formatted text in a single executable document.
George Bishop (view profile)
BTW, substitute Eq.5 from paper by Zaker, J Comp Phys 1969
https://doi.org/10.1016/00219991(69)900114
in Faddeeva definition w(z)=exp(z^2)*(1erf(i*z)) and isolate real and imaginary terms "to derive" in one line those approximations for smaller z. Too trivial actually...
George Bishop (view profile)
I analysed this code and it's appeared to be nothing but a modified Poppe and Wijers algorithm. Only a tiny bit of origin's area for smaller z is computed by Salzer's approximations of the Faddeeva function. Rest is Poppe and Wijers, I guess. Aren't so, Steven?
Kyle Thackston (view profile)
Much faster than the symbolic erf function. Worked for me well on MATLAB 2018a in Windows 10 after following the instructions Andrew posted about MinGW.
sreeraj t (view profile)
So, how do you run this one?
Chaud Solaire (view profile)
I cannot run C++ files via Matlab! I followed all recommended instructions and could not do anything because of permanent errors. It's a headache and terrible frustration to adjust a computer to run these nonnative C++ files from the Matlab environment. Could you please translate these C++ files to Matlab instead?
junye li (view profile)
Faddeeva_erfi(1+1i)
ans =
0.1905 + 1.3162i
which is different from
>> double(erf(sym(1+1i)))
ans =
1.3162 + 0.1905i
Why? I make mex function using Mac.
swank (view profile)
AriessBB (view profile)
Run with Matlab 2017a and compiled with MinGW64 Compiler (C++), it gets
Faddeeva.cc: In function 'cmplx Faddeeva::w(cmplx, double)':
Faddeeva.cc:970:54: error: '_copysign' was not declared in this scope (0.5*c)*copysign(sum5sum4, creal(z)));
Faddeeva.cc:179:26: note: in definition of macro 'C' # define C(a,b) cmplx(a,b)
Is that problem with the MinGW64 Compiler or else?
Andrew (view profile)
I found another Faddeeva / Voigt function that is just as good (maybe better?), and already in MATLAB format (no need for compilers):
https://www.mathworks.com/matlabcentral/fileexchange/47801thevoigtcomplexerrorfunctionsecondversion
Andrew (view profile)
@Andrew
I have discovered that MinGW is not compatible with MATLAB R2015a.
Instead I was able to use SDK 7.1 on Windows 10.
To help others do this, If you install SDK 7.1 from this link:
http://www.microsoft.com/enus/download/details.aspx?id=8279
and still have problems, such as typing into MATLAB
mex setup
and then receiving the error:
No supported SDK or compiler was found on this computer.
For a list of supported compilers, see
<http://www.mathworks.com/support/compilers/R2012a/win64.html>
This means that you need to follow the instructions on this page:
https://www.mathworks.com/matlabcentral/answers/95039whydoesthesdk71installationfailwithaninstallationfailedmessageonmywindowssystem
This worked for me with Windows 10 MATLAB 2015a.
Andrew (view profile)
Has anyone got a C++ compiler for this to work on:
MATLAB R2015a (8.5.0.197613)
64bit (win46)
Windows 10
?
If so, please tell me step by step instructions (I have spent a day trying to make this work, and really need to get fitting data for a publication that is under a review deadline.)
Janos (view profile)
This is a fantastic implementation. This code works about 2000x faster for me (when tested with large multidimensional arrays) than the builtin Matlab erfi function.
Karan Gill (view profile)
Alternatively, the Symbolic Math Toolbox provides the error and dawson functions for complex inputs.
Take the error function for example:
>> double(erf(sym(1+1i)))
ans =
1.3162 + 0.1905i
You could define an anonymous function to make it easier:
>> erfCmplx = @(x) double(erf(sym(x)))
erfCmplx =
@(x)double(erf(sym(x)))
>> erfCmplx(1+1i)
ans =
1.3162 + 0.1905i
The Symbolic Math Toolbox functions are:
erf: http://www.mathworks.com/help/symbolic/erf.html
erfc: http://www.mathworks.com/help/symbolic/erfc.html
erfcinv: http://www.mathworks.com/help/symbolic/erfcinv.html
erfi: http://www.mathworks.com/help/symbolic/erfi.html
erfinv: http://www.mathworks.com/help/symbolic/erfinv.html
dawson: http://www.mathworks.com/help/symbolic/dawson.html
Brian Hannan (view profile)
Ian (view profile)
One update I would really like to see in this package is the derivative functions, at least for W(z). As pointed out in Zaghloul and Ali, the derivative functions of W(z) (equations 2123) become numerically unstable near dV/dx = 0 (V=real(W(z)), at the peak of the Voigt function. This can cause problems when trying to compute analytical Jacobians for doing nonlinear fits of the Voigt function to optical spectra.
Since you are already using the Zaghloul and Ali algorithm in this region, it would be helpful to also use their method to output a function, say Faddeeva_dw(z) = dV/dx + i*dL/dx (L=imag(W(z)). The corresponding y derivatives can then be trivially computed.
Ian (view profile)
Works as advertised and is extremely fast.
Niklas (view profile)
Steven G. Johnson (view profile)
It won't work with lcc, since that is a C compiler, not a C++ compiler. The failure to compile with Visual C++ is a bug, which is fixed in the latest release on my web site (and which should appear on Matlab Central shortly).
alan (view profile)
Hello, Love this program. I have it up and running on my mac to simulate voigt broadening. Im trying to add this functionality to a lab computer running win32 matlab 7.12.0 (R2011a), but I cannot successfully compile with mex.
With Lccwin32 C 2.4.1 in C:\PROGRA~1\MATLAB\R2011a\sys\lcc: lcc preprocessor error: .\Faddeeva.hh:30 .\Faddeeva_mex.cc:35 Faddeeva_w_mex.cc:3 Could not find include file complex
full verbose: http://pastebin.com/YyLYYe8C
With Microsoft Visual C++ 2010 Express in C:\Program Files\Microsoft Visual Studio 10.0: Faddeeva.cc Faddeeva.cc(184) : error C2124: divide or mod by zero Faddeeva.cc(822) : error C3861: 'copysign': identifier not found C:\PROGRA~1\MATLAB\R2011A\BIN\MEX.PL: Error: Compile of 'Faddeeva.cc' failed.
full verbose: http://pastebin.com/K4LPd9FN
The only instances of complex.h are in a pythonwx folder not related to matlab. I see that it was released for R2012a. Does the package need a specific C++ compiler?