summaryrefslogtreecommitdiff
path: root/langevin.cpp
diff options
context:
space:
mode:
Diffstat (limited to 'langevin.cpp')
-rw-r--r--langevin.cpp8
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);