Return to mi.c CVS log

Up to [local] / OpenXM_contrib2 / asir2000 / engine

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

1.1       noro        1: /* $OpenXM: OpenXM/src/asir99/engine/mi.c,v 1.1.1.1 1999/11/10 08:12:26 noro Exp $ */
                      2: #include "ca.h"
                      3:
                      4: extern int current_mod;
                      5:
                      6: void addmi(a,b,c)
                      7: MQ a,b;
                      8: MQ *c;
                      9: {
                     10:        if ( !current_mod )
                     11:                error("addmi : invalid modulus");
                     12:        if ( (a && NID(a) != N_M) || (b && NID(b) != N_M) )
                     13:                error("addmi : invalid argument");
                     14:        addmq(current_mod,a,b,c);
                     15: }
                     16:
                     17: void submi(a,b,c)
                     18: MQ a,b;
                     19: MQ *c;
                     20: {
                     21:        if ( !current_mod )
                     22:                error("submi : invalid modulus");
                     23:        if ( (a && NID(a) != N_M) || (b && NID(b) != N_M) )
                     24:                error("submi : invalid argument");
                     25:        submq(current_mod,a,b,c);
                     26: }
                     27:
                     28: void mulmi(a,b,c)
                     29: MQ a,b;
                     30: MQ *c;
                     31: {
                     32:        if ( !current_mod )
                     33:                error("mulmi : invalid modulus");
                     34:        if ( (a && NID(a) != N_M) || (b && NID(b) != N_M) )
                     35:                error("mulmi : invalid argument");
                     36:        mulmq(current_mod,a,b,c);
                     37: }
                     38:
                     39: void divmi(a,b,c)
                     40: MQ a,b;
                     41: MQ *c;
                     42: {
                     43:        if ( !current_mod )
                     44:                error("divmi : invalid modulus");
                     45:        if ( (a && NID(a) != N_M) || (b && NID(b) != N_M) )
                     46:                error("divmi : invalid argument");
                     47:        divmq(current_mod,a,b,c);
                     48: }
                     49:
                     50: void chsgnmi(a,c)
                     51: MQ a,*c;
                     52: {
                     53:        if ( !current_mod )
                     54:                error("chsgnmi : invalid modulus");
                     55:        if ( a && NID(a) != N_M )
                     56:                error("chsgnmi : invalid argument");
                     57:        submq(current_mod,0,a,c);
                     58: }
                     59:
                     60: void pwrmi(a,b,c)
                     61: MQ a;
                     62: Q b;
                     63: MQ *c;
                     64: {
                     65:        if ( !current_mod )
                     66:                error("pwrmi : invalid modulus");
                     67:        if ( a && NID(a) != N_M )
                     68:                error("pwrmi : invalid argument");
                     69:        pwrmq(current_mod,a,b,c);
                     70: }
                     71:
                     72: int cmpmi(a,b)
                     73: MQ a,b;
                     74: {
                     75:        int c;
                     76:
                     77:        if ( !current_mod )
                     78:                error("cmpmi : invalid modulus");
                     79:        if ( (a && NID(a) != N_M) || (b && NID(b) != N_M) )
                     80:                error("cmpmi : invalid argument");
                     81:        if ( !a )
                     82:                if ( !b )
                     83:                        return 0;
                     84:                else
                     85:                        return CONT(b) < 0;
                     86:        else
                     87:                if ( !b )
                     88:                        return CONT(a) > 0;
                     89:                else {
                     90:                        c = CONT(a)-CONT(b);
                     91:                        if ( c % current_mod )
                     92:                                return c > 0;
                     93:                        else
                     94:                                return 0;
                     95:                }
                     96: }