Commit e83671e1 authored by Michel Juillard's avatar Michel Juillard
Browse files

debugging get_epsilon

parent 3306088f
......@@ -50,6 +50,9 @@ function get_epsilonh!(epsilon::AbstractVector{U}, H::AbstractMatrix{U},
tmp2::AbstractVector{U}) where U <: AbstractFloat
copy!(tmp1, iFv)
mul!(tmp1, transpose(K), r, -1.0, 1.0)
@show epsilon
@show H
@show tmp1
mul!(epsilon, H, tmp1)
Supports Markdown
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment