fun_ieee.c File Reference

#include "calc.h"

Include dependency graph for fun_ieee.c:

Go to the source code of this file.

Macros
#define	NEXT(n, i)   (((n) + (i)/(n)) >> 1)

Functions
static double	validate_rad2angle (double a)
static double	validate_angle2rad (calc_number_t *c)
void	apply_int_mask (calc_number_t *r)
double	asinh (double x)
double	acosh (double x)
double	atanh (double x)
void	rpn_sin (calc_number_t *c)
void	rpn_cos (calc_number_t *c)
void	rpn_tan (calc_number_t *c)
void	rpn_asin (calc_number_t *c)
void	rpn_acos (calc_number_t *c)
void	rpn_atan (calc_number_t *c)
void	rpn_sinh (calc_number_t *c)
void	rpn_cosh (calc_number_t *c)
void	rpn_tanh (calc_number_t *c)
void	rpn_asinh (calc_number_t *c)
void	rpn_acosh (calc_number_t *c)
void	rpn_atanh (calc_number_t *c)
void	rpn_int (calc_number_t *c)
void	rpn_frac (calc_number_t *c)
void	rpn_reci (calc_number_t *c)
void	rpn_fact (calc_number_t *c)
__int64	logic_dbl2int (calc_number_t *a)
double	logic_int2dbl (calc_number_t *a)
void	rpn_not (calc_number_t *c)
void	rpn_pi (calc_number_t *c)
void	rpn_2pi (calc_number_t *c)
void	rpn_sign (calc_number_t *c)
void	rpn_exp2 (calc_number_t *c)
void	rpn_exp3 (calc_number_t *c)
static __int64	myabs64 (__int64 number)
static unsigned __int64	sqrti (unsigned __int64 number)
void	rpn_sqrt (calc_number_t *c)
static __int64	cbrti (__int64 x)
void	rpn_cbrt (calc_number_t *c)
void	rpn_exp (calc_number_t *c)
void	rpn_exp10 (calc_number_t *c)
void	rpn_ln (calc_number_t *c)
void	rpn_log (calc_number_t *c)
static double	stat_sum (void)
static double	stat_sum2 (void)
void	rpn_ave (calc_number_t *c)
void	rpn_ave2 (calc_number_t *c)
void	rpn_sum (calc_number_t *c)
void	rpn_sum2 (calc_number_t *c)
static void	rpn_s_ex (calc_number_t *c, int pop_type)
void	rpn_s (calc_number_t *c)
void	rpn_s_m1 (calc_number_t *c)
void	rpn_dms2dec (calc_number_t *c)
void	rpn_dec2dms (calc_number_t *c)
void	rpn_zero (calc_number_t *c)
void	rpn_copy (calc_number_t *dst, calc_number_t *src)
int	rpn_is_zero (calc_number_t *c)
void	rpn_alloc (calc_number_t *c)
void	rpn_free (calc_number_t *c)

Macro Definition Documentation

NEXT

#define NEXT	(		n,
			i
	)		(((n) + (i)/(n)) >> 1)

Function Documentation

acosh()

double acosh

(

double

x

)

Definition at line 54 of file fun_ieee.c.

{
    // must be x>=1, if not return Nan (Not a Number)
    if(!(x>=1.0)) return sqrt(-1.0);
 
    // return only the positive result (as sqrt does).
    return log(x+sqrt(x*x-1.0));
}

Referenced by rpn_acosh().

apply_int_mask()

void apply_int_mask

(

calc_number_t *

r

)

Definition at line 26 of file fun_ieee.c.

{
    unsigned __int64 mask;
 
    switch (calc.size) {
    case IDC_RADIO_QWORD:
        mask = _UI64_MAX;
        break;
    case IDC_RADIO_DWORD:
        mask = ULONG_MAX;
        break;
    case IDC_RADIO_WORD:
        mask = USHRT_MAX;
        break;
    case IDC_RADIO_BYTE:
        mask = UCHAR_MAX;
        break;
    default:
        mask = (unsigned __int64)-1;
    }
    r->i &= mask;
}

asinh()

double asinh

(

double

x

)

Definition at line 49 of file fun_ieee.c.

{
    return log(x+sqrt(x*x+1));
}

Referenced by rpn_asinh().

atanh()

double atanh

(

double

x

)

Definition at line 63 of file fun_ieee.c.

{
    // must be x>-1, x<1, if not return Nan (Not a Number)
    if(!(x>-1.0 && x<1.0)) return sqrt(-1.0);
 
    return log((1.0+x)/(1.0-x))/2.0;
}

Referenced by rpn_atanh().

cbrti()

static __int64 cbrti ( __int64  x )

static

Definition at line 370 of file fun_ieee.c.

                                {
   __int64 s, y, b;
 
   s = 60;
   y = 0;
   while(s >= 0) {
      y = 2*y;
      b = (3*y*(y + 1) + 1) << s;
      s = s - 3;
      if (x >= b) {
         x = x - b;
         y = y + 1;
      }
   }
   return y;
}

Referenced by rpn_cbrt().

logic_dbl2int()

__int64 logic_dbl2int

(

calc_number_t *

a

)

Definition at line 261 of file fun_ieee.c.

{
    double   int_part;
    int      width;
 
    modf(a->f, &int_part);
    width = (int_part==0) ? 1 : (int)log10(fabs(int_part))+1;
    if (width > 63) {
        calc.is_nan = TRUE;
        return 0;
    }
    return (__int64)int_part;
}

Referenced by rpn_and_f(), rpn_not(), rpn_or_f(), rpn_xor_f(), and run_operator().

logic_int2dbl()

double logic_int2dbl

(

calc_number_t *

a

)

Definition at line 275 of file fun_ieee.c.

{
    return (double)a->i;
}

Referenced by run_operator().

myabs64()

static __int64 myabs64 ( __int64  number )

static

Definition at line 328 of file fun_ieee.c.

{
    return (number < 0) ? 0-number : number;
}

Referenced by sqrti().

rpn_2pi()

void rpn_2pi

(

calc_number_t *

c

)

Definition at line 295 of file fun_ieee.c.

{
    c->f = CALC_PI*2;
}

rpn_acos()

void rpn_acos

(

calc_number_t *

c

)

Definition at line 150 of file fun_ieee.c.

{
    c->f = validate_rad2angle(acos(c->f));
    if (_isnan(c->f))
        calc.is_nan = TRUE;
}

rpn_acosh()

void rpn_acosh

(

calc_number_t *

c

)

Definition at line 188 of file fun_ieee.c.

{
    c->f = acosh(c->f);
    if (_isnan(c->f))
        calc.is_nan = TRUE;
}

rpn_alloc()

void rpn_alloc

(

calc_number_t *

c

)

Definition at line 601 of file fun_ieee.c.

{
}

rpn_asin()

void rpn_asin

(

calc_number_t *

c

)

Definition at line 144 of file fun_ieee.c.

{
    c->f = validate_rad2angle(asin(c->f));
    if (_isnan(c->f))
        calc.is_nan = TRUE;
}

rpn_asinh()

void rpn_asinh

(

calc_number_t *

c

)

Definition at line 182 of file fun_ieee.c.

{
    c->f = asinh(c->f);
    if (_isnan(c->f))
        calc.is_nan = TRUE;
}

rpn_atan()

void rpn_atan

(

calc_number_t *

c

)

Definition at line 156 of file fun_ieee.c.

{
    c->f = validate_rad2angle(atan(c->f));
    if (_isnan(c->f))
        calc.is_nan = TRUE;
}

rpn_atanh()

void rpn_atanh

(

calc_number_t *

c

)

Definition at line 194 of file fun_ieee.c.

{
    c->f = atanh(c->f);
    if (_isnan(c->f))
        calc.is_nan = TRUE;
}

rpn_ave()

void rpn_ave

(

calc_number_t *

c

)

Definition at line 467 of file fun_ieee.c.

{
    double       ave = 0;
    int          n;
 
    ave = stat_sum();
    n = SendDlgItemMessage(calc.hStatWnd, IDC_LIST_STAT, LB_GETCOUNT, 0, 0);
 
    if (n)
        ave = ave / (double)n;
    if (calc.base == IDC_RADIO_DEC)
        c->f = ave;
    else
        c->i = (__int64)ave;
}

rpn_ave2()

void rpn_ave2

(

calc_number_t *

c

)

Definition at line 483 of file fun_ieee.c.

{
    double       ave = 0;
    int          n;
 
    ave = stat_sum2();
    n = SendDlgItemMessage(calc.hStatWnd, IDC_LIST_STAT, LB_GETCOUNT, 0, 0);
 
    if (n)
        ave = ave / (double)n;
    if (calc.base == IDC_RADIO_DEC)
        c->f = ave;
    else
        c->i = (__int64)ave;
}

rpn_cbrt()

void rpn_cbrt

(

calc_number_t *

c

)

Definition at line 387 of file fun_ieee.c.

{
    if (calc.base == IDC_RADIO_DEC)
#if defined(__GNUC__) && !defined(__REACTOS__)
        c->f = cbrt(c->f);
#else
        c->f = pow(c->f,1./3.);
#endif
    else {
        c->i = cbrti(c->i);
    }
}

rpn_copy()

void rpn_copy	(	calc_number_t *	dst,
		calc_number_t *	src
	)

Definition at line 591 of file fun_ieee.c.

{
    *dst = *src;
}

rpn_cos()

void rpn_cos

(

calc_number_t *

c

)

Definition at line 116 of file fun_ieee.c.

{
    double angle = validate_angle2rad(c);
 
    if (angle == CALC_PI_2 || angle == CALC_3_PI_2)
        c->f = 0;
    else
    if (angle == CALC_PI)
        c->f = -1;
    else
    if (angle == CALC_2_PI)
        c->f = 1;
    else
        c->f = cos(angle);
}

rpn_cosh()

void rpn_cosh

(

calc_number_t *

c

)

Definition at line 169 of file fun_ieee.c.

{
    c->f = cosh(c->f);
    if (_isnan(c->f))
        calc.is_nan = TRUE;
}

rpn_dec2dms()

void rpn_dec2dms

(

calc_number_t *

c

)

Definition at line 577 of file fun_ieee.c.

{
    double d, m, s;
 
    m = modf(c->f, &d) * 60;
    s = ceil(modf(m, &m) * 60);
    c->f = d + m/100. + s/10000.;
}

rpn_dms2dec()

void rpn_dms2dec

(

calc_number_t *

c

)

Definition at line 563 of file fun_ieee.c.

{
    double d, m, s;
 
    m = modf(c->f, &d) * 100;
    s = (modf(m, &m) * 100)+.5;
    modf(s, &s);
 
    m = m/60;
    s = s/3600;
 
    c->f = d + m + s;
}

rpn_exp()

void rpn_exp

(

calc_number_t *

c

)

Definition at line 400 of file fun_ieee.c.

{
    c->f = exp(c->f);
    if (_finite(c->f) == 0)
        calc.is_nan = TRUE;
}

rpn_exp10()

void rpn_exp10

(

calc_number_t *

c

)

Definition at line 407 of file fun_ieee.c.

{
    double int_part;
 
    modf(c->f, &int_part);
    if (fmod(int_part, 2.) == 0.)
        calc.is_nan = TRUE;
    else {
        c->f = pow(10., c->f);
        if (_finite(c->f) == 0)
            calc.is_nan = TRUE;
    }
}

rpn_exp2()

void rpn_exp2

(

calc_number_t *

c

)

Definition at line 308 of file fun_ieee.c.

{
    if (calc.base == IDC_RADIO_DEC) {
        c->f *= c->f;
        if (_finite(c->f) == 0)
            calc.is_nan = TRUE;
    } else
        c->i *= c->i;
}

rpn_exp3()

void rpn_exp3

(

calc_number_t *

c

)

Definition at line 318 of file fun_ieee.c.

{
    if (calc.base == IDC_RADIO_DEC) {
        c->f = pow(c->f, 3.);
        if (_finite(c->f) == 0)
            calc.is_nan = TRUE;
    } else
        c->i *= (c->i*c->i);
}

rpn_fact()

void rpn_fact

(

calc_number_t *

c

)

Definition at line 224 of file fun_ieee.c.

{
    double fact, mult, num;
 
    if (calc.base == IDC_RADIO_DEC)
        num = c->f;
    else
        num = (double)c->i;
    if (num > 1000) {
        calc.is_nan = TRUE;
        return;
    }
    if (num < 0) {
        calc.is_nan = TRUE;
        return;
    } else
    if (num == 0)
        fact = 1;
    else {
        rpn_int(c);
        fact = 1;
        mult = 2;
        while (mult <= num) {
            fact *= mult;
            mult++;
        }
        c->f = fact;
    }
    if (_finite(fact) == 0)
        calc.is_nan = TRUE;
    else
    if (calc.base == IDC_RADIO_DEC)
        c->f = fact;
    else
        c->i = (__int64)fact;
}

rpn_frac()

void rpn_frac

(

calc_number_t *

c

)

Definition at line 209 of file fun_ieee.c.

{
    double int_part;
 
    c->f = modf(calc.code.f, &int_part);
}

rpn_free()

void rpn_free

(

calc_number_t *

c

)

Definition at line 605 of file fun_ieee.c.

{
}

rpn_int()

void rpn_int

(

calc_number_t *

c

)

Definition at line 201 of file fun_ieee.c.

{
    double int_part;
 
    modf(calc.code.f, &int_part);
    c->f = int_part;
}

Referenced by rpn_fact().

rpn_is_zero()

int rpn_is_zero

(

calc_number_t *

c

)

Definition at line 596 of file fun_ieee.c.

{
    return (c->f == 0);
}

rpn_ln()

void rpn_ln

(

calc_number_t *

c

)

Definition at line 421 of file fun_ieee.c.

{
    if (c->f <= 0)
        calc.is_nan = TRUE;
    else
        c->f = log(c->f);
}

rpn_log()

void rpn_log

(

calc_number_t *

c

)

Definition at line 429 of file fun_ieee.c.

{
    if (c->f <= 0)
        calc.is_nan = TRUE;
    else
        c->f = log10(c->f);
}

rpn_not()

void rpn_not

(

calc_number_t *

c

)

Definition at line 280 of file fun_ieee.c.

{
    if (calc.base == IDC_RADIO_DEC) {
        calc_number_t n;
        n.i = logic_dbl2int(c);
        c->f = (long double)(~n.i);
    } else
        c->i = ~c->i;
}

rpn_pi()

void rpn_pi

(

calc_number_t *

c

)

Definition at line 290 of file fun_ieee.c.

{
    c->f = CALC_PI;
}

rpn_reci()

void rpn_reci

(

calc_number_t *

c

)

Definition at line 216 of file fun_ieee.c.

{
    if (c->f == 0)
        calc.is_nan = TRUE;
    else
        c->f = 1./c->f;
}

rpn_s()

void rpn_s

(

calc_number_t *

c

)

Definition at line 553 of file fun_ieee.c.

{
    rpn_s_ex(c, 0);
}

rpn_s_ex()

static void rpn_s_ex ( calc_number_t *  c, int  pop_type  )

static

Definition at line 519 of file fun_ieee.c.

{
    double       ave = 0;
    double       n = 0;
    double       dev = 0;
    double       num = 0;
    statistic_t *p = calc.stat;
 
    ave = stat_sum();
    n = (double)SendDlgItemMessage(calc.hStatWnd, IDC_LIST_STAT, LB_GETCOUNT, 0, 0);
 
    if (n == 0) {
        c->f = 0;
        return;
    }
    ave = ave / n;
 
    dev = 0;
    p = calc.stat;
    while (p != NULL) {
        if (p->base == IDC_RADIO_DEC)
            num = p->num.f;
        else
            num = (double)p->num.i;
        dev += pow(num-ave, 2.);
        p = (statistic_t *)(p->next);
    }
    dev = sqrt(dev/(pop_type ? n-1 : n));
    if (calc.base == IDC_RADIO_DEC)
        c->f = dev;
    else
        c->i = (__int64)dev;
}

Referenced by rpn_s(), and rpn_s_m1().

rpn_s_m1()

void rpn_s_m1

(

calc_number_t *

c

)

Definition at line 558 of file fun_ieee.c.

{
    rpn_s_ex(c, 1);
}

rpn_sign()

void rpn_sign

(

calc_number_t *

c

)

Definition at line 300 of file fun_ieee.c.

{
    if (calc.base == IDC_RADIO_DEC)
        c->f = 0-c->f;
    else
        c->i = 0-c->i;
}

rpn_sin()

void rpn_sin

(

calc_number_t *

c

)

Definition at line 101 of file fun_ieee.c.

{
    double angle = validate_angle2rad(c);
 
    if (angle == 0 || angle == CALC_PI)
        c->f = 0;
    else
    if (angle == CALC_3_PI_2)
        c->f = -1;
    else
    if (angle == CALC_2_PI)
        c->f = 1;
    else
        c->f = sin(angle);
}

rpn_sinh()

void rpn_sinh

(

calc_number_t *

c

)

Definition at line 163 of file fun_ieee.c.

{
    c->f = sinh(c->f);
    if (_isnan(c->f))
        calc.is_nan = TRUE;
}

rpn_sqrt()

void rpn_sqrt

(

calc_number_t *

c

)

Definition at line 358 of file fun_ieee.c.

{
    if (calc.base == IDC_RADIO_DEC) {
        if (c->f < 0)
            calc.is_nan = TRUE;
        else
            c->f = sqrt(c->f);
    } else {
        c->i = sqrti(c->i);
    }
}

rpn_sum()

void rpn_sum

(

calc_number_t *

c

)

Definition at line 499 of file fun_ieee.c.

{
    double sum = stat_sum();
 
    if (calc.base == IDC_RADIO_DEC)
        c->f = sum;
    else
        c->i = (__int64)sum;
}

rpn_sum2()

void rpn_sum2

(

calc_number_t *

c

)

Definition at line 509 of file fun_ieee.c.

{
    double sum = stat_sum2();
 
    if (calc.base == IDC_RADIO_DEC)
        c->f = sum;
    else
        c->i = (__int64)sum;
}

rpn_tan()

void rpn_tan

(

calc_number_t *

c

)

Definition at line 131 of file fun_ieee.c.

{
    double angle = validate_angle2rad(c);
 
    if (angle == CALC_PI_2 || angle == CALC_3_PI_2)
        calc.is_nan = TRUE;
    else
    if (angle == CALC_PI || angle == CALC_2_PI)
        c->f = 0;
    else
        c->f = tan(angle);
}

rpn_tanh()

void rpn_tanh

(

calc_number_t *

c

)

Definition at line 175 of file fun_ieee.c.

{
    c->f = tanh(c->f);
    if (_isnan(c->f))
        calc.is_nan = TRUE;
}

rpn_zero()

void rpn_zero

(

calc_number_t *

c

)

Definition at line 586 of file fun_ieee.c.

{
    c->f = 0;
}

sqrti()

static unsigned __int64 sqrti ( unsigned __int64  number )

static

Definition at line 333 of file fun_ieee.c.

{
/* modified form of Newton's method for approximating roots */
#define NEXT(n, i)  (((n) + (i)/(n)) >> 1)
    unsigned __int64 n, n1;
 
#ifdef __GNUC__
    if (number == 0xffffffffffffffffULL)
#else
    if (number == 0xffffffffffffffffUI64)
#endif
        return 0xffffffff;
 
    n  = 1;
    n1 = NEXT(n, number);
    while (myabs64(n1 - n) > 1) {
        n  = n1;
        n1 = NEXT(n, number);
    }
    while((n1*n1) > number)
        n1--;
    return n1;
#undef NEXT
}

Referenced by rpn_sqrt().

stat_sum()

static double stat_sum ( void  )

static

Definition at line 437 of file fun_ieee.c.

{
    double       sum = 0;
    statistic_t *p = calc.stat;
 
    while (p != NULL) {
        if (p->base == IDC_RADIO_DEC)
            sum += p->num.f;
        else
            sum += p->num.i;
        p = (statistic_t *)(p->next);
    }
    return sum;
}

Referenced by rpn_ave(), rpn_s_ex(), and rpn_sum().

stat_sum2()

static double stat_sum2 ( void  )

static

Definition at line 452 of file fun_ieee.c.

{
    double       sum = 0;
    statistic_t *p = calc.stat;
 
    while (p != NULL) {
        if (p->base == IDC_RADIO_DEC)
            sum += p->num.f * p->num.f;
        else
            sum += (double)p->num.i * (double)p->num.i;
        p = (statistic_t *)(p->next);
    }
    return sum;
}

Referenced by rpn_ave2(), and rpn_sum2().

validate_angle2rad()

static double validate_angle2rad ( calc_number_t *  c )

static

Definition at line 86 of file fun_ieee.c.

{
    switch (calc.degr) {
    case IDC_RADIO_DEG:
        c->f = c->f * (CALC_PI/180.0);
        break;
    case IDC_RADIO_RAD:
        break;
    case IDC_RADIO_GRAD:
        c->f = c->f * (CALC_PI/200.0);
        break;
    }
    return c->f;
}

Referenced by rpn_cos(), rpn_sin(), and rpn_tan().

validate_rad2angle()

static double validate_rad2angle ( double  a )

static

Definition at line 71 of file fun_ieee.c.

{
    switch (calc.degr) {
    case IDC_RADIO_DEG:
        a = a * (180.0/CALC_PI);
        break;
    case IDC_RADIO_RAD:
        break;
    case IDC_RADIO_GRAD:
        a = a * (200.0/CALC_PI);
        break;
    }
    return a;
}

Referenced by rpn_acos(), rpn_asin(), and rpn_atan().