diff --git a/src/periodicorbit/PeriodicOrbitCollocation.jl b/src/periodicorbit/PeriodicOrbitCollocation.jl
index 0f53966f..8640052e 100644
--- a/src/periodicorbit/PeriodicOrbitCollocation.jl
+++ b/src/periodicorbit/PeriodicOrbitCollocation.jl
@@ -403,7 +403,7 @@ $(SIGNATURES)
     n, m, Ntst = size(pb)
     L, ∂L = get_Ls(pb.mesh_cache)
     ω = pb.mesh_cache.gauss_weight
-    mesh = pb.mesh_cache.mesh
+    mesh = pb.mesh_cache.τs
 
     guj = zeros(Ty, n, m)
     uj  = zeros(Ty, n, m+1)
@@ -1203,11 +1203,11 @@ end
     rg = 1:nbcoll
     for _ in 1:Ntst
         F = lu(J[rg, rg .+ N])
-        P[rg, rg] .= (F.P \ F.L)
+        𝐅𝐬[rg, rg] .= (F.P \ F.L)
         rg = rg .+ nbcoll
     end
 
-    Fₚ = lu(P)
+    Fₚ = lu(𝐅𝐬)
     Jcond = Fₚ \ J
     rhs = Fₚ \ rhs0
 
diff --git a/test/stuartLandauCollocation.jl b/test/stuartLandauCollocation.jl
index feaf6235..a64c417b 100644
--- a/test/stuartLandauCollocation.jl
+++ b/test/stuartLandauCollocation.jl
@@ -42,7 +42,7 @@ length(prob_col)
 BK.get_times(prob_col)
 BK.get_max_time_step(prob_col)
 size(prob_col.mesh_cache)
-BK.update_mesh!(prob_col, prob_col.mesh_cache.mesh)
+BK.update_mesh!(prob_col, prob_col.mesh_cache.τs)
 PeriodicOrbitOCollProblem(10, 2) |> BK.get_mesh_size
 BK.get_Ls(prob_col)
 show(prob_col)