Merge pull request #291 from sandialabs/bugfix_eigensolver

Change eigensolver initial iterate (bugfix)
sandialabs · Oct 17, 2022 · 3f6cc9a · 3f6cc9a
2 parents ae7dedf + c9effa6
commit 3f6cc9a
Show file tree

Hide file tree

Showing 2 changed files with 2 additions and 2 deletions.
diff --git a/cmake/Compadre_Version.txt b/cmake/Compadre_Version.txt
@@ -1 +1 @@
-1.5.2
+1.5.3
diff --git a/src/Compadre_LinearAlgebra_Definitions.hpp b/src/Compadre_LinearAlgebra_Definitions.hpp
@@ -18,7 +18,7 @@ void largestTwoEigenvectorsThreeByThreeSymmetric(const member_type& teamMember,
         double eigenvalue_relative_tolerance = 1e-6; // TODO: use something smaller, but really anything close is acceptable for this manifold
 
 
-        double v[3] = {1,1,1};
+        double v[3] = {rand_gen.drand(maxRange),rand_gen.drand(maxRange),rand_gen.drand(maxRange)};
         double v_old[3] = {v[0], v[1], v[2]};
 
         double error = 1;