Skip to content

Commit

Permalink
chore: removed debugging statements from _construct_bound_mapping
Browse files Browse the repository at this point in the history
  • Loading branch information
phajy committed May 4, 2024
1 parent d17d644 commit b92a056
Showing 1 changed file with 0 additions and 12 deletions.
12 changes: 0 additions & 12 deletions src/fitting/binding.jl
Original file line number Diff line number Diff line change
@@ -1,52 +1,40 @@
_sort_binding!(binding) = sort!(binding, by = i -> i[1])

function _construct_bound_mapping(bindings, parameter_count)
println("Bindings are ", bindings)
println("Parameter count is ", parameter_count, " with length ", length(parameter_count))
remove = Int[]

parameter_mapping = map((1:length(parameter_count)...,)) do i
# indices of the full parameter array for the ith model
collect(range(_get_range(parameter_count, i)...))
end
println("Initial parameter mapping is ", parameter_mapping)

for binding in bindings
println(" Looking at binding ", binding)
# the model we use as a reference to bind *to*
reference = binding[1]
println(" Reference to bind to is ", reference)
for b in @views binding[2:end]
println(" Now considering binding ", b, " which has elements ", b[1], " and ", b[2])
# get the number of the parameter that we are binding
parameter_number = parameter_mapping[b[1]][b[2]]
parameter_mapping[b[1]][b[2]] = reference[2]
println(" Reference[2] is ", reference[2])
println(" Parameter mapping is ", parameter_mapping)

# mark for removal: find the parameter index in the global array
N = sum(length(parameter_mapping[q]) for q = 1:b[1]-1)
index = N + b[2]
push!(remove, index)
println(" Remove updated to be ", remove, " having added ", index, " with value ", parameter_number)

# need to now shuffle all the indices greater than this one down by 1
for k = b[2]+1:length(parameter_mapping[b[1]])
if (parameter_mapping[b[1]][k] > parameter_number)
parameter_mapping[b[1]][k] -= 1
println(" Changed parameter_mapping[", b[1], "][", k, "] to ", parameter_mapping[b[1]][k])
end
end
# and for subsequent models
for j = b[1]+1:length(parameter_count)
for k = 1:length(parameter_mapping[j])
if parameter_mapping[j][k] > parameter_number
parameter_mapping[j][k] -= 1
println(" Also changed mapping[", j, "][", k, "] to ", parameter_mapping[j][k])
end
end
end
println(" Parameter mapping is ", parameter_mapping)
end
end

Expand Down

0 comments on commit b92a056

Please sign in to comment.