Merge pull request #226 from ReactionMechanismGenerator/fix_transitorysens_interfaces

fix bug with normalization of transitory sensitivities with interfaces

Author: mjohnson541

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