def evolve(self, initial_conditions, times): a0, a_dot0, b0, b_dot0, c0, c_dot0 = initial_conditions V0 = a0*b0*c0 Ha0, Hb0, Hc0 = a_dot0/a0, b_dot0/b0, c_dot0/c0 Htot0 = Ha0 + Hb0 + Hc0 I0 = Hb0*Hc0 + Ha0*Hc0 + Ha0*Hb0 self.const = I0/3.0 - (2*self.k)/V0**(2.0/3.0) + (self.k/Htot0)*(Ha0/a0**2 + Hb0/b0**2 + Hc0/c0**2) AnisotropicUniverse.evolve(self, initial_conditions, times)
def evolve(self, initial_conditions, times):
a0, a_dot0, b0, b_dot0, c0, c_dot0 = initial_conditions
V0 = a0 * b0 * c0
Ha0, Hb0, Hc0 = a_dot0 / a0, b_dot0 / b0, c_dot0 / c0
Htot0 = Ha0 + Hb0 + Hc0
I0 = Hb0 * Hc0 + Ha0 * Hc0 + Ha0 * Hb0
self.const = I0 / 3.0 - (2 * self.k) / V0**(2.0 / 3.0) + (
self.k / Htot0) * (Ha0 / a0**2 + Hb0 / b0**2 + Hc0 / c0**2)
AnisotropicUniverse.evolve(self, initial_conditions, times)
def __init__(self, shape_or_k=0):
AnisotropicUniverse.__init__(self)
CurvedUniverse.__init__(self, shape_or_k)
def evolve(self, initial_conditions, times):
a0, a_dot0, b0, b_dot0 = initial_conditions[:-2]
Ha0, Hb0 = a_dot0 / a0, b_dot0 / b0
self.const = -(self.k / a0**2.0 + 2.0 * Ha0 * Hb0 + Ha0**2.0)
AnisotropicUniverse.evolve(self, initial_conditions, times)
def __init__(self, shape_or_k=0):
AnisotropicUniverse.__init__(self)
self.shapes = {'o': -1, 'f': 0, 'c': 1}
self.set_shape_or_k(shape_or_k)
def evolve(self, initial_conditions, times): a0, a_dot0, b0, b_dot0 = initial_conditions[:-2] Ha0, Hb0 = a_dot0/a0, b_dot0/b0 self.const = -(self.k/a0**2.0 + 2.0*Ha0*Hb0 + Ha0**2.0) AnisotropicUniverse.evolve(self, initial_conditions, times)
def __init__(self, shape_or_k = 0): AnisotropicUniverse.__init__(self) CurvedUniverse.__init__(self, shape_or_k)
def __init__(self, shape_or_k = 0):
AnisotropicUniverse.__init__(self)
self.shapes = {'o': -1, 'f': 0, 'c': 1}
self.set_shape_or_k(shape_or_k)
