OpenXM_contrib2/asir2018/builtin/dp-supp.c - diff

Return to dp-supp.c CVS log

Up to [local] / OpenXM_contrib2 / asir2018 / builtin

Diff for /OpenXM_contrib2/asir2018/builtin/dp-supp.c between version 1.13 and 1.16

-version 1.13, 2019/12/27 08:13:59
+version 1.16, 2020/12/05 03:27:20
 Line 45
 Line 45
 Line 45
   * 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.12 2019/12/12 04:44: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"
 Line 315  int _dl_redble(DL d1,DL d2,int nvar)
 Line 315  int _dl_redble(DL d1,DL d2,int nvar)
 Line 315  int _dl_redble(DL d1,DL d2,int nvar)
      return 1;
  }
+ int _dl_redble_ext(DL d1,DL d2,DL quo,int nvar)
+ {
+   int i;
+   if ( (quo->td = d2->td-d1->td) < 0 )
+     return 0;
+   for ( i = 0; i < nvar; i++ )
+     if ( (quo->d[i] = d2->d[i]-d1->d[i]) < 0 )
+       break;
+   if ( i < nvar )
+     return 0;
+   else
+     return 1;
+ }
  void insert_to_node(DL d,NODE *n,int nvar)
  {
    DL d1;
-Line 3323  void dp_mod(DP p,int mod,NODE subst,DP *rp)
+Line 3338  void dp_mod(DP p,int mod,NODE subst,DP *rp)
 Line 3323  void dp_mod(DP p,int mod,NODE subst,DP *rp)
 Line 3338  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;
-Line 3622  void dpm_rest(DPM p,DPM *rp)
+Line 3660  void dpm_rest(DPM p,DPM *rp)
 Line 3622  void dpm_rest(DPM p,DPM *rp)
 Line 3660  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)

FreeBSD-CVSweb <freebsd-cvsweb@FreeBSD.org>