A direct and local method for computing polynomial Pythagorean-normal patches with global G1 continuity. (September 2018)