]> https://gitweb.dealii.org/ - dealii-svn.git/commitdiff
Fixed an issue with define flag for Xgetri.
authorkronbichler <kronbichler@0785d39b-7218-0410-832d-ea1e28bc413d>
Wed, 19 Nov 2008 16:55:17 +0000 (16:55 +0000)
committerkronbichler <kronbichler@0785d39b-7218-0410-832d-ea1e28bc413d>
Wed, 19 Nov 2008 16:55:17 +0000 (16:55 +0000)
git-svn-id: https://svn.dealii.org/trunk@17644 0785d39b-7218-0410-832d-ea1e28bc413d

deal.II/base/include/base/config.h.in
deal.II/contrib/blastemplates/lapack_templates.h.in
deal.II/contrib/blastemplates/templates.pl
deal.II/lac/include/lac/full_matrix.h
deal.II/lac/include/lac/lapack_templates.h
deal.II/lac/source/full_matrix.cc

index 7e10dd2e47d441d9e5cd155e1566671c89ae6557..684d4b1f237842730a09b2fd8f6d7540dc48c840 100644 (file)
    linking. For the details, look at aclocal.m4 in the top-level directory. */
 #undef DEAL_II_WEAK_LINKAGE_BUG
 
+/* Define to 1 if you have the <Amesos.h> header file. */
+#undef HAVE_AMESOS_H
+
+/* Define to 1 if you have the <AztecOO.h> header file. */
+#undef HAVE_AZTECOO_H
+
+/* Define to 1 if you have the <AztecOO_Operator.h> header file. */
+#undef HAVE_AZTECOO_OPERATOR_H
+
 /* Define if the compiler provides __builtin_expect */
 #undef HAVE_BUILTIN_EXPECT
 
 /* Define to 1 if you have the `dgetrf_' function. */
 #undef HAVE_DGETRF_
 
+/* Define to 1 if you have the `dgetri_' function. */
+#undef HAVE_DGETRI_
+
 /* Define to 1 if you have the `dgetrs_' function. */
 #undef HAVE_DGETRS_
 
 /* Define to 1 if you have the `dstev_' function. */
 #undef HAVE_DSTEV_
 
+/* Define to 1 if you have the <Epetra_CrsGraph.h> header file. */
+#undef HAVE_EPETRA_CRSGRAPH_H
+
+/* Define to 1 if you have the <Epetra_CrsMatrix.h> header file. */
+#undef HAVE_EPETRA_CRSMATRIX_H
+
+/* Define to 1 if you have the <Epetra_Import.h> header file. */
+#undef HAVE_EPETRA_IMPORT_H
+
+/* Define to 1 if you have the <Epetra_LinearProblem.h> header file. */
+#undef HAVE_EPETRA_LINEARPROBLEM_H
+
+/* Define to 1 if you have the <Epetra_Map.h> header file. */
+#undef HAVE_EPETRA_MAP_H
+
+/* Define to 1 if you have the <Epetra_MultiVector.h> header file. */
+#undef HAVE_EPETRA_MULTIVECTOR_H
+
+/* Define to 1 if you have the <Epetra_Operator.h> header file. */
+#undef HAVE_EPETRA_OPERATOR_H
+
+/* Define to 1 if you have the <Epetra_SerialComm.h> header file. */
+#undef HAVE_EPETRA_SERIALCOMM_H
+
+/* Define to 1 if you have the <Epetra_Vector.h> header file. */
+#undef HAVE_EPETRA_VECTOR_H
+
 /* Define to 1 if you have the `gethostname' function. */
 #undef HAVE_GETHOSTNAME
 
 /* Availability of the MA47 algorithm from HSL */
 #undef HAVE_HSL_MA47
 
+/* Define to 1 if you have the <Ifpack.h> header file. */
+#undef HAVE_IFPACK_H
+
 /* Define to 1 if you have the <inttypes.h> header file. */
 #undef HAVE_INTTYPES_H
 
-/* Define to 1 if you have the `amd' library (-lamd). */
-#undef HAVE_LIBAMD
-
 /* Defined if deal.II was configured with BLAS support */
 #undef HAVE_LIBBLAS
 
    libstdc++ interface. */
 #undef HAVE_LIBSTDCXX_DEMANGLER
 
-/* UMFPACK is $withval */
+/* UMFPACK is */
 #undef HAVE_LIBUMFPACK
 
 /* Define to 1 if you have the `z' library (-lz). */
 /* Define to 1 if you have the <memory.h> header file. */
 #undef HAVE_MEMORY_H
 
+/* Define to 1 if you have the <ml_MultiLevelPreconditioner.h> header file. */
+#undef HAVE_ML_MULTILEVELPRECONDITIONER_H
+
 /* Define if you have the rand_r function */
 #undef HAVE_RAND_R
 
+/* Define to 1 if you have the <Sacado.hpp> header file. */
+#undef HAVE_SACADO_HPP
+
 /* Define to 1 if you have the `sgeevx_' function. */
 #undef HAVE_SGEEVX_
 
 /* Define to 1 if you have the `sgetrf_' function. */
 #undef HAVE_SGETRF_
 
+/* Define to 1 if you have the `sgetri_' function. */
+#undef HAVE_SGETRI_
+
 /* Define to 1 if you have the `sgetrs_' function. */
 #undef HAVE_SGETRS_
 
 /* Define to 1 if you have the <sys/types.h> header file. */
 #undef HAVE_SYS_TYPES_H
 
+/* Define to 1 if you have the <Teuchos_ParameterList.hpp> header file. */
+#undef HAVE_TEUCHOS_PARAMETERLIST_HPP
+
+/* Define to 1 if you have the <Teuchos_RCP.hpp> header file. */
+#undef HAVE_TEUCHOS_RCP_HPP
+
+/* Define to 1 if you have the <Teuchos_RefCountPtr.hpp> header file. */
+#undef HAVE_TEUCHOS_REFCOUNTPTR_HPP
+
+/* Define to 1 if you have the <Thyra_AztecOOLinearOpWithSolveFactory.hpp>
+   header file. */
+#undef HAVE_THYRA_AZTECOOLINEAROPWITHSOLVEFACTORY_HPP
+
+/* Define to 1 if you have the <Thyra_DefaultBlockedLinearOpDecl.hpp> header
+   file. */
+#undef HAVE_THYRA_DEFAULTBLOCKEDLINEAROPDECL_HPP
+
+/* Define to 1 if you have the <Thyra_DefaultBlockedLinearOp.hpp> header file.
+   */
+#undef HAVE_THYRA_DEFAULTBLOCKEDLINEAROP_HPP
+
+/* Define to 1 if you have the <Thyra_DefaultInverseLinearOp.hpp> header file.
+   */
+#undef HAVE_THYRA_DEFAULTINVERSELINEAROP_HPP
+
+/* Define to 1 if you have the <Thyra_EpetraLinearOp.hpp> header file. */
+#undef HAVE_THYRA_EPETRALINEAROP_HPP
+
+/* Define to 1 if you have the <Thyra_EpetraThyraWrappers.hpp> header file. */
+#undef HAVE_THYRA_EPETRATHYRAWRAPPERS_HPP
+
+/* Define to 1 if you have the <Thyra_InverseLinearOperator.hpp> header file.
+   */
+#undef HAVE_THYRA_INVERSELINEAROPERATOR_HPP
+
+/* Define to 1 if you have the <Thyra_LinearOperatorDecl.hpp> header file. */
+#undef HAVE_THYRA_LINEAROPERATORDECL_HPP
+
+/* Define to 1 if you have the <Thyra_LinearOperatorImpl.hpp> header file. */
+#undef HAVE_THYRA_LINEAROPERATORIMPL_HPP
+
+/* Define to 1 if you have the <Thyra_LinearOpWithSolveFactoryHelpers.hpp>
+   header file. */
+#undef HAVE_THYRA_LINEAROPWITHSOLVEFACTORYHELPERS_HPP
+
+/* Define to 1 if you have the <Thyra_MultiVectorBase.hpp> header file. */
+#undef HAVE_THYRA_MULTIVECTORBASE_HPP
+
+/* Define to 1 if you have the <Thyra_MultiVectorDefaultBase.hpp> header file.
+   */
+#undef HAVE_THYRA_MULTIVECTORDEFAULTBASE_HPP
+
+/* Define to 1 if you have the <Thyra_VectorDecl.hpp> header file. */
+#undef HAVE_THYRA_VECTORDECL_HPP
+
+/* Define to 1 if you have the <Thyra_VectorImpl.hpp> header file. */
+#undef HAVE_THYRA_VECTORIMPL_HPP
+
+/* Define to 1 if you have the <Thyra_VectorSpaceImpl.hpp> header file. */
+#undef HAVE_THYRA_VECTORSPACEIMPL_HPP
+
 /* Define to 1 if you have the <unistd.h> header file. */
 #undef HAVE_UNISTD_H
 
index 7800f15b1c6129c34d22ff5d4072ed8960e68fac..c795262b03d5b4ca1da828b0402380e0f1b8d280 100644 (file)
@@ -15,6 +15,9 @@ void dgemm_ (const char* transa, const char* transb,
 // Compute LU factorization
 void dgetrf_ (const int* m, const int* n, double* A,
              const int* lda, int* ipiv, int* info);
+// Invert matrix from LU factorization
+void dgetri_ (const int* n, double* A, const int* lda, 
+             int* ipiv, double* inv_work, const int* lwork, int* info);
 // Apply forward/backward substitution to LU factorization
 void dgetrs_ (const char* trans, const int* n, const int* nrhs,
              const double* A, const int* lda, const int* ipiv,
index 855076e9a3a13ff9f9618549d8ee3781099de9bc..2ddbb78eec92e071794a78379aa6f94829e96ebb 100644 (file)
@@ -2,7 +2,7 @@
 #    $Id$
 #    Version: $Name$
 #
-#    Copyright (C) 2005, 2006, 2007 by the deal authors
+#    Copyright (C) 2005, 2006, 2007, 2008 by the deal authors
 #
 #    This file is subject to QPL and may not be  distributed
 #    without copyright and license information. Please refer
@@ -23,7 +23,7 @@ print << 'EOT'
 //    This file was automatically generated from blas.h.in
 //    See blastemplates in the deal.II contrib directory
 //
-//    Copyright (C) 2005, 2006, 2007 by the deal authors
+//    Copyright (C) 2005, 2006, 2007, 2008 by the deal authors
 //
 //    This file is subject to QPL and may not be  distributed
 //    without copyright and license information. Please refer
index ffbebca87ca238d6e044cd58135cbaf2a16878b2..7a27ce59dc849ab4f8a5908016e5029a68d749a0 100644 (file)
@@ -1156,7 +1156,7 @@ class FullMatrix : public Table<2,number>
 
   private:
 
-#ifdef HAVE_LIBLAPACK
+#if defined(HAVE_DGETRF_) && defined (HAVE_SGETRF_) && defined(HAVE_DGETRI_) && defined (HAVE_SGETRI_)
                                      /**
                                      * The vector storing the
                                      * permutations applied for
index 89fcebf77fc62fdbc89dc1d32344dfafa3c6201a..fafd383f88c612f9aed4e5e8c2061d4424b0f125 100644 (file)
@@ -218,7 +218,7 @@ getrf (const int*, const int*, float*, const int*, int*, int*)
 #endif
 
 
-#ifdef HAVE_DGETRF_
+#ifdef HAVE_DGETRI_
 inline void
 getri (const int* n, double* A, const int* lda, int* ipiv, double* inv_work, const int* lwork, int* info)
 {
@@ -233,7 +233,7 @@ getri (const int*, double*, const int*, int*, double*, const int*, int*)
 #endif
 
 
-#ifdef HAVE_SGETRF_
+#ifdef HAVE_SGETRI_
 inline void
 getri (const int* n, float* A, const int* lda, int* ipiv, float* inv_work, const int* lwork, int* info)
 {
index 30ec55d91d62eef679bb7a339cece89683fb88a7..bf7d74b8b4630d9e724bad16ef1c7e7e9f689ad0 100644 (file)
@@ -23,7 +23,7 @@ DEAL_II_NAMESPACE_OPEN
                                    // inversion since it only works with
                                    // floats and doubles in case LAPACK was
                                    // detected by configure.
-#if defined(HAVE_DGETRF_) && defined (HAVE_SGETRF_) 
+#if defined(HAVE_DGETRF_) && defined (HAVE_SGETRF_) && defined(HAVE_DGETRI_) && defined (HAVE_SGETRI_)
 
 template <>
 void

In the beginning the Universe was created. This has made a lot of people very angry and has been widely regarded as a bad move.

Douglas Adams


Typeset in Trocchi and Trocchi Bold Sans Serif.