aboutsummaryrefslogtreecommitdiff
diff options
context:
space:
mode:
authorjthorn <jthorn@f88db872-0e4f-0410-b76b-b9085cfa78c5>2003-01-22 20:22:54 +0000
committerjthorn <jthorn@f88db872-0e4f-0410-b76b-b9085cfa78c5>2003-01-22 20:22:54 +0000
commit2e821f2b71325aa1b56baedbf9c3a30c4af413b9 (patch)
tree243c71a80c7cd2a46367e61ae788d3f87e930fd4
parenta7bef848ed4af0592373acbcf26616a4ef71d3a3 (diff)
only call the single precision lapack/ilucg routines if #ifdef FP_IS_FLOAT,
only call the double precision lapack/ilucg routines if #ifdef FP_IS_DOUBLE this way we can link ok even on systems that have (eg) only double precision LAPACK, so long as that's allwe're using git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@936 f88db872-0e4f-0410-b76b-b9085cfa78c5
-rw-r--r--src/elliptic/ilucg_wrapper.F778
-rw-r--r--src/elliptic/lapack_wrapper.F774
2 files changed, 10 insertions, 2 deletions
diff --git a/src/elliptic/ilucg_wrapper.F77 b/src/elliptic/ilucg_wrapper.F77
index d1ab803..03aa5ca 100644
--- a/src/elliptic/ilucg_wrapper.F77
+++ b/src/elliptic/ilucg_wrapper.F77
@@ -16,6 +16,7 @@ c
#include "config.h"
#ifdef HAVE_ROW_SPARSE_JACOBIAN__ILUCG
+ #ifdef FP_IS_FLOAT
subroutine silucg_wrapper(n,
$ ia, ja, a,
$ b, x,
@@ -41,9 +42,11 @@ c
$ istatus)) ierror = 1
return
end
-#endif // HAVE_ROW_SPARSE_JACOBIAN__ILUCG
+ #endif /* FP_IS_FLOAT */
+#endif /* HAVE_ROW_SPARSE_JACOBIAN__ILUCG */
#ifdef HAVE_ROW_SPARSE_JACOBIAN__ILUCG
+ #ifdef FP_IS_DOUBLE
subroutine dilucg_wrapper(n,
$ ia, ja, a,
$ b, x,
@@ -69,4 +72,5 @@ c
$ istatus)) ierror = 1
return
end
-#endif // HAVE_ROW_SPARSE_JACOBIAN__ILUCG
+ #endif /* FP_IS_DOUBLE */
+#endif /* HAVE_ROW_SPARSE_JACOBIAN__ILUCG */
diff --git a/src/elliptic/lapack_wrapper.F77 b/src/elliptic/lapack_wrapper.F77
index 25796b3..c9f2fa2 100644
--- a/src/elliptic/lapack_wrapper.F77
+++ b/src/elliptic/lapack_wrapper.F77
@@ -15,6 +15,7 @@ c
#include "config.h"
#ifdef HAVE_DENSE_JACOBIAN__LAPACK
+ #ifdef FP_IS_FLOAT
subroutine sgecon_wrapper(norm_int,
$ N, A, LDA, anorm, rcond,
$ WORK, IWORK, info)
@@ -34,9 +35,11 @@ c
end if
return
end
+ #endif /* FP_IS_FLOAT */
#endif /* HAVE_DENSE_JACOBIAN__LAPACK */
#ifdef HAVE_DENSE_JACOBIAN__LAPACK
+ #ifdef FP_IS_DOUBLE
subroutine dgecon_wrapper(norm_int,
$ N, A, LDA, anorm, rcond,
$ WORK, IWORK, info)
@@ -56,4 +59,5 @@ c
end if
return
end
+ #endif /* FP_IS_DOUBLE */
#endif /* HAVE_DENSE_JACOBIAN__LAPACK */