Annotation of OpenXM_contrib2/asir2018/builtin/reduct.c, Revision 1.1
1.1 ! noro 1: /*
! 2: * Copyright (c) 1994-2000 FUJITSU LABORATORIES LIMITED
! 3: * All rights reserved.
! 4: *
! 5: * FUJITSU LABORATORIES LIMITED ("FLL") hereby grants you a limited,
! 6: * non-exclusive and royalty-free license to use, copy, modify and
! 7: * redistribute, solely for non-commercial and non-profit purposes, the
! 8: * computer program, "Risa/Asir" ("SOFTWARE"), subject to the terms and
! 9: * conditions of this Agreement. For the avoidance of doubt, you acquire
! 10: * only a limited right to use the SOFTWARE hereunder, and FLL or any
! 11: * third party developer retains all rights, including but not limited to
! 12: * copyrights, in and to the SOFTWARE.
! 13: *
! 14: * (1) FLL does not grant you a license in any way for commercial
! 15: * purposes. You may use the SOFTWARE only for non-commercial and
! 16: * non-profit purposes only, such as academic, research and internal
! 17: * business use.
! 18: * (2) The SOFTWARE is protected by the Copyright Law of Japan and
! 19: * international copyright treaties. If you make copies of the SOFTWARE,
! 20: * with or without modification, as permitted hereunder, you shall affix
! 21: * to all such copies of the SOFTWARE the above copyright notice.
! 22: * (3) An explicit reference to this SOFTWARE and its copyright owner
! 23: * shall be made on your publication or presentation in any form of the
! 24: * results obtained by use of the SOFTWARE.
! 25: * (4) In the event that you modify the SOFTWARE, you shall notify FLL by
! 26: * e-mail at risa-admin@sec.flab.fujitsu.co.jp of the detailed specification
! 27: * for such modification or the source code of the modified part of the
! 28: * SOFTWARE.
! 29: *
! 30: * THE SOFTWARE IS PROVIDED AS IS WITHOUT ANY WARRANTY OF ANY KIND. FLL
! 31: * MAKES ABSOLUTELY NO WARRANTIES, EXPRESSED, IMPLIED OR STATUTORY, AND
! 32: * EXPRESSLY DISCLAIMS ANY IMPLIED WARRANTY OF MERCHANTABILITY, FITNESS
! 33: * FOR A PARTICULAR PURPOSE OR NONINFRINGEMENT OF THIRD PARTIES'
! 34: * RIGHTS. NO FLL DEALER, AGENT, EMPLOYEES IS AUTHORIZED TO MAKE ANY
! 35: * MODIFICATIONS, EXTENSIONS, OR ADDITIONS TO THIS WARRANTY.
! 36: * UNDER NO CIRCUMSTANCES AND UNDER NO LEGAL THEORY, TORT, CONTRACT,
! 37: * OR OTHERWISE, SHALL FLL BE LIABLE TO YOU OR ANY OTHER PERSON FOR ANY
! 38: * DIRECT, INDIRECT, SPECIAL, INCIDENTAL, PUNITIVE OR CONSEQUENTIAL
! 39: * DAMAGES OF ANY CHARACTER, INCLUDING, WITHOUT LIMITATION, DAMAGES
! 40: * ARISING OUT OF OR RELATING TO THE SOFTWARE OR THIS AGREEMENT, DAMAGES
! 41: * FOR LOSS OF GOODWILL, WORK STOPPAGE, OR LOSS OF DATA, OR FOR ANY
! 42: * DAMAGES, EVEN IF FLL SHALL HAVE BEEN INFORMED OF THE POSSIBILITY OF
! 43: * SUCH DAMAGES, OR FOR ANY CLAIM BY ANY OTHER PARTY. EVEN IF A PART
! 44: * OF THE SOFTWARE HAS BEEN DEVELOPED BY A THIRD PARTY, THE THIRD PARTY
! 45: * DEVELOPER SHALL HAVE NO LIABILITY IN CONNECTION WITH THE USE,
! 46: * PERFORMANCE OR NON-PERFORMANCE OF THE SOFTWARE.
! 47: *
! 48: * $OpenXM$
! 49: */
! 50: #include "ca.h"
! 51: #include "parse.h"
! 52:
! 53: void Pred(), Predc(), Pprim();
! 54:
! 55: struct ftab reduct_tab[] = {
! 56: {"red",Pred,1},
! 57: {"redc",Predc,2},
! 58: {"prim",Pprim,-2},
! 59: {0,0,0},
! 60: };
! 61:
! 62: void Pred(NODE arg,Obj *rp)
! 63: {
! 64: Obj a,b,t;
! 65: LIST l;
! 66: V v;
! 67: int row,col,len;
! 68: VECT vect;
! 69: MAT mat;
! 70: int i,j;
! 71: NODE n0,n,nd;
! 72: struct oNODE arg0;
! 73: MP m,mp,mp0;
! 74: DP d;
! 75:
! 76: a = (Obj)ARG0(arg);
! 77: if ( !a ) {
! 78: *rp = 0;
! 79: return;
! 80: }
! 81: switch ( OID(a) ) {
! 82: case O_N: case O_P: case O_R:
! 83: reductr(CO,(Obj)ARG0(arg),rp);
! 84: break;
! 85: case O_LIST:
! 86: n0 = 0;
! 87: for ( nd = BDY((LIST)a); nd; nd = NEXT(nd) ) {
! 88: NEXTNODE(n0,n);
! 89: arg0.body = (pointer)BDY(nd);
! 90: arg0.next = 0;
! 91: Pred(&arg0,&b);
! 92: BDY(n) = (pointer)b;
! 93: }
! 94: if ( n0 )
! 95: NEXT(n) = 0;
! 96: MKLIST(l,n0);
! 97: *rp = (Obj)l;
! 98: break;
! 99: case O_VECT:
! 100: len = ((VECT)a)->len;
! 101: MKVECT(vect,len);
! 102: for ( i = 0; i < len; i++ ) {
! 103: arg0.body = (pointer)BDY((VECT)a)[i];
! 104: arg0.next = 0;
! 105: Pred(&arg0,&b);
! 106: BDY(vect)[i] = (pointer)b;
! 107: }
! 108: *rp = (Obj)vect;
! 109: break;
! 110: case O_MAT:
! 111: row = ((MAT)a)->row;
! 112: col = ((MAT)a)->col;
! 113: MKMAT(mat,row,col);
! 114: for ( i = 0; i < row; i++ )
! 115: for ( j = 0; j < col; j++ ) {
! 116: arg0.body = (pointer)BDY((MAT)a)[i][j];
! 117: arg0.next = 0;
! 118: Pred(&arg0,&b);
! 119: BDY(mat)[i][j] = (pointer)b;
! 120: }
! 121: *rp = (Obj)mat;
! 122: break;
! 123: case O_DP:
! 124: mp0 = 0;
! 125: for ( m = BDY((DP)a); m; m = NEXT(m) ) {
! 126: arg0.body = (pointer)C(m);
! 127: arg0.next = 0;
! 128: Pred(&arg0,&b);
! 129: if ( b ) {
! 130: NEXTMP(mp0,mp);
! 131: C(mp) = b;
! 132: mp->dl = m->dl;
! 133: }
! 134: }
! 135: if ( mp0 ) {
! 136: MKDP(NV((DP)a),mp0,d);
! 137: d->sugar = ((DP)a)->sugar;
! 138: *rp = (Obj)d;
! 139: } else
! 140: *rp = 0;
! 141:
! 142: break;
! 143: default:
! 144: error("red : invalid argument");
! 145: }
! 146: }
! 147:
! 148: void Predc(NODE arg,P *rp)
! 149: {
! 150: asir_assert(ARG0(arg),O_P,"redc");
! 151: asir_assert(ARG1(arg),O_P,"redc");
! 152: remsdcp(CO,(P)ARG0(arg),(P)ARG1(arg),rp);
! 153: }
! 154:
! 155: void Pprim(NODE arg,P *rp)
! 156: {
! 157: P t,p,p1,r;
! 158: V v;
! 159: VL vl;
! 160:
! 161: asir_assert(ARG0(arg),O_P,"prim");
! 162: p = (P)ARG0(arg);
! 163: if ( NUM(p) )
! 164: *rp = (P)ONE;
! 165: else {
! 166: if ( argc(arg) == 2 ) {
! 167: v = VR((P)ARG1(arg));
! 168: change_mvar(CO,p,v,&p1);
! 169: if ( VR(p1) != v ) {
! 170: *rp = (P)ONE; return;
! 171: } else {
! 172: reordvar(CO,v,&vl); pcp(vl,p1,&r,&t);
! 173: restore_mvar(CO,r,v,rp);
! 174: }
! 175: } else
! 176: pcp(CO,p,rp,&t);
! 177: }
! 178: }
FreeBSD-CVSweb <freebsd-cvsweb@FreeBSD.org>