Return to fctrdata CVS log

Up to [local] / OpenXM_contrib2 / asir2000 / lib

Annotation of OpenXM_contrib2/asir2000/lib/fctrdata, Revision 1.3

1.2       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
1.3     ! noro       26:  * e-mail at risa-admin@sec.flab.fujitsu.co.jp of the detailed specification
1.2       noro       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:  *
1.3     ! noro       48:  * $OpenXM: OpenXM_contrib2/asir2000/lib/fctrdata,v 1.2 2000/08/21 08:31:41 noro Exp $
1.2       noro       49: */
1.1       noro       50: /* examples for factorization
                     51:  * Wang[1] - Wang[19] : factor.tst of REDUCE
                     52:  * Sq[1] - Sq[10] : examples for square free factorization
                     53:  * Alg[1] - Alg[8] : examples for algebraic factorization
                     54:  */
                     55:
                     56: Wang = newvect(20)$
                     57: Sq = newvect(20)$
                     58: Alg = newvect(20)$
                     59:
                     60: A = x*y+z+10$ B = x*z+y+30$ C = x+y*z+20$
                     61: Wang[1] = A*B*C$
                     62:
                     63: A=x^3*(z+y)+z-11$ B=x^2*(z^2+y^2)+y+90$
                     64: Wang[2] = (A*B)$
                     65:
                     66: A = x^3*y^2+x*z^4+x+z$ B = x^3+x*y*z+y^2+y*z^3$
                     67: Wang[3] = (A*B) $
                     68:
                     69: A = x^2*z+y^4*z^2+5$ B = x*y^3+z^2$ C = -x^3*y+z^2+3$ D = x^3*y^4+z^2$
                     70: Wang[4] = (A*B*C*D) $
                     71:
                     72: A = 3*u^2*x^3*y^4*z+x*z^2+y^2*z^2+19*y^2$ B = u^2*y^4*z^2+x^2*z+5$ C = u^2+x^3*y^4+z^2$
                     73: Wang[5] = (A*B*C) $
                     74:
                     75: A = w^4*x^5*y^6-w^4*z^3+w^2*x^3*y+x*y^2*z^2$
                     76: B = w^4*z^6-w^3*x^3*y-w^2*x^2*y^2*z^2+x^5*z-x^4*y^2+y^2*z^3$
                     77: C = -x^5*z^3+x^2*y^3+y*z$
                     78: Wang[6] = (A*B*C) $
                     79:
                     80: A = x+y+z-2$ B = x+y+z-2$ C = x+y+z-3$ D = x+y+z-3$ E = x+y+z-3$
                     81: Wang[7] = (A*B*C*D*E) $
                     82:
                     83: A = -z^31-w^12*z^20+y^18-y^14+x^2*y^2+x^21+w^2$
                     84: B = -15*y^2*z^16+29*w^4*x^12*z^3+21*x^3*z^2+3*w^15*y^20$
                     85: Wang[8] = (A*B) $
                     86:
                     87: A = 18*u^2*w^3*x*z^2+10*u^2*w*x*y^3+15*u*z^2+6*w^2*y^3*z^2$
                     88: B = u^4*x*z^2$
                     89: C = 25*u^2*w^3*y*z^4+32*u^2*w^4*y^4*z^3-48*u^2*x^2*y^3*z^3-
                     90: 2*u^2*w*x^2*y^2+44*u*w*x*y^4*z^4-8*u*w*x^3*z^4+4*w^2*x+
                     91: 11*w^2*x^3*y+12*y^3*z^2$
                     92: Wang[9] = (A*B*C) $
                     93:
                     94: A = 31*u^2*x*z+35*w^2*y^2+40*w*x^2+6*x*y$
                     95: B = 42*u^2*w^2*y^2+47*u^2*w^2*z+22*u^2*w^2+9*u^2*w*x^2+21
                     96: *u^2*w*x*y*z+37*u^2*y^2*z+u^2*w^2*x*y^2*z^2+8*u^2*w^2
                     97: *z^2+24*u^2*w*x*y^2*z^2+24*u^2*x^2*y*z^2+12*u^2*x*y^2
                     98: *z^2+13*u*w^2*x^2*y^2+27*u*w^2*x^2*y+39*u*w*x*z+43*u*
                     99: x^2*y+44*u*w^2* z^2+37*w^2*x*y+29*w^2*y^2+31*w^2*y*z^2
                    100: +12*w*x^2*y*z+43*w*x*y*z^2+22*x*y^2+23*x*y*z+24*x*y+41*y^2*z$
                    101: Wang[10] = (A*B) $
                    102:
                    103: A = -36*u^2*w^3*x*y*z^3-31*u^2*w^3*y^2+20*u^2*w^2*x^2*y^2
                    104: *z^2-36*u^2*w*x*y^3*z+46*u^2*w*x+9*u^2*y^2-36*u*w^2*y^3
                    105: +9*u*w*y^3-5*u*w*x^2*y^3+48*u*w*x^3*y^2*z+23*u*w*x^3*y^2
                    106: -43*u*x^3*y^3*z^3-46*u*x^3*y^2+29*w^3*x*y^3*z^2-
                    107: 14*w^3*x^3*y^3*z^2-45*x^3-8*x*y^2$
                    108: B = 13*u^3*w^2*x*y*z^3-4*u*x*y^2-w^3*z^3-47*x*y$
                    109: Wang[11] = (A*B) $
                    110:
                    111: A = x+y+z-3$ B = x+y+z-3$ C = x+y+z-3$
                    112: Wang[12] = (A*B*C) $
                    113:
                    114: A = 2*w*z+45*x^3-9*y^3-y^2+3*z^3$ B = w^2*z^3-w^2+47*x*y$
                    115: Wang[13] = (A*B) $
                    116:
                    117: A = 18*x^4*y^5+41*x^4*y^2-37*x^4+26*x^3*y^4+38*x^2*y^4-29*x^2*y^3-22*y^5$
                    118: B = 33*x^5*y^6-22*x^4+35*x^3*y+11*y^2$
                    119: Wang[14] = (A*B) $
                    120:
                    121: A = 12*w^2*x*y*z^3-w^2*z^3+w^2-29*x-3*x*y^2$
                    122: B = 14*w^2*y^2+2*w*z+18*x^3*y-8*x*y^2-y^2+3*z^3$
                    123: C = x^6*y^3*z^2$
                    124: Wang[15] = (A*B*C) $
                    125:
                    126: A1 = 4096*x^10+8192*x^9-3008*x^8-30848*x^7+21056*x^6+146496*x^5
                    127:        -221360*x^4+1232*x^3+144464*x^2-78488*x+11993$
                    128: A2 = 4096*x^10+8192*x^9+1600*x^8-20608*x^7+20032*x^6+87360*x^5
                    129:        -105904*x^4+18544*x^3+11888*x^2-3416*x+41$
                    130: A3 = 8192*x^10+12288*x^9+66560*x^8-22528*x^7-138240*x^6+572928*x^5
                    131:        -90496*x^4-356032*x^3+113032*x^2+23420*x-8179$
                    132: A4 = 8192*x^10+20480*x^9+58368*x^8-161792*x^7+198656*x^6+199680*x^5
                    133:        -414848*x^4-4160*x^3+171816*x^2-48556*x+469$
                    134: Wang[16] = (A1*A2*A3*A4)$
                    135:
                    136: A1 = x^25-25*x^20-3500*x^15-57500*x^10+21875*x^5-3125$
                    137: Wang[17] = (A1)$
                    138:
                    139: A1 = x^18+9*x^17+45*x^16+126*x^15+189*x^14+27*x^13-540*x^12-1215*x^11
                    140:        +1377*x^10+15444*x^9+46899*x^8+90153*x^7+133893*x^6+125388*x^5
                    141:        +29160*x^4-32076*x^3+26244*x^2-8748*x+2916$
                    142: Wang[18] = (A1)$
                    143:
                    144: A1 = x^16+4*x^12-16*x^11+80*x^9+2*x^8+160*x^7+128*x^6-160*x^5+28*x^4-48*x^3+128*x^2-16*x+1$
                    145: A2 = x^16+4*x^12+16*x^11-80*x^9+2*x^8-160*x^7+128*x^6+160*x^5+28*x^4+48*x^3+128*x^2+16*x+1$
                    146: Wang[19] = (A1*A2)$
                    147:
                    148: Sq[1]= ((3*x^2+x+1)*y^2+2*x*y+x^2+x)$
                    149: Sq[2]= ((y^4+x^3)*(y^3+x^2)^2*(y^2+x)^3)$
                    150: Sq[3]= ((z+y+x+1)*(z-y+x+4)^2*(z-2*y+x+7)^3)$
                    151: Sq[4]= ((y+x^2+5)*(6*y+2*x^3+31)^3*(x*y^3+8*y-x)^5)$
                    152: Sq[5]= ((z^2+x*y*z+x^2)^2*(3*z^2+(y^2+x)*z-4)^3*(z^3+z^2+(x-1)*y)^4)$
                    153: Sq[6]= ((3*y^4+x+5)*(6*y^2+2*x+31)^2*((x^2+x+1)*y^3-y+x+8)^4*(x*y+8*y+x)^10)$
                    154: Sq[7]= ((x^4+351*x^3+27*x^2-31*x+1)^2*(6*x^5+251*x^3-372*x^2+15*x-323)^6)$
                    155: Sq[8]= ((x^4+351*x^3+27*x^2-31*x+1)^3*(6*x^5+251*x^3-372*x^2+15*x-323)^9)$
                    156: Sq[9]= ((x^3+75*x^2+68*x+1)^3*(x^2+15*x+35)^6*(7*x^2+750*x+137)^9*(x+75)^12)$
                    157: Sq[10]= ((x^3+75*x^2+68*x+1)^5*(x^2+15*x+35)^10*(7*x^2+750*x+137)^15*(x+75)^20)$
                    158:
                    159: Alg[1]= (x^2-x+3)$
                    160: Alg[2]= (x^3+2)$
                    161: Alg[3]= (x^4-x+1)$
                    162: Alg[4]= (x^6+3*x^5+6*x^4+x^3-3*x^2+12*x+16)$
                    163: Alg[5]= (x^9-15*x^6-87*x^3-125)$
                    164: Alg[6]= (x^9-54)$
                    165: Alg[7]= (x^3-19)$
                    166: Alg[8]= (x^2+x+7)$
                    167:
                    168: end$