Mercurial > octave

changeset 25546:93b564c789aa

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
maint: merge stable to default.
author Rik <rik@octave.org>
date Mon, 02 Jul 2018 14:52:28 -0700
parents 15460ae5462d (current diff) d068c71d4eed (diff)
children 2b9a30925a9c
files doc/interpreter/module.mk
diffstat 3 files changed, 8 insertions(+), 8 deletions(-) [+]
line wrap: on
line diff
--- a/libinterp/corefcn/__betainc__.cc	Mon Jul 02 22:27:00 2018 +0200
+++ b/libinterp/corefcn/__betainc__.cc	Mon Jul 02 14:52:28 2018 -0700
@@ -77,7 +77,7 @@
         b = args(2).float_array_value ();
 
       // Initialize variables used in algorithm
-      static const float tiny = pow (2, -50);
+      static const float tiny = octave::math::exp2 (-50.0f);
       static const float eps = std::numeric_limits<float>::epsilon ();
       float xj, x2, y, Cj, Dj, aj, bj, Deltaj, alpha_j, beta_j;
       int j, maxit;
@@ -148,7 +148,7 @@
         b = args(2).array_value ();
 
       // Initialize variables used in algorithm
-      static const double tiny = pow (2, -100);
+      static const double tiny = octave::math::exp2 (-100.0);
       static const double eps = std::numeric_limits<double>::epsilon ();
       double xj, x2, y, Cj, Dj, aj, bj, Deltaj, alpha_j, beta_j;
       int j, maxit;
--- a/libinterp/corefcn/__expint__.cc	Mon Jul 02 22:27:00 2018 +0200
+++ b/libinterp/corefcn/__expint__.cc	Mon Jul 02 14:52:28 2018 -0700
@@ -61,7 +61,7 @@
         x = args(0).float_complex_array_value ();
 
       // Initialize variables used in algorithm
-      static const FloatComplex tiny = pow (2, -50);
+      static const FloatComplex tiny = octave::math::exp2 (-50.0f);
       static const float eps = std::numeric_limits<float>::epsilon ();
       const FloatComplex cone (1.0, 0.0);
       const FloatComplex czero (0.0, 0.0);
@@ -95,7 +95,7 @@
               Dj = cone / Dj;
               Deltaj = Cj * Dj;
               y *= Deltaj;
-              alpha_j = floor ((j + 1) / 2);
+              alpha_j = (j + 1) / 2;
               if ((j % 2) == 0)
                 beta_j = xj;
               else
@@ -119,7 +119,7 @@
         x = args(0).complex_array_value ();
 
       // Initialize variables used in algorithm
-      static const Complex tiny = pow (2, -100);
+      static const Complex tiny = octave::math::exp2 (-100.0);
       static const double eps = std::numeric_limits<double>::epsilon ();
       const Complex cone (1.0, 0.0);
       const Complex czero (0.0, 0.0);
@@ -153,7 +153,7 @@
               Dj = cone / Dj;
               Deltaj = Cj * Dj;
               y *= Deltaj;
-              alpha_j = floor ((j + 1) / 2);
+              alpha_j = (j + 1) / 2;
               if ((j % 2) == 0)
                 beta_j = xj;
               else
--- a/libinterp/corefcn/__gammainc__.cc	Mon Jul 02 22:27:00 2018 +0200
+++ b/libinterp/corefcn/__gammainc__.cc	Mon Jul 02 14:52:28 2018 -0700
@@ -69,7 +69,7 @@
         a = args(1).float_array_value ();
 
       // Initialize variables used in algorithm
-      static const float tiny = pow (2, -50);
+      static const float tiny = octave::math::exp2 (-50.0f);
       static const float eps = std::numeric_limits<float>::epsilon();
       float y, Cj, Dj, bj, aj, Deltaj;
       int j, maxit;
@@ -124,7 +124,7 @@
         a = args(1).array_value ();
 
       // Initialize variables used in algorithm
-      static const double tiny = pow (2, -100);
+      static const double tiny = octave::math::exp2 (-100.0);
       static const double eps = std::numeric_limits<double>::epsilon();
       double y, Cj, Dj, bj, aj, Deltaj;
       int j, maxit;