Polynomial-based non-uniform interpolatory subdivision with features control