Currently, the autodiff backend in internal optimization solver can only be AutoFiniteDiff(), because the value update in EvalSol after collocation equation evaluation in loss function __mirk_loss!, however, such operation prevent proper sparsity detection with SparseConnectivityTracer.jl. We need to restruct the loss funtion being passed to optimization solver so that more autodiff backends can be used.

This issue can affect #374