def format_dipole_derivatives(property):
    """
    Reshape the extracted dipole derivatives to the correct
    format. Format is Natoms x (dx dy dz) x (property x y z)

    Args:
        property (numpy.array):

    Returns:
        numpy.array: Reshaped array
    """
    N, _ = property.shape
    N = N // 3
    property = property.reshape(N, 3, 3)
    return property