The singularity at the origin for this function has eight tangent directions. The function is extremely flat near the singularity, so this model is a good test-case for singularity-detection.
This function can also be defined in the polar form as r=a sin(2 theta)