endfunction
[Xr,Yr]=erf_ref();
y=erf(Xr);
-assert_checkalmostequal(Yr,y);
+assert_checkalmostequal(Yr,y,sqrt(%eps),sqrt(%eps));
// Tests are from faddeeva code:
assert_checkalmostequal(erf(1 + 2*%i), -0.5366435657785650339917955593141927494421 + -5.049143703447034669543036958614140565553*%i, 1.e-10);
assert_checkalmostequal(erf(-1 + 2*%i), 0.5366435657785650339917955593141927494421 + -5.049143703447034669543036958614140565553*%i, 1.e-10);