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; |