solver.py 文件源码

def solve_dynamic_eigen_model(system):
    from scipy import linalg as sl
    if not system._is_inited:
        system.calc_KG()
    system.calc_deleted_KG_matrix()
    system.calc_MG()
    system.calc_deleted_MG_matrix()

    w1,system.model = sl.eigh(system.KG_keeped,system.MG_keeped)
    system.w = np.sqrt(w1)
    T = 2*np.pi/system.w
    system.freq = 1/T