コード例 #1
0
ファイル: err.py プロジェクト: zimoun/mtf
    }
    # },
    #   { 'name': 'B',
    #     'phys': 1,
    #     'union': 2,
    #     }
    # },
    #   { 'name': 'C',
    #     'phys': 4,
    #     'union': 3,
    # }
]

mtf = MultiTrace(kRef, meshname, dd)

Aw = mtf.A_weak_form()

At, X, J, iJ = mtf.tolinop()

shape = mtf.shape

A = 2.0 * At
A2 = A * iJ * A

Ce = 0.5 * J - At
Ci = 0.5 * J + At

Ce2 = Ce * iJ * Ce
Ci2 = Ci * iJ * Ci

x = np.random.rand(shape[0]) + 1j * np.random.rand(shape[0])