fmm bug

[dialuser]

I got the following RunTimeError when trying to run GeoNet fastmarching on Ubuntu Linux. I wonder if anyone has seen this before.

The printout is pasted below;

cc@geoml:/data/GeoNet$ python pygeonet_fast_marching.py
Average Computation Time (0.7 GB DEM): 5 - 7 min
Computing percentage drainage area of each indexed basin
Calc Time: 1.3560327330196742
[[ 13 33 113]
[ 66 46 0]]
Preparing to calculate cost function
normalizing curvature
Curvature min: 0.0 exp(min): 1.0
Curvature max: 1.0 exp(max): 20.085536923187668
Calculating local costs
Evaluating local cost func.
1/cost min: 1.0
1/cost max: 467626.4802614126
('writing geotiff', 'skunk_costfunction.tif')
Performing fast marching
start point : [13 66]
RAM usage before FMM 0: svmem(total=16818544640, available=16060346368, percent=4.5, used=412332032, free=7899185152, active=2472325120, inactive=5654888448, buffers=386818048, cached=8120209408, shared=741376, slab=722386944)
start point : [33 46]
RAM usage before FMM 1: svmem(total=16818544640, available=16026443776, percent=4.7, used=446234624, free=7865282560, active=2506403840, inactive=5654888448, buffers=386818048, cached=8120209408, shared=741376, slab=722386944)
Traceback (most recent call last):
File "pygeonet_fast_marching.py", line 220, in
main()
File "pygeonet_fast_marching.py", line 212, in main
geodesicDistanceArray = Fast_Marching(fastMarchingStartPointListFMM, basinIndexArray, flowArray, reciprocalLocalCostArray)
File "pygeonet_fast_marching.py", line 138, in Fast_Marching
travelTimearray = skfmm.travel_time(phi, speed, dx=.01)
File "/home/cc/anaconda3/envs/geonet/lib/python3.8/site-packages/skfmm/pfmm.py", line 172, in travel_time
t = cFastMarcher(phi, dx, flag, speed, ext_mask,
RuntimeError: Negative discriminant in time marcher quadratic.