More comments

Author: odow

src/Nonlinear/operators.jl

@@ -625,6 +625,7 @@ end
 
 _nan_to_zero(x) = isnan(x) ? 0.0 : x
 
+# No docstring because this function is still a WIP.
 function eval_multivariate_hessian(
     registry::OperatorRegistry,
     op::Symbol,
@@ -657,8 +658,8 @@ function eval_multivariate_hessian(
         # ∇f(x)  = x[2]*x[1]^(x[2]-1)
         #          x[1]^x[2]*log(x[1])
         #
-        # ∇²f(x) = x[2]*(x[2]-1)*x[1]^(x[2]-2) x[1]^(x[2]-1)*(x[2]*log(x[1])+1)
-        #                      .               x[1]^x[2]*log(x[1])^2
+        # ∇²f(x) = x[2]*(x[2]-1)*x[1]^(x[2]-2)
+        #          x[1]^(x[2]-1)*(x[2]*log(x[1])+1) x[1]^x[2]*log(x[1])^2
         ln = x[1] > 0 ? log(x[1]) : NaN
         if x[2] == one(T)
             H[2, 1] = _nan_to_zero(ln + one(T))
@@ -678,8 +679,8 @@ function eval_multivariate_hessian(
         # ∇f(x) = 1/x[2]
         #         -x[1]/x[2]^2
         #
-        # ∇²(x) = 0.0          -1/x[2]^2
-        #          .           2x[1]/x[2]^3
+        # ∇²(x) = 0.0
+        #         -1/x[2]^2 2x[1]/x[2]^3
         d = 1 / x[2]^2
         H[2, 1] = -d
         H[2, 2] = 2 * x[1] * d / x[2]

src/Nonlinear/types.jl

@@ -151,6 +151,8 @@ mutable struct NonlinearData <: MOI.AbstractNLPEvaluator
     last_constraint_index::Int64
     # Fields for initialize
     julia_expressions::Vector{Any}  # Any because expressions may be constants
+    # ordered_constraints is needed because `OrderedDict` doesn't support
+    # looking up a key by the linear index.
     ordered_constraints::Vector{ConstraintIndex}
     constraint_dual::Vector{Float64}
     # Timers