Skip to content

Commit

Permalink
correction of linear solver in NS MinAug
Browse files Browse the repository at this point in the history
remove updatesection! in init of Fold PO Coll
  • Loading branch information
rveltz committed Aug 24, 2024
1 parent e6482d6 commit 8705731
Show file tree
Hide file tree
Showing 2 changed files with 2 additions and 2 deletions.
2 changes: 1 addition & 1 deletion src/periodicorbit/codim2/MinAugNS.jl
Original file line number Diff line number Diff line change
Expand Up @@ -126,7 +126,7 @@ function NSMALinearSolver(x, p::𝒯, Ο‰::𝒯, 𝐍𝐒::NeimarkSackerProblemMi
~cv && @debug "[codim2 NS] Linear solver for N did not converge."

# we solve Nα΅—[w, Οƒ2] = [0, 1]
w, Οƒ2, cv, itw = nstest(JNSβ˜…, b, a, zero(𝒯), 𝐍𝐒.zero, one(𝒯); lsbd = 𝐍𝐒.linbdsolver)
w, Οƒ2, cv, itw = nstest(JNSβ˜…, b, a, zero(𝒯), 𝐍𝐒.zero, one(𝒯); lsbd = 𝐍𝐒.linbdsolverAdjoint)
~cv && @debug "[codim2 NS] Linear solver for Nα΅— did not converge."

Ξ΄ = getdelta(POWrap)
Expand Down
2 changes: 1 addition & 1 deletion src/periodicorbit/codim2/PeriodicOrbitCollocation.jl
Original file line number Diff line number Diff line change
Expand Up @@ -128,7 +128,7 @@ function continuation_coll_fold(br::AbstractResult{Tkind, Tprob},

# update section
# THIS IS A HACK, SHOULD BE SAVED FOR PROPER BRANCHING ETC
updatesection!(coll, Ο•, nothing)
# updatesection!(coll, Ο•, nothing)

_finsol = modify_po_finalise(FoldMAProblem(FoldProblemMinimallyAugmented(WrapPOColl(coll)), lens2), kwargs, coll.update_section_every_step)

Expand Down

0 comments on commit 8705731

Please sign in to comment.