| version 1.3, 2000/08/22 05:04:05 |
version 1.4, 2000/12/22 10:03:28 |
|
|
| * 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.3 2000/08/22 05:04:05 noro Exp $ |
| */ |
*/ |
| #include "ca.h" |
#include "ca.h" |
| #include "base.h" |
#include "base.h" |
|
|
| { |
{ |
| 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); |
|
|
| |
|
| 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; |