(sx1, sx2, atol)
| 114 | |
| 115 | |
| 116 | def same_geometry(sx1, sx2, atol): |
| 117 | geom1 = sx1.ReadGeometry() |
| 118 | geom2 = sx2.ReadGeometry() |
| 119 | if geom1.shape[0] != geom2.shape[0]: |
| 120 | return False |
| 121 | else: |
| 122 | return np.all(np.isclose(geom1, geom2, rtol=1e-3, atol=atol)) |
| 123 | |
| 124 | |
| 125 | def compute_areas(geom, connect): |
no outgoing calls
no test coverage detected