|
a/LMathLib.mod |
|
b/LMathLib.mod |
|
... |
|
... |
31 |
(* eingefuegt *)
|
31 |
(* eingefuegt *)
|
32 |
(* 21.04.17, MRi: RationalApprox und gcd eingefuegt *)
|
32 |
(* 21.04.17, MRi: RationalApprox und gcd eingefuegt *)
|
33 |
(* 15.01.18, MRi: CalcPrimes eingefuegt *)
|
33 |
(* 15.01.18, MRi: CalcPrimes eingefuegt *)
|
34 |
(* 03.03.18, MRi: Log2 eingefuegt (von P. Moylan) *)
|
34 |
(* 03.03.18, MRi: Log2 eingefuegt (von P. Moylan) *)
|
35 |
(* 11.03.18, MRi: ceil & entier hierher verschoben *)
|
35 |
(* 11.03.18, MRi: ceil & entier hierher verschoben *)
|
|
|
36 |
(* 26.07.18, MRi: SmallR4 eingefuegt *)
|
36 |
(*------------------------------------------------------------------------*)
|
37 |
(*------------------------------------------------------------------------*)
|
37 |
(* Implementation : Michael Riedl *)
|
38 |
(* Implementation : Michael Riedl *)
|
38 |
(* Licence : GNU Lesser General Public License (LGPL) *)
|
39 |
(* Licence : GNU Lesser General Public License (LGPL) *)
|
39 |
(*------------------------------------------------------------------------*)
|
40 |
(*------------------------------------------------------------------------*)
|
40 |
|
41 |
|
|
... |
|
... |
664 |
EinsPlusEps:=1.0+eps;
|
665 |
EinsPlusEps:=1.0+eps;
|
665 |
END;
|
666 |
END;
|
666 |
RETURN 2.0*eps;
|
667 |
RETURN 2.0*eps;
|
667 |
END CalcMachEps;
|
668 |
END CalcMachEps;
|
668 |
|
669 |
|
|
|
670 |
PROCEDURE SmallR4() : LONGREAL;
|
|
|
671 |
|
|
|
672 |
VAR s,klein : REAL;
|
|
|
673 |
BEGIN
|
|
|
674 |
s:=1.0;
|
|
|
675 |
REPEAT
|
|
|
676 |
klein := s;
|
|
|
677 |
s := s / 16.0;
|
|
|
678 |
UNTIL (s = 0.0);
|
|
|
679 |
klein := 1.0*(klein / MachEpsR4);
|
|
|
680 |
RETURN VAL(LONGREAL,klein);
|
|
|
681 |
END SmallR4;
|
|
|
682 |
|
669 |
BEGIN
|
683 |
BEGIN
|
670 |
IF (MachEps # CalcMachEps()) THEN
|
684 |
IF (MachEps # CalcMachEps()) THEN
|
671 |
ErrOut("Fehler in LMathLib CalcMachEps() # MachEps, bitte pruefen ");
|
685 |
ErrOut("Fehler in LMathLib CalcMachEps() # MachEps, bitte pruefen ");
|
672 |
END;
|
686 |
END;
|
673 |
(*
|
687 |
(*
|