316
316
# ## Other Tests ###
317
317
318
318
# Checks that identifiability can be assessed for coupled CRN/DAE systems.
319
+ # `remove_conserved = false` is used to remove info print statement from log.
319
320
let
320
321
rs = @reaction_network begin
321
322
@parameters k c1 c2
@@ -329,19 +330,20 @@ let
329
330
@unpack p, d, k, c1, c2 = rs
330
331
331
332
# Tests identifiability assessment when all unknowns are measured.
332
- gi_1 = assess_identifiability (rs; measured_quantities = [:X , :V , :C ], loglevel)
333
- li_1 = assess_local_identifiability (rs; measured_quantities = [:X , :V , :C ], loglevel)
334
- ifs_1 = find_identifiable_functions (rs; measured_quantities = [:X , :V , :C ], loglevel)
333
+ remove_conserved = false
334
+ gi_1 = assess_identifiability (rs; measured_quantities = [:X , :V , :C ], loglevel, remove_conserved)
335
+ li_1 = assess_local_identifiability (rs; measured_quantities = [:X , :V , :C ], loglevel, remove_conserved)
336
+ ifs_1 = find_identifiable_functions (rs; measured_quantities = [:X , :V , :C ], loglevel, remove_conserved)
335
337
@test sym_dict (gi_1) == Dict ([:X => :globally , :C => :globally , :V => :globally , :k => :globally ,
336
338
:c1 => :nonidentifiable , :c2 => :nonidentifiable , :p => :globally , :d => :globally ])
337
339
@test sym_dict (li_1) == Dict ([:X => 1 , :C => 1 , :V => 1 , :k => 1 , :c1 => 0 , :c2 => 0 , :p => 1 , :d => 1 ])
338
340
@test issetequal (ifs_1, [d, p, k, c1 + c2])
339
341
340
342
# Tests identifiability assessment when only variables are measured.
341
343
# Checks that a parameter in an equation can be set as known.
342
- gi_2 = assess_identifiability (rs; measured_quantities = [:V , :C ], known_p = [:c1 ], loglevel)
343
- li_2 = assess_local_identifiability (rs; measured_quantities = [:V , :C ], known_p = [:c1 ], loglevel)
344
- ifs_2 = find_identifiable_functions (rs; measured_quantities = [:V , :C ], known_p = [:c1 ], loglevel)
344
+ gi_2 = assess_identifiability (rs; measured_quantities = [:V , :C ], known_p = [:c1 ], loglevel, remove_conserved )
345
+ li_2 = assess_local_identifiability (rs; measured_quantities = [:V , :C ], known_p = [:c1 ], loglevel, remove_conserved )
346
+ ifs_2 = find_identifiable_functions (rs; measured_quantities = [:V , :C ], known_p = [:c1 ], loglevel, remove_conserved )
345
347
@test sym_dict (gi_2) == Dict ([:X => :nonidentifiable , :C => :globally , :V => :globally , :k => :nonidentifiable ,
346
348
:c1 => :globally , :c2 => :nonidentifiable , :p => :nonidentifiable , :d => :globally ])
347
349
@test sym_dict (li_2) == Dict ([:X => 0 , :C => 1 , :V => 1 , :k => 0 , :c1 => 1 , :c2 => 0 , :p => 0 , :d => 1 ])
0 commit comments