From fc9b890744c5e16099f971ebec45b2458f2541ca Mon Sep 17 00:00:00 2001 From: Frondelius Date: Tue, 22 Nov 2016 15:31:23 +0200 Subject: [PATCH] fine tuning --- geometry/3d_beam_nx/beam_3d.jl | 14 +++++--------- 1 file changed, 5 insertions(+), 9 deletions(-) diff --git a/geometry/3d_beam_nx/beam_3d.jl b/geometry/3d_beam_nx/beam_3d.jl index efac829..53b92b4 100644 --- a/geometry/3d_beam_nx/beam_3d.jl +++ b/geometry/3d_beam_nx/beam_3d.jl @@ -136,16 +136,12 @@ function lsq_fit(elements, field) SparseArrays.droptol!(A, 1.0e-6) nz = get_nonzero_rows(A) dropzeros!(A) - #try - F = ldltfact(A) + F = ldltfact(A) - x = zeros(size(b)...) - x[nz, :] = F \ b[nz, :] - nodal_values = Dict(i => vec(x[i,:]) for i=1:size(x,1)) - #catch - err("Problem with zeros or tolerances") - dump(nz) - #end + x = zeros(size(b)...) + x[nz, :] = F \ b[nz, :] + + nodal_values = Dict(i => vec(x[i,:]) for i=1:size(x,1)) return nodal_values end