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)
Example #2
0
 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)
Example #3
0
 def __init__(self, shape_or_k=0):
     AnisotropicUniverse.__init__(self)
     CurvedUniverse.__init__(self, shape_or_k)
Example #4
0
 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)
Example #5
0
 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)
Example #6
0
	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)
Example #7
0
	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)