Return to lineq-4.rtf CVS log

Up to [local] / OpenXM / src / cfep / tests

Annotation of OpenXM/src/cfep/tests/lineq-4.rtf, Revision 1.1

1.1     ! takayama    1: {\rtf1\mac\ansicpg10001\cocoartf824\cocoasubrtf330
        !             2: {\fonttbl\f0\fswiss\fcharset77 Helvetica;\f1\fnil\fcharset78 HiraKakuPro-W3;\f2\fnil\fcharset77 Monaco;
        !             3: }
        !             4: {\colortbl;\red255\green255\blue255;\red106\green117\blue255;}
        !             5: \paperw11900\paperh16840\margl1440\margr1440\vieww14600\viewh10420\viewkind0
        !             6: \pard\tx560\tx1120\tx1680\tx2240\tx2800\tx3360\tx3920\tx4480\tx5040\tx5600\tx6160\tx6720\ql\qnatural\pardirnatural
        !             7:
        !             8: \f0\fs24 \cf0 /*
        !             9: \f1 \'98\'41\'97\'a7\'95\'fb\'92\'f6\'8e\'ae\'82\'c9\'82\'c2\'82\'a2\'82\'c4\'82\'cc\'83\'76\'83\'8d\'83\'4f\'83\'89\'83\'80\'82\'f0\'8f\'91\'82\'a2\'82\'c4\'83\'65\'83\'58\'83\'67\'82\'c9\'82\'b7\'82\'e9
        !            10: \f0 . \
        !            11: */\
        !            12: \pard\tx560\tx1120\tx1680\tx2240\tx2800\tx3360\tx3920\tx4480\tx5040\tx5600\tx6160\tx6720\ql\qnatural\pardirnatural
        !            13:
        !            14: \f2\fs20 \cf0 def gaussE(A,Y) \{\
        !            15:   N = size(A)[0];\
        !            16:   for (K=1; K<=N-1; K++) \{\
        !            17:     for (I=K+1; I<=N; I++) \{\
        !            18:       P = A[I-1][K-1]/A[K-1][K-1];\
        !            19:       for (J=K; J<=N; J++) \{\
        !            20:         A[I-1][J-1] = A[I-1][J-1]-P*A[K-1][J-1];\
        !            21:       \}\
        !            22:       Y[I-1] = Y[I-1]-P*Y[K-1];\
        !            23:     \} \
        !            24:     print(A); print(Y); print("----------------");\
        !            25:   \}\
        !            26: \}\
        !            27: \
        !            28: def gtest1() \{\
        !            29:   A=newmat(3,3,\cf2 [[1,-1,2],[-1,2,-3],[3,1,1]]\cf0 );\
        !            30:   Y=newvect(3,[5,-6,8]);\
        !            31:   print(A); print(Y); print("----------------");\
        !            32:   gaussE(A,Y);\
        !            33:   return [A,Y];\
        !            34: \}\
        !            35: \
        !            36: def laplacian(N) \{\
        !            37:    A = newmat(N,N);\
        !            38:    for (I=0; I<N; I++) \{\
        !            39:         if (I > 0) A[I][I-1] = 1;\
        !            40:        A[I][I] = -2;\
        !            41:           if  (l < N-1) A[I][I+1] = 1;\
        !            42:    \}\
        !            43:    return A;\
        !            44: \}\
        !            45: [gtest1(), laplacian(3)];\
        !            46: }