Merge branch 'main' into plug_gaps

Author: Quafadas

build.mill

@@ -39,8 +39,8 @@ object V:
   val scalaTags: Dep = mvn"com.lihaoyi::scalatags::0.13.1"
   val scalaVersion = "3.7.2"
   val munitVersion = "1.1.1"
-  val blas: Dep = mvn"dev.ludovic.netlib:blas:3.1.1"
-  val lapack: Dep = mvn"dev.ludovic.netlib:lapack:3.1.1"
+  val blas: Dep = mvn"dev.ludovic.netlib:blas:3.2.0"
+  val lapack: Dep = mvn"dev.ludovic.netlib:lapack:3.2.0"
   val scalaJavaTime: Dep = mvn"io.github.cquiroz::scala-java-time::2.6.0"
   val catsVersion = "2.13.0"
   val disciplineVersion = "2.0.0"

vecxt/src-jvm/svd.scala

@@ -157,7 +157,10 @@ object Svd:
     val ldu = math.max(1, m)
     val ldvt = math.max(1, vtRows)
 
-    val uArr = Array.ofDim[Double](ldu * uCols)
+    // Allocate ldu * m (not ldu * uCols) because netlib validates u.length >= ldu * m
+    // regardless of ReducedSVD mode (changed in netlib 3.2.0). Only the first uCols
+    // columns contain meaningful data; the Matrix wrapper uses uCols as the column count.
+    val uArr = Array.ofDim[Double](ldu * m)
     val vtArr = Array.ofDim[Double](ldvt * n)
 
     val workQuery = Array.ofDim[Double](1)