//! @file a68g-double.h
//! @author J. Marcel van der Veer

//! @section Copyright
//!
//! This file is part of Algol68G - an Algol 68 compiler-interpreter.
//! Copyright 2001-2024 J. Marcel van der Veer [algol68g@xs4all.nl].

//! @section License
//!
//! This program is free software; you can redistribute it and/or modify it 
//! under the terms of the GNU General Public License as published by the 
//! Free Software Foundation; either version 3 of the License, or 
//! (at your option) any later version.
//!
//! This program is distributed in the hope that it will be useful, but 
//! WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY 
//! or FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License for 
//! more details. You should have received a copy of the GNU General Public 
//! License along with this program. If not, see [http://www.gnu.org/licenses/].

//! @section Synopsis
//!
//! LONG REAL definitions.

#if !defined (__A68G_DOUBLE_H__)
#define __A68G_DOUBLE_H__

#if (A68_LEVEL >= 3)

#define MODCHK(p, m, c) (!(MODULAR_MATH (p) && (m == M_LONG_BITS)) && (c))

#if defined (HAVE_IEEE_754)
#define CHECK_DOUBLE_REAL(p, u) PRELUDE_ERROR (!finite_double (u), p, ERROR_INFINITE, M_LONG_REAL)
#define CHECK_DOUBLE_COMPLEX(p, u, v)\
  PRELUDE_ERROR (isinf_double (u), p, ERROR_INFINITE, M_LONG_REAL);\
  PRELUDE_ERROR (isinf_double (v), p, ERROR_INFINITE, M_LONG_REAL);
#else
#define CHECK_DOUBLE_REAL(p, u) {;}
#define CHECK_DOUBLE_COMPLEX(p, u, v) {;}
#endif

#define LONG_INT_BASE (9223372036854775808.0q)
#define HW(z) ((z).u[1])
#define LW(z) ((z).u[0])
#define D_NEG(d) ((HW(d) & D_SIGN) != 0)
#define D_LT(u, v) (HW (u) < HW (v) ? A68_TRUE : (HW (u) == HW (v) ? LW (u) < LW (v) : A68_FALSE))

#define RADIX (65536)
#define RADIX_Q (65536.0q)
#define CONST_2_UP_112_Q (5192296858534827628530496329220096.0q)

#define IS_ZERO(u) (HW (u) == 0 && LW (u) == 0)
#define EQ(u, v) (HW (u) == HW (v) && LW (u) == LW (v))
#define GT(u, v) (HW (u) != HW (v) ? HW (u) > HW (v) : LW (u) > LW (v))
#define GE(u, v) (HW (u) != HW (v) ? HW (u) >= HW (v) : LW (u) >= LW (v))

#define acos_double acosq
#define acosh_double acoshq
#define asin_double asinq
#define asinh_double asinhq
#define atan2_double atan2q
#define atan_double atanq
#define atanh_double atanhq
#define cacos_double cacosq
#define cacosh_double cacoshq
#define casin_double casinq
#define casinh_double casinhq
#define catan_double catanq
#define catanh_double catanhq
#define cbrt_double cbrtq
#define ccos_double ccosq
#define ccosh_double ccoshq
#define cexp_double cexpq
#define cimag_double cimagq
#define clog_double clogq
#define cos_double cosq
#define cosh_double coshq
#define creal_double crealq
#define csin_double csinq
#define csinh_double csinhq
#define csqrt_double csqrtq
#define ctan_double ctanq
#define ctanh_double ctanhq
#define erfc_double erfcq
#define erf_double erfq
#define exp_double expq
#define fabs_double fabsq
#define finite_double finiteq
#define floor_double floorq
#define fmod_double fmodq
#define isinf_double isinfq
#define lgamma_double lgammaq
#define log10_double log10q
#define log_double logq
#define pow_double powq
#define sin_double sinq
#define sinh_double sinhq
#define sqrt_double sqrtq
#define tan_double tanq
#define tanh_double tanhq
#define tgamma_double tgammaq
#define trunc_double truncq

#define DBLEQ(z) ((dble_double (A68 (f_entry), (z))).f)

#define ABSQ(n) ((n) >= 0.0q ? (n) : -(n))

#define POP_LONG_COMPLEX(p, re, im) {\
  POP_OBJECT (p, im, A68_LONG_REAL);\
  POP_OBJECT (p, re, A68_LONG_REAL);\
  }

#define set_lw(z, k) {LW(z) = k; HW(z) = 0;}
#define set_hw(z, k) {LW(z) = 0; HW(z) = k;}
#define set_hwlw(z, h, l) {LW (z) = l; HW (z) = h;}
#define D_ZERO(z) (HW (z) == 0 && LW (z) == 0)

#define add_double(p, m, w, u, v) {\
    DOUBLE_NUM_T _ww_;\
    LW (_ww_) = LW (u) + LW (v);\
    HW (_ww_) = HW (u) + HW (v);\
    PRELUDE_ERROR (MODCHK (p, m, HW (_ww_) < HW (v)), p, ERROR_MATH, (m));\
    if (LW (_ww_) < LW (v)) {\
      HW (_ww_)++;\
      PRELUDE_ERROR (MODCHK (p, m, HW (_ww_) < 1), p, ERROR_MATH, (m));\
    }\
    w = _ww_;\
  }

#define sub_double(p, m, w, u, v) {\
    DOUBLE_NUM_T _ww_;\
    LW (_ww_) = LW (u) - LW (v);\
    HW (_ww_) = HW (u) - HW (v);\
    PRELUDE_ERROR (MODCHK (p, m, HW (_ww_) > HW (u)), p, ERROR_MATH, (m));\
    if (LW (_ww_) > LW (u)) {\
      PRELUDE_ERROR (MODCHK (p, m, HW (_ww_) == 0), p, ERROR_MATH, (m));\
      HW (_ww_)--;\
    }\
    w = _ww_;\
  }

static inline DOUBLE_NUM_T dble (DOUBLE_T x)
{
  DOUBLE_NUM_T w;
  w.f = x;
  return w;
}

static inline int sign_double_int (DOUBLE_NUM_T w)
{
  if (D_NEG (w)) {
    return -1;
  } else if (D_ZERO (w)) {
    return 0;
  } else {
    return 1;
  }
}

static inline int sign_double (DOUBLE_NUM_T w)
{
  if (w.f < 0.0q) {
    return -1;
  } else if (w.f == 0.0q) {
    return 0;
  } else {
    return 1;
  }
}

static inline DOUBLE_NUM_T abs_double_int (DOUBLE_NUM_T z)
{
  DOUBLE_NUM_T w;
  LW (w) = LW (z);
  HW (w) = HW (z) & (~D_SIGN);
  return w;
}

static inline DOUBLE_NUM_T neg_double_int (DOUBLE_NUM_T z)
{
  DOUBLE_NUM_T w;
  LW (w) = LW (z);
  if (D_NEG (z)) {
    HW (w) = HW (z) & (~D_SIGN);
  } else {
    HW (w) = HW (z) | D_SIGN;
  }
  return w;
}

void m64to128 (DOUBLE_NUM_T *, UNSIGNED_T, UNSIGNED_T);
void m128to128 (NODE_T *, MOID_T *, DOUBLE_NUM_T *, DOUBLE_NUM_T, DOUBLE_NUM_T);
DOUBLE_NUM_T double_udiv (NODE_T *, MOID_T *, DOUBLE_NUM_T, DOUBLE_NUM_T, int);
DOUBLE_NUM_T double_uadd (NODE_T *, MOID_T *, DOUBLE_NUM_T, DOUBLE_NUM_T);
DOUBLE_NUM_T double_usub (NODE_T *, MOID_T *, DOUBLE_NUM_T, DOUBLE_NUM_T);
DOUBLE_NUM_T double_umul (NODE_T *, MOID_T *, DOUBLE_NUM_T, DOUBLE_NUM_T);
DOUBLE_NUM_T double_sadd (NODE_T *, DOUBLE_NUM_T, DOUBLE_NUM_T);
DOUBLE_NUM_T double_ssub (NODE_T *, DOUBLE_NUM_T, DOUBLE_NUM_T);
DOUBLE_NUM_T double_smul (NODE_T *, DOUBLE_NUM_T, DOUBLE_NUM_T);
DOUBLE_NUM_T double_sdiv (NODE_T *, DOUBLE_NUM_T, DOUBLE_NUM_T, int);

int sign_double_int (DOUBLE_NUM_T);
int sign_double (DOUBLE_NUM_T);
int string_to_double_int (NODE_T *, A68_LONG_INT *, char *);
DOUBLE_T a68_hypot_double (DOUBLE_T, DOUBLE_T);
DOUBLE_T string_to_double (char *, char **);
DOUBLE_T inverf_double (DOUBLE_T);
DOUBLE_NUM_T abs_double_int (DOUBLE_NUM_T);
DOUBLE_NUM_T bits_to_double_int (NODE_T *, char *);
DOUBLE_NUM_T dble_double (NODE_T *, REAL_T);
DOUBLE_NUM_T double_int_to_double (NODE_T *, DOUBLE_NUM_T);
DOUBLE_NUM_T double_strtou (NODE_T *, char *);
DOUBLE_NUM_T double_udiv (NODE_T *, MOID_T *, DOUBLE_NUM_T, DOUBLE_NUM_T, int);
DOUBLE_T a68_neginf_double (void);
DOUBLE_T a68_posinf_double (void);
void deltagammainc_double (DOUBLE_T *, DOUBLE_T *, DOUBLE_T, DOUBLE_T, DOUBLE_T, DOUBLE_T);

GPROC genie_infinity_double;
GPROC genie_minus_infinity_double;
GPROC genie_gamma_inc_g_double;
GPROC genie_gamma_inc_f_double;
GPROC genie_gamma_inc_h_double;
GPROC genie_gamma_inc_gf_double;
GPROC genie_abs_double_compl;
GPROC genie_abs_double_int;
GPROC genie_abs_double;
GPROC genie_acos_double_compl;
GPROC genie_acosdg_double;
GPROC genie_acosh_double_compl;
GPROC genie_acosh_double;
GPROC genie_acos_double;
GPROC genie_acotdg_double;
GPROC genie_acot_double;
GPROC genie_asec_double;
GPROC genie_asecdg_double;
GPROC genie_acsc_double;
GPROC genie_acscdg_double;
GPROC genie_add_double_compl;
GPROC genie_add_double_bits;
GPROC genie_add_double_int;
GPROC genie_add_double;
GPROC genie_add_double;
GPROC genie_and_double_bits;
GPROC genie_arg_double_compl;
GPROC genie_asin_double_compl;
GPROC genie_asindg_double;
GPROC genie_asindg_double;
GPROC genie_asinh_double_compl;
GPROC genie_asinh_double;
GPROC genie_asin_double;
GPROC genie_atan2dg_double;
GPROC genie_atan2_double;
GPROC genie_atan_double_compl;
GPROC genie_atandg_double;
GPROC genie_atanh_double_compl;
GPROC genie_atanh_double;
GPROC genie_atan_double;
GPROC genie_bin_double_int;
GPROC genie_clear_double_bits;
GPROC genie_conj_double_compl;
GPROC genie_cas_double;
GPROC genie_cos_double_compl;
GPROC genie_cosdg_double;
GPROC genie_cosdg_double;
GPROC genie_cosh_double_compl;
GPROC genie_cosh_double;
GPROC genie_cospi_double;
GPROC genie_cos_double;
GPROC genie_cotdg_double;
GPROC genie_cotpi_double;
GPROC genie_cot_double;
GPROC genie_sec_double;
GPROC genie_secdg_double;
GPROC genie_csc_double;
GPROC genie_cscdg_double;
GPROC genie_curt_double;
GPROC genie_divab_double_compl;
GPROC genie_divab_double;
GPROC genie_divab_double;
GPROC genie_div_double_compl;
GPROC genie_div_double_int;
GPROC genie_double_bits_pack;
GPROC genie_double_max_bits;
GPROC genie_double_max_int;
GPROC genie_double_max_real;
GPROC genie_double_min_real;
GPROC genie_double_small_real;
GPROC genie_double_zeroin;
GPROC genie_elem_double_bits;
GPROC genie_entier_double;
GPROC genie_eq_double_compl;
GPROC genie_eq_double_bits;
GPROC genie_eq_double_int;
GPROC genie_eq_double_int;
GPROC genie_eq_double;
GPROC genie_eq_double;
GPROC genie_eq_double;
GPROC genie_eq_double;
GPROC genie_erfc_double;
GPROC genie_erf_double;
GPROC genie_exp_double_compl;
GPROC genie_exp_double;
GPROC genie_gamma_double;
GPROC genie_ge_double_bits;
GPROC genie_ge_double_int;
GPROC genie_ge_double_int;
GPROC genie_ge_double;
GPROC genie_ge_double;
GPROC genie_ge_double;
GPROC genie_ge_double;
GPROC genie_gt_double_bits;
GPROC genie_gt_double_int;
GPROC genie_gt_double_int;
GPROC genie_gt_double;
GPROC genie_gt_double;
GPROC genie_gt_double;
GPROC genie_gt_double;
GPROC genie_i_double_compl;
GPROC genie_i_int_double_compl;
GPROC genie_im_double_compl;
GPROC genie_inverfc_double;
GPROC genie_inverf_double;
GPROC genie_le_double_bits;
GPROC genie_le_double_int;
GPROC genie_le_double_int;
GPROC genie_lengthen_bits_to_double_bits;
GPROC genie_lengthen_double_compl_to_long_mp_complex;
GPROC genie_lengthen_complex_to_double_compl;
GPROC genie_lengthen_double_int_to_mp;
GPROC genie_lengthen_int_to_double_int;
GPROC genie_lengthen_double_to_mp;
GPROC genie_lengthen_real_to_double;
GPROC genie_le_double;
GPROC genie_le_double;
GPROC genie_le_double;
GPROC genie_le_double;
GPROC genie_ln_double_compl;
GPROC genie_lngamma_double;
GPROC genie_ln_double;
GPROC genie_log_double;
GPROC genie_lt_double_bits;
GPROC genie_lt_double_int;
GPROC genie_lt_double_int;
GPROC genie_lt_double;
GPROC genie_lt_double;
GPROC genie_lt_double;
GPROC genie_lt_double;
GPROC genie_minusab_double_compl;
GPROC genie_minusab_double_int;
GPROC genie_minusab_double_int;
GPROC genie_minusab_double;
GPROC genie_minusab_double;
GPROC genie_minus_double_compl;
GPROC genie_minus_double_int;
GPROC genie_minus_double;
GPROC genie_modab_double_int;
GPROC genie_modab_double_int;
GPROC genie_mod_double_bits;
GPROC genie_mod_double_int;
GPROC genie_mul_double_compl;
GPROC genie_mul_double_int;
GPROC genie_mul_double;
GPROC genie_mul_double;
GPROC genie_ne_double_compl;
GPROC genie_ne_double_bits;
GPROC genie_ne_double_int;
GPROC genie_ne_double_int;
GPROC genie_ne_double_int;
GPROC genie_ne_double_int;
GPROC genie_ne_double;
GPROC genie_ne_double;
GPROC genie_ne_double;
GPROC genie_ne_double;
GPROC genie_ne_double;
GPROC genie_ne_double;
GPROC genie_ne_double;
GPROC genie_ne_double;
GPROC genie_next_random_double;
GPROC genie_not_double_bits;
GPROC genie_odd_double_int;
GPROC genie_or_double_bits;
GPROC genie_overab_double_int;
GPROC genie_overab_double_int;
GPROC genie_over_double_bits;
GPROC genie_over_double_int;
GPROC genie_over_double;
GPROC genie_over_double;
GPROC genie_pi_double;
GPROC genie_plusab_double_compl;
GPROC genie_plusab_double_int;
GPROC genie_plusab_double_int;
GPROC genie_plusab_double;
GPROC genie_pow_double_compl_int;
GPROC genie_pow_double_int_int;
GPROC genie_pow_double;
GPROC genie_pow_double_int;
GPROC genie_re_double_compl;
GPROC genie_rol_double_bits;
GPROC genie_ror_double_bits;
GPROC genie_round_double;
GPROC genie_set_double_bits;
GPROC genie_shl_double_bits;
GPROC genie_shorten_double_compl_to_complex;
GPROC genie_shorten_double_bits_to_bits;
GPROC genie_shorten_long_int_to_int;
GPROC genie_shorten_long_mp_complex_to_double_compl;
GPROC genie_shorten_mp_to_double_int;
GPROC genie_shorten_mp_to_double;
GPROC genie_shorten_double_to_real;
GPROC genie_shr_double_bits;
GPROC genie_sign_double_int;
GPROC genie_sign_double;
GPROC genie_sin_double_compl;
GPROC genie_sindg_double;
GPROC genie_sinh_double_compl;
GPROC genie_sinh_double;
GPROC genie_sinpi_double;
GPROC genie_sin_double;
GPROC genie_sqrt_double_compl;
GPROC genie_sqrt_double;
GPROC genie_sqrt_double;
GPROC genie_sqrt_double;
GPROC genie_sub_double_compl;
GPROC genie_sub_double_bits;
GPROC genie_sub_double_int;
GPROC genie_sub_double;
GPROC genie_sub_double;
GPROC genie_tan_double_compl;
GPROC genie_tandg_double;
GPROC genie_tanh_double_compl;
GPROC genie_tanh_double;
GPROC genie_tanpi_double;
GPROC genie_tan_double;
GPROC genie_timesab_double_compl;
GPROC genie_timesab_double_int;
GPROC genie_timesab_double_int;
GPROC genie_timesab_double;
GPROC genie_timesab_double;
GPROC genie_times_double_bits;
GPROC genie_widen_double_int_to_double;
GPROC genie_xor_double_bits;
GPROC genie_beta_inc_cf_double;
GPROC genie_beta_double;
GPROC genie_ln_beta_double;
GPROC genie_gamma_inc_double;
GPROC genie_zero_double_int;

#endif

#endif
