From 302832a1a29cd9af1ef01d48446acb08ac9588b3 Mon Sep 17 00:00:00 2001
From: ErikQQY <2283984853@qq.com>
Date: Thu, 15 Aug 2024 21:25:55 +0800
Subject: [PATCH] Fix VectorOfArray lastindex

Signed-off-by: ErikQQY <2283984853@qq.com>
---
 src/BoundaryValueDiffEq.jl | 8 ++++----
 1 file changed, 4 insertions(+), 4 deletions(-)

diff --git a/src/BoundaryValueDiffEq.jl b/src/BoundaryValueDiffEq.jl
index ea942f84..9dcb6f22 100644
--- a/src/BoundaryValueDiffEq.jl
+++ b/src/BoundaryValueDiffEq.jl
@@ -57,10 +57,10 @@ end
 
     bc1! = @closure (residual, u, p, t) -> begin
         residual[1] = u[:, 1][1] - 5
-        residual[2] = u[:, lastindex(u)][1]
+        residual[2] = u[:, end][1]
     end
 
-    bc1 = @closure (u, p, t) -> [u[:, 1][1] - 5, u[:, lastindex(u)][1]]
+    bc1 = @closure (u, p, t) -> [u[:, 1][1] - 5, u[:, end][1]]
 
     bc1_a! = @closure (residual, ua, p) -> (residual[1] = ua[1] - 5)
     bc1_b! = @closure (residual, ub, p) -> (residual[1] = ub[1])
@@ -102,14 +102,14 @@ end
 
     bc1_nlls! = @closure (resid, sol, p, t) -> begin
         solₜ₁ = sol[:, 1]
-        solₜ₂ = sol[:, lastindex(sol)]
+        solₜ₂ = sol[:, end]
         resid[1] = solₜ₁[1]
         resid[2] = solₜ₂[1] - 1
         resid[3] = solₜ₂[2] + 1.729109
         return nothing
     end
     bc1_nlls = @closure (sol, p, t) -> [
-        sol[:, 1][1], sol[:, lastindex(sol)][1] - 1, sol[:, lastindex(sol)][2] + 1.729109]
+        sol[:, 1][1], sol[:, end][1] - 1, sol[:, end][2] + 1.729109]
 
     bc1_nlls_a! = @closure (resid, ua, p) -> (resid[1] = ua[1])
     bc1_nlls_b! = @closure (resid, ub, p) -> (resid[1] = ub[1] - 1;