void c_qd_selfadd_d(double a, double *b) { qd_real bb(b); bb += a; TO_DOUBLE_PTR(bb, b); }
void c_qd_pi(double *a) { TO_DOUBLE_PTR(qd_real::_pi, a); }
void c_dd_log10(const double *a, double *b) { dd_real bb; bb = log10(dd_real(a)); TO_DOUBLE_PTR(bb, b); }
void c_dd_atan2(const double *a, const double *b, double *c) { dd_real cc; cc = atan2(dd_real(a), dd_real(b)); TO_DOUBLE_PTR(cc, c); }
void c_dd_nroot(const double *a, int n, double *b) { dd_real bb; bb = nroot(dd_real(a), n); TO_DOUBLE_PTR(bb, b); }
void c_dd_floor(const double *a, double *b) { dd_real bb; bb = floor(dd_real(a)); TO_DOUBLE_PTR(bb, b); }
void c_qd_sqrt(const double *a, double *b) { qd_real bb; bb = sqrt(qd_real(a)); TO_DOUBLE_PTR(bb, b); }
void c_qd_mul_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); }
void c_qd_selfdiv_dd(const double *a, double *b) { qd_real bb(b); bb /= dd_real(a); TO_DOUBLE_PTR(bb, b); }
void c_qd_selfdiv_d(double a, double *b) { qd_real bb(b); bb /= a; TO_DOUBLE_PTR(bb, b); }
void c_qd_selfmul_d(double a, double *b) { qd_real bb(b); bb *= a; TO_DOUBLE_PTR(bb, b); }
/* selfmul */ void c_qd_selfmul(const double *a, double *b) { qd_real bb(b); bb *= qd_real(a); TO_DOUBLE_PTR(bb, b); }
void c_qd_selfsub_d(double a, double *b) { qd_real bb(b); bb -= a; TO_DOUBLE_PTR(bb, b); }
void c_qd_sub_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_npwr(const double *a, int n, double *b) { qd_real bb; bb = npwr(qd_real(a), n); TO_DOUBLE_PTR(bb, b); }
void c_qd_sub_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_cos(const double *a, double *b) { qd_real bb; bb = cos(qd_real(a)); TO_DOUBLE_PTR(bb, b); }
void c_qd_mul_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_atanh(const double *a, double *b) { qd_real bb; bb = atanh(qd_real(a)); TO_DOUBLE_PTR(bb, b); }
void c_dd_nint(const double *a, double *b) { dd_real bb; bb = nint(dd_real(a)); TO_DOUBLE_PTR(bb, b); }
void c_qd_sincosh(const double *a, double *s, double *c) { qd_real ss, cc; sincosh(qd_real(a), ss, cc); TO_DOUBLE_PTR(cc, c); TO_DOUBLE_PTR(ss, s); }
void c_dd_ceil(const double *a, double *b) { dd_real bb; bb = ceil(dd_real(a)); TO_DOUBLE_PTR(bb, b); }
void c_qd_read(const char *s, double *a) { qd_real aa(s); TO_DOUBLE_PTR(aa, a); }
void c_dd_tan(const double *a, double *b) { dd_real bb; bb = tan(dd_real(a)); TO_DOUBLE_PTR(bb, b); }
void c_qd_rand(double *a) { qd_real aa; aa = qdrand(); TO_DOUBLE_PTR(aa, a); }
void c_dd_asinh(const double *a, double *b) { dd_real bb; bb = asinh(dd_real(a)); TO_DOUBLE_PTR(bb, b); }
/* div */ void c_qd_div(const double *a, const double *b, double *c) { qd_real cc; cc = qd_real(a) / qd_real(b); TO_DOUBLE_PTR(cc, c); }