From d34b71e60ecb7475d65e9b7024e580b6c899ed7f Mon Sep 17 00:00:00 2001 From: VvanderHeiden Date: Fri, 6 Sep 2024 11:35:18 +0200 Subject: [PATCH] Make stiffness matrix square when elastic input is depth dependent. --- src/calculator/FaultSlip.m | 6 +++--- 1 file changed, 3 insertions(+), 3 deletions(-) diff --git a/src/calculator/FaultSlip.m b/src/calculator/FaultSlip.m index c0cfd7c..525d4b2 100644 --- a/src/calculator/FaultSlip.m +++ b/src/calculator/FaultSlip.m @@ -219,9 +219,9 @@ nx = length(L); Kline = -nx/(2*pi*(L(1)-L(end))) ./ ( [0:nx-1]'.^2-0.25 ); Ko = toeplitz(Kline); - % K = mu_II.*Ko; - K = Ko*mu_II; - K(and(K<0,K>(min(min(K)/10000)))) = 0; % set very small changes to 0 to avoid continued interactions of the two peaks to + % K = Ko*mu_II; % if stiffness not depth dependent + K = Ko .* mu_II'; % take depth dependent stiffness into account + K(and(K<0,K>(min(min(K)/10000)))) = 0; % set very small changes to 0 to avoid continued interactions of the two peaks to end