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 ] )