[cig-commits] r13135 - in short/3D/PyLith/trunk/libsrc: feassemble utils
brad at geodynamics.org
brad at geodynamics.org
Sun Oct 26 16:25:41 PDT 2008
Author: brad
Date: 2008-10-26 16:25:41 -0700 (Sun, 26 Oct 2008)
New Revision: 13135
Modified:
short/3D/PyLith/trunk/libsrc/feassemble/ElasticityImplicit.cc
short/3D/PyLith/trunk/libsrc/utils/lapack.h
Log:
Revert to non-PETSc interface to dgesvd.
Modified: short/3D/PyLith/trunk/libsrc/feassemble/ElasticityImplicit.cc
===================================================================
--- short/3D/PyLith/trunk/libsrc/feassemble/ElasticityImplicit.cc 2008-10-26 23:19:09 UTC (rev 13134)
+++ short/3D/PyLith/trunk/libsrc/feassemble/ElasticityImplicit.cc 2008-10-26 23:25:41 UTC (rev 13135)
@@ -421,9 +421,9 @@
const int n2 = n*n;
for (int i = 0; i < n2; ++i)
elemMat[i] = _cellMatrix[i];
- LAPACKdgesvd("N", "N", &n, &n, elemMat, &n, svalues,
- &sdummy, &idummy, &sdummy, &idummy, work,
- &lwork, &lierr);
+ lapack_dgesvd("N", "N", &n, &n, elemMat, &n, svalues,
+ &sdummy, &idummy, &sdummy, &idummy, work,
+ &lwork, &lierr);
if (lierr)
throw std::runtime_error("Lapack SVD failed");
minSV = svalues[n-1];
Modified: short/3D/PyLith/trunk/libsrc/utils/lapack.h
===================================================================
--- short/3D/PyLith/trunk/libsrc/utils/lapack.h 2008-10-26 23:19:09 UTC (rev 13134)
+++ short/3D/PyLith/trunk/libsrc/utils/lapack.h 2008-10-26 23:25:41 UTC (rev 13135)
@@ -20,21 +20,23 @@
#if !defined(pylith_utils_lapack_h)
#define pylith_utils_lapack_h
+#define lapack_dgesvd dgesvd_
+
extern "C" {
- EXTERN void LAPACKdgesvd(const char*,
- const char*,
- PetscBLASInt*,
- PetscBLASInt*,
- PetscScalar*,
- PetscBLASInt*,
- PetscReal*,
- PetscScalar*,
- PetscBLASInt*,
- PetscScalar*,
- PetscBLASInt*,
- PetscScalar*,
- PetscBLASInt*,
- PetscBLASInt*);
+ EXTERN void lapack_dgesvd(const char*,
+ const char*,
+ PetscBLASInt*,
+ PetscBLASInt*,
+ PetscScalar*,
+ PetscBLASInt*,
+ PetscReal*,
+ PetscScalar*,
+ PetscBLASInt*,
+ PetscScalar*,
+ PetscBLASInt*,
+ PetscScalar*,
+ PetscBLASInt*,
+ PetscBLASInt*);
}
#endif // pylith_utils_lapack_h
More information about the CIG-COMMITS
mailing list