Skip to content

Commit

Permalink
Tweak variable names in new_branches!
Browse files Browse the repository at this point in the history
  • Loading branch information
tkf committed Jun 18, 2018
1 parent 1e90247 commit 230ba04
Showing 1 changed file with 6 additions and 6 deletions.
12 changes: 6 additions & 6 deletions src/continuations/branching.jl
Original file line number Diff line number Diff line change
Expand Up @@ -96,21 +96,21 @@ end


function new_branches!(cache, opts, sbint::SimpleBifurcationInterval)
u0, tJ, L, Q = find_simple_bifurcation!(cache, opts, sbint)
ub0, tJ, L, Q = find_simple_bifurcation!(cache, opts, sbint)
# TODO: handle not-found case

tv1, tv2 = solve_simple_bifurcation!(cache, opts, u0, tJ, L, Q)
tv1, tv2 = solve_simple_bifurcation!(cache, opts, ub0, tJ, L, Q)
tv0 = sbint.u1 .- sbint.u0
tv0 /= norm(tv0)
@assert abs(tv1 tv0) > abs(tv2 tv0) # TODO: just warn
tv = tv2

args = (sbint.h, sbint.direction)
u1, h1 = predictor_corrector_step!(cache, opts, u0, tv, args...)
u2, h2 = predictor_corrector_step!(cache, opts, u0, -tv, args...)
ub1, h1 = predictor_corrector_step!(cache, opts, ub0, tv, args...)
ub2, h2 = predictor_corrector_step!(cache, opts, ub0, -tv, args...)

return [
(u0, u1, sbint.direction, h1),
(u0, u2, sbint.direction, h2),
(ub0, ub1, sbint.direction, h1),
(ub0, ub2, sbint.direction, h2),
]
end

0 comments on commit 230ba04

Please sign in to comment.