version 1.3, 2000/08/22 05:04:05 |
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.2 2000/08/21 08:31:27 noro 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 PARI |
|
#include "genpari.h" |
|
void patori(GEN,Obj *); |
|
void patori_i(GEN,N *); |
|
void ritopa(Obj,GEN *); |
|
void ritopa_i(N,int,GEN *); |
|
#endif |
|
|
|
void toreim(a,rp,ip) |
void toreim(a,rp,ip) |
Num a; |
Num a; |
|
|
{ |
{ |
if ( !a ) |
if ( !a ) |
*rp = *ip = 0; |
*rp = *ip = 0; |
|
#if defined(INTERVAL) |
|
else if ( NID(a) <= N_PRE_C ) { |
|
#else |
else if ( NID(a) <= N_B ) { |
else if ( NID(a) <= N_B ) { |
|
#endif |
*rp = a; *ip = 0; |
*rp = a; *ip = 0; |
} else { |
} else { |
*rp = ((C)a)->r; *ip = ((C)a)->i; |
*rp = ((C)a)->r; *ip = ((C)a)->i; |
|
|
*c = b; |
*c = b; |
else if ( !b ) |
else if ( !b ) |
*c = a; |
*c = a; |
|
#if defined(INTERVAL) |
|
else if ( (NID(a) <= N_PRE_C) && (NID(b) <= N_PRE_C ) ) |
|
#else |
else if ( (NID(a) <= N_B) && (NID(b) <= N_B ) ) |
else if ( (NID(a) <= N_B) && (NID(b) <= N_B ) ) |
|
#endif |
addnum(0,a,b,c); |
addnum(0,a,b,c); |
else { |
else { |
toreim(a,&ar,&ai); toreim(b,&br,&bi); |
toreim(a,&ar,&ai); toreim(b,&br,&bi); |
|
|
chsgnnum(b,c); |
chsgnnum(b,c); |
else if ( !b ) |
else if ( !b ) |
*c = a; |
*c = a; |
|
#if defined(INTERVAL) |
|
else if ( (NID(a) <= N_PRE_C) && (NID(b) <= N_PRE_C ) ) |
|
#else |
else if ( (NID(a) <= N_B) && (NID(b) <= N_B ) ) |
else if ( (NID(a) <= N_B) && (NID(b) <= N_B ) ) |
|
#endif |
subnum(0,a,b,c); |
subnum(0,a,b,c); |
else { |
else { |
toreim(a,&ar,&ai); toreim(b,&br,&bi); |
toreim(a,&ar,&ai); toreim(b,&br,&bi); |
|
|
|
|
if ( !a || !b ) |
if ( !a || !b ) |
*c = 0; |
*c = 0; |
|
#if defined(INTERVAL) |
|
else if ( (NID(a) <= N_PRE_C) && (NID(b) <= N_PRE_C ) ) |
|
#else |
else if ( (NID(a) <= N_B) && (NID(b) <= N_B ) ) |
else if ( (NID(a) <= N_B) && (NID(b) <= N_B ) ) |
|
#endif |
mulnum(0,a,b,c); |
mulnum(0,a,b,c); |
else { |
else { |
toreim(a,&ar,&ai); toreim(b,&br,&bi); |
toreim(a,&ar,&ai); toreim(b,&br,&bi); |
|
|
error("divcplx : division by 0"); |
error("divcplx : division by 0"); |
else if ( !a ) |
else if ( !a ) |
*c = 0; |
*c = 0; |
|
#if defined(INTERVAL) |
|
else if ( (NID(a) <= N_PRE_C) && (NID(b) <= N_PRE_C ) ) |
|
#else |
else if ( (NID(a) <= N_B) && (NID(b) <= N_B ) ) |
else if ( (NID(a) <= N_B) && (NID(b) <= N_B ) ) |
|
#endif |
divnum(0,a,b,c); |
divnum(0,a,b,c); |
else { |
else { |
toreim(a,&ar,&ai); toreim(b,&br,&bi); |
toreim(a,&ar,&ai); toreim(b,&br,&bi); |
|
|
{ |
{ |
int ei; |
int ei; |
Num t; |
Num t; |
extern long prec; |
|
|
|
if ( !e ) |
if ( !e ) |
*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 PARI |
error("pwrcplx : not implemented (use eval())."); |
GEN pa,pe,z; |
else { |
int ltop,lbot; |
|
|
|
ltop = avma; ritopa((Obj)a,&pa); ritopa((Obj)e,&pe); lbot = avma; |
|
z = gerepile(ltop,lbot,gpui(pa,pe,prec)); |
|
patori(z,(Obj *)c); cgiv(z); |
|
#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 ) |
|
|
|
|
if ( !a ) |
if ( !a ) |
*c = 0; |
*c = 0; |
|
#if defined(INTERVAL) |
|
else if ( NID(a) <= N_PRE_C ) |
|
#else |
else if ( NID(a) <= N_B ) |
else if ( NID(a) <= N_B ) |
|
#endif |
chsgnnum(a,c); |
chsgnnum(a,c); |
else { |
else { |
chsgnnum(((C)a)->r,&r); chsgnnum(((C)a)->i,&i); |
chsgnnum(((C)a)->r,&r); chsgnnum(((C)a)->i,&i); |
|
|
int s; |
int s; |
|
|
if ( !a ) { |
if ( !a ) { |
|
#if defined(INTERVAL) |
|
if ( !b || (NID(b)<=N_PRE_C) ) |
|
#else |
if ( !b || (NID(b)<=N_B) ) |
if ( !b || (NID(b)<=N_B) ) |
|
#endif |
return compnum(0,a,b); |
return compnum(0,a,b); |
else |
else |
return -1; |
return -1; |
} else if ( !b ) { |
} else if ( !b ) { |
|
#if defined(INTERVAL) |
|
if ( !a || (NID(a)<=N_PRE_C) ) |
|
#else |
if ( !a || (NID(a)<=N_B) ) |
if ( !a || (NID(a)<=N_B) ) |
|
#endif |
return compnum(0,a,b); |
return compnum(0,a,b); |
else |
else |
return 1; |
return 1; |