diff options
| -rw-r--r-- | langevin.cpp | 8 |
1 files changed, 4 insertions, 4 deletions
diff --git a/langevin.cpp b/langevin.cpp index 7590d2b..6f26278 100644 --- a/langevin.cpp +++ b/langevin.cpp @@ -40,7 +40,7 @@ Vector initializeVector(unsigned N, Rng& r) { return z; } -std::tuple<std::complex<double>, Vector, Matrix> hamGradHess(const Tensor& J, const Vector& z) { +std::tuple<Scalar, Vector, Matrix> hamGradHess(const Tensor& J, const Vector& z) { Matrix Jz = contractDown(J, z); Vector Jzz = Jz * z; @@ -139,13 +139,13 @@ int main(int argc, char* argv[]) { } } - std::complex<double> κ(Rκ, Iκ); + Scalar κ(Rκ, Iκ); double σ = sqrt(factorial(p) / (2.0 * pow(N, p - 1))); Rng r; complex_normal_distribution<> d(0, σ, κ); - Tensor J = generateCouplings<std::complex<double>, p>(N, d, r.engine()); + Tensor J = generateCouplings<Scalar, p>(N, d, r.engine()); Vector z = initializeVector(N, r); std::function<bool(double, unsigned)> findSaddle = [δ](double W, unsigned) { @@ -155,7 +155,7 @@ int main(int argc, char* argv[]) { Vector zm = langevin(J, z, T, γ, findSaddle, r); - std::complex<double> H; + Scalar H; Vector dH; std::tie(H, dH, std::ignore) = hamGradHess(J, zm); |