« previous next »
Pages: [1]   Go Down

Author Topic: intel decimal math library  (Read 6441 times)

0 Members and 1 Guest are viewing this topic.

Offline jack

  • Seasoned Forum Regular
  • Posts: 408
    • View Profile
intel decimal math library
« on: April 03, 2022, 11:31:04 am »
here's a small test of the Intel decimal math library https://www.intel.com/content/www/us/en/developer/articles/tool/intel-decimal-floating-point-math-library.html, if there's enough interest I may expand the function declarations
Code: QB64: [Select]
  1.  
  3.     Function __bid_strtod64~&& (s$, Byval ptr&)
  4.     Function add64~&& Alias "__bid64_add" (ByVal x~&&, Byval y~&&, Byval rnd_mode&, flag&)
  5.     Function sub64~&& Alias "__bid64_sub" (ByVal x~&&, Byval y~&&, Byval rnd_mode&, flag&)
  6.     Function mul64~&& Alias "__bid64_mul" (ByVal x~&&, Byval y&&, Byval rnd_mode&, flag&)
  7.     Function div64~&& Alias "__bid64_div" (ByVal x~&&, Byval y~&&, Byval rnd_mode&, flag&)
  8.  
  9.     Function is_equal64& Alias "__bid64_quiet_equal" (ByVal x~&&, Byval y~&&, flag&)
  10.     Function is_less64& Alias "__bid64_quiet_less" (ByVal x~&&, Byval y~&&, flag&)
  11.     Function is_greater64& Alias "__bid64_quiet_greater" (ByVal x~&&, Byval y~&&, flag&)
  12.     Function is_less_equal64& Alias "__bid64_quiet_less_equal" (ByVal x~&&, Byval y~&&, flag&)
  13.     Function is_greater_equal64& Alias "__bid64_quiet_greater_equal" (ByVal x~&&, Byval y~&&, flag&)
  14.     Function is_unequal64& Alias "__bid64_quiet_not_equal" (ByVal x~&&, Byval y~&&, flag&)
  15.     Function isZero64& Alias "__bid64_isZero" (ByVal x~&&)
  16.     Sub __bid64_to_string (ps$, Byval x~&&, pf&)
  17.     Function copy64~&& Alias "__bid64_copy" (ByVal x~&&)
  18.     Function negate64~&& Alias "__bid64_negate" (ByVal x~&&)
  19.     Function abs64~&& Alias "__bid64_abs" (ByVal x~&&)
  20.     Function sqrt64~&& Alias "__bid64_sqrt" (ByVal x~&&, Byval rnd_mode&, flag&)
  21.     Function cbrt64~&& Alias "__bid64_cbrt" (ByVal x~&&, Byval rnd_mode&, flag&)
  22.     Function sin64~&& Alias "__bid64_sin" (ByVal x~&&, Byval rnd_mode&, flag&)
  23.     Function cos64~&& Alias "__bid64_cos" (ByVal x~&&, Byval rnd_mode&, flag&)
  24.     Function tan64~&& Alias "__bid64_tan" (ByVal x~&&, Byval rnd_mode&, flag&)
  25.     Function asin64~&& Alias "__bid64_asin" (ByVal x~&&, Byval rnd_mode&, flag&)
  26.     Function acos64~&& Alias "__bid64_acos" (ByVal x~&&, Byval rnd_mode&, flag&)
  27.     Function atan64~&& Alias "__bid64_atan" (ByVal x~&&, Byval rnd_mode&, flag&)
  28.     Function sinh64~&& Alias "__bid64_sinh" (ByVal x~&&, Byval rnd_mode&, flag&)
  29.     Function cosh64~&& Alias "__bid64_cosh" (ByVal x~&&, Byval rnd_mode&, flag&)
  30.     Function tanh64~&& Alias "__bid64_tanh" (ByVal x~&&, Byval rnd_mode&, flag&)
  31.     Function asinh64~&& Alias "__bid64_asinh" (ByVal x~&&, Byval rnd_mode&, flag&)
  32.     Function acosh64~&& Alias "__bid64_acosh" (ByVal x~&&, Byval rnd_mode&, flag&)
  33.     Function atanh64~&& Alias "__bid64_atanh" (ByVal x~&&, Byval rnd_mode&, flag&)
  34.     Function atan264~&& Alias "__bid64_atan2" (ByVal x~&&, Byval y~&&, Byval rnd_mode&, flag&)
  35.     Function hypot64~&& Alias "__bid64_hypot" (ByVal x~&&, Byval y~&&, Byval rnd_mode&, flag&)
  36.     Function log1064~&& Alias "__bid64_log10" (ByVal x~&&, Byval rnd_mode&, flag&)
  37.     Function exp1064~&& Alias "__bid64_exp10" (ByVal x~&&, Byval rnd_mode&, flag&)
  38.     Function log64~&& Alias "__bid64_log" (ByVal x~&&, Byval rnd_mode&, flag&)
  39.     Function exp64~&& Alias "__bid64_exp" (ByVal x~&&, Byval rnd_mode&, flag&)
  40.     Function pow64~&& Alias "__bid64_pow" (ByVal x~&&, Byval y~&&, Byval rnd_mode&, flag&)
  41.     Function log264~&& Alias "__bid64_log2" (ByVal x~&&, Byval rnd_mode&, flag&)
  42.     Function exp264~&& Alias "__bid64_exp2" (ByVal x~&&, Byval rnd_mode&, flag&)
  43.     Function erf64~&& Alias "__bid64_erf" (ByVal x~&&, Byval rnd_mode&, flag&)
  44.     Function erfc64~&& Alias "__bid64_erfc" (ByVal x~&&, Byval rnd_mode&, flag&)
  45.     Function tgamma64~&& Alias "__bid64_tgamma" (ByVal x~&&, Byval rnd_mode&, flag&)
  46.     Function lgamma64~&& Alias "__bid64_lgamma" (ByVal x~&&, Byval rnd_mode&, flag&)
  47.     Function floor64l& Alias "__bid64_to_int32_floor" (ByVal x~&&, flag&)
  48.     Function ceil64l& Alias "__bid64_to_int32_ceil" (ByVal x~&&, flag&)
  49.     Function floor64ul~& Alias "__bid64_to_uint32_floor" (ByVal x~&&, flag&)
  50.     Function ceil64ul~& Alias "__bid64_to_uint32_ceil" (ByVal x~&&, flag&)
  51.     Function floor64ll~&& Alias "__bid64_to_int64_floor" (ByVal x~&&, flag&)
  52.     Function ceil64ll~&& Alias "__bid64_to_int64_ceil" (ByVal x~&&, flag&)
  53.     Function floor64ull~&& Alias "__bid64_to_uint64_floor" (ByVal x~&&, flag&)
  54.     Function ceil64ull~&& Alias "__bid64_to_uint64_ceil" (ByVal x~&&, flag&)
  55.     Function tosingle64! Alias "__bid64_to_binary32" (ByVal x~&&, Byval rnd_mode&, flag&)
  56.     Function todouble# Alias "__bid64_to_binary64" (ByVal x~&&, Byval rnd_mode&, flag&)
  57.     Function singleto64~&& Alias "__binary32_to_bid64" (ByVal x As Single, Byval rnd_mode&, flag&)
  58.     Function doubleto64~&& Alias "__binary64_to_bid64" (ByVal x As Double, Byval rnd_mode&, flag&)
  60.  
  62. Dim As Long flag, rm
  63.  
  64. rm = 0
  65. i = strto64("-1")
  66. y = strto64("1")
  67. z = strto64(".1")
  68. While is_less_equal64(i, y, flag)
  69.     Print tostring(i)
  70.     i = add64(i, z, rm, flag)
  72.  
  73. Function strto64~&& (s As String)
  74.     Dim As String sn
  75.     sn = s + Chr$(0)
  76.     strto64~&& = __bid_strtod64(sn, 0)
  78.  
  80.     Dim As Long ex, ex1, i, ln
  81.     Dim As String s, sd, sdl, sdr, se, sign
  82.     Dim As _Unsigned _Integer64 one, x
  83.     Dim As Long flag, rm
  84.     If isZero64(z) Then
  85.         tostring = " 0"
  86.         Exit Function
  87.     End If
  88.     rm = 0
  89.     x = z
  90.     one = __bid_strtod64("1.00000000000000000000", 0)
  91.     s = Space$(256)
  92.     __bid64_to_string s, mul64(x, one, 0, flag), flag
  93.     sd = _Trim$(s)
  94.     ln = Len(sd)
  95.     i = InStr(sd, "E")
  96.     sdl = Left$(sd, i - 1)
  97.     ex = Val(Mid$(sd, i + 1))
  98.     ln = Len(sdl)
  99.     If ex = 0 Then
  100.         sd = sdl
  101.     ElseIf ex > 0 Or ex < 0 Then
  102.         ex1 = ln + ex - 2
  103.         If ex1 >= 0 And ex1 < (ln - 1) Then
  104.             sd = Left$(sdl, ex1 + 2) + "." + Mid$(sdl, ex1 + 3)
  105.         ElseIf ex1 < 0 And ex1 > -5 Then
  106.             sign = Left$(sd, 1)
  107.             sd = sign + "0." + String$(Abs(ex1) - 1, "0") + Mid$(sdl, 2)
  108.         Else
  109.             se = _Trim$(Str$(ex1))
  110.             If ex1 > 0 Then
  111.                 se = "+" + se
  112.             End If
  113.             sd = sdl
  114.             sdl = Left$(sd, 2)
  115.             sdr = Mid$(sd, 3)
  116.             sd = sdl + "." + sdr + "E" + se
  117.         End If
  118.     End If
  119.     If Left$(sd, 1) = "+" Then
  120.         sd = " " + Mid$(sd, 2)
  121.     End If
  122.     tostring = sd
  124.  
here's the 64-bit dll
* libbid.zip (Filesize: 3.1 MB, Downloads: 373)
« Last Edit: April 03, 2022, 01:18:31 pm by jack »
Logged

Offline jack

  • Seasoned Forum Regular
  • Posts: 408
    • View Profile
Re: intel decimal math library
« Reply #1 on: April 03, 2022, 09:33:25 pm »
the Intel math library has 3 types decimal32, decimal64 and decimal128
the decimal64 is easy to use in QB64 but for the decimal128 I had make a wrapper
Code: QB64: [Select]
  1. Type BID_UINT128
  5.  
  7.     Sub binary64_to_bid128 (result As BID_UINT128, x As Double, Byval rnd_mode~&, pfpsf~&)
  8.     Sub binary80_to_bid128 (result As BID_UINT128, x As _Float, Byval rnd_mode~&, pfpsf~&)
  9.     Sub bid128_add (result As BID_UINT128, x As BID_UINT128, y As BID_UINT128, Byval rnd_mode~&, pfpsf~&)
  10.     Sub bid128_sub (result As BID_UINT128, x As BID_UINT128, y As BID_UINT128, Byval rnd_mode~&, pfpsf~&)
  11.     Sub bid128_mul (result As BID_UINT128, x As BID_UINT128, y As BID_UINT128, Byval rnd_mode~&, pfpsf~&)
  12.     Sub bid128_div (result As BID_UINT128, x As BID_UINT128, y As BID_UINT128, Byval rnd_mode~&, pfpsf~&)
  13.     Sub bid128_sqrt (result As BID_UINT128, x As BID_UINT128, Byval rnd_mode~&, pfpsf~&)
  14.     Sub bid128_cbrt (result As BID_UINT128, x As BID_UINT128, Byval rnd_mode~&, pfpsf~&)
  15.     Sub bid128_exp (result As BID_UINT128, x As BID_UINT128, Byval rnd_mode~&, pfpsf~&)
  16.     Sub bid128_log (result As BID_UINT128, x As BID_UINT128, Byval rnd_mode~&, pfpsf~&)
  17.     Sub bid128_pow (result As BID_UINT128, x As BID_UINT128, y As BID_UINT128, Byval rnd_mode~&, pfpsf~&)
  18.     Sub bid128_atan2 (result As BID_UINT128, x As BID_UINT128, y As BID_UINT128, Byval rnd_mode~&, pfpsf~&)
  19.     Sub bid128_fmod (result As BID_UINT128, x As BID_UINT128, y As BID_UINT128, pfpsf~&)
  20.     Sub bid128_modf (result As BID_UINT128, x As BID_UINT128, y As BID_UINT128, pfpsf~&)
  21.     Sub bid128_hypot (result As BID_UINT128, x As BID_UINT128, y As BID_UINT128, Byval rnd_mode~&, pfpsf~&)
  22.     Sub bid128_sin (result As BID_UINT128, x As BID_UINT128, Byval rnd_mode~&, pfpsf~&)
  23.     Sub bid128_cos (result As BID_UINT128, x As BID_UINT128, Byval rnd_mode~&, pfpsf~&)
  24.     Sub bid128_tan (result As BID_UINT128, x As BID_UINT128, Byval rnd_mode~&, pfpsf~&)
  25.     Sub bid128_asin (result As BID_UINT128, x As BID_UINT128, Byval rnd_mode~&, pfpsf~&)
  26.     Sub bid128_acos (result As BID_UINT128, x As BID_UINT128, Byval rnd_mode~&, pfpsf~&)
  27.     Sub bid128_atan (result As BID_UINT128, x As BID_UINT128, Byval rnd_mode~&, pfpsf~&)
  28.     Sub bid128_sinh (result As BID_UINT128, x As BID_UINT128, Byval rnd_mode~&, pfpsf~&)
  29.     Sub bid128_cosh (result As BID_UINT128, x As BID_UINT128, Byval rnd_mode~&, pfpsf~&)
  30.     Sub bid128_tanh (result As BID_UINT128, x As BID_UINT128, Byval rnd_mode~&, pfpsf~&)
  31.     Sub bid128_asinh (result As BID_UINT128, x As BID_UINT128, Byval rnd_mode~&, pfpsf~&)
  32.     Sub bid128_acosh (result As BID_UINT128, x As BID_UINT128, Byval rnd_mode~&, pfpsf~&)
  33.     Sub bid128_atanh (result As BID_UINT128, x As BID_UINT128, Byval rnd_mode~&, pfpsf~&)
  34.     Sub bid128_log1p (result As BID_UINT128, x As BID_UINT128, Byval rnd_mode~&, pfpsf~&)
  35.     Sub bid128_expm1 (result As BID_UINT128, x As BID_UINT128, Byval rnd_mode~&, pfpsf~&)
  36.     Sub bid128_log10 (result As BID_UINT128, x As BID_UINT128, Byval rnd_mode~&, pfpsf~&)
  37.     Sub bid128_log2 (result As BID_UINT128, x As BID_UINT128, Byval rnd_mode~&, pfpsf~&)
  38.     Sub bid128_exp2 (result As BID_UINT128, x As BID_UINT128, Byval rnd_mode~&, pfpsf~&)
  39.     Sub bid128_exp10 (result As BID_UINT128, x As BID_UINT128, Byval rnd_mode~&, pfpsf~&)
  40.     Sub bid128_erf (result As BID_UINT128, x As BID_UINT128, Byval rnd_mode~&, pfpsf~&)
  41.     Sub bid128_erfc (result As BID_UINT128, x As BID_UINT128, Byval rnd_mode~&, pfpsf~&)
  42.     Sub bid128_tgamma (result As BID_UINT128, x As BID_UINT128, Byval rnd_mode~&, pfpsf~&)
  43.     Sub bid128_lgamma (result As BID_UINT128, x As BID_UINT128, Byval rnd_mode~&, pfpsf~&)
  44.     Function bid128_quiet_equal& (x As BID_UINT128, y As BID_UINT128, pfpsf~&)
  45.     Function bid128_quiet_greater& (x As BID_UINT128, y As BID_UINT128, pfpsf~&)
  46.     Function bid128_quiet_greater_equal& (x As BID_UINT128, y As BID_UINT128, pfpsf~&)
  47.     Function bid128_quiet_greater_unordered& (x As BID_UINT128, y As BID_UINT128, pfpsf~&)
  48.     Function bid128_quiet_less& (x As BID_UINT128, y As BID_UINT128, pfpsf~&)
  49.     Function bid128_quiet_less_equal& (x As BID_UINT128, y As BID_UINT128, pfpsf~&)
  50.     Function bid128_quiet_less_unordered& (x As BID_UINT128, y As BID_UINT128, pfpsf~&)
  51.     Function bid128_quiet_not_equal& (x As BID_UINT128, y As BID_UINT128, pfpsf~&)
  52.     Function bid128_quiet_not_greater& (x As BID_UINT128, y As BID_UINT128, pfpsf~&)
  53.     Function bid128_quiet_not_less& (x As BID_UINT128, y As BID_UINT128, pfpsf~&)
  54.     Function bid128_quiet_ordered& (x As BID_UINT128, y As BID_UINT128, pfpsf~&)
  55.     Function bid128_quiet_unordered& (x As BID_UINT128, y As BID_UINT128, pfpsf~&)
  56.     Function bid128_signaling_greater& (x As BID_UINT128, y As BID_UINT128, pfpsf~&)
  57.     Function bid128_signaling_greater_equal& (x As BID_UINT128, y As BID_UINT128, pfpsf~&)
  58.     Function bid128_signaling_greater_unordered& (x As BID_UINT128, y As BID_UINT128, pfpsf~&)
  59.     Function bid128_signaling_less& (x As BID_UINT128, y As BID_UINT128, pfpsf~&)
  60.     Function bid128_signaling_less_equal& (x As BID_UINT128, y As BID_UINT128, pfpsf~&)
  61.     Function bid128_signaling_less_unordered& (x As BID_UINT128, y As BID_UINT128, pfpsf~&)
  62.     Function bid128_signaling_not_greater& (x As BID_UINT128, y As BID_UINT128, pfpsf~&)
  63.     Function bid128_signaling_not_less& (x As BID_UINT128, y As BID_UINT128, pfpsf~&)
  64.  
  65.     Sub bid128_copy (result As BID_UINT128, x As BID_UINT128)
  66.     Sub bid128_negate (result As BID_UINT128, x As BID_UINT128)
  67.     Sub bid128_abs (result As BID_UINT128, x As BID_UINT128)
  68.     Sub bid128_to_string (S As String, x As BID_UINT128, pfpsf~&)
  69.     Function bid128_to_binary32! (x As BID_UINT128, Byval rnd_mode~&, pfpsf~&)
  70.     Function bid128_to_binary64# (x As BID_UINT128, Byval rnd_mode~&, pfpsf~&)
  71.     Function bid128_to_binary80## (x As BID_UINT128, Byval rnd_mode~&, pfpsf~&)
  72.     Sub bid_strtod128 (result As BID_UINT128, ps_in As String, Byval endptr&)
  73.  
  74.     Function bid128_to_int8_rnint%% (x As BID_UINT128, pfpsf~&)
  75.     Function bid128_to_int8_xrnint%% (x As BID_UINT128, pfpsf~&)
  76.     Function bid128_to_int8_rninta%% (x As BID_UINT128, pfpsf~&)
  77.     Function bid128_to_int8_xrninta%% (x As BID_UINT128, pfpsf~&)
  78.     Function bid128_to_int8_int%% (x As BID_UINT128, pfpsf~&)
  79.     Function bid128_to_int8_xint%% (x As BID_UINT128, pfpsf~&)
  80.     Function bid128_to_int8_floor%% (x As BID_UINT128, pfpsf~&)
  81.     Function bid128_to_int8_xfloor%% (x As BID_UINT128, pfpsf~&)
  82.     Function bid128_to_int8_ceil%% (x As BID_UINT128, pfpsf~&)
  83.     Function bid128_to_int8_xceil%% (x As BID_UINT128, pfpsf~&)
  84.  
  85.     Function bid128_to_int16_rnint% (x As BID_UINT128, pfpsf~&)
  86.     Function bid128_to_int16_xrnint% (x As BID_UINT128, pfpsf~&)
  87.     Function bid128_to_int16_rninta% (x As BID_UINT128, pfpsf~&)
  88.     Function bid128_to_int16_xrninta% (x As BID_UINT128, pfpsf~&)
  89.     Function bid128_to_int16_int% (x As BID_UINT128, pfpsf~&)
  90.     Function bid128_to_int16_xint% (x As BID_UINT128, pfpsf~&)
  91.     Function bid128_to_int16_floor% (x As BID_UINT128, pfpsf~&)
  92.     Function bid128_to_int16_xfloor% (x As BID_UINT128, pfpsf~&)
  93.     Function bid128_to_int16_ceil% (x As BID_UINT128, pfpsf~&)
  94.     Function bid128_to_int16_xceil% (x As BID_UINT128, pfpsf~&)
  95.  
  96.     Function bid128_to_uint8_rnint~%% (x As BID_UINT128, pfpsf~&)
  97.     Function bid128_to_uint8_xrnint~%% (x As BID_UINT128, pfpsf~&)
  98.     Function bid128_to_uint8_rninta~%% (x As BID_UINT128, pfpsf~&)
  99.     Function bid128_to_uint8_xrninta~%% (x As BID_UINT128, pfpsf~&)
  100.     Function bid128_to_uint8_int~%% (x As BID_UINT128, pfpsf~&)
  101.     Function bid128_to_uint8_xint~%% (x As BID_UINT128, pfpsf~&)
  102.     Function bid128_to_uint8_floor~%% (x As BID_UINT128, pfpsf~&)
  103.     Function bid128_to_uint8_xfloor~%% (x As BID_UINT128, pfpsf~&)
  104.     Function bid128_to_uint8_ceil~%% (x As BID_UINT128, pfpsf~&)
  105.     Function bid128_to_uint8_xceil~%% (x As BID_UINT128, pfpsf~&)
  106.  
  107.     Function bid128_to_uint16_rnint~% (x As BID_UINT128, pfpsf~&)
  108.     Function bid128_to_uint16_xrnint~% (x As BID_UINT128, pfpsf~&)
  109.     Function bid128_to_uint16_rninta~% (x As BID_UINT128, pfpsf~&)
  110.     Function bid128_to_uint16_xrninta~% (x As BID_UINT128, pfpsf~&)
  111.     Function bid128_to_uint16_int~% (x As BID_UINT128, pfpsf~&)
  112.     Function bid128_to_uint16_xint~% (x As BID_UINT128, pfpsf~&)
  113.     Function bid128_to_uint16_floor~% (x As BID_UINT128, pfpsf~&)
  114.     Function bid128_to_uint16_xfloor~% (x As BID_UINT128, pfpsf~&)
  115.     Function bid128_to_uint16_ceil~% (x As BID_UINT128, pfpsf~&)
  116.     Function bid128_to_uint16_xceil~% (x As BID_UINT128, pfpsf~&)
  117.  
  118.     Function bid128_to_int32_rnint& (x As BID_UINT128, pfpsf~&)
  119.     Function bid128_to_int32_xrnint& (x As BID_UINT128, pfpsf~&)
  120.     Function bid128_to_int32_rninta& (x As BID_UINT128, pfpsf~&)
  121.     Function bid128_to_int32_xrninta& (x As BID_UINT128, pfpsf~&)
  122.     Function bid128_to_int32_int& (x As BID_UINT128, pfpsf~&)
  123.     Function bid128_to_int32_xint& (x As BID_UINT128, pfpsf~&)
  124.     Function bid128_to_int32_floor& (x As BID_UINT128, pfpsf~&)
  125.     Function bid128_to_int32_xfloor& (x As BID_UINT128, pfpsf~&)
  126.     Function bid128_to_int32_ceil& (x As BID_UINT128, pfpsf~&)
  127.     Function bid128_to_int32_xceil& (x As BID_UINT128, pfpsf~&)
  128.  
  129.     Function bid128_to_uint32_rnint~& (x As BID_UINT128, pfpsf~&)
  130.     Function bid128_to_uint32_xrnint~& (x As BID_UINT128, pfpsf~&)
  131.     Function bid128_to_uint32_rninta~& (x As BID_UINT128, pfpsf~&)
  132.     Function bid128_to_uint32_xrninta~& (x As BID_UINT128, pfpsf~&)
  133.     Function bid128_to_uint32_int~& (x As BID_UINT128, pfpsf~&)
  134.     Function bid128_to_uint32_xint~& (x As BID_UINT128, pfpsf~&)
  135.     Function bid128_to_uint32_floor~& (x As BID_UINT128, pfpsf~&)
  136.     Function bid128_to_uint32_xfloor~& (x As BID_UINT128, pfpsf~&)
  137.     Function bid128_to_uint32_ceil~& (x As BID_UINT128, pfpsf~&)
  138.     Function bid128_to_uint32_xceil~& (x As BID_UINT128, pfpsf~&)
  139.  
  140.     Function bid128_to_int64_rnint&& (x As BID_UINT128, pfpsf~&)
  141.     Function bid128_to_int64_xrnint&& (x As BID_UINT128, pfpsf~&)
  142.     Function bid128_to_int64_rninta&& (x As BID_UINT128, pfpsf~&)
  143.     Function bid128_to_int64_xrninta&& (x As BID_UINT128, pfpsf~&)
  144.     Function bid128_to_int64_int&& (x As BID_UINT128, pfpsf~&)
  145.     Function bid128_to_int64_xint&& (x As BID_UINT128, pfpsf~&)
  146.     Function bid128_to_int64_floor&& (x As BID_UINT128, pfpsf~&)
  147.     Function bid128_to_int64_xfloor&& (x As BID_UINT128, pfpsf~&)
  148.     Function bid128_to_int64_ceil&& (x As BID_UINT128, pfpsf~&)
  149.     Function bid128_to_int64_xceil&& (x As BID_UINT128, pfpsf~&)
  150.  
  151.     Function bid128_to_uint64_rnint~&& (x As BID_UINT128, pfpsf~&)
  152.     Function bid128_to_uint64_xrnint~&& (x As BID_UINT128, pfpsf~&)
  153.     Function bid128_to_uint64_rninta~&& (x As BID_UINT128, pfpsf~&)
  154.     Function bid128_to_uint64_xrninta~&& (x As BID_UINT128, pfpsf~&)
  155.     Function bid128_to_uint64_int~&& (x As BID_UINT128, pfpsf~&)
  156.     Function bid128_to_uint64_xint~&& (x As BID_UINT128, pfpsf~&)
  157.     Function bid128_to_uint64_floor~&& (x As BID_UINT128, pfpsf~&)
  158.     Function bid128_to_uint64_xfloor~&& (x As BID_UINT128, pfpsf~&)
  159.     Function bid128_to_uint64_ceil~&& (x As BID_UINT128, pfpsf~&)
  160.     Function bid128_to_uint64_xceil~&& (x As BID_UINT128, pfpsf~&)
  161.  
  162.     Sub bid128_round_integral_exact (result As BID_UINT128, x As BID_UINT128, Byval rnd_mode~&, pfpsf~&)
  163.     Sub bid128_round_integral_nearest_even (result As BID_UINT128, x As BID_UINT128, pfpsf~&)
  164.     Sub bid128_round_integral_negative (result As BID_UINT128, x As BID_UINT128, pfpsf~&)
  165.     Sub bid128_round_integral_positive (result As BID_UINT128, x As BID_UINT128, pfpsf~&)
  166.     Sub bid128_round_integral_zero (result As BID_UINT128, x As BID_UINT128, pfpsf~&)
  167.     Sub bid128_round_integral_nearest_away (result As BID_UINT128, x As BID_UINT128, pfpsf~&)
  168.     Sub bid128_nextup (result As BID_UINT128, x As BID_UINT128, pfpsf~&)
  169.     Sub bid128_nextdown (result As BID_UINT128, x As BID_UINT128, pfpsf~&)
  170.     Sub bid128_nextafter (result As BID_UINT128, x As BID_UINT128, y As BID_UINT128, pfpsf~&)
  171.     Sub bid128_minnum (result As BID_UINT128, x As BID_UINT128, y As BID_UINT128, pfpsf~&)
  172.     Sub bid128_minnum_mag (result As BID_UINT128, x As BID_UINT128, y As BID_UINT128, pfpsf~&)
  173.     Sub bid128_maxnum (result As BID_UINT128, x As BID_UINT128, y As BID_UINT128, pfpsf~&)
  174.     Sub bid128_maxnum_mag (result As BID_UINT128, x As BID_UINT128, y As BID_UINT128, pfpsf~&)
  175.     Sub bid128_from_int32 (result As BID_UINT128, Byval x&)
  176.     Sub bid128_from_uint32 (result As BID_UINT128, Byval x~&)
  177.     Sub bid128_from_int64 (result As BID_UINT128, Byval x&&)
  178.     Sub bid128_from_uint64 (result As BID_UINT128, Byval x~&&)
  179.     Function bid128_isSigned& (x As BID_UINT128)
  180.     Function bid128_isNormal& (x As BID_UINT128)
  181.     Function bid128_isSubnormal& (x As BID_UINT128)
  182.     Function bid128_isFinite& (x As BID_UINT128)
  183.     Function bid128_isZero& (x As BID_UINT128)
  184.     Function bid128_isInf& (x As BID_UINT128)
  185.     Function bid128_isSignaling& (x As BID_UINT128)
  186.     Function bid128_isCanonical& (x As BID_UINT128)
  187.     Function bid128_isNaN& (x As BID_UINT128)
  188.  
  190.  
  193.  
  194. Dim As BID_UINT128 i, y, z, x
  195. Dim As Long flag, rm
  196.  
  197. rm = 0
  198. bid128_from_int32 i, -1 'strto128 i, "1"
  199. bid128_from_int32 y, 1 'strto128 y, "20"
  200. strto128 z, ".1"
  201. While bid128_quiet_less_equal(i, y, flag)
  202.     'bid128_sin x, i, rm, flag
  203.     'bid128_asin x, x, rm, flag
  204.     'bid128_log10 x, i, rm, flag
  205.     'Print d128tostring$(x)
  206.     Print d128tostring$(i)
  207.     bid128_add i, i, z, rm, flag
  209.  
  210. Sub strto128 (result As BID_UINT128, s As String)
  211.     Dim As String sn
  212.     sn = s + Chr$(0)
  213.     bid_strtod128 result, sn, 0
  215.  
  216. Function d128tostring$ (z As BID_UINT128)
  217.     Dim As Long ex, ex1, i, ln
  218.     Dim As String s, sd, sdl, sdr, se, sign
  219.     Dim As BID_UINT128 one, x
  220.     Dim As Long flag, rm
  221.     If bid128_isZero(z) Then
  222.         d128tostring = " 0"
  223.         Exit Function
  224.     End If
  225.     rm = 0
  226.     x = z
  227.     bid128_copy x, z
  228.     bid_strtod128 one, "1.00000000000000000000000000000000000", 0
  229.     s = Space$(256)
  230.     bid128_mul x, x, one, 0, flag
  231.     bid128_to_string s, x, flag
  232.     sd = _Trim$(s)
  233.     ln = Len(sd)
  234.     i = InStr(sd, "E")
  235.     sdl = Left$(sd, i - 1)
  236.     ex = Val(Mid$(sd, i + 1))
  237.     ln = Len(sdl)
  238.     If ex = 0 Then
  239.         sd = sdl
  240.     ElseIf ex > 0 Or ex < 0 Then
  241.         ex1 = ln + ex - 2
  242.         If ex1 >= 0 And ex1 < (ln - 1) Then
  243.             sd = Left$(sdl, ex1 + 2) + "." + Mid$(sdl, ex1 + 3)
  244.         ElseIf ex1 < 0 And ex1 > -5 Then
  245.             sign = Left$(sd, 1)
  246.             sd = sign + "0." + String$(Abs(ex1) - 1, "0") + Mid$(sdl, 2)
  247.         Else
  248.             se = _Trim$(Str$(ex1))
  249.             If ex1 > 0 Then
  250.                 se = "+" + se
  251.             End If
  252.             sd = sdl
  253.             sdl = Left$(sd, 2)
  254.             sdr = Mid$(sd, 3)
  255.             sd = sdl + "." + sdr + "E" + se
  256.         End If
  257.     End If
  258.     If Left$(sd, 1) = "+" Then
  259.         sd = " " + Mid$(sd, 2)
  260.     End If
  261.     d128tostring = sd
  263.  

the 64-bit dll

output
Code: [Select]
-1.000000000000000000000000000000000
-0.9000000000000000000000000000000000
-0.8000000000000000000000000000000000
-0.7000000000000000000000000000000000
-0.6000000000000000000000000000000000
-0.5000000000000000000000000000000000
-0.4000000000000000000000000000000000
-0.3000000000000000000000000000000000
-0.2000000000000000000000000000000000
-0.1000000000000000000000000000000000
 0
 0.1000000000000000000000000000000000
 0.2000000000000000000000000000000000
 0.3000000000000000000000000000000000
 0.4000000000000000000000000000000000
 0.5000000000000000000000000000000000
 0.6000000000000000000000000000000000
 0.7000000000000000000000000000000000
 0.8000000000000000000000000000000000
 0.9000000000000000000000000000000000
 1.000000000000000000000000000000000

Press any key to continue
« Last Edit: April 07, 2022, 01:55:16 pm by jack »
Logged

Offline jack

  • Seasoned Forum Regular
  • Posts: 408
    • View Profile
Re: intel decimal math library
« Reply #2 on: April 05, 2022, 09:55:53 pm »
updated the code and dll, many functions added
Logged
Pages: [1]   Go Up
« previous next »