def test_split_compression_zone(compr_zone, neutral_axis, A_gross, desired):

    # ----- Setup --------
    # Turn tuple of desired polygons(s) into collection (MultiPolygon object)
    desired = MultiPolygon(desired)

    # ----- Exercise -----
    actual = su.split_compression_zone(compr_zone, neutral_axis, A_gross)

    # Turn tuple of actual polygon(s) into collection (MultiPolygon object)
    actual = MultiPolygon(actual)

    # ----- Verify -------
    # Compare by shapely built-in method
    # FIXME Find out if this is actually a correct way to check if they are equal
    assert actual.equals(desired)