We model a cubic-to-tetragonal martensitic transition by a Ginzburg-Landau free energy in the symmetric strain tensor. We show in three dimensions (3D) that solving the St. Venant compatibility relations for strain, treated as independent field equations, generates three anisotropic long-range potentials between the two order parameter components. These potentials encode 3D discrete symmetries, express the energetics of lattice integrity, and determine 3D textures. Simulation predictions include twins with temperature-varying orientation, helical twins, competing metastable states, and compatibility-induced elastic frustration. Our approach also applies to improper ferroelastics.