diff --git a/p_space.jl b/p_space.jl index 2e5e2de..dfb23b8 100644 --- a/p_space.jl +++ b/p_space.jl @@ -26,9 +26,9 @@ end get_V_matrix(V_pq, p, w) = V_pq.(p, transpose(p)) .* transpose(w) -get_K_matrix(p, μ) = collect(Diagonal(p.*p ./ (2*μ))) +get_T_matrix(p, μ) = collect(Diagonal(p.*p ./ (2*μ))) -get_H_matrix(V_pq, p, w, μ=0.5) = get_K_matrix(p, μ) + get_V_matrix(V_pq, p, w) +get_H_matrix(V_pq, p, w, μ=0.5) = get_T_matrix(p, μ) + get_V_matrix(V_pq, p, w) function identify_pole_i(p, evals, μ=0.5) mesh_Es = (p.*p) ./ (2*μ)