def METHOD_NAME(self):
    L_actual = mods.trapezoids_and_rectangles(self.net)
    L_desired = np.array([[0.25, 0.3, 0.45],
                          [0.45, 0.2, 0.35]])
    assert_allclose(L_actual, L_desired)
    # Overlapping pores
    self.net["pore.diameter"][0] = 1.2
    L_actual = mods.trapezoids_and_rectangles(self.net)
    L_desired = np.array([[5.7875e-01, 1.0000e-15, 4.2125e-01],
                          [4.5000e-01, 2.0000e-01, 3.5000e-01]])
    assert_allclose(L_actual, L_desired)
    self.net["pore.diameter"][0] = 0.5