rlatvec.append(
            np.array(2 * math.pi * np.cross(latvec[0, :], latvec[1, :]) / volume)