From ce468e8ed61667cadfb0c05dd8755318b37fa8be Mon Sep 17 00:00:00 2001 From: Jukka Aho Date: Sun, 18 Mar 2018 17:37:17 +0700 Subject: [PATCH] Fix bug in eigenvalue solver (#191) This closes issue #190. --- src/solvers_modal.jl | 5 +++-- 1 file changed, 3 insertions(+), 2 deletions(-) diff --git a/src/solvers_modal.jl b/src/solvers_modal.jl index 450274e..06a47b0 100644 --- a/src/solvers_modal.jl +++ b/src/solvers_modal.jl @@ -12,6 +12,7 @@ julia> push!(solver, problems...) julia> solver() """ type Modal <: AbstractSolver + time :: Float64 geometric_stiffness :: Bool eigvals :: Vector eigvecs :: Matrix @@ -27,7 +28,7 @@ type Modal <: AbstractSolver end function Modal(nev=10, which=:SM) - solver = Modal(false, [], Matrix{Float64}(0,0), nev, which, + solver = Modal(0.0, false, [], Matrix{Float64}(0,0), nev, which, false, [], true, true, false, false, 0.0) end @@ -313,7 +314,7 @@ function update_xdmf!(solver::Solver{Modal}) xdmf = first(results_writers) - @timeit "fetch geometry" X_ = solver("geometry", solver.time) + @timeit "fetch geometry" X_ = solver("geometry", solver.properties.time) node_ids = keys(X_) @timeit "create node permutation" P = Dict(j=>i for (i, j) in enumerate(node_ids)) nnodes = length(X_)