void c_qd_asinh(const double *a, double *b) { qd_real bb; bb = asinh(qd_real(a)); TO_DOUBLE_PTR(bb, b); }
void c_qd_tan(const double *a, double *b) { qd_real bb; bb = tan(qd_real(a)); TO_DOUBLE_PTR(bb, b); }
void c_qd_atan2(const double *a, const double *b, double *c) { qd_real cc; cc = atan2(qd_real(a), qd_real(b)); TO_DOUBLE_PTR(cc, c); }
void c_qd_log10(const double *a, double *b) { qd_real bb; bb = log10(qd_real(a)); TO_DOUBLE_PTR(bb, b); }
void c_qd_add_qd_dd(const double *a, const double *b, double *c) { qd_real cc; cc = qd_real(a) + dd_real(b); TO_DOUBLE_PTR(cc, c); }
void c_qd_floor(const double *a, double *b) { qd_real bb; bb = floor(qd_real(a)); TO_DOUBLE_PTR(bb, b); }
void c_qd_ceil(const double *a, double *b) { qd_real bb; bb = ceil(qd_real(a)); TO_DOUBLE_PTR(bb, b); }
void c_qd_nroot(const double *a, int n, double *b) { qd_real bb; bb = nroot(qd_real(a), n); TO_DOUBLE_PTR(bb, b); }
void c_qd_nint(const double *a, double *b) { qd_real bb; bb = nint(qd_real(a)); TO_DOUBLE_PTR(bb, b); }
/* selfdiv */ void c_qd_selfdiv(const double *a, double *b) { qd_real bb(b); bb /= qd_real(a); TO_DOUBLE_PTR(bb, b); }
void c_qd_div_d_qd(double a, const double *b, double *c) { qd_real cc; cc = a / qd_real(b); TO_DOUBLE_PTR(cc, c); }
void c_qd_div_qd_d(const double *a, double b, double *c) { qd_real cc; cc = qd_real(a) / b; TO_DOUBLE_PTR(cc, c); }
void c_qd_div_dd_qd(const double *a, const double *b, double *c) { qd_real cc; cc = dd_real(a) / qd_real(b); TO_DOUBLE_PTR(cc, c); }
scalar_pt<qd_real>() : x(qd_real(0)), y(qd_real(0)) {}