QHull error with newer scipy versions when running the multispectral part of SICOR.
This seems to come from a QHull call:
[...]
- r143311 tested simplicialtop
vertices: p3548(v3549) p3798(v3069) p4073(v2896) p1547(v1820)
between f378675 and f247
- r143310 tested simplicialtop
vertices: p3548(v3549) p4073(v2896) p3573(v2367) p1547(v1820)
between f378672 and f247
- r141959 tested simplicialtop simplicialbot
vertices: p2599(v3522) p2124(v3459) p3349(v1492) p848(v332)
between f247 and f379330
- r139275
vertices: p699(v3462) p623(v3374) p398(v2411) p199(v559)
between f379723 and f247
- r135507
vertices: p623(v3374) p398(v2411) p199(v559) p174(v231)
between f379723 and f247
- r139273 nonconvex
vertices: p699(v3462) p623(v3374) p1149(v1074) p199(v559)
between f247 and f379734
- r139274 nonconvex simplicialbot
vertices: p699(v3462) p623(v3374) p1199(v2944) p1149(v1074)
between f247 and f379734
- r135540 nonconvex
vertices: p623(v3374) p1199(v2944) p1424(v2426) p1149(v1074)
between f379734 and f247
- r135497 tested simplicialtop
vertices: p623(v3374) p649(v3026) p149(v1184) p174(v231)
between f379765 and f247
- r135508 tested nonconvex simplicialbot
vertices: p623(v3374) p1149(v1074) p199(v559) p174(v231)
between f247 and f379737
ERRONEOUS OTHER FACET:
- f379734
- flags: bottom newfacet keepcentrum newmerge
- merges: 2
- normal: 0.9101 0.004034 7.039e-05 2.129e-05 -0.4144
- offset: -573.2704
- maxoutside: 1.713715e-09
- vertices: p2469(v3561) p699(v3462) p623(v3374) p1199(v2944) p1424(v2426) p1149(v1074) p199(v559)
- neighboring facets: f247 f318154 f464 f379737 f379707 f379723
- ridges:
- r139274 nonconvex simplicialbot
vertices: p699(v3462) p623(v3374) p1199(v2944) p1149(v1074)
between f247 and f379734
- r139273 nonconvex
vertices: p699(v3462) p623(v3374) p1149(v1074) p199(v559)
between f247 and f379734
- r135540 nonconvex
vertices: p623(v3374) p1199(v2944) p1424(v2426) p1149(v1074)
between f379734 and f247
- r143986
vertices: p2469(v3561) p623(v3374) p1424(v2426) p1149(v1074)
between f379734 and f318154
- r143862 simplicialbot
vertices: p2469(v3561) p699(v3462) p1199(v2944) p1149(v1074)
between f464 and f379734
- r143853
vertices: p2469(v3561) p699(v3462) p1149(v1074) p199(v559)
between f464 and f379734
- r143866
vertices: p2469(v3561) p1199(v2944) p1424(v2426) p1149(v1074)
between f464 and f379734
- r144069 simplicialbot
vertices: p2469(v3561) p623(v3374) p1149(v1074) p199(v559)
between f379734 and f379737
- r144053 simplicialbot
vertices: p2469(v3561) p699(v3462) p623(v3374) p1199(v2944)
between f379707 and f379734
- r144058
vertices: p2469(v3561) p623(v3374) p1199(v2944) p1424(v2426)
between f379734 and f379707
- r144036
vertices: p2469(v3561) p699(v3462) p623(v3374) p199(v559)
between f379734 and f379723
While executing: | qhull d Qc Qbb Qt Q12 Qz
Options selected for Qhull 2019.1.r 2019/06/21:
run-id 380203081 delaunay Qcoplanar-keep Qbbound-last Qtriangulate
Q12-allow-wide Qz-infinity-point _pre-merge _zero-centrum Qinterior-keep
Pgood _max-width 5.3e+02 Error-roundoff 2.6e-12 _one-merge 2.9e-11
Visible-distance 1.6e-11 U-max-coplanar 1.6e-11 Width-outside 3.1e-11
_wide-facet 9.4e-11 _maxoutside 3.8e-08
Last point added to hull was p2469. Last merge was #45690.
At error exit:
Delaunay triangulation by the convex hull of 5001 points in 5-d:
Number of input sites and at-infinity: 3560
Total number of deleted points due to merging: 1
Number of nearly incident points: 148
Number of Delaunay regions: 83584
Number of non-simplicial Delaunay regions: 2013
Statistics for: | qhull d Qc Qbb Qt Q12 Qz
Number of points processed: 3561
Number of hyperplanes created: 339676
Number of facets in hull: 84760
Number of distance tests for qhull: 25946002
Number of distance tests for merging: 962588942
Number of distance tests for checking: 0
Number of merged facets: 46712
Maximum distance of point above facet: 3.8e-08 (1231.8x)
Maximum distance of vertex below facet: -0.0018 (58807128.2x)
precision problems (corrected unless 'Q0' or an error)
108 flipped facets
4402 coplanar horizon facets for new vertices
148 coplanar points during partitioning
295 nearly singular or axis-parallel hyperplanes
295 zero divisors during back substitute
319 zero divisors during gaussian elimination
7 dupridges with multiple neighbors
10 dupridges with flip facet into good neighbor
A wide merge error has occurred. Qhull has produced a wide facet due to facet merges and vertex merges.
This usually occurs when the input is nearly degenerate and substantial merging has occurred.
See http://www.qhull.org/html/qh-impre.htm#limit
As suggested here, downgrading scipy to version 1.1.0 fixes the issue. However, this should at least be reflected in the setup.py and the requirements files. The better solution would be to find the actual cause of the error - pinning scipy might cause incompatibilities in future Python environments.