Local point constraints used during fitting.
These constraints are compiled into rows of Aeq * xi = beq or Aineq * xi <= bineq by evaluating the spline basis or its derivatives at specified points.
Aeq * xi = beq
Aineq * xi <= bineq