version 1.7, 2000/12/08 02:39:04 |
version 1.9, 2000/12/08 08:26:08 |
|
|
* 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/builtin/dp-supp.c,v 1.6 2000/12/05 08:29:43 noro Exp $ |
* $OpenXM: OpenXM_contrib2/asir2000/builtin/dp-supp.c,v 1.8 2000/12/08 06:43:09 noro Exp $ |
*/ |
*/ |
#include "ca.h" |
#include "ca.h" |
#include "base.h" |
#include "base.h" |
|
|
} |
} |
} |
} |
|
|
void dp_ptozp_d(dist,ndist,p,rp) |
extern int mpi_mag; |
NODE dist; |
extern int PCoeffs; |
int ndist; |
|
|
void dp_ptozp_d(p,rp) |
DP p,*rp; |
DP p,*rp; |
{ |
{ |
int i,j,k,l,n,nsep; |
int i,j,k,l,n,nsep; |
|
|
N qn,gn; |
N qn,gn; |
double get_rtime(); |
double get_rtime(); |
int blen; |
int blen; |
|
NODE dist; |
|
int ndist; |
double t0; |
double t0; |
double t_e,t_d,t_d1,t_c; |
double t_e,t_d,t_d1,t_c; |
|
extern int DP_NFStat; |
|
extern LIST Dist; |
|
|
if ( !p ) |
if ( !p ) |
*rp = 0; |
*rp = 0; |
else { |
else { |
|
if ( PCoeffs ) { |
|
dp_ptozp(p,rp); return; |
|
} |
|
if ( !Dist || p_mag(BDY(p)->c) <= mpi_mag ) { |
|
dist = 0; ndist = 0; |
|
if ( DP_NFStat ) fprintf(asir_out,"L"); |
|
} else { |
|
dist = BDY(Dist); ndist = length(dist); |
|
if ( DP_NFStat ) fprintf(asir_out,"D"); |
|
} |
for ( m = BDY(p), n = 0; m; m = NEXT(m), n++ ); |
for ( m = BDY(p), n = 0; m; m = NEXT(m), n++ ); |
nsep = ndist + 1; |
nsep = ndist + 1; |
if ( n <= nsep ) { |
if ( n <= nsep ) { |
|
|
} |
} |
} |
} |
|
|
void dp_ptozp2_d(dist,ndist,p0,p1,hp,rp) |
void dp_ptozp2_d(p0,p1,hp,rp) |
NODE dist; |
|
int ndist; |
|
DP p0,p1; |
DP p0,p1; |
DP *hp,*rp; |
DP *hp,*rp; |
{ |
{ |
DP t,s,h,r; |
DP t,s,h,r; |
MP m,mr,mr0,m0; |
MP m,mr,mr0,m0; |
|
|
addd(CO,p0,p1,&t); dp_ptozp_d(dist,ndist,t,&s); |
addd(CO,p0,p1,&t); dp_ptozp_d(t,&s); |
if ( !p0 ) { |
if ( !p0 ) { |
h = 0; r = s; |
h = 0; r = s; |
} else if ( !p1 ) { |
} else if ( !p1 ) { |