Skip to content

Commit d9d06fb

Browse files
committed
remove curvature computation with BK
1 parent 51e0992 commit d9d06fb

File tree

2 files changed

+3
-17
lines changed

2 files changed

+3
-17
lines changed

ext/BifurcationKitExt/problems.jl

Lines changed: 0 additions & 17 deletions
Original file line numberDiff line numberDiff line change
@@ -40,23 +40,6 @@ function correct_guess(cache, options::BK.NewtonPar)
4040
return sol.u
4141
end
4242

43-
function get_curvature(prob, c::Chart, par)
44-
u0 = c.u
45-
Φ = c.Φ
46-
n, m = size(prob)
47-
d = n-m
48-
J = jacobian(prob, u0, par)
49-
_A = vcat(J, Φ')
50-
A = zeros(d, n, d)
51-
for i = 1:d
52-
Φi = Φ[:, i]
53-
for j = 1:d
54-
Φj = Φ[:, j]
55-
d2 = d2F(prob, u0, par, Φi, Φj)
56-
A[i, :, j] .= _A \ vcat(d2, zeros(d))
57-
end
58-
end
59-
norm(A)
6043
end
6144

6245
##############################################################################################################

test/plane.jl

Lines changed: 3 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -52,3 +52,6 @@ S = continuation(prob,
5252

5353
@test S[1] isa MPC.Chart
5454
@test length(S) == 109
55+
56+
show(S)
57+
show(S[1])

0 commit comments

Comments
 (0)