|
118 | 118 |
|
119 | 119 | # creates a *sparse* reaction complex matrix
|
120 | 120 | function reactioncomplexes(::Type{SparseMatrixCSC{Int, Int}}, rn::ReactionSystem,
|
121 |
| - complextorxsmap) |
| 121 | + complextorxsmap) |
122 | 122 | # computes the I, J, and V vectors used for the sparse matrix (read about sparse matrix
|
123 | 123 | # representation for more information)
|
124 | 124 | complexes = collect(keys(complextorxsmap))
|
|
480 | 480 | function subnetworkmapping(linkageclass, allrxs, complextorxsmap, p)
|
481 | 481 | # Finds the reactions that are part of teh sub-reaction network.
|
482 | 482 | rxinds = sort!(collect(Set(rxidx for rcidx in linkageclass
|
483 |
| - for rxidx in complextorxsmap[rcidx]))) |
| 483 | + for rxidx in complextorxsmap[rcidx]))) |
484 | 484 | rxs = allrxs[rxinds]
|
485 | 485 | specset = Set(s for rx in rxs for s in rx.substrates if !isconstant(s))
|
486 | 486 | for rx in rxs
|
@@ -531,7 +531,7 @@ function subnetworks(rs::ReactionSystem)
|
531 | 531 | newrxs, newspecs, newps = subnetworkmapping(lcs[i], rxs, complextorxsmap, p)
|
532 | 532 | newname = Symbol(nameof(rs), "_", i)
|
533 | 533 | push!(subnetworks,
|
534 |
| - ReactionSystem(reacs, t, specs, newps; name = newname, spatial_ivs)) |
| 534 | + ReactionSystem(reacs, t, specs, newps; name = newname, spatial_ivs)) |
535 | 535 | end
|
536 | 536 | subnetworks
|
537 | 537 | end
|
@@ -745,8 +745,8 @@ function cache_conservationlaw_eqs!(rn::ReactionSystem, N::AbstractMatrix, col_o
|
745 | 745 | # finds the coefficient (in the conservation law) of the species that is eliminated
|
746 | 746 | # by this conservation law
|
747 | 747 | scaleby = (N[i, depidx] != 1) ? N[i, depidx] : one(eltype(N))
|
748 |
| - (scaleby != 0) || error("Error, found a zero in the conservation law matrix where " |
749 |
| - * "one was not expected.") |
| 748 | + (scaleby != 0) || |
| 749 | + error("Error, found a zero in the conservation law matrix where one was not expected.") |
750 | 750 |
|
751 | 751 | # creates, for this conservation law, the sum of all independent species (weighted by
|
752 | 752 | # the ratio between the coefficient of the species and the species which is elimianted
|
|
0 commit comments