math/blas, math/cblas, math/lapack, math/lapacke: upgrade to 3.12.1

Release notes at https://github.com/Reference-LAPACK/lapack/releases/tag/v3.12.1

Remark: many patches are needed, see
https://github.com/Reference-LAPACK/lapack/pull/1093

13 files changed, 257 insertions, 10 deletions

diff --git a/math/lapack/files/patch-CMakeLists.txt b/math/lapack/files/patch-CMakeLists.txt
deleted file mode 100644
index db50ddd14eee..000000000000
--- a/math/lapack/files/patch-CMakeLists.txt
+++ /dev/null
@@ -1,10 +0,0 @@
---- CMakeLists.txt.orig	2023-11-24 20:41:15 UTC
-+++ CMakeLists.txt
-@@ -223,6 +223,7 @@ if(BLAS_LIBRARIES)
- 
- # Check the usage of the user provided BLAS libraries
- if(BLAS_LIBRARIES)
-+  enable_language(Fortran)
-   include(CheckFortranFunctionExists)
-   set(CMAKE_REQUIRED_LIBRARIES ${BLAS_LIBRARIES})
-   CHECK_FORTRAN_FUNCTION_EXISTS("dgemm" BLAS_FOUND)
diff --git a/math/lapack/files/patch-SRC_DEPRECATED_cgelqs.f b/math/lapack/files/patch-SRC_DEPRECATED_cgelqs.f
new file mode 100644
index 000000000000..1d842f5ddc7e
--- /dev/null
+++ b/math/lapack/files/patch-SRC_DEPRECATED_cgelqs.f
@@ -0,0 +1,24 @@
+--- SRC/DEPRECATED/cgelqs.f.orig	2024-12-03 11:39:11 UTC
++++ SRC/DEPRECATED/cgelqs.f
+@@ -15,8 +15,9 @@
+ *       INTEGER            INFO, LDA, LDB, LWORK, M, N, NRHS
+ *       ..
+ *       .. Array Arguments ..
++* Patch from https://github.com/Reference-LAPACK/lapack/pull/1094
+ *       COMPLEX            A( LDA, * ), B( LDB, * ), TAU( * ),
+-*      $                   WORK( LWORK )
++*      $                   WORK( * )
+ *       ..
+ *
+ *
+@@ -127,8 +128,9 @@
+       INTEGER            INFO, LDA, LDB, LWORK, M, N, NRHS
+ *     ..
+ *     .. Array Arguments ..
++* Patch from https://github.com/Reference-LAPACK/lapack/pull/1094
+       COMPLEX            A( LDA, * ), B( LDB, * ), TAU( * ),
+-     $                   WORK( LWORK )
++     $                   WORK( * )
+ *     ..
+ *
+ *  =====================================================================
diff --git a/math/lapack/files/patch-SRC_DEPRECATED_cgelsx.f b/math/lapack/files/patch-SRC_DEPRECATED_cgelsx.f
new file mode 100644
index 000000000000..512fc7cc5349
--- /dev/null
+++ b/math/lapack/files/patch-SRC_DEPRECATED_cgelsx.f
@@ -0,0 +1,13 @@
+--- SRC/DEPRECATED/cgelsx.f.orig	2024-12-03 11:39:11 UTC
++++ SRC/DEPRECATED/cgelsx.f
+@@ -364,8 +364,8 @@
+ *
+ *     B(1:M,1:NRHS) := Q**H * B(1:M,1:NRHS)
+ *
+-      CALL CUNM2R( 'Left', 'Conjugate transpose', M, NRHS, MN, A, LDA,
+-     $             WORK( 1 ), B, LDB, WORK( 2*MN+1 ), INFO )
++      CALL CUNM2R( 'Left', 'Conjugate transpose', M, NRHS, MN, A,
++     $              LDA, WORK( 1 ), B, LDB, WORK( 2*MN+1 ), INFO )
+ *
+ *     workspace NRHS
+ *
diff --git a/math/lapack/files/patch-SRC_DEPRECATED_cgeqrs.f b/math/lapack/files/patch-SRC_DEPRECATED_cgeqrs.f
new file mode 100644
index 000000000000..8f68964a6fa3
--- /dev/null
+++ b/math/lapack/files/patch-SRC_DEPRECATED_cgeqrs.f
@@ -0,0 +1,24 @@
+--- SRC/DEPRECATED/cgeqrs.f.orig	2024-12-03 11:39:11 UTC
++++ SRC/DEPRECATED/cgeqrs.f
+@@ -15,8 +15,9 @@
+ *       INTEGER            INFO, LDA, LDB, LWORK, M, N, NRHS
+ *       ..
+ *       .. Array Arguments ..
++* Patch from https://github.com/Reference-LAPACK/lapack/pull/1094
+ *       COMPLEX            A( LDA, * ), B( LDB, * ), TAU( * ),
+-*      $                   WORK( LWORK )
++*      $                   WORK( * )
+ *       ..
+ *
+ *
+@@ -127,8 +128,9 @@
+       INTEGER            INFO, LDA, LDB, LWORK, M, N, NRHS
+ *     ..
+ *     .. Array Arguments ..
++* Patch from https://github.com/Reference-LAPACK/lapack/pull/1094
+       COMPLEX            A( LDA, * ), B( LDB, * ), TAU( * ),
+-     $                   WORK( LWORK )
++     $                   WORK( * )
+ *     ..
+ *
+ *  =====================================================================
diff --git a/math/lapack/files/patch-SRC_DEPRECATED_dgelqs.f b/math/lapack/files/patch-SRC_DEPRECATED_dgelqs.f
new file mode 100644
index 000000000000..840d365db804
--- /dev/null
+++ b/math/lapack/files/patch-SRC_DEPRECATED_dgelqs.f
@@ -0,0 +1,24 @@
+--- SRC/DEPRECATED/dgelqs.f.orig	2024-12-03 11:39:11 UTC
++++ SRC/DEPRECATED/dgelqs.f
+@@ -15,8 +15,9 @@
+ *       INTEGER            INFO, LDA, LDB, LWORK, M, N, NRHS
+ *       ..
+ *       .. Array Arguments ..
++* Patch from https://github.com/Reference-LAPACK/lapack/pull/1094
+ *       DOUBLE PRECISION   A( LDA, * ), B( LDB, * ), TAU( * ),
+-*      $                   WORK( LWORK )
++*      $                   WORK( * )
+ *       ..
+ *
+ *
+@@ -127,8 +128,9 @@
+       INTEGER            INFO, LDA, LDB, LWORK, M, N, NRHS
+ *     ..
+ *     .. Array Arguments ..
++* Patch from https://github.com/Reference-LAPACK/lapack/pull/1094
+       DOUBLE PRECISION   A( LDA, * ), B( LDB, * ), TAU( * ),
+-     $                   WORK( LWORK )
++     $                   WORK( * )
+ *     ..
+ *
+ *  =====================================================================
diff --git a/math/lapack/files/patch-SRC_DEPRECATED_dgeqpf.f b/math/lapack/files/patch-SRC_DEPRECATED_dgeqpf.f
new file mode 100644
index 000000000000..0677f69762ed
--- /dev/null
+++ b/math/lapack/files/patch-SRC_DEPRECATED_dgeqpf.f
@@ -0,0 +1,13 @@
+--- SRC/DEPRECATED/dgeqpf.f.orig	2024-12-03 11:39:11 UTC
++++ SRC/DEPRECATED/dgeqpf.f
+@@ -218,8 +218,8 @@
+          MA = MIN( ITEMP, M )
+          CALL DGEQR2( M, MA, A, LDA, TAU, WORK, INFO )
+          IF( MA.LT.N ) THEN
+-            CALL DORM2R( 'Left', 'Transpose', M, N-MA, MA, A, LDA, TAU,
+-     $                   A( 1, MA+1 ), LDA, WORK, INFO )
++            CALL DORM2R( 'Left', 'Transpose', M, N-MA, MA, A, LDA,
++     $                    TAU, A( 1, MA+1 ), LDA, WORK, INFO )
+          END IF
+       END IF
+ *
diff --git a/math/lapack/files/patch-SRC_DEPRECATED_dgeqrs.f b/math/lapack/files/patch-SRC_DEPRECATED_dgeqrs.f
new file mode 100644
index 000000000000..f48a6bbcf312
--- /dev/null
+++ b/math/lapack/files/patch-SRC_DEPRECATED_dgeqrs.f
@@ -0,0 +1,24 @@
+--- SRC/DEPRECATED/dgeqrs.f.orig	2024-12-03 11:39:11 UTC
++++ SRC/DEPRECATED/dgeqrs.f
+@@ -15,8 +15,9 @@
+ *       INTEGER            INFO, LDA, LDB, LWORK, M, N, NRHS
+ *       ..
+ *       .. Array Arguments ..
++* Patch from https://github.com/Reference-LAPACK/lapack/pull/1094
+ *       DOUBLE PRECISION   A( LDA, * ), B( LDB, * ), TAU( * ),
+-*      $                   WORK( LWORK )
++*      $                   WORK( * )
+ *       ..
+ *
+ *
+@@ -127,8 +128,9 @@
+       INTEGER            INFO, LDA, LDB, LWORK, M, N, NRHS
+ *     ..
+ *     .. Array Arguments ..
++* Patch from https://github.com/Reference-LAPACK/lapack/pull/1094
+       DOUBLE PRECISION   A( LDA, * ), B( LDB, * ), TAU( * ),
+-     $                   WORK( LWORK )
++     $                   WORK( * )
+ *     ..
+ *
+ *  =====================================================================
diff --git a/math/lapack/files/patch-SRC_DEPRECATED_dlahrd.f b/math/lapack/files/patch-SRC_DEPRECATED_dlahrd.f
new file mode 100644
index 000000000000..291efc2ead01
--- /dev/null
+++ b/math/lapack/files/patch-SRC_DEPRECATED_dlahrd.f
@@ -0,0 +1,26 @@
+--- SRC/DEPRECATED/dlahrd.f.orig	2024-12-03 11:39:11 UTC
++++ SRC/DEPRECATED/dlahrd.f
+@@ -231,8 +231,9 @@
+ *
+ *           w := T**T *w
+ *
+-            CALL DTRMV( 'Upper', 'Transpose', 'Non-unit', I-1, T, LDT,
+-     $                  T( 1, NB ), 1 )
++* Patch from https://github.com/Reference-LAPACK/lapack/pull/1093
++            CALL DTRMV( 'Upper', 'Transpose', 'Non-unit', I-1,
++     $                  T, LDT, T( 1, NB ), 1 )
+ *
+ *           b2 := b2 - V2*w
+ *
+@@ -270,8 +271,9 @@
+ *        Compute T(1:i,i)
+ *
+          CALL DSCAL( I-1, -TAU( I ), T( 1, I ), 1 )
+-         CALL DTRMV( 'Upper', 'No transpose', 'Non-unit', I-1, T, LDT,
+-     $               T( 1, I ), 1 )
++* Patch from https://github.com/Reference-LAPACK/lapack/pull/1093
++         CALL DTRMV( 'Upper', 'No transpose', 'Non-unit', I-1,
++     $               T, LDT, T( 1, I ), 1 )
+          T( I, I ) = TAU( I )
+ *
+    10 CONTINUE
diff --git a/math/lapack/files/patch-SRC_DEPRECATED_sgelqs.f b/math/lapack/files/patch-SRC_DEPRECATED_sgelqs.f
new file mode 100644
index 000000000000..6f54e915d803
--- /dev/null
+++ b/math/lapack/files/patch-SRC_DEPRECATED_sgelqs.f
@@ -0,0 +1,24 @@
+--- SRC/DEPRECATED/sgelqs.f.orig	2024-12-03 11:39:11 UTC
++++ SRC/DEPRECATED/sgelqs.f
+@@ -15,8 +15,9 @@
+ *       INTEGER            INFO, LDA, LDB, LWORK, M, N, NRHS
+ *       ..
+ *       .. Array Arguments ..
++* Patch from https://github.com/Reference-LAPACK/lapack/pull/1094
+ *       REAL               A( LDA, * ), B( LDB, * ), TAU( * ),
+-*      $                   WORK( LWORK )
++*      $                   WORK( * )
+ *       ..
+ *
+ *
+@@ -127,8 +128,9 @@
+       INTEGER            INFO, LDA, LDB, LWORK, M, N, NRHS
+ *     ..
+ *     .. Array Arguments ..
++* Patch from https://github.com/Reference-LAPACK/lapack/pull/1094
+       REAL               A( LDA, * ), B( LDB, * ), TAU( * ),
+-     $                   WORK( LWORK )
++     $                   WORK( * )
+ *     ..
+ *
+ *  =====================================================================
diff --git a/math/lapack/files/patch-SRC_DEPRECATED_sgeqpf.f b/math/lapack/files/patch-SRC_DEPRECATED_sgeqpf.f
new file mode 100644
index 000000000000..b8c568852a78
--- /dev/null
+++ b/math/lapack/files/patch-SRC_DEPRECATED_sgeqpf.f
@@ -0,0 +1,13 @@
+--- SRC/DEPRECATED/sgeqpf.f.orig	2024-12-03 11:39:11 UTC
++++ SRC/DEPRECATED/sgeqpf.f
+@@ -218,8 +218,8 @@
+          MA = MIN( ITEMP, M )
+          CALL SGEQR2( M, MA, A, LDA, TAU, WORK, INFO )
+          IF( MA.LT.N ) THEN
+-            CALL SORM2R( 'Left', 'Transpose', M, N-MA, MA, A, LDA, TAU,
+-     $                   A( 1, MA+1 ), LDA, WORK, INFO )
++            CALL SORM2R( 'Left', 'Transpose', M, N-MA, MA, A, LDA,
++     $                    TAU, A( 1, MA+1 ), LDA, WORK, INFO )
+          END IF
+       END IF
+ *
diff --git a/math/lapack/files/patch-SRC_DEPRECATED_sgeqrs.f b/math/lapack/files/patch-SRC_DEPRECATED_sgeqrs.f
new file mode 100644
index 000000000000..8d7017109436
--- /dev/null
+++ b/math/lapack/files/patch-SRC_DEPRECATED_sgeqrs.f
@@ -0,0 +1,24 @@
+--- SRC/DEPRECATED/sgeqrs.f.orig	2024-12-03 11:39:11 UTC
++++ SRC/DEPRECATED/sgeqrs.f
+@@ -15,8 +15,9 @@
+ *       INTEGER            INFO, LDA, LDB, LWORK, M, N, NRHS
+ *       ..
+ *       .. Array Arguments ..
++* Patch from https://github.com/Reference-LAPACK/lapack/pull/1094
+ *       REAL               A( LDA, * ), B( LDB, * ), TAU( * ),
+-*      $                   WORK( LWORK )
++*      $                   WORK( * )
+ *       ..
+ *
+ *
+@@ -127,8 +128,9 @@
+       INTEGER            INFO, LDA, LDB, LWORK, M, N, NRHS
+ *     ..
+ *     .. Array Arguments ..
++* Patch from https://github.com/Reference-LAPACK/lapack/pull/1094
+       REAL               A( LDA, * ), B( LDB, * ), TAU( * ),
+-     $                   WORK( LWORK )
++     $                   WORK( * )
+ *     ..
+ *
+ *  =====================================================================
diff --git a/math/lapack/files/patch-SRC_DEPRECATED_zgelqs.f b/math/lapack/files/patch-SRC_DEPRECATED_zgelqs.f
new file mode 100644
index 000000000000..9967747bbeeb
--- /dev/null
+++ b/math/lapack/files/patch-SRC_DEPRECATED_zgelqs.f
@@ -0,0 +1,24 @@
+--- SRC/DEPRECATED/zgelqs.f.orig	2024-12-03 11:39:11 UTC
++++ SRC/DEPRECATED/zgelqs.f
+@@ -15,8 +15,9 @@
+ *       INTEGER            INFO, LDA, LDB, LWORK, M, N, NRHS
+ *       ..
+ *       .. Array Arguments ..
++* Patch from https://github.com/Reference-LAPACK/lapack/pull/1094
+ *       COMPLEX*16         A( LDA, * ), B( LDB, * ), TAU( * ),
+-*      $                   WORK( LWORK )
++*      $                   WORK( * )
+ *       ..
+ *
+ *
+@@ -127,8 +128,9 @@
+       INTEGER            INFO, LDA, LDB, LWORK, M, N, NRHS
+ *     ..
+ *     .. Array Arguments ..
++* Patch from https://github.com/Reference-LAPACK/lapack/pull/1094
+       COMPLEX*16         A( LDA, * ), B( LDB, * ), TAU( * ),
+-     $                   WORK( LWORK )
++     $                   WORK( * )
+ *     ..
+ *
+ *  =====================================================================
diff --git a/math/lapack/files/patch-SRC_DEPRECATED_zgeqrs.f b/math/lapack/files/patch-SRC_DEPRECATED_zgeqrs.f
new file mode 100644
index 000000000000..8bbac0dd3e47
--- /dev/null
+++ b/math/lapack/files/patch-SRC_DEPRECATED_zgeqrs.f
@@ -0,0 +1,24 @@
+--- SRC/DEPRECATED/zgeqrs.f.orig	2024-12-03 11:39:11 UTC
++++ SRC/DEPRECATED/zgeqrs.f
+@@ -15,8 +15,9 @@
+ *       INTEGER            INFO, LDA, LDB, LWORK, M, N, NRHS
+ *       ..
+ *       .. Array Arguments ..
++* Patch from https://github.com/Reference-LAPACK/lapack/pull/1094
+ *       COMPLEX*16         A( LDA, * ), B( LDB, * ), TAU( * ),
+-*      $                   WORK( LWORK )
++*      $                   WORK( * )
+ *       ..
+ *
+ *
+@@ -127,8 +128,9 @@
+       INTEGER            INFO, LDA, LDB, LWORK, M, N, NRHS
+ *     ..
+ *     .. Array Arguments ..
++* Patch from https://github.com/Reference-LAPACK/lapack/pull/1094
+       COMPLEX*16         A( LDA, * ), B( LDB, * ), TAU( * ),
+-     $                   WORK( LWORK )
++     $                   WORK( * )
+ *     ..
+ *
+ *  =====================================================================