Purity calculation