vishapoberon/compiler
1
0
Fork 0
mirror of https://github.com/vishapoberon/compiler.git synced 2026-04-05 23:22:25 +00:00
Code Issues Projects Releases Packages Wiki Activity Actions

Math, MathL: add procedure fcmp

Browse source
This commit is contained in:
Carlos Une 2025-05-08 07:41:02 -03:00 • committed by Norayr Chilingarian
parent 94683df1d8
commit 28c327be28
2 changed files with 52 additions and 0 deletions
Download patch file Download diff file Expand all files Collapse all files

26
src/runtime/Math.Mod
View file

@ -741,6 +741,32 @@ BEGIN t := ABS(x);
RETURN arcsinh(x/sqrt(ONE - x * x)) RETURN arcsinh(x/sqrt(ONE - x * x))
END arctanh; END arctanh;
PROCEDURE fcmp* (x, y, epsilon: REAL): INTEGER;
(* fcmp: this procedure determines whether `x` and `y` are approximately equal
to a relative accuracy `epsilon`.
References:
The implementation is based on the GNU Scientific Library (GSL).
https://www.gnu.org/software/gsl/doc/html/math.html#approximate-comparison-of-floating-point-numbers
*)
VAR max, exponent0, delta, difference: REAL;
BEGIN
IF ABS(x) > ABS(y) THEN
max := x;
ELSE
max := y;
END;
exponent0 := exponent(max);
delta := 2.0*epsilon*power(2.0, exponent0);
difference := x - y;
IF difference > delta THEN
RETURN 1;
ELSIF difference < -delta THEN
RETURN -1;
ELSE
RETURN 0; (* approximately equal *)
END;
END fcmp;
PROCEDURE ToREAL(h: HUGEINT): REAL; PROCEDURE ToREAL(h: HUGEINT): REAL;
BEGIN RETURN SYSTEM.VAL(REAL, h) BEGIN RETURN SYSTEM.VAL(REAL, h)
END ToREAL; END ToREAL;

26
src/runtime/MathL.Mod
View file

@ -665,6 +665,32 @@ BEGIN t := ABS(x);
RETURN arcsinh(x/sqrt(ONE-x*x)) RETURN arcsinh(x/sqrt(ONE-x*x))
END arctanh; END arctanh;
PROCEDURE fcmp* (x, y, epsilon: LONGREAL): INTEGER;
(* fcmp: this procedure determines whether `x` and `y` are approximately equal
to a relative accuracy `epsilon`.
References:
The implementation is based on the GNU Scientific Library (GSL).
https://www.gnu.org/software/gsl/doc/html/math.html#approximate-comparison-of-floating-point-numbers
*)
VAR max, exponent0, delta, difference: LONGREAL;
BEGIN
IF ABS(x) > ABS(y) THEN
max := x;
ELSE
max := y;
END;
exponent0 := exponent(max);
delta := 2.0D0*epsilon*power(2.0D0, exponent0);
difference := x - y;
IF difference > delta THEN
RETURN 1;
ELSIF difference < -delta THEN
RETURN -1;
ELSE
RETURN 0; (* approximately equal *)
END;
END fcmp;
PROCEDURE ToLONGREAL(h: HUGEINT): LONGREAL; PROCEDURE ToLONGREAL(h: HUGEINT): LONGREAL;
BEGIN RETURN SYSTEM.VAL(LONGREAL, h) BEGIN RETURN SYSTEM.VAL(LONGREAL, h)
END ToLONGREAL; END ToLONGREAL;