version 1.3, 2000/05/30 01:35:12 |
version 1.4, 2000/07/13 05:09:01 |
|
|
/* $OpenXM: OpenXM_contrib2/asir2000/engine/distm.c,v 1.2 2000/05/29 08:54:46 noro Exp $ */ |
/* $OpenXM: OpenXM_contrib2/asir2000/engine/distm.c,v 1.3 2000/05/30 01:35:12 noro Exp $ */ |
#include "ca.h" |
#include "ca.h" |
#include "inline.h" |
#include "inline.h" |
|
|
|
|
/* offset of h-degree */ |
/* offset of h-degree */ |
NEWDL(d,n); |
NEWDL(d,n); |
d->td = d->d[n-1] = d0->d[n-1]+d1->d[n-1]; |
d->td = d->d[n-1] = d0->d[n-1]+d1->d[n-1]; |
NEWMP(mr); mr->c = (P)ONEM; mr->dl = d; |
NEWMP(mr); mr->c = (P)ONEM; mr->dl = d; NEXT(mr) = 0; |
MKDP(n,mr,r); r->sugar = d->td; |
MKDP(n,mr,r); r->sugar = d->td; |
} else |
} else |
r = (DP)ONEM; |
r = (DP)ONEM; |
|
|
d->td = d->d[n-1] = d0->d[n-1]+d1->d[n-1]; |
d->td = d->d[n-1] = d0->d[n-1]+d1->d[n-1]; |
else |
else |
d->td = 0; |
d->td = 0; |
NEWMP(mr); mr->c = STOI(c); mr->dl = d; |
NEWMP(mr); mr->c = STOI(c); mr->dl = d; NEXT(mr) = 0; |
MKDP(n,mr,r); r->sugar = d->td; |
MKDP(n,mr,r); r->sugar = d->td; |
|
|
/* homogenized computation; dx-xd=h^2 */ |
/* homogenized computation; dx-xd=h^2 */ |