Most sophisticated branching possible implemented where we compute the
cv function only on those variables not totally smoothed out.
Finalizing branching, needed to retain W to get the proper polynomial
but remoive the x's for which the h's/lambda's are at their upper
bound to avoid computing kernels for irrelevant variables Also fixed a
glitch with the lower bound for cateorical predictors