From bcb39756d18e9574d1c14f09bbe646d2b6b69729 Mon Sep 17 00:00:00 2001
From: Omar Ashour <omar.ashour22@gmail.com>
Date: Wed, 24 Feb 2021 21:41:59 -0800
Subject: [PATCH] Fix CMP 4 Coefficients

---
 src/CubicSolvers.jl | 8 ++++----
 1 file changed, 4 insertions(+), 4 deletions(-)

diff --git a/src/CubicSolvers.jl b/src/CubicSolvers.jl
index 94071a6..cfcf661 100644
--- a/src/CubicSolvers.jl
+++ b/src/CubicSolvers.jl
@@ -327,10 +327,10 @@ See also: [`solve!`](@ref), [`Operators`](@ref)
 function T4A_CMP!(ψₒ, ψᵢ, dx, ops)
     T2A!(ψₒ, ψᵢ, dx/2, ops)
     T2A!(ψₒ, ψₒ, dx/2, ops)
-    ψₒ .*= -1/3
+    ψₒ .*= 4/3
 
     T2A!(ops.ψ₁, ψᵢ, dx,ops)
-    ops.ψ₁ .*= 4/3
+    ops.ψ₁ .*= -1/3
 
     ψₒ .+= ops.ψ₁
 end
@@ -345,10 +345,10 @@ See also: [`solve!`](@ref), [`Operators`](@ref)
 function T4B_CMP!(ψₒ, ψᵢ, dx, ops)
     T2B!(ψₒ, ψᵢ, dx/2, ops)
     T2B!(ψₒ, ψₒ, dx/2, ops)
-    ψₒ .*= -1/3
+    ψₒ .*= 4/3
 
     T2B!(ops.ψ₁, ψᵢ, dx,ops)
-    ops.ψ₁ .*= 4/3
+    ops.ψ₁ .*= -1/3
 
     ψₒ .+= ops.ψ₁
 end