From 8ed9b01f00d638709d905c105d8929af89bea135 Mon Sep 17 00:00:00 2001
From: Ian Bell <ian.bell@nist.gov>
Date: Fri, 4 Nov 2022 12:27:01 -0400
Subject: [PATCH] Fix bug when alignment vector is provided but empty
---
include/teqp/algorithms/critical_tracing.hpp | 2 +-
1 file changed, 1 insertion(+), 1 deletion(-)
diff --git a/include/teqp/algorithms/critical_tracing.hpp b/include/teqp/algorithms/critical_tracing.hpp
index bae1a54..752fc3d 100644
--- a/include/teqp/algorithms/critical_tracing.hpp
+++ b/include/teqp/algorithms/critical_tracing.hpp
@@ -109,7 +109,7 @@ struct CriticalTracing {
else {
throw std::invalid_argument("More than one non-zero concentration value found; not currently supported");
}
- if (alignment_v0 && ed.eigenvectorscols.col(0).matrix().dot(alignment_v0.value().matrix()) < 0) {
+ if (alignment_v0 && alignment_v0.value().size() > 0 && ed.eigenvectorscols.col(0).matrix().dot(alignment_v0.value().matrix()) < 0) {
ed.eigenvectorscols.col(0) *= -1;
}
--
GitLab