Force OpenBLAS Haswell kernels for COLMAP

colmap_sfm.py

@@ -75,6 +75,12 @@ def run_sfm(
             "Upload several overlapping photos of the same scene."
         )
     log(f"[colmap] {len(names)} images prepared at <= {max_image_size}px")
+    try:
+        import os as _os
+        cpu = next((l.split(":", 1)[1].strip() for l in open("/proc/cpuinfo") if l.startswith("model name")), "?")
+        log(f"[colmap] cpu={cpu} ncpu={_os.cpu_count()} OPENBLAS_CORETYPE={_os.environ.get('OPENBLAS_CORETYPE')}")
+    except Exception:
+        pass
 
     database_path = workdir / "colmap.db"
     if database_path.exists():
@@ -201,6 +207,9 @@ def run_sfm_subprocess(image_paths, workdir, max_image_size=1280, matcher="exhau
     env = {k: v for k, v in os.environ.items() if k not in ("LD_PRELOAD",)}
     env["LD_LIBRARY_PATH"] = os.environ.get("_ORIG_LD_LIBRARY_PATH", "")
     env["OPENBLAS_NUM_THREADS"] = "1"
+    # Force a well-tested OpenBLAS kernel set: auto-detected kernels on some Space CPUs
+    # give garbage results (every image fails to register after the initial pair).
+    env["OPENBLAS_CORETYPE"] = "Haswell"
     env.pop("CUDA_VISIBLE_DEVICES", None)
     proc = subprocess.run([sys.executable, os.path.abspath(__file__), str(spec)], env=env,
                           capture_output=True, text=True, timeout=timeout)