from bosonstar.ComplexBosonStar import Complex_Boson_Star


# =====================
#  All imporntnat definitions
# =====================

# Physics defintions
phi0 = 0.40         # centeral phi
D = 5.0             # Dimension (total not only spacial)
Lambda = -0.2       # Cosmological constant
# Solver definitions
Rstart = 3
Rend = 50.00
deltaR = 1
N = 100000
e_pow_minus_delta_guess = 0.4999
verbose = 2
eps = 1e-10  # Small epsilon to avoid r \neq 0

# ====================================
#   Main routine
# ====================================

pewpew = Complex_Boson_Star(e_pow_minus_delta_guess, phi0, D, Lambda, verbose)

pewpew.print_parameters()

alpha0 = pewpew.radial_walker(Rstart, Rend, deltaR, N, eps)

# =====================================
#   Output and plotting
# =====================================
soldict = pewpew.get_solution()

# Makes sure that lapse goes asymptotically to 1
# (Not an essential step, but recommended)
pewpew.normalise_edelta()

pewpew.check_Einstein_equation()

# ===============================
path = pewpew.get_path()
pewpew.plot_solution()
pewpew.print_solution()

[((25, 9, 25, 78), 'bosonstar.ComplexBosonStar.Complex_Boson_Star', 'Complex_Boson_Star', ({(25, 28, 25, 51): 'e_pow_minus_delta_guess', (25, 53, 25, 57): 'phi0', (25, 59, 25, 60): 'D', (25, 62, 25, 68): 'Lambda', (25, 70, 25, 77): 'verbose'}, {}), '(e_pow_minus_delta_guess, phi0, D, Lambda, verbose)', False, 'from bosonstar.ComplexBosonStar import Complex_Boson_Star\n')]