Return to vect.c CVS log

Up to [local] / OpenXM_contrib2 / asir2000 / engine

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

1.1       noro        1: /* $OpenXM: OpenXM/src/asir99/engine/vect.c,v 1.1.1.1 1999/11/10 08:12:26 noro Exp $ */
                      2: #include "ca.h"
                      3:
                      4: void addvect(vl,a,b,c)
                      5: VL vl;
                      6: VECT a,b,*c;
                      7: {
                      8:        int len,i;
                      9:        VECT t;
                     10:        pointer *ab,*bb,*tb;
                     11:
                     12:        if ( !a )
                     13:                *c = b;
                     14:        else if ( !b )
                     15:                *c = a;
                     16:        else if ( a->len != b->len ) {
                     17:                *c = 0; error("addvect : size mismatch");
                     18:        } else {
                     19:                len = a->len;
                     20:                MKVECT(t,len);
                     21:                for ( i = 0, ab = BDY(a), bb = BDY(b), tb = BDY(t); i < len; i++ )
                     22:                        addr(vl,(Obj)ab[i],(Obj)bb[i],(Obj *)&tb[i]);
                     23:                *c = t;
                     24:        }
                     25: }
                     26:
                     27: void subvect(vl,a,b,c)
                     28: VL vl;
                     29: VECT a,b,*c;
                     30: {
                     31:        int len,i;
                     32:        VECT t;
                     33:        pointer *ab,*bb,*tb;
                     34:
                     35:        if ( !a )
                     36:                chsgnvect(b,c);
                     37:        else if ( !b )
                     38:                *c = a;
                     39:        else if (a->len != b->len ) {
                     40:                *c = 0; error("subvect : size mismatch");
                     41:        } else {
                     42:                len = a->len;
                     43:                MKVECT(t,len);
                     44:                for ( i = 0, ab = BDY(a), bb = BDY(b), tb = BDY(t);
                     45:                                i < len; i++ )
                     46:                                subr(vl,(Obj)ab[i],(Obj)bb[i],(Obj *)&tb[i]);
                     47:                *c = t;
                     48:        }
                     49: }
                     50:
                     51: void mulvect(vl,a,b,c)
                     52: VL vl;
                     53: Obj a,b,*c;
                     54: {
                     55:        if ( !a || !b )
                     56:                *c = 0;
                     57:        else if ( OID(a) <= O_R )
                     58:                mulrvect(vl,a,(VECT)b,(VECT *)c);
                     59:        else if ( OID(b) <= O_R )
                     60:                mulrvect(vl,b,(VECT)a,(VECT *)c);
                     61:        else
                     62:                notdef(vl,a,b,c);
                     63: }
                     64:
                     65: void divvect(vl,a,b,c)
                     66: VL vl;
                     67: Obj a,b,*c;
                     68: {
                     69:        Obj t;
                     70:
                     71:        if ( !b )
                     72:                error("divvect : division by 0");
                     73:        else if ( !a )
                     74:                *c = 0;
                     75:        else if ( OID(b) > O_R )
                     76:                notdef(vl,a,b,c);
                     77:        else {
                     78:                divr(vl,(Obj)ONE,b,&t); mulrvect(vl,(Obj)t,(VECT)a,(VECT *)c);
                     79:        }
                     80: }
                     81:
                     82: void chsgnvect(a,b)
                     83: VECT a,*b;
                     84: {
                     85:        VECT t;
                     86:        int len,i;
                     87:        pointer *ab,*tb;
                     88:
                     89:        if ( !a )
                     90:                *b = 0;
                     91:        else {
                     92:                len = a->len;
                     93:                MKVECT(t,len);
                     94:                for ( i = 0, ab = BDY(a), tb = BDY(t);
                     95:                        i < len; i++ )
                     96:                        chsgnr((Obj)ab[i],(Obj *)&tb[i]);
                     97:                *b = t;
                     98:        }
                     99: }
                    100:
                    101: void mulrvect(vl,a,b,c)
                    102: VL vl;
                    103: Obj a;
                    104: VECT b,*c;
                    105: {
                    106:        int len,i;
                    107:        VECT t;
                    108:        pointer *bb,*tb;
                    109:
                    110:        if ( !a || !b )
                    111:                *c = 0;
                    112:        else {
                    113:                len = b->len;
                    114:                MKVECT(t,len);
                    115:                for ( i = 0, bb = BDY(b), tb = BDY(t); i < len; i++ )
                    116:                        mulr(vl,a,(Obj)bb[i],(Obj *)&tb[i]);
                    117:                *c = t;
                    118:        }
                    119: }
                    120:
                    121: int compvect(vl,a,b)
                    122: VL vl;
                    123: VECT a,b;
                    124: {
                    125:        int i,len,t;
                    126:
                    127:        if ( !a )
                    128:                return b?-1:0;
                    129:        else if ( !b )
                    130:                return 1;
                    131:        else if ( a->len != b->len )
                    132:                return a->len>b->len ? 1 : -1;
                    133:        else {
                    134:                for ( i = 0, len = a->len; i < len; i++ )
                    135:                        if ( t = compr(vl,(Obj)BDY(a)[i],(Obj)BDY(b)[i]) )
                    136:                                return t;
                    137:                return 0;
                    138:        }
                    139: }