**sf** fails with Overlay-NG (GEOS 3.9.0)

[rsbivand]

testthat/test_geos.R fails at lines 213, 214, because st_difference() output vertices are ordered differently with Overlay-NG:

> out1[[2]][[1]]
     [,1] [,2]
[1,]  0.0  0.5
[2,]  1.0  1.5
[3,]  2.0  0.5
[4,]  1.5  0.5
[5,]  1.0  1.0
[6,]  0.5  0.5
[7,]  0.0  0.5
> correct_geom[[2]][[1]]
     [,1] [,2]
[1,]  0.5  0.5
[2,]  0.0  0.5
[3,]  1.0  1.5
[4,]  2.0  0.5
[5,]  1.5  0.5
[6,]  1.0  1.0
[7,]  0.5  0.5

[rsbivand]

Do you use CPL_geos_normalize() always in returning GEOS geometry output? See https://lists.osgeo.org/pipermail/geos-devel/2020-October/009758.html