Return to ddN-27.c CVS log

Up to [local] / OpenXM_contrib2 / asir2000 / engine-27

Annotation of OpenXM_contrib2/asir2000/engine-27/ddN-27.c, Revision 1.1

1.1     ! noro        1: /* $OpenXM: OpenXM/src/asir99/engine-27/ddN-27.c,v 1.1.1.1 1999/11/10 08:12:27 noro Exp $ */
        !             2: #include "ca-27.h"
        !             3: #include "base.h"
        !             4: #include "inline.h"
        !             5:
        !             6: void divn_27(n1,n2,nq,nr)
        !             7: N n1,n2,*nq,*nr;
        !             8: {
        !             9:        int tmp;
        !            10:        int br;
        !            11:        int j,d,d1,d2,dd;
        !            12:        int *m1,*m1e,*m2,*mr;
        !            13:        N q,r;
        !            14:
        !            15:        if ( !n2 )
        !            16:                error("divn: division by 0");
        !            17:        else if ( !n1 ) {
        !            18:                *nq = 0; *nr = 0;
        !            19:        } else if ( UNIN(n2) ) {
        !            20:                *nq = n1; *nr = 0;
        !            21:        } else if ( (d1 = PL(n1)) < (d2 = PL(n2)) ) {
        !            22:                *nr = n1; *nq = 0;
        !            23:        } else if ( d2 == 1 ) {
        !            24:                if ( d1 == 1 ) {
        !            25:                        DQR(BD(n1)[0],BD(n2)[0],br,tmp)
        !            26:                        STON(br,*nq); STON(tmp,*nr);
        !            27:                } else {
        !            28:                        tmp = divin_27(n1,BD(n2)[0],nq); STON(tmp,*nr);
        !            29:                }
        !            30:                return;
        !            31:        } else if ( (d1 == d2) && ((tmp = cmpn(n1,n2)) <= 0) )
        !            32:                if ( !tmp ) {
        !            33:                        *nr = 0; *nq = ONEN;
        !            34:                } else {
        !            35:                        *nr = n1; *nq = 0;
        !            36:                }
        !            37:        else {
        !            38:                d = BASE27/(BD(n2)[d2-1] + 1); dd = d1-d2;
        !            39:                m1 = W_ALLOC(d1+1); m2 = W_ALLOC(d2+1);
        !            40:                mulin_27(n1,d,m1); mulin_27(n2,d,m2);
        !            41:                *nq = q = NALLOC(dd+1);
        !            42:                divnmain_27(d1,d2,m1,m2,BD(q));
        !            43:                PL(q) = (BD(q)[dd]?dd+1:dd);
        !            44:                for ( j = d2-1; (j >= 0) && (m1[j] == 0); j--);
        !            45:                if ( j == -1 )
        !            46:                        *nr = 0;
        !            47:                else {
        !            48:                        *nr = r = NALLOC(j+1);
        !            49:                        for ( br = 0, m1e = m1+j, mr = BD(r)+j;
        !            50:                                m1e >= m1; m1e--, mr-- ) {
        !            51:                                DSAB27(d,br,*m1e,*mr,br)
        !            52:                        }
        !            53:                        PL(r) = (BD(r)[j]?j+1:j);
        !            54:                }
        !            55:        }
        !            56: }
        !            57:
        !            58: void divnmain_27(d1,d2,m1,m2,q)
        !            59: int d1,d2;
        !            60: int *m1,*m2,*q;
        !            61: {
        !            62:        unsigned int tmp;
        !            63:        int *n1,*n2;
        !            64:        int i,j;
        !            65:        int l,r;
        !            66:        int u,qhat;
        !            67:        int v1,v2;
        !            68:
        !            69:        v1 = m2[d2-1]; v2 = m2[d2-2];
        !            70:        for ( j = d1-d2, m1 += j, q += j; j >= 0; j--, q--, m1-- ) {
        !            71:                n1 = m1+d2; n2 = m1+d2-1;
        !            72:                if ( *n1 == v1 ) {
        !            73:                        qhat = BASE27 - 1; r = *n1 + *n2;
        !            74:                } else {
        !            75:                        DSAB27(v1,*n1,*n2,qhat,r)
        !            76:                }
        !            77:                DM27(v2,qhat,u,l)
        !            78:                while ( 1 ) {
        !            79:                        if ((r > u) || ((r == u) && (*(n1-2) >= l)))
        !            80:                                break;
        !            81:                        if ( ( l -= v2 ) < 0 ) {
        !            82:                                l += BASE27; u--;
        !            83:                        }
        !            84:                r += v1; qhat--;
        !            85:                }
        !            86:                if ( qhat ) {
        !            87:                        for ( i = d2, u = 0, n1 = m1, n2 = m2; i > 0; i--, n1++, n2++ ) {
        !            88:                                DMA27(*n2,qhat,u,u,l)
        !            89:                                if ( (l = *n1 - l) < 0 ) {
        !            90:                                        u++; *n1 = l + BASE27;
        !            91:                                } else
        !            92:                                        *n1 = l;
        !            93:                        }
        !            94:                        if ( *n1 < u ) {
        !            95:                                for ( i = d2, qhat--, r = 0, n1 = m1, n2 = m2;
        !            96:                                        i > 0; i--, n1++, n2++ ) {
        !            97:                                        l = *n1 + *n2 + r; *n1 = l % BASE27; r = l / BASE27;
        !            98:                                }
        !            99:                        }
        !           100:                        *n1 = 0;
        !           101:                }
        !           102:                *q = qhat;
        !           103:        }
        !           104: }
        !           105:
        !           106: void mulin_27(n,d,p)
        !           107: N n;
        !           108: int d;
        !           109: int *p;
        !           110: {
        !           111:        int carry;
        !           112:        int *m,*me;
        !           113:
        !           114:        bzero((char *)p,(int)((PL(n)+1)*sizeof(int)));
        !           115:        for ( carry = 0, m = BD(n), me = m+PL(n); m < me; p++, m++ ) {
        !           116:                carry += *p;
        !           117:                DMA27(*m,d,carry,carry,*p)
        !           118:        }
        !           119:        *p = carry;
        !           120: }
        !           121:
        !           122: int divin_27(n,dvr,q)
        !           123: N n;
        !           124: int dvr;
        !           125: N *q;
        !           126: {
        !           127:        int up,d;
        !           128:        int *m,*mq,*md;
        !           129:        N nq;
        !           130:
        !           131:        d = PL(n); m = BD(n);
        !           132:        if ( ( d == 1 ) && ( dvr > *m ) ) {
        !           133:                *q = 0;
        !           134:                return( *m );
        !           135:        }
        !           136:        *q = nq = NALLOC(d);
        !           137:        for ( md = m+d-1, mq = BD(nq)+d-1, up = 0; md >= m; md--, mq-- ) {
        !           138:                DSAB27(dvr,up,*md,*mq,up)
        !           139:        }
        !           140:        PL(nq) = (BD(nq)[d-1]?d:d-1);
        !           141:        return ( up );
        !           142: }