Return to Ebug.c CVS log

Up to [local] / OpenXM_contrib2 / asir2000 / engine

Annotation of OpenXM_contrib2/asir2000/engine/Ebug.c, Revision 1.1.1.1

1.1       noro        1: /* $OpenXM: OpenXM/src/asir99/engine/Ebug.c,v 1.1.1.1 1999/11/10 08:12:26 noro Exp $ */
                      2: #include "ca.h"
                      3:
                      4: void substvp(vl,f,vn,g)
                      5: VL vl;
                      6: P f;
                      7: VN vn;
                      8: P *g;
                      9: {
                     10:        V v;
                     11:        int i;
                     12:        P h,h1;
                     13:        Q t;
                     14:
                     15:        h = f;
                     16:        for ( i = 0; v = vn[i].v; i++ ) {
                     17:                STOQ(vn[i].n,t);
                     18:                substp(vl,h,v,(P)t,&h1); h = h1;
                     19:        }
                     20:        *g = h;
                     21: }
                     22:
                     23: void affine(vl,f,vn,fr)
                     24: VL vl;
                     25: P f;
                     26: VN vn;
                     27: P *fr;
                     28: {
                     29:        int i,j,n;
                     30:        P vv,g,g1,t,u;
                     31:        Q s;
                     32:        int *dlist;
                     33:        P **plist;
                     34:
                     35:        for ( n = 0; vn[n].v; n++);
                     36:        dlist = (int *)ALLOCA((n+1)*sizeof(int));
                     37:        plist = (P **)ALLOCA((n+1)*sizeof(P *));
                     38:        for ( i = 0; vn[i].v; i++ ) {
                     39:                if ( !vn[i].n )
                     40:                        continue;
                     41:                dlist[i] = getdeg(vn[i].v,f);
                     42:                plist[i] = (P *)ALLOCA((dlist[i]+1)*sizeof(P));
                     43:
                     44:                MKV(vn[i].v,t);
                     45:                if ( vn[i].n ) {
                     46:                        STOQ(vn[i].n,s); addp(vl,t,(P)s,&vv);
                     47:                } else
                     48:                        vv = t;
                     49:
                     50:                for ( j = 0, t = (P)ONE; j < dlist[i]; j++ ) {
                     51:                        plist[i][j] = t;
                     52:                        mulp(vl,t,vv,&u);
                     53:                        t = u;
                     54:                }
                     55:                plist[i][j] = t;
                     56:        }
                     57:
                     58:        g = f;
                     59:        for ( i = 0; vn[i].v; i++ ) {
                     60:                if ( !vn[i].n )
                     61:                        continue;
                     62:                affinemain(vl,g,vn[i].v,dlist[i],plist[i],&g1); g = g1;
                     63:        }
                     64:        *fr = g;
                     65: }
                     66: