hcnm function

Hierarchical Constrained Newton method