Skip to content

Commit

Permalink
fixing some issues on Julia 1.7-rc1
Browse files Browse the repository at this point in the history
  • Loading branch information
mateuszbaran committed Sep 21, 2021
1 parent cf00e22 commit c723345
Show file tree
Hide file tree
Showing 3 changed files with 3 additions and 1 deletion.
1 change: 1 addition & 0 deletions src/nlsolve.jl
Original file line number Diff line number Diff line change
Expand Up @@ -55,6 +55,7 @@ function _inverse_retract_nlsolve(
retract!(M, F, p, project(M, p, X), retraction; kwargs...)
project_point && project!(M, q, q)
F .-= q
project_point && project!(M, F, F)

This comment has been minimized.

Copy link
@kellertuer

kellertuer Sep 21, 2021

Member

Why is this necessary? F has the same values as q and q was just projected in the line before?

This comment has been minimized.

Copy link
@mateuszbaran

mateuszbaran Sep 21, 2021

Author Member

That's because F may be outside the manifold after subtracting q. The respective test was failing because the returned point was off manifold by about 1e-8 IIRC.

This comment has been minimized.

Copy link
@kellertuer

kellertuer Sep 21, 2021

Member

Ah I did somehow miss the minus when I wrote the question.

This comment has been minimized.

Copy link
@mateuszbaran

mateuszbaran Sep 21, 2021

Author Member

No problem :). Probably the thing we actually don't need is projecting q, so it looks like it might be good to pay a bit more attention to this code.

return F
end
isdefined(Manifolds, :NLsolve) ||
Expand Down
1 change: 1 addition & 0 deletions test/groups/special_euclidean.jl
Original file line number Diff line number Diff line change
Expand Up @@ -78,6 +78,7 @@ Random.seed!(10)
X_pts;
test_diff=true,
diff_convs=[(), (LeftAction(),), (RightAction(),)],
atol=1e-9,
)
end
end
Expand Down
2 changes: 1 addition & 1 deletion test/manifolds/rotations.jl
Original file line number Diff line number Diff line change
Expand Up @@ -116,7 +116,7 @@ include("../utils.jl")
point_distributions=[ptd],
tvector_distributions=[tvd],
basis_types_to_from=basis_types,
exp_log_atol_multiplier=20,
exp_log_atol_multiplier=80,
retraction_atol_multiplier=12,
test_inplace=true,
)
Expand Down

0 comments on commit c723345

Please sign in to comment.