diff --git a/src/pathways.jl b/src/pathways.jl index 1a9730bad8fefce3bb92c59f9163298577957bc5..aee7119080fbd6ae600cf55906ecb2774076bb27 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