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