- Added the <~> approximative comparison operator (alias for `appro…

…x_cmp`). Example (solution to `x` for x**3 = 10): say bsearch_inverse(100, {|x| x**3 <~> 10 }) #=> 2.15443469003188372...
trizen · Dec 21, 2019 · 7f9680f · 7f9680f
1 parent 972ef55
commit 7f9680f
Show file tree

Hide file tree

Showing 3 changed files with 22 additions and 15 deletions.
diff --git a/lib/Sidef/Parser.pm b/lib/Sidef/Parser.pm
@@ -304,6 +304,7 @@ package Sidef::Parser {
 
                   %% ≅
                   ~~ !~
+                  <~>
                   <=> =~=
                   <<= >>=
                   << >>

diff --git a/lib/Sidef/Types/Number/Number.pm b/lib/Sidef/Types/Number/Number.pm
@@ -9368,6 +9368,9 @@ package Sidef::Types::Number::Number {
             $max = 11465;
         }
 
+        require Memoize;
+        Memoize::memoize('_prime_count');
+
         my $k = 0;
 
         while (1) {
@@ -9387,6 +9390,8 @@ package Sidef::Types::Number::Number {
             }
         }
 
+        Memoize::unmemoize('_prime_count');
+
         while (!Math::Prime::Util::GMP::is_semiprime($k)) {
             --$k;
         }
@@ -13810,6 +13815,7 @@ package Sidef::Types::Number::Number {
         *{__PACKAGE__ . '::' . 'δ'}   = \&kronecker_delta;
         *{__PACKAGE__ . '::' . '=~='} = \&approx_eq;
         *{__PACKAGE__ . '::' . '≅'}   = \&approx_eq;
+        *{__PACKAGE__ . '::' . '<~>'} = \&approx_cmp;
     }
 }
 

diff --git a/lib/Sidef/Types/Number/Number.pod b/lib/Sidef/Types/Number/Number.pod
@@ -228,6 +228,21 @@ Aliases: I<cmp>
 
 =cut
 
+=head2 approx_cmp
+
+    approx_cmp(a, b)
+    approx_cmp(a, b, k)
+
+Approximate comparison of C<a> and C<b>.
+
+Equivalent with:
+
+    a.round(k) <=> b.round(k)
+
+When C<k> is omitted, it uses the default floating-point precision to deduce C<k>.
+
+=cut
+
 =head2 ==
 
     a == b
@@ -564,21 +579,6 @@ Return the
 
 =cut
 
-=head2 approx_cmp
-
-    approx_cmp(a, b)
-    approx_cmp(a, b, k)
-
-Approximate comparison of C<a> and C<b>.
-
-Equivalent with:
-
-    a.round(k) <=> b.round(k)
-
-When C<k> is omitted, it uses the default floating-point precision to deduce C<k>.
-
-=cut
-
 =head2 approx_ge
 
     approx_ge(a, b)
-Original file line number
+Diff line change
@@ Expand Up / @@ -304,6 +304,7 @@ package Sidef::Parser { @@
                       %% ≅
                       ~~ !~
+                      <~>
                       <=> =~=
                       <<= >>=
                       << >>
@@ Expand Down @@