From 2e5393a4b692341e9c68c39aa25904133c549d18 Mon Sep 17 00:00:00 2001
From: Ryan Latture <ryan.latture@gmail.com>
Date: Tue, 5 Jul 2016 11:17:33 -0700
Subject: [PATCH] Added prune step before compressing global stiffness matrix.

---
 src/threed_beam_fea.cpp | 3 ++-
 1 file changed, 2 insertions(+), 1 deletion(-)

diff --git a/src/threed_beam_fea.cpp b/src/threed_beam_fea.cpp
index 2fa1158..5b33675 100644
--- a/src/threed_beam_fea.cpp
+++ b/src/threed_beam_fea.cpp
@@ -423,8 +423,9 @@ namespace fea {
         if (forces.size() > 0) {
             loadForces(force_vec, forces);
         }
-
+ 
         // compress global stiffness matrix since all non-zero values have been added.
+        Kg.prune(1.e-14);
         Kg.makeCompressed();
 
         // initialize solver based on whether MKL should be used