fix zeroing constant comps

This commit is contained in:
Gaspard Jankowiak 2026-01-12 09:04:17 +01:00
commit ca567122ff

View file

@ -28,9 +28,8 @@ true
function check(f, Jf, x, constant_components::Vector{Int}=Int[]; taylortestplot::Bool=false, taylortestdirection=nothing, f_kwargs...)
if isnothing(taylortestdirection)
direction = 2 * (rand(size(x)...) .- 0.5)
for cst in constant_components
direction[cst] = 0
end
direction[constant_components] .= 0
if direction isa Array
direction ./= norm(direction)