axiom-developer
[Top][All Lists]
Advanced

[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

[Axiom-developer] [#138 asin(complex(1.0,0.0))]


From: unknown
Subject: [Axiom-developer] [#138 asin(complex(1.0,0.0))]
Date: Fri, 01 Apr 2005 16:59:20 -0600

Changes http://page.axiom-developer.org/zope/mathaction/138AsinComplex1000/diff
--
This and the log(i)-log(-i) bug can be fixed by using GCL's asin and log here.  
I tried 

     if $ has Field and $ has sqrt: $ -> $ then
--       asin x == atan(x/sqrt(1-x**2))
       asin x == ASIN(COMPLEX(real x,imag x)$Lisp)$Lisp

but the form of x when complex in axiom is not that of lisp, i.e. asin gets 
passed a non-number:

   >> System error:
   ((1 . 0) 0 . 0) is not of type NUMBER.


In general, I think it a good idea to avoid having to duplicate the 
functionality of the standard lisp functions in axiom, but there may be some 
philosophical reason for doing this due to the domain structures or some such.  
If one needs the correct logic for a duplicate implementation, one can look at 
the gcl_numlib.lsp file in the lsp/ subdirectory.




--
forwarded from http://page.axiom-developer.org/zope/mathaction/address@hidden




reply via email to

[Prev in Thread] Current Thread [Next in Thread]