Abstract: A standing wave pattern appears on the free surface of a fluid layer when it is subjected to vertical oscillation of sufficiently high amplitude. Like Taylor-Couette flow (TC) and Rayleigh-Benard convection (RB), the Faraday instability is one of the archetypical pattern-forming systems. Unlike TC and RB, the wavelength is controlled by the forcing frequency rather than by the fluid depth, making it easy to destabilize multiple wavelengths everywhere simultaneously. Starting in the 1990s, experimental realizations using this technique produced fascinating phenomena such as quasipatterns and superlattices which in turn led to new mathematical theories of pattern formation. Another difference is that the Faraday instability has been the subject of surprisingly little numerical study, lagging behind TC and RB by several decades. The first 3D simulation reproduced hexagonal standing waves, which were succeeded by long-time recurrent alternation between quasi-hexagonal and beaded striped patterns, interconnected by spatio-temporal symmetries. In a large domain, a supersquare is observed in which diagonal subsquares are synchronized. A liquid drop subjected to an oscillatory radial force comprises a spherical version of the Faraday instability. Simulations show Platonic solids alternating with their duals while drifting.