Fix bugs in BVP test problems

Signed-off-by: ErikQQY <2283984853@qq.com>

Author: ErikQQY

lib/BVProblemLibrary/Project.toml

@@ -1,6 +1,6 @@
 name = "BVProblemLibrary"
 uuid = "ded0fc24-dfea-4565-b1d9-79c027d14d84"
-version = "0.1.3"
+version = "0.1.4"
 
 [deps]
 DiffEqBase = "2b5f629d-d688-5b77-993f-72d75c75574e"

lib/BVProblemLibrary/src/BVProblemLibrary.jl

@@ -24,6 +24,8 @@ function flat_moon_bca!(res_a, u_a, p)
     res_a[4] = u_a[4]
 end
 function flat_moon_bcb!(res_b, u_b, p)
+    h = 185.2
+    Vc = 1627
     res_b[1] = u_b[5] - h
     res_b[2] = u_b[6] - Vc
     res_b[3] = u_b[7]
@@ -98,7 +100,7 @@ function flat_earth_f!(du, u, p, t)
     Vc = sqrt(398600.4 / (6378.14 + 300)) * 1000
     h = 300000
     g = 9.80665
-    thrust2Weight = 3
+    Thrust2Weight = 3
     acc = Thrust2Weight * g
 
     du[1] = (u[3] * (Vc / h))

lib/BVProblemLibrary/src/linear.jl

@@ -454,10 +454,10 @@ function prob_bvp_linear_8_f!(du, u, p, t)
     du[1] = u[2]
     du[2] = -1 / p * u[2]
 end
-function prob_bvp_linear_8_bca!(res_a!, u_a, p)
+function prob_bvp_linear_8_bca!(res_a, u_a, p)
     res_a[1] = u_a[1] - 1
 end
-function prob_bvp_linear_8_bcb!(res_b!, u_b, p)
+function prob_bvp_linear_8_bcb!(res_b, u_b, p)
     res_b[1] = u_b[1] - 2
 end
 prob_bvp_linear_8_function = BVPFunction(prob_bvp_linear_8_f!, (prob_bvp_linear_8_bca!, prob_bvp_linear_8_bcb!),

lib/BVProblemLibrary/src/nonlinear.jl

@@ -407,10 +407,10 @@ function prob_bvp_nonlinear_8_f!(du, u, p, t)
     du[1] = u[2]
     du[2] = 1 / p * prob_bvp_nonlinear_8_f(u[1], u[2])
 end
-function prob_bvp_nonlinear_8_bca!(resid_a, u_a, p)
+function prob_bvp_nonlinear_8_bca!(res_a, u_a, p)
     res_a[1] = u_a[1] - 1
 end
-function prob_bvp_nonlinear_8_bcb!(resid_b, u_b, p)
+function prob_bvp_nonlinear_8_bcb!(res_b, u_b, p)
     res_b[1] = u_b[1] + 1 / 3
 end
 prob_bvp_nonlinear_8_function = BVPFunction(prob_bvp_nonlinear_8_f!, (prob_bvp_nonlinear_8_bca!, prob_bvp_nonlinear_8_bcb!), bcresid_prototype = (zeros(1), zeros(1)), twopoint = Val(true))
@@ -615,10 +615,10 @@ function prob_bvp_nonlinear_12_f!(du, u, p, t)
     du[1] = u[2]
     du[2] = 1 / p * prob_bvp_nonlinear_12_f(u[1], u[2])
 end
-function prob_bvp_nonlinear_12_bca!(resid_a, u_a, p)
+function prob_bvp_nonlinear_12_bca!(res_a, u_a, p)
     res_a[1] = u_a[1] + 7 / 6
 end
-function prob_bvp_nonlinear_12_bcb!(resid_b, u_b, p)
+function prob_bvp_nonlinear_12_bcb!(res_b, u_b, p)
     res_b[1] = u_b[1] - 3 / 2
 end
 prob_bvp_nonlinear_12_function = BVPFunction(prob_bvp_nonlinear_12_f!, (prob_bvp_nonlinear_12_bca!, prob_bvp_nonlinear_12_bcb!), bcresid_prototype = (zeros(1), zeros(1)), twopoint = Val(true))
@@ -733,7 +733,7 @@ function prob_bvp_nonlinear_14_f!(du, u, p, t)
     du[3] = u[4]
     du[4] = prob_bvp_nonlinear_14_f(u[1], u[2], u[3], u[4], p)
 end
-function prob_bvp_nonlinear_14_bca!(res_a, u_b, p)
+function prob_bvp_nonlinear_14_bca!(res_a, u_a, p)
     res_a[1] = u_a[1]
     res_a[2] = u_a[2]
 end