changed: better jacobian and hessian defaults for the sparse case

Author: franckgaga

src/controller/nonlinmpc.jl

@@ -5,7 +5,7 @@ const DEFAULT_NONLINMPC_JACDENSE  = AutoForwardDiff()
 const DEFAULT_NONLINMPC_JACSPARSE = AutoSparse(
     AutoForwardDiff();
     sparsity_detector=TracerSparsityDetector(),
-    coloring_algorithm=GreedyColoringAlgorithm(),
+    coloring_algorithm=GreedyColoringAlgorithm(ALL_COLORING_ORDERS),
 )
 const DEFAULT_NONLINMPC_HESSIAN = DEFAULT_NONLINMPC_JACSPARSE
 
@@ -291,10 +291,11 @@ NonLinMPC controller with a sample time Ts = 10.0 s:
     AutoSparse(
         AutoForwardDiff(); 
         sparsity_detector  = TracerSparsityDetector(), 
-        coloring_algorithm = GreedyColoringAlgorithm()
+        coloring_algorithm = GreedyColoringAlgorithm(ALL_COLORING_ORDERS)
     )
     ```
-    This is also the sparse backend selected for the Hessian of the Lagrangian function if 
+    that is, it will test all coloring orders at preparation and keep the best. This is
+    also the sparse backend selected for the Hessian of the Lagrangian function if 
     `oracle=true` and `hessian=true`, which is the second exception. Second order 
     derivatives are only supported with `oracle=true` option.
     

src/estimator/mhe/construct.jl

@@ -5,7 +5,7 @@ const DEFAULT_NONLINMHE_JACOBIAN  = AutoForwardDiff()
 const DEFAULT_NONLINMHE_HESSIAN   = AutoSparse(
     AutoForwardDiff();
     sparsity_detector=TracerSparsityDetector(),
-    coloring_algorithm=GreedyColoringAlgorithm(),
+    coloring_algorithm=GreedyColoringAlgorithm(ALL_COLORING_ORDERS),
 ) 
 
 @doc raw"""
@@ -382,9 +382,11 @@ MovingHorizonEstimator estimator with a sample time Ts = 10.0 s:
     AutoSparse(
         AutoForwardDiff(); 
         sparsity_detector  = TracerSparsityDetector(), 
-        coloring_algorithm = GreedyColoringAlgorithm()
+        coloring_algorithm = GreedyColoringAlgorithm(ALL_COLORING_ORDERS)
     )
     ```
+    that is, it will test all coloring orders at preparation and keep the best. 
+    
     The slack variable ``ε`` relaxes the constraints if enabled, see [`setconstraint!`](@ref). 
     It is disabled by default for the MHE (from `Cwt=Inf`) but it should be activated for
     problems with two or more types of bounds, to ensure feasibility (e.g. on the estimated

src/general.jl

@@ -6,6 +6,15 @@ const DEFAULT_LWT = 0.0
 const DEFAULT_CWT = 1e5
 const DEFAULT_EWT = 0.0
 
+"All deterministic algorithms for matrix coloring order in `SparseMatrixColoring.jl`."
+const ALL_COLORING_ORDERS = (
+    NaturalOrder(),
+    LargestFirst(),
+    SmallestLast(),
+    IncidenceDegree(),
+    DynamicLargestFirst(),
+)
+
 "Termination status that means 'no solution available'."
 const ERROR_STATUSES = (
     JuMP.INFEASIBLE, JuMP.DUAL_INFEASIBLE, JuMP.LOCALLY_INFEASIBLE,