[gdal-dev] Could vector reproject be affected if alternative BLAS/LAPACK implementations are in use?

[Chris Toney]

Hi,
Is it possible that vector reproject (via `GDALVectorTranslate()`)
could be affected if alternative BLAS/LAPACK implementations are
present (i.e, from Intel MKL or OpenBLAS)?

"ERROR 1: Reprojection failed" may originate from
https://github.com/OSGeo/gdal/blob/e48cc4f31f43cf872c036d2c4e65e7368166c34a/apps/ogr2ogr_lib.cpp#L6378.
This is followed by seg fault.

The error is reported with GDAL 3.10.3, GEOS 3.13.0, PROJ 9.6.2 on
Fedora 42. The question is motivated by the issue at:
https://github.com/USDAForestService/gdalraster/issues/769.

The same code/inputs test without issue on a wide range of different
platforms in several CI and package management scenarios. The error
message and seg fault appear to only occur (sometimes not always) when
the alternative BLAS/LAPACK are present.

The code that generates the error is at:
https://github.com/USDAForestService/gdalraster/blob/d7c620b74a987cbc3f96dec245d37809e09f858a/src/gdal_exp.cpp#L1791

Input and output datasets are GPKG and reprojection is from NAD83 GCS
(EPSG:4269) to EPSG:32100.

Thanks.

Chris Toney