Changeset 257 for Include/basic/function.sbp

Timestamp:

May 19, 2007, 3:29:33 PM (18 years ago)

Author:

イグトランス (egtra)

Message:

VersionTest追加、Log1p追加

File:

• Include/basic/function.sbp (modified) (10 diffs)

Include/basic/function.sbp

@@ -9 +9 @@
 Const _System_LOG2 = 0.6931471805599453094172321214581765680755
 Const _System_SQRT2 = 1.41421356237309504880168872421
+Const _System_Log_N = 7 As Long
 
 
@@ -21 +22 @@
         Exit Function
     End If
-    Dim pSrc As *QWord, pDest As *QWord
-    pSrc = VarPtr(x) As *QWord
-    pDest = VarPtr(ldexp) As *QWord
-
+    Dim pSrc = VarPtr(x) As *QWord
+    Dim pDest = VarPtr(ldexp) As *QWord
     n += (pSrc[0] >> 52) As DWord And &h7FF
-
     pDest[0] = n << 52 Or (pSrc[0] And &h800FFFFFFFFFFFFF)
 End Function
@@ -37 +35 @@
     End If
 
-    Dim pSrc As *QWord, pDest As *QWord
-    pSrc = VarPtr(x) As *QWord
-    pDest = VarPtr(frexp) As *QWord
+    Dim pSrc = VarPtr(x) As *QWord
+    Dim pDest = VarPtr(frexp) As *QWord
     n = ((pSrc[0] >> 52) As DWord And &h7FF) - 1022
     pDest[0] = (pSrc[0] And &h800FFFFFFFFFFFFF) Or &h3FE0000000000000
@@ -100 +97 @@
 End Function
 
-Function _System_GetInf(sign As BOOL) As Double
+Function _System_GetInf(sign As Boolean) As Double
     Dim s = 0 As QWord
     If sign Then s = 1 << 63
@@ -115 +112 @@
 End Function
 
-Function _System_GetInf(sign As BOOL) As Double
+Function _System_GetInf(sign As Boolean) As Double
     Dim s = 0 As DWord
     If sign Then s = (1 As DWord) << 31
@@ -240 +237 @@
 End Function
 
+Function Log1p(x As Double) As Double
+    If x < -1 Or IsNaN(x) Then
+        Log1p = _System_GetNaN()
+    ElseIf x = 0 Then
+        x = 0
+    ElseIf IsInf(x) Then
+        Log1p = x
+    Else
+        Log1p = _System_Log1p(x)
+    End If
+End Function
+
+Function _System_Log1p(x As Double) As Double
+    Dim s = 0 As Double
+    Dim i = 7 As Long
+    While i >= 1
+        Dim t = (i * x) As Double
+        s = t / (2 + t / (2 * i + 1 + s))
+        i--
+    Wend
+    Return x / (1 + s)
+End Function
+
 Function Sgn(number As Double) As Long
     Sgn = Math.Sign(number)
@@ -293 +313 @@
     Dim p As *DWord
     p = VarPtr(x) As *DWord
-    IsNaN = FALSE
+    IsNaN = False
     If (p[1] And &H7FF00000) = &H7FF00000 Then
         If (p[0] <> 0) Or ((p[1] And &HFFFFF) <> 0) Then
@@ -299 +319 @@
         End If
     End If
-
-'   IsNaN=FALSE
 End Function
 
@@ -307 +325 @@
     p = VarPtr(x) As *DWord
     p[1] And= &h7fffffff
-    nan = _System_GetInf(FALSE)
+    nan = _System_GetInf(False)
     IsInf = (memcmp(p As *Byte, VarPtr(nan), SizeOf (Double)) = 0)
 End Function
@@ -321 +339 @@
     p[1] And= &H7FF00000
     p[0] = 0
-    nan = _System_GetInf(/*x,*/ FALSE)
+    nan = _System_GetInf(/*x,*/ False)
     IsFinite = (memcmp(p As BytePtr, VarPtr(nan), SizeOf (Double)) = 0)
 End Function