Surfaces containing two circles through each point and decomposition of quaternionic matrices

We find all analytic surfaces in space R^3 such that through each point of the surface one can draw two circular arcs fully contained in the surface. The proof uses a new decomposition technique for quaternionic matrices.