+ fcts re_wofz, im_wofz

pub/lib/fbase.cpp

@@ -711,6 +711,16 @@ double func_voigt_hwhm(double s, double g)
     return voigt_hwhm(s, g);
 }
 
+double func_re_wofz(double x, double y)
+{
+    return creal(w_of_z({x, y}));
+}
+
+double func_im_wofz(double x, double y)
+{
+    return cimag(w_of_z({x, y}));
+}
+
 double func_zorn(double v, double vavg, double s)
 {
     if (s <= 0 || s >= 1)
@@ -1063,6 +1073,10 @@ void fbase_initialize()
     G->register_fct_d_dd(m, func_heat_sphere);
     m = {"voigt_hwhm", 3, "(sigma,gamma): hwhm of voigt(x,sigma,gamma)"};
     G->register_fct_d_dd(m, func_voigt_hwhm);
+    m = {"im_wofz", 3, "(x,y): imaginary part of Faddeeva's function"};
+    G->register_fct_d_dd(m, func_im_wofz);
+    m = {"re_wofz", 3, "(x,y): real part of Faddeeva's function"};
+    G->register_fct_d_dd(m, func_re_wofz);
 
     // f(3 args)
     m = {"rehavneg", 3, "(x,y,z): real part of the Havriliak-Negami function"};