Return to alpi CVS log

Up to [local] / OpenXM_contrib2 / asir2000 / lib

Annotation of OpenXM_contrib2/asir2000/lib/alpi, Revision 1.2

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
        !            26:  * e-mail at risa-admin@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: OpenXM_contrib2/asir2000/lib/alpi,v 1.1.1.1 1999/12/03 07:39:11 noro Exp $
        !            49: */
1.1       noro       50: A4 = [
                     51: x+y+z+t,
                     52: x*y+y*z+z*t+t*x,
                     53: x*y*z+y*z*t+z*t*x+t*x*y,
                     54: x*y*z*t-1
                     55: ]$
                     56:
                     57: A5 = [
                     58: x+y+z+t+u,
                     59: x*y+y*z+z*t+t*u+u*x,
                     60: x*y*z+y*z*t+z*t*u+t*u*x+u*x*y,
                     61: x*y*z*t+y*z*t*u+z*t*u*x+t*u*x*y+u*x*y*z,
                     62: x*y*z*t*u-1
                     63: ]$
                     64:
                     65: MA = [
                     66: a+b+c+d+e,
                     67: a*b+b*c+c*d+d*e+e*a,
                     68: a*b*c+b*c*d+c*d*e+d*e*a+e*a*b,
                     69: b*c*d+b*c*d*e+c*d*e*a+d*e*a*b+e*a*b*c,
                     70: a*b*c*d*e-1
                     71: ]$
                     72:
                     73: La = [
                     74: a^2*b*c+a*b^2*c+a*b*c^2+a*b*c+a*b+a*c+b*c,
                     75: a^2*b^2*c+a*b^2*c^2+a^2*b*c+a*b*c+b*c+a+c,
                     76: a^2*b^2*c^2+a^2*b^2*c+a*b^2*c+a*b*c+a*c+c+1
                     77: ]$
                     78:
                     79: Mo = [x*b-y*a,(x-1)*d-y*(c-1),b^2+a^2-r^2,(c-1)^2+d^2-s^2,(a-c)^2+(b-d)^2-t^2]$
                     80:
                     81: Ca = [
                     82: y^2*z+2*x*y*t-2*x-z,
                     83: -x^3*z+4*x*y^2*z+4*x^2*y*t+2*y^3*t+4*x^2-10*y^2+4*x*z-10*y*t+2,
                     84: 2*y*z*t+x*t^2-x-2*z,
                     85: -x*z^3+4*y*z^2*t+4*x*z*t^2+2*y*t^3+4*x*z+4*z^2-10*y*t-10*t^2+2
                     86: ]$
                     87:
                     88: P4 = [x*(y+z+t)-a,y*(x+z+t)-b,z*(x+y+t)-c,t*(x+y+z)-d]$
                     89:
                     90: M1 = [
                     91: x^2+y^2+z^2,
                     92: z*y*x,
                     93: y^2*x^2+z^2*x^2+z^2*y^2,
                     94: u^2+1/3*t^2,
                     95: u^3-t^2*u,
                     96: 2*x^2*u-y^2*u-z^2*u+y^2*t-z^2*t,
                     97: -y^2*x^2*u-z^2*x^2*u+2*z^2*y^2*u-y^2*x^2*t+z^2*x^2*t,
                     98: 2*x^2*u^2-y^2*u^2-z^2*u^2-2*y^2*t*u+2*z^2*t*u-2/3*x^2*t^2+1/3*y^2*t^2
                     99: +1/3*z^2*t^2,
                    100: -y^2*x^2*u^2-z^2*x^2*u^2+2*z^2*y^2*u^2+2*y^2*x^2*t*u-2*z^2*x^2*t*u
                    101:     +1/3*y^2*x^2*t^2+1/3*z^2*x^2*t^2-2/3*z^2*y^2*t^2,
                    102: -3*y^2*x^4*t*u^2+3*z^2*x^4*t*u^2+3*y^4*x^2*t*u^2-3*z^4*x^2*t*u^2
                    103:      -3*z^2*y^4*t*u^2+3*z^4*y^2*t*u^2+1/3*y^2*x^4*t^3-1/3*z^2*x^4*t^3
                    104:      -1/3*y^4*x^2*t^3+1/3*z^4*x^2*t^3+1/3*z^2*y^4*t^3-1/3*z^4*y^2*t^3
                    105: ]$
                    106:
                    107: M2 = [
                    108: y+u+v-1, z+t+2*u-3, y+t+2*v-1, x-y-z-t-u-v,
                    109: t*u*x^2-1569/31250*y*z^3, z*v-587/15625*y*t
                    110: ]$
                    111:
                    112: M3=[
                    113: x^2-x+2*y^2+2*z^2+2*t^2, 2*x*y+2*y*z+2*z*t-y,
                    114: 2*x*z+y^2+2*y*t-z, x+2*y+2*z+2*t-1
                    115: ]$
                    116:
                    117: Rb = [ x^31-x^6-x-y, x^8-z, x^10-t]$
                    118:
                    119: Va = [
                    120: -z^2*u+2*y*z*v-x*v^2-y^2*a+x*u*a+aA^3,
                    121: -z*t*u+y*t*v+y*z*w-x*v*w-y^2*b+x*u*b+bB^3,
                    122: -t^2*u+2*y*t*w-x*w^2-y^2*c+x*u*c+cC^3,
                    123: -z*t*v+z^2*w+y*t*a-x*w*a-y*z*b+x*v*b+dD^3,
                    124: -t^2*v+z*t*w+y*t*b-x*w*b-y*z*c+x*v*c+eE^3,
                    125: -t*v^2+z*v*w+t*u*a-y*w*a-z*u*b+y*v*b+fF^3,
                    126: -t*v*w+z*w^2+t*u*b-y*w*b-z*u*c+y*v*c+gG^3,
                    127: -t^2*a+2*z*t*b-x*b^2-z^2*c+x*a*c+hH^3,
                    128: -t*w*a+t*v*b+z*w*b-y*b^2-z*v*c+y*a*c+iI^3,
                    129: -w^2*a+2*v*w*b-u*b^2-v^2*c+u*a*c+jJ^3
                    130: ]$
                    131:
                    132: C2=[
                    133: x^3*y^2+4*x^2*y^2*z-x^2*y*z^2+288*x^2*y^2+207*x^2*y*z+1152*x*y^2*z+156*
                    134:     x*y*z^2+x*z^3-3456*x^2*y+20736*x*y^2+19008*x*y*z+82944*y^2*z+432*x*z^2-
                    135:     497664*x*y+62208*x*z+2985984*x,
                    136: y^3*t^3+4*y^3*t^2-y^2*z*t^2+4*y^2*t^3-48*y^2*t^2-5*y*z*t^2+108*y*z*t+
                    137:     z^2*t+144*z*t-1728*z,
                    138: x^2*z^2*t-4*x*z^2*t^2-z^3*t^2-x^3*z-156*x^2*z*t-207*x*z^2*t-1152*x*z*
                    139:     t^2-288*z^2*t^2-432*x^2*z-19008*x*z*t+3456*z^2*t-82944*x*t^2-20736*z*t^2-
                    140:     62208*x*z+497664*z*t-2985984*z,
                    141: y^3*t^3-x*y^2*t^2+4*y^3*t^2+4*y^2*t^3-5*x*y^2*t-48*y^2*t^2+x^2*y+108*x*
                    142:     y*t+144*x*y-1728*x
                    143: ]$
                    144:
                    145: A=nm(z^3+(432+156*y-x*y)*z^2+(4*x*y^2+207*x*y+1152*y^2+19008*y+62208+82944*y^2/x)*z
                    146: +(x*y+144*y-1728)^2)$
                    147: B=nm(red(t^3-(x-4*y+48)*t^2/(y+4)+(-5*y+108)*x*t/(y*(y+4))+x*(x*y+144*y-1728)/(y^2*(y+4))))$
                    148: C=subst(A,x,xx,y,yy,z,zz,xx,z,yy,t,zz,x)$
                    149: D=subst(B,x,xx,y,yy,t,tt,xx,z,yy,t,tt,y)$ Co=[A,B,C,D]$
                    150: OCo=[x,y,z,t]$
                    151:
                    152: IP1=[x^2*y*z^4-t,x^5*y^7-z^2*u,-x^3*z*v+y^2,-z^5*w+x*y^3]$
                    153: IP2=[-y^82*a+x^32*z^23,x^45-y^13*z^21*b,y^33*z^12-x^41*c,
                    154: -y^33*z^12*d+x^22,x^5*y^17*z^22*e-1,x*y*z*tT-1]$
                    155: OIP2=[tT,x,y,z,a,b,c,d,e]$
                    156: end$