diff --git a/include/teqp/algorithms/critical_tracing.hpp b/include/teqp/algorithms/critical_tracing.hpp index bae1a547962e7f264f0b9d8e85014238c1514a57..752fc3dd73ae55b1c090145b58fa427706b633bd 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; }