Changeset b859805f in OpenModelica

Timestamp:

2019-09-25T11:31:20+02:00 (5 years ago)

Author:

Julien Schueller <schueller@…>

Parents:

77fc13e9

git-author:

Julien Schueller <schueller@…> (09/25/19 11:31:17)

git-committer:

Julien Schueller <schueller@…> (09/25/19 11:31:20)

Message:

Drop deprecated LAPACK apis

dgelsx, dgegv, dgeqpf are deprecated since lapack 3.1.1 (2007), remove them

they were also dropped in MSL:
​https://github.com/modelica/ModelicaStandardLibrary/pull/2857

File:

• OMCompiler/Compiler/runtime/lapackimpl.c (modified) (3 diffs)

OMCompiler/Compiler/runtime/lapackimpl.c

@@ -356 +356 @@
     void **outWORK, int *INFO)
 {
-#ifndef NO_LAPACK
-  integer n, lda, ldb, ldvl, ldvr, lwork, info = 0;
-  double *a, *b, *work, *alphar, *alphai, *beta, *vl, *vr;
-
-  n = N;
-  lda = LDA;
-  ldb = LDB;
-  ldvl = LDVL;
-  ldvr = LDVR;
-  lwork = LWORK;
-
-  a = alloc_real_matrix(lda, n, A);
-  b = alloc_real_matrix(ldb, n, B);
-  alphar = alloc_zeroed_real_vector(n);
-  alphai = alloc_zeroed_real_vector(n);
-  beta = alloc_zeroed_real_vector(n);
-  vl = alloc_zeroed_real_matrix(ldvl, n);
-  vr = alloc_zeroed_real_matrix(ldvl, n);
-  work = alloc_real_vector(lwork, inWORK);
-
-  dgegv_(&*jobvl, &*jobvr, &n, a, &lda, b, &ldb, alphar, alphai, beta, vl,
-    &ldvl, vr, &ldvr, work, &lwork, &info);
-
-  *ALPHAR = mk_rml_real_vector(n, alphar);
-  *ALPHAI = mk_rml_real_vector(n, alphai);
-  *BETA = mk_rml_real_vector(n, beta);
-  *VL = mk_rml_real_matrix(ldvl, n, vl);
-  *VR = mk_rml_real_matrix(ldvl, n, vr);
-  *outWORK = mk_rml_real_vector(lwork, work);
-  *INFO = info;
-
-  free(a);
-  free(b);
-  free(alphar);
-  free(alphai);
-  free(beta);
-  free(vl);
-  free(vr);
-  free(work);
-#else
-  MMC_THROW();
-#endif
+  MMC_THROW();
 }
 
@@ -438 +397 @@
     void **outA, void **outB, void **outJPVT, int *RANK, int *INFO)
 {
-#ifndef NO_LAPACK
-  integer m, n, nrhs, lda, ldb, rank = 0, info = 0, lwork;
-  double *a, *b, *work;
-  integer *jpvt;
-
-  m = M;
-  n = N;
-  nrhs = NRHS;
-  lda = LDA;
-  ldb = LDB;
-  lwork = (integer)fmax(fmin(M, N) + 3*N, 2*fmin(M, N) + nrhs);
-
-  a = alloc_real_matrix(lda, n, inA);
-  b = alloc_real_matrix(ldb, nrhs, inB);
-  work = alloc_real_vector(lwork, WORK);
-  jpvt = alloc_int_vector(n, inJPVT);
-
-  dgelsx_(&m, &n, &nrhs, a, &lda, b, &ldb, jpvt, &rcond, &rank, work, &info);
-
-  *outA = mk_rml_real_matrix(lda, n, a);
-  *outB = mk_rml_real_matrix(lda, nrhs, b);
-  *outJPVT = mk_rml_int_vector(n, jpvt);
-  *RANK = rank;
-  *INFO = info;
-
-  free(a);
-  free(b);
-  free(work);
-  free(jpvt);
-#else
-  MMC_THROW();
-#endif
+  MMC_THROW();
 }
 
@@ -793 +721 @@
     void *WORK, void **outA, void **outJPVT, void **TAU, int *INFO)
 {
-#ifndef NO_LAPACK
-  integer m, n, lda, lwork, ldtau, info = 0;
-  double *a, *tau, *work;
-  integer *jpvt;
-
-  m = M;
-  n = N;
-  lda = LDA;
-  lwork = 3 * n;
-  ldtau = (m < n ? m : n);
-
-  a = alloc_real_matrix(lda, n, inA);
-  jpvt = alloc_int_vector(n, inJPVT);
-  tau = alloc_zeroed_real_vector(ldtau);
-  work = alloc_real_vector(lwork, WORK);
-
-  dgeqpf_(&m, &n, a, &lda, jpvt, tau, work, &info);
-
-  *outA = mk_rml_real_matrix(lda, n, a);
-  *outJPVT = mk_rml_int_vector(n, jpvt);
-  *TAU = mk_rml_real_vector(ldtau, tau);
-  *INFO = info;
-
-  free(a);
-  free(jpvt);
-  free(tau);
-  free(work);
-#else
-  MMC_THROW();
-#endif
+  MMC_THROW();
 }