diff options
author | Jaron Kent-Dobias <jaron@kent-dobias.com> | 2021-01-05 11:14:16 +0100 |
---|---|---|
committer | Jaron Kent-Dobias <jaron@kent-dobias.com> | 2021-01-05 11:14:16 +0100 |
commit | 3541fa50708470789592ce9ac1f3baf1ee35943f (patch) | |
tree | 297f5e0cebb25906458e0560103654ffb99d1e15 | |
parent | 1c79cff86e5cfae48e00184842101a9678b89903 (diff) | |
download | code-3541fa50708470789592ce9ac1f3baf1ee35943f.tar.gz code-3541fa50708470789592ce9ac1f3baf1ee35943f.tar.bz2 code-3541fa50708470789592ce9ac1f3baf1ee35943f.zip |
Cleaned up some notition.
-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); |