def dy_dt(t, QC_ins, QC, dQC_ins):
    Cs = QC[:(- 1)]
    dydt_cstr_no_rxn_fixed_aer(QC_ins, dQC_ins, V_arr, Q_e_arr, _dstate, Cs)
    _dstate[:(- 1)] += r(*Cs)
    _update_dstate()
def dy_dt(t, QC_ins, QC, dQC_ins):
    Cs = QC[:(- 1)]
    dydt_cstr_no_rxn_fixed_aer(QC_ins, dQC_ins, V_arr, Q_e_arr, _dstate, Cs)
    _dstate[:(- 1)] += r(*Cs)
    _update_dstate()