login  home  contents  what's new  discussion  bug reports     help  links  subscribe  changes  refresh  edit

Edit detail for #333 asin 1.0 over Complex Float fails revision 1 of 2

1 2
Editor:
Time: 2007/11/17 22:28:24 GMT-8
Note: patch possible.

changed:
-
The first real calculus is right, the second fails.
\begin{axiom}
asin 1.0
asin (1.0+0*%i)
\end{axiom}

The problem is perhaps line 162 of trigcat.spad.
There is no test for the ... / sqrt (1-x^2) in asin definition.

asin x == atan(x/sqrt(1-x**2))


From FrancoisMaltey Fri Jan 5 07:57:48 -0600 2007
From: Francois Maltey
Date: Fri, 05 Jan 2007 07:57:48 -0600
Subject: patch possible.
Message-ID: <20070105075748-0600@wiki.axiom-developer.org>

I propose to remplace in asin x == atan (x/sqrt(1-x**2)) in trigcat.spad line 162 by

       asin x ==
         x=1  => 2*atan 1
         x=-1 => -2*atan 1
         atan(x/sqrt(1-x**2))

The result seems right...

Submitted by : (unknown) at: 2007-11-17T22:28:24-08:00 (16 years ago)
Name :
Axiom Version :
Category : Severity : Status :
Optional subject :  
Optional comment :

The first real calculus is right, the second fails.

axiom
asin 1.0
LatexWiki Image(1)
Type: Float
axiom
asin (1.0+0*%i) >> Error detected within library code: catdef: division by zero

The problem is perhaps line 162 of trigcat.spad. There is no test for the ... / sqrt (1-x^2) in asin definition.

asin x == atan(x/sqrt(1-x**2))

patch possible. --Francois Maltey, Fri, 05 Jan 2007 07:57:48 -0600 reply
I propose to remplace in asin x == atan (x/sqrt(1-x**2)) in trigcat.spad line 162 by

asin x == x=1 => 2atan 1 x=-1 => -2atan 1 atan(x/sqrt(1-x**2))

The result seems right...