| version 1.6, 2009/03/27 14:42:29 |
version 1.7, 2015/08/20 08:42:07 |
|
|
| * DEVELOPER SHALL HAVE NO LIABILITY IN CONNECTION WITH THE USE, |
* DEVELOPER SHALL HAVE NO LIABILITY IN CONNECTION WITH THE USE, |
| * PERFORMANCE OR NON-PERFORMANCE OF THE SOFTWARE. |
* PERFORMANCE OR NON-PERFORMANCE OF THE SOFTWARE. |
| * |
* |
| * $OpenXM: OpenXM_contrib2/asir2000/engine/cplx.c,v 1.5 2003/02/14 22:29:08 ohara Exp $ |
* $OpenXM: OpenXM_contrib2/asir2000/engine/cplx.c,v 1.6 2009/03/27 14:42:29 ohara Exp $ |
| */ |
*/ |
| #include "ca.h" |
#include "ca.h" |
| #include "base.h" |
#include "base.h" |
| #if defined(PARI) |
|
| #include "genpari.h" |
|
| #endif |
|
| |
|
| void toreim(a,rp,ip) |
void toreim(a,rp,ip) |
| Num a; |
Num a; |
|
|
| *c = (Num)ONE; |
*c = (Num)ONE; |
| else if ( !a ) |
else if ( !a ) |
| *c = 0; |
*c = 0; |
| else if ( !INT(e) ) { |
else if ( !INT(e) ) |
| #if defined(PARI) |
error("pwrcplx : not implemented (use eval())."); |
| gpui_ri((Obj)a,(Obj)c,(Obj *)c); |
else { |
| #else |
|
| error("pwrcplx : can't calculate a fractional power"); |
|
| #endif |
|
| } else { |
|
| ei = SGN((Q)e)*QTOS((Q)e); |
ei = SGN((Q)e)*QTOS((Q)e); |
| pwrcplx0(a,ei,&t); |
pwrcplx0(a,ei,&t); |
| if ( SGN((Q)e) < 0 ) |
if ( SGN((Q)e) < 0 ) |
![[BACK]](/icons/cvsweb/back.gif){kind=icon}
Return to cplx.c CVS log ![[TXT]](/icons/cvsweb/text.gif){kind=icon}
![[DIR]](/icons/cvsweb/dir.gif){kind=icon}
Up to [local] / OpenXM_contrib2 / asir2000 / engine