Почему `eigs_gen` от броненосца дает разные фазы при наименьшем собственном значении при многократном запуске ⇐ C++

[Anonymous]

Я использую функцию eigs_gen из библиотеки C++ Armadillo, чтобы найти основное состояние квантово-физического гамильтониана, хранящегося в виде сложной разреженной матрицы sp_cx_mat H;. Это работает, но результаты различаются случайным образом каждый раз, когда я запускаю eigs_gen. Рассмотрим этот пример, где я вычисляю собственные значения H дважды:

#include #include #include #include #include #include #include использование пространства имен std; использование пространства имен Arma; интервал основной() { std::cout