lr = LRSplineSurface([2,2], [3,3]); lr.refine(1:size(lr.elements,1), 'elements'); lr.refine(1:size(lr.elements,1), 'elements'); does not give uniform refinement even if it should.