Basically I think all you need to do to increase numerical precision is add a conditional here, where for `continuous` you just do `normal = center / norm(center)`, or (for `fs`) execute the next six lines. Right?