-
Notifications
You must be signed in to change notification settings - Fork 0
/
explore.py
executable file
·30 lines (25 loc) · 1017 Bytes
/
explore.py
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
#!/usr/bin/env python3
# -*- coding: utf-8 -*-
"""
explore
Used to explore initial data
"""
from run import params, initial_data
from evolver.initialize import unpack
from math import pi
from evolver.eoms import compute_hartree, compute_rho, compute_deltarho2
unpacked_data = unpack(initial_data, params.total_wavenumbers)
a, adot, phi0, phi0dot, phiA, phidotA, psiA, phiB, phidotB, psiB = unpacked_data
H = adot/a
phi2pt, phi2ptdt, phi2ptgrad = compute_hartree(phiA, phidotA, phiB, phidotB, params)
rho = compute_rho(phi0, phi0dot, params.model)
deltarho2 = compute_deltarho2(a, phi0, phi2pt, phi2ptdt, phi2ptgrad, params.model)
print("rho:", rho)
print("deltarho2:", deltarho2)
print("deltarho2/rho:", deltarho2/rho)
print("kappa^4/(2 pi^2):", params.kappa**4/(2*pi**2))
print("deltarho2 / (kappa^4/(2 pi^2)):", deltarho2/(params.kappa**4/(2*pi**2)))
print()
print("phi2pt:", phi2pt)
print("kappa^2/(4 pi^2):", params.kappa**2/4/pi**2)
print("phi2pt / (kappa^2/(4 pi^2)):", phi2pt / (params.kappa**2/4/pi**2))