=================================================================== RCS file: /home/cvs/OpenXM_contrib2/asir2018/builtin/dp-supp.c,v retrieving revision 1.14 retrieving revision 1.16 diff -u -p -r1.14 -r1.16 --- OpenXM_contrib2/asir2018/builtin/dp-supp.c 2020/07/02 09:24:16 1.14 +++ OpenXM_contrib2/asir2018/builtin/dp-supp.c 2020/12/05 03:27:20 1.16 @@ -45,7 +45,7 @@ * DEVELOPER SHALL HAVE NO LIABILITY IN CONNECTION WITH THE USE, * PERFORMANCE OR NON-PERFORMANCE OF THE SOFTWARE. * - * $OpenXM: OpenXM_contrib2/asir2018/builtin/dp-supp.c,v 1.13 2019/12/27 08:13:59 noro Exp $ + * $OpenXM: OpenXM_contrib2/asir2018/builtin/dp-supp.c,v 1.15 2020/10/26 02:41:05 noro Exp $ */ #include "ca.h" #include "base.h" @@ -3338,6 +3338,29 @@ void dp_mod(DP p,int mod,NODE subst,DP *rp) } } +void dpm_mod(DPM p,int mod,DPM *rp) +{ + DMM m,mr,mr0; + P t; + V v; + NODE tn; + + if ( !p ) + *rp = 0; + else { + for ( mr0 = 0, m = BDY(p); m; m = NEXT(m) ) { + ptomp(mod,(P)m->c,&t); + if ( t ) { + NEXTDMM(mr0,mr); mr->c = (Obj)t; mr->dl = m->dl; + } + } + if ( mr0 ) { + NEXT(mr) = 0; MKDPM(p->nv,mr0,*rp); (*rp)->sugar = p->sugar; + } else + *rp = 0; + } +} + void dp_rat(DP p,DP *rp) { MP m,mr,mr0; @@ -3637,6 +3660,23 @@ void dpm_rest(DPM p,DPM *rp) if ( *rp ) (*rp)->sugar = p->sugar; } +} + +int dp_getdeg(DP p) +{ + int max,n,i; + MP m; + int *d; + + if ( !p ) return 0; + n = p->nv; + max = 0; + for ( m = BDY(p); m; m = NEXT(m) ) { + d = m->dl->d; + for ( i = 0; i < n; i++ ) + if ( d[i] > max ) max = d[i]; + } + return max; } int dpm_getdeg(DPM p,int *r)