paweld/fpc
1
0
Fork 0
mirror of https://gitlab.com/freepascal.org/fpc/source.git synced 2025-08-17 16:29:24 +02:00
Code Issues Packages Projects Releases Wiki Activity

* fixed power() in genmath.inc (code duplication from math.pp for **

Browse Source 
support!)
  * fixed power() in math.pp to give an error from 0^0
This commit is contained in:
Jonas Maebe 2004-12-05 16:43:57 +00:00
parent 80606b3569
commit 08da4e9278
2 changed files with 44 additions and 7 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

37
rtl/inc/genmath.inc
View File

@ -1198,7 +1198,30 @@ type
{$endif} {$endif}
{$ifndef FPC_SYSTEM_HAS_POWER} {$ifndef FPC_SYSTEM_HAS_POWER}
function intpower(base : real;const exponent: longint) : real;
var
i : longint;
begin
i:=abs(exponent);
intpower:=1.0;
while i>0 do
begin
while (i and 1)=0 do
begin
i:=i shr 1;
base:=sqr(base);
end;
dec(i);
intpower:=intpower*base;
end;
if exponent<0 then
intpower:=1.0/intpower;
end;
function power(bas,expo : real) : real; function power(bas,expo : real) : real;
var
sign: real;
begin begin
if bas=0.0 then if bas=0.0 then
begin begin
@ -1209,10 +1232,11 @@ type
end end
else if expo=0.0 then else if expo=0.0 then
power:=1 power:=1
else else if (abs(expo)<=high(longint)) and (frac(expo)=0.0) then
power := intpower(bas,trunc(expo))
{ bas < 0 is not allowed } { bas < 0 is not allowed }
if bas<0.0 then else if bas<0.0 then
handleerror(207) HandleError(207)
else else
power:=exp(ln(bas)*expo); power:=exp(ln(bas)*expo);
end; end;
@ -1307,7 +1331,12 @@ function fpc_int64_to_double(i : int64): double; compilerproc;
{ {
$Log$ $Log$
Revision 1.29 2004-11-21 15:35:23 peter Revision 1.30 2004-12-05 16:43:57 jonas
* fixed power() in genmath.inc (code duplication from math.pp for **
support!)
* fixed power() in math.pp to give an error from 0^0
Revision 1.29 2004/11/21 15:35:23 peter
* float routines all use internproc and compilerproc helpers * float routines all use internproc and compilerproc helpers
Revision 1.28 2004/11/20 15:49:21 jonas Revision 1.28 2004/11/20 15:49:21 jonas

14
rtl/objpas/math.pp
View File

@ -620,7 +620,10 @@ function power(base,exponent : float) : float;
begin begin
if Exponent=0.0 then if Exponent=0.0 then
result:=1.0 if base <> 0.0 then
result:=1.0
else
InvalidArgument
else if (base=0.0) and (exponent>0.0) then else if (base=0.0) and (exponent>0.0) then
result:=0.0 result:=0.0
else if (abs(exponent)<=maxint) and (frac(exponent)=0.0) then else if (abs(exponent)<=maxint) and (frac(exponent)=0.0) then
@ -1339,7 +1342,12 @@ end;
end. end.
{ {
$Log$ $Log$
Revision 1.24 2004-12-04 23:38:59 florian Revision 1.25 2004-12-05 16:43:57 jonas
* fixed power() in genmath.inc (code duplication from math.pp for **
support!)
* fixed power() in math.pp to give an error from 0^0
Revision 1.24 2004/12/04 23:38:59 florian
* fixed power(float,float) for negative exponents * fixed power(float,float) for negative exponents
Revision 1.23 2004/07/25 16:46:08 michael Revision 1.23 2004/07/25 16:46:08 michael
@ -1394,4 +1402,4 @@ end.
Revision 1.7 2002/09/07 16:01:22 peter Revision 1.7 2002/09/07 16:01:22 peter
* old logs removed and tabs fixed * old logs removed and tabs fixed
} }