From: Ramakrishnan Muthukrishnan <vu3rdd@gmail.com>
To: guile-devel <guile-devel@gnu.org>
Subject: Re: fix for expt bug
Date: Tue, 2 Nov 2010 11:06:06 +0530 [thread overview]
Message-ID: <AANLkTim5SwwEtAArLdd5_BnoNAtQE79rrhLBBny=Ky1+@mail.gmail.com> (raw)
In-Reply-To: <AANLkTim=2gFn2U8ca8ZDxOcy19k+v_dGBy8q3e0Z6tLF@mail.gmail.com>
[-- Attachment #1: Type: text/plain, Size: 48 bytes --]
Here is the updated patch.
thanks
Ramakrishnan
[-- Attachment #2: expt-fix.patch --]
[-- Type: text/x-patch, Size: 3066 bytes --]
From e320d79c8f3cd8b7ddcb9c2d13356e34a3346cfe Mon Sep 17 00:00:00 2001
From: Ramakrishnan Muthukrishnan <vu3rdd@gmail.com>
Date: Sun, 31 Oct 2010 23:22:52 +0530
Subject: [PATCH] Fix for bug #31464. expt needs to treat negative bases specially.
* libguile/numbers.c: If base is negative, expt needs to find
-x^n = (-1^n) * (|x|^n). We find x^n and then if n is odd, we
also multiply the result with -1. These changes apply only for
cases where n is an integer.
* test-suite/tests/numbers.test: Two new test cases for expt. For
cases where the base is negative and the power to be raised is
not an integer, the result should be a complex number.
---
libguile/numbers.c | 13 ++++++++++++-
test-suite/tests/numbers.test | 15 ++++++++++++++-
2 files changed, 26 insertions(+), 2 deletions(-)
diff --git a/libguile/numbers.c b/libguile/numbers.c
index fbc6cc8..4b81d98 100644
--- a/libguile/numbers.c
+++ b/libguile/numbers.c
@@ -5445,12 +5445,23 @@ SCM_DEFINE (scm_expt, "expt", 2, 0, 0,
"Return @var{x} raised to the power of @var{y}.")
#define FUNC_NAME s_scm_expt
{
- if (scm_is_true (scm_exact_p (x)) && scm_is_integer (y))
+ if (scm_is_true (scm_exact_p (y)) && scm_is_integer (y))
return scm_integer_expt (x, y);
else if (scm_is_real (x) && scm_is_real (y) && scm_to_double (x) >= 0.0)
{
return scm_from_double (pow (scm_to_double (x), scm_to_double (y)));
}
+ else if (scm_is_real (x) && scm_is_integer (y) && (scm_to_double (x) < 0))
+ {
+ double result;
+
+ result = pow (-scm_to_double (x), scm_to_double (y));
+
+ if (scm_is_true (scm_odd_p (y)))
+ return scm_from_double (-result);
+ else
+ return scm_from_double (result);
+ }
else
return scm_exp (scm_product (scm_log (x), y));
}
diff --git a/test-suite/tests/numbers.test b/test-suite/tests/numbers.test
index 3c3e14f..637449f 100644
--- a/test-suite/tests/numbers.test
+++ b/test-suite/tests/numbers.test
@@ -2892,7 +2892,20 @@
(pass-if "(= 1 (expt 0 0))" (= 1 (expt 0 0)))
(pass-if "(= 1 (expt 0 0.0))" (= 1 (expt 0 0.0)))
(pass-if "(= 1 (expt 0.0 0))" (= 1 (expt 0.0 0)))
- (pass-if "(= 1 (expt 0.0 0.0))" (= 1 (expt 0.0 0.0))))
+ (pass-if "(= 1 (expt 0.0 0.0))" (= 1 (expt 0.0 0.0)))
+ ;; tests for non zero values of base and exponent.
+ (pass-if "(eqv-loosely? -2742638075.5 (expt -2742638075.5 1))"
+ (eqv-loosely? -2742638075.5 (expt -2742638075.5 1)))
+ (pass-if "(eqv-loosely? 4.0 (expt -2.0 2.0))"
+ (eqv-loosely? 4.0 (expt -2.0 2.0)))
+ (pass-if "(eqv-loosely? -0.125 (expt -2.0 -3.0))"
+ (eqv-loosely? -0.125 (expt -2.0 -3.0)))
+ (pass-if "(eqv-loosely? -0.125 (expt -2 -3.0))"
+ (eqv-loosely? -0.125 (expt -2 -3.0)))
+ (pass-if "(eqv-loosely? 0.25 (expt 2.0 -2.0))"
+ (eqv-loosely? 0.25 (expt 2.0 -2.0)))
+ (pass-if "(eqv-loosely? 1.0000000000000002+1.7320508075688772i (expt -8 1/3))"
+ (eqv-loosely? 1.0000000000000002+1.7320508075688772i (expt -8 1/3))))
;;;
;;; asinh
--
1.7.2.3
next prev parent reply other threads:[~2010-11-02 5:36 UTC|newest]
Thread overview: 24+ messages / expand[flat|nested] mbox.gz Atom feed top
2010-10-31 18:12 fix for expt bug Ramakrishnan Muthukrishnan
2010-11-01 0:03 ` Mark H Weaver
2010-11-01 4:03 ` Ramakrishnan Muthukrishnan
2010-11-01 18:23 ` Ramakrishnan Muthukrishnan
2010-11-02 4:55 ` Mark H Weaver
2010-11-02 5:22 ` Ramakrishnan Muthukrishnan
2010-11-02 5:36 ` Ramakrishnan Muthukrishnan [this message]
2010-11-02 19:19 ` Mark H Weaver
2010-11-02 23:00 ` Ludovic Courtès
2010-11-03 2:10 ` Mark H Weaver
2010-11-03 9:41 ` Ramakrishnan Muthukrishnan
2010-11-03 15:32 ` Mark H Weaver
2010-11-03 16:41 ` Ramakrishnan Muthukrishnan
2010-11-03 16:22 ` Mark H Weaver
2010-11-03 17:53 ` Ramakrishnan Muthukrishnan
2010-11-03 18:25 ` Ramakrishnan Muthukrishnan
2010-11-03 23:27 ` Ludovic Courtès
2010-11-04 17:27 ` Mark H Weaver
2010-11-08 21:08 ` Ludovic Courtès
2010-11-20 21:49 ` Andy Wingo
2010-11-21 22:34 ` Ludovic Courtès
2010-11-05 2:18 ` [PATCH] Fix bugs in expt and integer-expt Mark H Weaver
2010-11-10 22:01 ` Ludovic Courtès
2011-01-20 22:47 ` Andy Wingo
Reply instructions:
You may reply publicly to this message via plain-text email
using any one of the following methods:
* Save the following mbox file, import it into your mail client,
and reply-to-all from there: mbox
Avoid top-posting and favor interleaved quoting:
https://en.wikipedia.org/wiki/Posting_style#Interleaved_style
List information: https://www.gnu.org/software/guile/
* Reply using the --to, --cc, and --in-reply-to
switches of git-send-email(1):
git send-email \
--in-reply-to='AANLkTim5SwwEtAArLdd5_BnoNAtQE79rrhLBBny=Ky1+@mail.gmail.com' \
--to=vu3rdd@gmail.com \
--cc=guile-devel@gnu.org \
/path/to/YOUR_REPLY
https://kernel.org/pub/software/scm/git/docs/git-send-email.html
* If your mail client supports setting the In-Reply-To header
via mailto: links, try the mailto: link
Be sure your reply has a Subject: header at the top and a blank line
before the message body.
This is a public inbox, see mirroring instructions
for how to clone and mirror all data and code used for this inbox;
as well as URLs for read-only IMAP folder(s) and NNTP newsgroup(s).