From be4035546e3e148efb5123b19c1208afcdae885e Mon Sep 17 00:00:00 2001 From: Matt Johnson Date: Wed, 26 Jul 2023 12:42:54 -0700 Subject: [PATCH] fix bug with normalization of transitory sensitivities with interfaces --- src/TransitorySensitivities.jl | 10 ++++++++++ 1 file changed, 10 insertions(+) diff --git a/src/TransitorySensitivities.jl b/src/TransitorySensitivities.jl index d37f0f076..85476b814 100644 --- a/src/TransitorySensitivities.jl +++ b/src/TransitorySensitivities.jl @@ -163,6 +163,11 @@ function normalizefulltransitorysensitivities!(dSdt,ssys::SystemSimulation,t) @views ns = y[sim.domain.indexes[1]:sim.domain.indexes[2]] @views dSdt[sim.domain.indexes[1]:sim.domain.indexes[2],:] ./= ns end + for (i,inter) in enumerate(ssys.interfaces) + if isa(inter, AbstractReactiveInternalInterface) + @views dSdt[:,inter.parameterindexes[1]:inter.parameterindexes[2]] .*= ssys.p[inter.parameterindexes[1]:inter.parameterindexes[2]]' + end + end return dSdt end @@ -205,6 +210,11 @@ function normalizeadjointtransitorysensitivities!(dSdt,ssys::SystemSimulation,t, dSdt ./= y[ind] end end + for (i,inter) in enumerate(ssys.interfaces) + if isa(inter, AbstractReactiveInternalInterface) + @views dSdt[inter.parameterindexes[1]:inter.parameterindexes[2]] .*= ssys.p[inter.parameterindexes[1]:inter.parameterindexes[2]] + end + end return dSdt end