Source code for mpet.electrode.materials.testRS_ps

import numpy as np



[docs]
def testRS_ps(self, y, ybar, muR_ref, ISfuncs=None):
    muRtheta = -self.eokT*2.
    muRhomog = self.reg_sln(y, self.get_trode_param("Omega_a"), ISfuncs)
    muRnonHomog = self.general_non_homog(y, ybar)
    muR = muRhomog + muRnonHomog
    actR = np.exp(muR/self.T)
    muR += muRtheta + muR_ref
    return muR, actR