Return to cplxnum.c CVS log

Up to [local] / OpenXM_contrib2 / asir2000 / builtin

Annotation of OpenXM_contrib2/asir2000/builtin/cplxnum.c, Revision 1.1

1.1     ! noro        1: /* $OpenXM: OpenXM/src/asir99/builtin/cplxnum.c,v 1.1.1.1 1999/11/10 08:12:25 noro Exp $ */
        !             2: #include "ca.h"
        !             3: #include "parse.h"
        !             4:
        !             5: void Pconj(),Preal(),Pimag();
        !             6:
        !             7: #if defined(THINK_C)
        !             8: void cplx_conj(Obj,Obj *);
        !             9: void cplx_real(Obj,Obj *);
        !            10: void cplx_imag(Obj,Obj *);
        !            11: #else
        !            12: void cplx_conj();
        !            13: void cplx_real();
        !            14: void cplx_imag();
        !            15: #endif
        !            16:
        !            17: struct ftab cplx_tab[] = {
        !            18:        {"conj",Pconj,1},
        !            19:        {"real",Preal,1},
        !            20:        {"imag",Pimag,1},
        !            21:        {0,0,0},
        !            22: };
        !            23:
        !            24: void Pconj(arg,rp)
        !            25: NODE arg;
        !            26: Obj *rp;
        !            27: {
        !            28:        cplx_conj((Obj)ARG0(arg),rp);
        !            29: }
        !            30:
        !            31: void Preal(arg,rp)
        !            32: NODE arg;
        !            33: Obj *rp;
        !            34: {
        !            35:        cplx_real((Obj)ARG0(arg),rp);
        !            36: }
        !            37:
        !            38: void Pimag(arg,rp)
        !            39: NODE arg;
        !            40: Obj *rp;
        !            41: {
        !            42:        cplx_imag((Obj)ARG0(arg),rp);
        !            43: }
        !            44:
        !            45: void cplx_conj(p,r)
        !            46: Obj p;
        !            47: Obj *r;
        !            48: {
        !            49:        C c;
        !            50:        DCP dc,dcr,dcr0;
        !            51:        P t;
        !            52:
        !            53:        if (  !p )
        !            54:                *r = 0;
        !            55:        else
        !            56:                switch ( OID(p) ) {
        !            57:                        case O_N:
        !            58:                                if ( NID((Num)p) <= N_B )
        !            59:                                        *r = p;
        !            60:                                else {
        !            61:                                        NEWC(c); c->r = ((C)p)->r; chsgnnum(((C)p)->i,&c->i);
        !            62:                                        *r = (Obj)c;
        !            63:                                }
        !            64:                                break;
        !            65:                        case O_P:
        !            66:                                for ( dcr0 = 0, dc = DC((P)p); dc; dc = NEXT(dc) ) {
        !            67:                                        NEXTDC(dcr0,dcr); cplx_conj((Obj)COEF(dc),(Obj *)&COEF(dcr));
        !            68:                                        DEG(dcr) = DEG(dc);
        !            69:                                }
        !            70:                                NEXT(dcr) = 0; MKP(VR((P)p),dcr0,t); *r = (Obj)t;
        !            71:                                break;
        !            72:                        default:
        !            73:                                error("cplx_conj : not implemented"); break;
        !            74:                }
        !            75: }
        !            76:
        !            77: void cplx_real(p,r)
        !            78: Obj p;
        !            79: Obj *r;
        !            80: {
        !            81:        DCP dc,dcr,dcr0;
        !            82:        P t;
        !            83:
        !            84:        if (  !p )
        !            85:                *r = 0;
        !            86:        else
        !            87:                switch ( OID(p) ) {
        !            88:                        case O_N:
        !            89:                                if ( NID((Num)p) <= N_B )
        !            90:                                        *r = p;
        !            91:                                else
        !            92:                                        *r = (Obj)((C)p)->r;
        !            93:                                break;
        !            94:                        case O_P:
        !            95:                                for ( dcr0 = 0, dc = DC((P)p); dc; dc = NEXT(dc) ) {
        !            96:                                        cplx_real((Obj)COEF(dc),(Obj *)&t);
        !            97:                                        if ( t ) {
        !            98:                                                NEXTDC(dcr0,dcr); DEG(dcr) = DEG(dc); COEF(dcr) = t;
        !            99:                                        }
        !           100:                                }
        !           101:                                if ( !dcr0 )
        !           102:                                        *r = 0;
        !           103:                                else {
        !           104:                                        NEXT(dcr) = 0; MKP(VR((P)p),dcr0,t); *r = (Obj)t;
        !           105:                                }
        !           106:                                break;
        !           107:                        default:
        !           108:                                error("cplx_real : not implemented"); break;
        !           109:                }
        !           110: }
        !           111:
        !           112: void cplx_imag(p,r)
        !           113: Obj p;
        !           114: Obj *r;
        !           115: {
        !           116:        DCP dc,dcr,dcr0;
        !           117:        P t;
        !           118:
        !           119:        if (  !p )
        !           120:                *r = 0;
        !           121:        else
        !           122:                switch ( OID(p) ) {
        !           123:                        case O_N:
        !           124:                                if ( NID((Num)p) <= N_B )
        !           125:                                        *r = 0;
        !           126:                                else
        !           127:                                        *r = (Obj)((C)p)->i;
        !           128:                                break;
        !           129:                        case O_P:
        !           130:                                for ( dcr0 = 0, dc = DC((P)p); dc; dc = NEXT(dc) ) {
        !           131:                                        cplx_imag((Obj)COEF(dc),(Obj *)&t);
        !           132:                                        if ( t ) {
        !           133:                                                NEXTDC(dcr0,dcr); DEG(dcr) = DEG(dc); COEF(dcr) = t;
        !           134:                                        }
        !           135:                                }
        !           136:                                if ( !dcr0 )
        !           137:                                        *r = 0;
        !           138:                                else {
        !           139:                                        NEXT(dcr) = 0; MKP(VR((P)p),dcr0,t); *r = (Obj)t;
        !           140:                                }
        !           141:                                break;
        !           142:                        default:
        !           143:                                error("cplx_imag : not implemented"); break;
        !           144:                }
        !           145: }