Characterizing modulated structures with first-principles calculations: a unified superspace scheme of ordering in mullite