[Nonlinear.SymbolicAD] simplify quadratic functions if possible

[odow]

No description provided.

[odow]

Now this would help a bunch of #2553

julia> op(head, args...) = MOI.ScalarNonlinearFunction(head, Any[args...])
op (generic function with 1 method)

julia> f = op(:/, 
           op(:+, 
               op(:-, 
                   op(:+, 
                       op(:^, MOI.VariableIndex(9), (2)),
                       op(:^, MOI.VariableIndex(6), (2)),
                       op(:^, MOI.VariableIndex(7), (2)),
                   ), 
                   0.0,
               ), 
               op(:-, 
                   op(:+, 
                       op(:^, MOI.VariableIndex(9), (2)),
                       op(:^, MOI.VariableIndex(6), (2)),
                       op(:^, MOI.VariableIndex(7), (2)),
                   ),
                   0.0,
               ),
               op(:+, 
                   op(:-, 
                       op(:+, 
                           MOI.VariableIndex(6),
                       ),
                   ),
               ),
               op(:+, 
                   op(:-, 
                       op(:+, 
                           MOI.VariableIndex(7),
                       ),
                   ),
               ), 
               op(:+, 
                   op(:+, 
                       MOI.VariableIndex(6),
                   ),
               ), 
               op(:+, 
                   op(:+, 
                       MOI.VariableIndex(7),
                   ),
               ),
           ),
           2.0,
       )
/(+(-(+(^(MOI.VariableIndex(9), (2)), ^(MOI.VariableIndex(6), (2)), ^(MOI.VariableIndex(7), (2))), 0.0), -(+(^(MOI.VariableIndex(9), (2)), ^(MOI.VariableIndex(6), (2)), ^(MOI.VariableIndex(7), (2))), 0.0), +(-(+(MOI.VariableIndex(6)))), +(-(+(MOI.VariableIndex(7)))), +(+(MOI.VariableIndex(6))), +(+(MOI.VariableIndex(7)))), 2.0)

julia> MOI.Nonlinear.SymbolicAD.simplify(f)
0.0 + 1.0 MOI.VariableIndex(6)² + 1.0 MOI.VariableIndex(7)² + 1.0 MOI.VariableIndex(9)²