Return to katsura CVS log

Up to [local] / OpenXM_contrib2 / asir2000 / lib

Annotation of OpenXM_contrib2/asir2000/lib/katsura, Revision 1.1.1.1

1.1       noro        1: /* $OpenXM: OpenXM/src/asir99/lib/katsura,v 1.1.1.1 1999/11/10 08:12:30 noro Exp $ */
                      2: def mku(L,N)
                      3: {
                      4:        if (L < 0)
                      5:                L = -L;
                      6:        if (L <= N)
                      7:                return strtov("u"+rtostr(L));
                      8:        return 0;
                      9: }
                     10:
                     11: /* generate katsura systems */
                     12:
                     13: def katsura(N)
                     14: {
                     15:        R = [];
                     16:        for(L = 0; L < N; L++) {
                     17:                for(A = -mku(L,N),I = -N; I <= N; I++)
                     18:                        A += mku(L-I,N)*mku(I,N);
                     19:                R = cons(A,R);
                     20:        }
                     21:        for(A = 0, L = -N; L <= N; L++)
                     22:                A += mku(L,N);
                     23:        R = cons(A-1,R);
                     24:        return R;
                     25: }
                     26:
                     27: /* generate homogenized katsura systems */
                     28:
                     29: def hkatsura(N)
                     30: {
                     31:        R = [];
                     32:        for(L = 0; L < N; L++) {
                     33:                for(A = -mku(L,N)*t,I = -N; I <= N; I++)
                     34:                        A += mku(L-I,N)*mku(I,N);
                     35:                R = cons(A,R);
                     36:        }
                     37:        for(A = 0, L = -N; L <= N; L++)
                     38:                A += mku(L,N);
                     39:        R = cons(A-t,R);
                     40:        return R;
                     41: }
                     42: end$