diff --git a/scri/asymptotic_bondi_data/from_initial_values.py b/scri/asymptotic_bondi_data/from_initial_values.py index 1d1041b4..2760c107 100644 --- a/scri/asymptotic_bondi_data/from_initial_values.py +++ b/scri/asymptotic_bondi_data/from_initial_values.py @@ -100,7 +100,7 @@ def asany_atleast2d_complex(a): ψ3_0 = -ð(conjugate(σ_1)) ψ3_1 = -2 * ð(conjugate(σ_2)) abd.psi3 = u * ψ3_1 + ψ3_0 - # ψ₂ = ∫ (-ðψ₃ + σψ₄) du + # ψ₂ = ∫ (ðψ₃ + σψ₄) du ψ2_0 = ( ModesTimeSeries(psi2, u, spin_weight=0, multiplication_truncator=max).real - 1j * (σ_0.bar.eth_GHP.eth_GHP + σ_0 * σ_1.bar).imag diff --git a/tests/test_abd_ivp.py b/tests/test_abd_ivp.py index 44afefa5..bc423a05 100644 --- a/tests/test_abd_ivp.py +++ b/tests/test_abd_ivp.py @@ -55,7 +55,7 @@ def validator(abd): def test1(): """Add ℓ=0 term to ψ₂ intial value - Ensures that first terms of ψ̇₁u = -ðψ₂ + 2σψ₃ and ψ̇₀ = -ðψ₁ + 3σψ₂ don't get excited + Ensures that first terms of ψ̇₁u = ðψ₂ + 2σψ₃ and ψ̇₀ = ðψ₁ + 3σψ₂ don't get excited """ @@ -85,7 +85,7 @@ def validator(abd): def test2(): """Add ℓ=1 term to ψ₂ intial value - Checks first term of ψ̇₁ = -ðψ₂ + 2σψ₃ + Checks first term of ψ̇₁ = ðψ₂ + 2σψ₃ """ @@ -115,7 +115,7 @@ def validator(abd): def test3(): """Add ℓ=2 term to ψ₂ intial value - Checks first term of ψ̇₀ = -ðψ₁ + 3σψ₂ + Checks first term of ψ̇₀ = ðψ₁ + 3σψ₂ """ @@ -147,7 +147,7 @@ def validator(abd): def test4(): """Add nonzero constant term to σ - Checks first term of ψ̇₁ = -ðψ₂ + 2σψ₃ and second term of ψ̇₀ = -ðψ₁ + 3σψ₂ + Checks first term of ψ̇₁ = ðψ₂ + 2σψ₃ and second term of ψ̇₀ = ðψ₁ + 3σψ₂ After satisfaction of the reality condition on the mass aspect, ψ₂ has nonzero modes in {(0, 0), (2, -2), (2, 2)}, so ψ₁ should have nonzero modes in {(2, -2), (2, 2)}. Here, @@ -182,7 +182,7 @@ def validator(abd): def test5(): """Add nonzero ℓ=2 term to σ̇ - Checks second term of ψ̇₁ = -ðψ₂ + 2σψ₃ + Checks second term of ψ̇₁ = ðψ₂ + 2σψ₃ """ @@ -212,7 +212,7 @@ def validator(abd): def test6(): """Add nonzero ℓ=2 term to σ̈ - Checks second term of ψ̇₂ = -ðψ₃ + σψ₄ + Checks second term of ψ̇₂ = ðψ₃ + σψ₄ """