def erf(self, o, fp):
f, = o.ufl_operands
return fp * (2.0 / sqrt(pi) * exp(-f**2))
def erf(self, o, fp):
f, = o.ufl_operands
return fp * (2.0 / sqrt(pi) * exp(-f**2))
def erf(self, o, a):
f, fp = a
o = self.reuse_if_possible(o, f)
op = fp*(2.0/sqrt(pi)*exp(-f**2))
return (o, op)
def erf(self, f):
"d/dx erf x = 2/sqrt(pi)*exp(-x^2)"
x, = f.operands()
return (2.0/sqrt(pi)*exp(-x**2),)
def Kinetic_oxy(Temperature):
"""Kinetic constant - mathematical expression"""
k_oxy = exp(alfa) * exp(-beta / Temperature)
return k_oxy
def erf(self, o, a):
f, fp = a
o = self.reuse_if_possible(o, f)
op = fp * (2.0 / sqrt(pi) * exp(-f**2))
return (o, op)
def Kinetic_oxy(Temperature):
"""mathematical expression of
kinetic temperature-dependent constant"""
k_oxy = exp(alfa) * exp(-beta / Temperature)
return k_oxy