Problem. Evaluation of Hoppe94 Loop Subdivision Surfaces

Given a point P(k, v, w) in the initial control mesh, where k is the triangle index, (v, w) is the corresponding barycentric coordiants(u+v+w=1). What is it posiition on the limint surface?

Papers give use the answer:
[1] J.Stam. Exact evaluation of catmull-clark subdivision surfaces at arbitrary parameter values. 1998.
[2] J.Stam. Evaluation of Loop Subdivision Surfaces.
I have implemented the [2], which is a trivial job, for most data (I use the lpdata50Nt.dat) can be downloaded from J.Stam's webpage.
However, if the loop masks are changed, for example, the h.hoppe97, then Stam98's evaluation method and data can not be used directly, because he did not take the feature types into consideration. So.... a big question, i'm still strugglling with it:-)
[3] Denis Zorin.Evaluation of Piecewise Smooth Subdivision Surfaces.2002