Equating two symbolically symmetric matrices in JuMP without redundant constraints

You can do

@constraint(test_model, vectorize(X - Y, SymmetricMatrixShape(n)) .== 0)
1 Like