Update

Author: odow

test/FileFormats/MPS/MPS.jl

@@ -1277,6 +1277,109 @@ function test_issue_2538()
     return
 end
 
+function test_qmatrix_objective()
+    file = """
+    NAME
+    ROWS
+        N  OBJ
+    COLUMNS
+        x         OBJ        1
+        y         OBJ        1
+    RHS
+    RANGES
+    BOUNDS
+        FR bounds    x
+        FR bounds    y
+    QMATRIX
+        x         x         10
+        x         y         2.0
+        y         x         2.0
+        y         y         2.0
+    ENDATA
+    """
+    io = IOBuffer()
+    print(io, file)
+    seekstart(io)
+    model = MPS.Model()
+    read!(io, model)
+    x, y = MOI.get.(model, MOI.VariableIndex, ["x", "y"])
+    F = MOI.ScalarQuadraticFunction{Float64}
+    @test isapprox(
+        MOI.get(model, MOI.ObjectiveFunction{F}()),
+        1.0 * x + 1.0 * y + 5.0 * x * x + 2.0 * x * y + 1.0 * y * y,
+    )
+    return
+end
+
+function test_qsection_objective()
+    file = """
+    NAME
+    ROWS
+        N  OBJ
+    COLUMNS
+        x         OBJ        1
+        y         OBJ        1
+    RHS
+    RANGES
+    BOUNDS
+        FR bounds    x
+        FR bounds    y
+    QSECTION      OBJ
+        x         x         10
+        x         y         2.0
+        y         y         2.0
+    ENDATA
+    """
+    io = IOBuffer()
+    print(io, file)
+    seekstart(io)
+    model = MPS.Model()
+    read!(io, model)
+    x, y = MOI.get.(model, MOI.VariableIndex, ["x", "y"])
+    F = MOI.ScalarQuadraticFunction{Float64}
+    @test isapprox(
+        MOI.get(model, MOI.ObjectiveFunction{F}()),
+        1.0 * x + 1.0 * y + 5.0 * x * x + 2.0 * x * y + 1.0 * y * y,
+    )
+    return
+end
+
+function test_qsection_row()
+    file = """
+    NAME
+    ROWS
+        N  OBJ
+        L  c1
+    COLUMNS
+        x         c1        1
+        y         c1        1
+    RHS
+        rhs       c1        1
+    RANGES
+    BOUNDS
+        FR bounds    x
+        FR bounds    y
+    QSECTION      c1
+        x         x         10
+        x         y         2.0
+        y         y         2.0
+    ENDATA
+    """
+    io = IOBuffer()
+    print(io, file)
+    seekstart(io)
+    model = MPS.Model()
+    read!(io, model)
+    x, y = MOI.get.(model, MOI.VariableIndex, ["x", "y"])
+    c1 = MOI.get(model, MOI.ConstraintIndex, "c1")
+    @test isapprox(
+        MOI.get(model, MOI.ConstraintFunction(), c1),
+        1.0 * x + 1.0 * y + 10.0 * x * x + 4.0 * x * y + 2.0 * y * y,
+    )
+    return
+end
+
+
 function test_qcmatrix_read_gurobi()
     file = """
     NAME