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