Example
using MadNLP, JuMP, IPMDiagnosis
model = Model(MadNLP.Optimizer)
set_attribute(model, "intermediate_callback", DiagnosisCallback())
I = 1:2
@variable(model, x[I])
@constraint(model, eq, x[1]^2 + x[2]^2 == 1)
@objective(model, Min, sum(x[i] for i in I))
optimize!(model)