Fixed bug in mixed-spin energy evalutaion.

1 files changed, 1 insertions, 1 deletions

diff --git a/p-spin.hpp b/p-spin.hpp
index a5f6c9a..7a49bca 100644
--- a/p-spin.hpp
+++ b/p-spin.hpp
@@ -46,7 +46,7 @@ std::tuple<Scalar, Vector<Scalar>, Matrix<Scalar>, Tensor<Scalar, 3>> hamGradHes
   Scalar H = Jzzz / pBang;
 
   if constexpr (sizeof...(Js) > 0) {
-    auto [Hs, dHs, ddHs, dddHs] = hamGradHess(z, Js...);
+    auto [Hs, dHs, ddHs, dddHs] = hamGradHessHelper(z, Js...);
     H += Hs;
     dH += dHs;
     ddH += ddHs;