def get_neqsys_chained_conditional(
self, rref_equil=False, rref_preserv=False, NumSys=NumSysLin, **kwargs
):
from pyneqsys import ConditionalNeqSys, ChainedNeqSys
def mk_factory(NS):
def factory(conds):
return self._SymbolicSys_from_NumSys(
NS, conds, rref_equil, rref_preserv, **kwargs
)
return factory
return ChainedNeqSys(
[
ConditionalNeqSys(
[
(self._fw_cond_factory(ri), self._bw_cond_factory(ri, NS.small))
for ri in self.phase_transfer_reaction_idxs()
],
mk_factory(NS),
)
for NS in NumSys
]
)
def get_neqsys_static_conditions(self, init_concs, rref_equil=False,
rref_preserv=False,
NumSys=NumSysLin, precipitates=None):
if precipitates is None:
precipitates = (False,)*len(self.phase_transfer_reaction_idxs())
from pyneqsys import ChainedNeqSys
return ChainedNeqSys([self._SymbolicSys_from_NumSys(
NS, precipitates, rref_equil, rref_preserv) for NS in NumSys])
def factory(conds):
return ChainedNeqSys(
[
self._SymbolicSys_from_NumSys(
NS, conds, rref_equil, rref_preserv, **kwargs
)
for NS in NumSys
]
)
