888     typedef struct { double real, imag; } __pyx_t_double_complex;  member
1064     #define __Pyx_CIMAG(z) ((z).imag())
1071     #define __Pyx_CIMAG(z) ((z).imag)
1076     #define __Pyx_SET_CIMAG(z,y) ((z).imag(y))
1141 static CYTHON_INLINE int __Pyx_PyErr_GivenExceptionMatches(PyObject *err, PyObject *type);
1142 static CYTHON_INLINE int __Pyx_PyErr_GivenExceptionMatches2(PyObject *err, PyObject *type1, PyObjec…
1145 #define __Pyx_PyErr_GivenExceptionMatches(err, type) PyErr_GivenExceptionMatches(err, type)  argument
1146 …PyErr_GivenExceptionMatches2(err, type1, type2) (PyErr_GivenExceptionMatches(err, type1) || PyErr_…  argument
1686   __pyx_v_Y = std::imag(__pyx_v_Z);  in __pyx_f_5qutip_2cy_4math_zerf()
2875       z.imag = y;  in __pyx_t_double_complex_from_parts()
2884        return (a.real == b.real) && (a.imag == b.imag);  in __Pyx_c_eq_double()
2889         z.imag = a.imag + b.imag;  in __Pyx_c_sum_double()
2895         z.imag = a.imag - b.imag;  in __Pyx_c_diff_double()
2900         z.real = a.real * b.real - a.imag * b.imag;  in __Pyx_c_prod_double()
2901         z.imag = a.real * b.imag + a.imag * b.real;  in __Pyx_c_prod_double()
2906         if (b.imag == 0) {  in __Pyx_c_quot_double()
2907             return __pyx_t_double_complex_from_parts(a.real / b.real, a.imag / b.real);  in __Pyx_c_quot_double()
2908         } else if (fabs(b.real) >= fabs(b.imag)) {  in __Pyx_c_quot_double()
2909             if (b.real == 0 && b.imag == 0) {  in __Pyx_c_quot_double()
2910                 return __pyx_t_double_complex_from_parts(a.real / b.real, a.imag / b.imag);  in __Pyx_c_quot_double()
2912                 double r = b.imag / b.real;  in __Pyx_c_quot_double()
2913                 double s = (double)(1.0) / (b.real + b.imag * r);  in __Pyx_c_quot_double()
2915                     (a.real + a.imag * r) * s, (a.imag - a.real * r) * s);  in __Pyx_c_quot_double()
2918             double r = b.real / b.imag;  in __Pyx_c_quot_double()
2919             double s = (double)(1.0) / (b.imag + b.real * r);  in __Pyx_c_quot_double()
2921                 (a.real * r + a.imag) * s, (a.imag * r - a.real) * s);  in __Pyx_c_quot_double()
2926         if (b.imag == 0) {  in __Pyx_c_quot_double()
2927             return __pyx_t_double_complex_from_parts(a.real / b.real, a.imag / b.real);  in __Pyx_c_quot_double()
2929             double denom = b.real * b.real + b.imag * b.imag;  in __Pyx_c_quot_double()
2931                 (a.real * b.real + a.imag * b.imag) / denom,  in __Pyx_c_quot_double()
2932                 (a.imag * b.real - a.real * b.imag) / denom);  in __Pyx_c_quot_double()
2939         z.imag = -a.imag;  in __Pyx_c_neg_double()
2943        return (a.real == 0) && (a.imag == 0);  in __Pyx_c_is_zero_double()
2948         z.imag = -a.imag;  in __Pyx_c_conj_double()
2954             return sqrt(z.real*z.real + z.imag*z.imag);  in __Pyx_c_abs_double()
2956             return hypot(z.real, z.imag);  in __Pyx_c_abs_double()
2962             if (b.imag == 0 && b.real == (int)b.real) {  in __Pyx_c_pow_double()
2964                     double denom = a.real * a.real + a.imag * a.imag;  in __Pyx_c_pow_double()
2966                     a.imag = -a.imag / denom;  in __Pyx_c_pow_double()
2972                         z.imag = 0;  in __Pyx_c_pow_double()
2986             if (a.imag == 0) {  in __Pyx_c_pow_double()
2989                 } else if (b.imag == 0) {  in __Pyx_c_pow_double()
2991                     z.imag = 0;  in __Pyx_c_pow_double()
3002                 theta = atan2(a.imag, a.real);  in __Pyx_c_pow_double()
3005             z_r = exp(lnr * b.real - theta * b.imag);  in __Pyx_c_pow_double()
3006             z_theta = theta * b.real + lnr * b.imag;  in __Pyx_c_pow_double()
3008             z.imag = z_r * sin(z_theta);  in __Pyx_c_pow_double()
3688 static int __Pyx_inner_PyErr_GivenExceptionMatches2(PyObject *err, PyObject* exc_type1, PyObject* e…  in __Pyx_inner_PyErr_GivenExceptionMatches2()  argument
3694     res = exc_type1 ? PyObject_IsSubclass(err, exc_type1) : 0;  in __Pyx_inner_PyErr_GivenExceptionMatches2()
3696         PyErr_WriteUnraisable(err);  in __Pyx_inner_PyErr_GivenExceptionMatches2()
3700         res = PyObject_IsSubclass(err, exc_type2);  in __Pyx_inner_PyErr_GivenExceptionMatches2()
3702             PyErr_WriteUnraisable(err);  in __Pyx_inner_PyErr_GivenExceptionMatches2()
3710 static CYTHON_INLINE int __Pyx_inner_PyErr_GivenExceptionMatches2(PyObject *err, PyObject* exc_type…  in __Pyx_inner_PyErr_GivenExceptionMatches2()  argument
3711     int res = exc_type1 ? __Pyx_IsSubtype((PyTypeObject*)err, (PyTypeObject*)exc_type1) : 0;  in __Pyx_inner_PyErr_GivenExceptionMatches2()
3713         res = __Pyx_IsSubtype((PyTypeObject*)err, (PyTypeObject*)exc_type2);  in __Pyx_inner_PyErr_GivenExceptionMatches2()
3739 static CYTHON_INLINE int __Pyx_PyErr_GivenExceptionMatches(PyObject *err, PyObject* exc_type) {  in __Pyx_PyErr_GivenExceptionMatches()  argument
3740     if (likely(err == exc_type)) return 1;  in __Pyx_PyErr_GivenExceptionMatches()
3741     if (likely(PyExceptionClass_Check(err))) {  in __Pyx_PyErr_GivenExceptionMatches()
3743             return __Pyx_inner_PyErr_GivenExceptionMatches2(err, NULL, exc_type);  in __Pyx_PyErr_GivenExceptionMatches()
3745             return __Pyx_PyErr_GivenExceptionMatchesTuple(err, exc_type);  in __Pyx_PyErr_GivenExceptionMatches()
3749     return PyErr_GivenExceptionMatches(err, exc_type);  in __Pyx_PyErr_GivenExceptionMatches()
3751 static CYTHON_INLINE int __Pyx_PyErr_GivenExceptionMatches2(PyObject *err, PyObject *exc_type1, PyO…  in __Pyx_PyErr_GivenExceptionMatches2()  argument
3754     if (likely(err == exc_type1 || err == exc_type2)) return 1;  in __Pyx_PyErr_GivenExceptionMatches2()
3755     if (likely(PyExceptionClass_Check(err))) {  in __Pyx_PyErr_GivenExceptionMatches2()
3756         return __Pyx_inner_PyErr_GivenExceptionMatches2(err, exc_type1, exc_type2);  in __Pyx_PyErr_GivenExceptionMatches2()
3758 …return (PyErr_GivenExceptionMatches(err, exc_type1) || PyErr_GivenExceptionMatches(err, exc_type2)…  in __Pyx_PyErr_GivenExceptionMatches2()