----------------
- ENH: Optimization: parallelized the bound constrained conjugate gradient without restart
- ENH: Optimization: implement bound constrained conjugate gradient without restart when the constrained set changes
- ENH: Optimization: implement bound constrained conjugate gradient with restart when the constrained set changes