def set_diff1rate(self):
mdiff = -self.ups0 * self.h * (1 - 2 * self.chi *
(1 - self.C * self.C))
dymdiff = cuda.dy4_m22(mdiff, self.dx())
dyC1 = cuda.dy4_m22((1 - self.C), self.dx())
dyyC1 = cuda.dyy4_m22((1 - self.C), self.dx2())
self.diff1rate = -(dymdiff * dyC1 + mdiff * dyyC1)
def set_diff2rate(self):
mdiff = -self.ups0 * self.h * (1 - 2 * self.chi *
(1 - self.C * self.C))
dymdiff = cuda.dy4_m22(mdiff, self.dx())
dyC2 = cuda.dy4_m22(self.C, self.dx())
dyyC2 = cuda.dyy4_m22(self.C, self.dx2())
self.diff2rate = -(dymdiff * dyC2 + mdiff * dyyC2)
def set_conv2rate(self):
if (not hasattr(self, "dyh")):
self.set_dyh()
if (not hasattr(self, "dyyyh")):
self.set_dyyyh()
if (not hasattr(self, "dyyyyh")):
self.set_dyyyyh()
if (not hasattr(self, "dfdh")):
self.set_dfdh()
dydfdh = cuda.dy4_m22(self.dfdh, self.dx())
dyydfdh = cuda.dyy4_m22(self.dfdh, self.dx2())
self.set_M2()
self.set_dyM2()
self.set_Gy()
self.set_dyGy()
self.conv2rate = -(self.dyM2 * (self.dyyyh - dydfdh) + self.M2 *
(self.dyyyyh - dyydfdh) - self.dyM2 * self.Gy -
self.M2 * self.dyGy + self.v * self.adv2)
def set_conv1rateComoving(self):
if (not hasattr(self, "dyh")):
self.set_dyh()
if (not hasattr(self, "dyyyh")):
self.set_dyyyh()
if (not hasattr(self, "dyyyyh")):
self.set_dyyyyh()
if (not hasattr(self, "dfdh")):
self.set_dfdh()
dydfdh = cuda.dy4_m22(self.dfdh, self.dx())
dyydfdh = cuda.dyy4_m22(self.dfdh, self.dx2())
self.set_M1()
self.set_dyM1()
self.set_Gy()
self.set_dyGy()
self.conv1rateComoving = -(self.dyM1 *
(self.dyyyh - dydfdh) + self.M1 *
(self.dyyyyh - dyydfdh) -
self.dyM1 * self.Gy - self.M1 * self.dyGy)
def set_diff2(self):
dyC2 = cuda.dy4_m22(self.C, self.dx())
self.diff2 = -self.ups0 * self.h * (1 - 2 * self.chi *
(1 - self.C) * self.C) * dyC2
def set_diff1(self):
dyC1 = cuda.dy4_m22((1 - self.C), self.dx())
self.diff1 = -self.ups0 * self.h * (1 - 2 * self.chi *
(1 - self.C) * self.C) * dyC1
def set_dyphi(self):
self.dyphi = cuda.dy4_m22(self.phi, self.dx())
def set_dyC2(self):
self.dyC2 = cuda.dy4_m22(self.C, self.dx())
def set_dyC1(self):
self.dyC1 = cuda.dy4_m22((1 - self.C), self.dx())