def METHOD_NAME(self):
    L_actual = mods.circles_and_rectangles(self.net)
    L_desired = np.array([[0.15      , 0.44688711, 0.40311289],
                          [0.40311289, 0.30965898, 0.28722813]])
    assert_allclose(L_actual, L_desired)
    # Incompatible data with model assumptions
    self.net["pore.diameter"][1] = 0.3
    with pytest.raises(Exception):
        L_actual = mods.circles_and_squares(self.net)
    self.net["pore.diameter"][1] = 0.9
    # Overlapping pores
    self.net["pore.diameter"][0] = 1.2
    L_actual = mods.circles_and_rectangles(self.net)
    L_desired = np.array([[5.78750000e-01, 1.00000000e-15, 4.21250000e-01],
                          [4.03112887e-01, 3.09658980e-01, 2.87228132e-01]])
    assert_allclose(L_actual, L_desired)
    self.net["pore.diameter"][0] = 0.5