fix bug with normalization of transitory sensitivities with interfaces

ReactionMechanismGenerator · Feb 4, 2024 · b3d12d8 · b3d12d8
1 parent a3b22d5
commit b3d12d8
Showing 1 changed file with 10 additions and 0 deletions.
diff --git 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