From 072e80f0d7181167a08f325c09b0d6dee5858e69 Mon Sep 17 00:00:00 2001
From: Michiel Cottaar <michiel.cottaar@ndcn.ox.ac.uk>
Date: Tue, 30 Jan 2024 15:03:36 +0000
Subject: [PATCH] use permutedims instead of adjoint
it does not try to take complex conjugates
---
src/pathways.jl | 3 ++-
1 file changed, 2 insertions(+), 1 deletion(-)
diff --git a/src/pathways.jl b/src/pathways.jl
index 1a9730b..aee7119 100644
--- a/src/pathways.jl
+++ b/src/pathways.jl
@@ -348,8 +348,9 @@ function update_gradient_tracker_till_time!(walker::PathwayWalker, key::Tuple, n
end
function update_gradient_tracker_till_time!(gradient_tracker::GradientTracker, new_time::VariableType)
+ @show gradient_tracker.qvec
gradient_tracker.bmat += (
- (gradient_tracker.qvec .* gradient_tracker.qvec') .*
+ (gradient_tracker.qvec .* permutedims(gradient_tracker.qvec)) .*
(new_time - gradient_tracker.last_gradient_time)
)
gradient_tracker.last_gradient_time = new_time
--
GitLab