the Compaq f90 compiler complains about empty (or all-comment)

files, so we still have to define empty subroutines even if they're
never called.  :( Oh well, memory's cheap...


git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@1010 f88db872-0e4f-0410-b76b-b9085cfa78c5

2 files changed, 28 insertions, 16 deletions

diff --git a/src/elliptic/ilucg_wrapper.F77 b/src/elliptic/ilucg_wrapper.F77
index 8e76936..32ee1ec 100644
--- a/src/elliptic/ilucg_wrapper.F77
+++ b/src/elliptic/ilucg_wrapper.F77
@@ -13,10 +13,14 @@ c ierror = (out) This is returned as 0 for a normal return,
 c                nonzero for an error return.
 c
 
+c
+c Note that the Compaq f90 compiler complains about empty (or all-comment)
+c files, so we still have to define an empty subroutine even if it's never
+c called.  :( Oh well, memory's cheap...
+c
+
 #include "config.h"
 
-#ifdef HAVE_ROW_SPARSE_JACOBIAN__ILUCG
-#ifdef FP_IS_FLOAT
 	subroutine silucg_wrapper(n,
      $				  ia, ja, a,
      $				  b, x,
@@ -30,6 +34,8 @@ c
 	real rtemp(*)
 	real eps
 	integer iter, istatus, ierror
+#ifdef HAVE_ROW_SPARSE_JACOBIAN__ILUCG
+#ifdef FP_IS_FLOAT
 c
 	logical silucg
 c
@@ -40,13 +46,11 @@ c
      $		   itemp, rtemp,
      $		   eps, iter,
      $		   istatus)) ierror = 1
-	return
-	end
 #endif /* FP_IS_FLOAT */
 #endif	/* HAVE_ROW_SPARSE_JACOBIAN__ILUCG */
+	return
+	end
 
-#ifdef HAVE_ROW_SPARSE_JACOBIAN__ILUCG
-#ifdef FP_IS_DOUBLE
 	subroutine dilucg_wrapper(n,
      $				  ia, ja, a,
      $				  b, x,
@@ -60,6 +64,8 @@ c
 	double precision rtemp(*)
 	double precision eps
 	integer iter, istatus, ierror
+#ifdef HAVE_ROW_SPARSE_JACOBIAN__ILUCG
+#ifdef FP_IS_DOUBLE
 c
 	logical dilucg
 c
@@ -70,7 +76,7 @@ c
      $		   itemp, rtemp,
      $		   eps, iter,
      $		   istatus)) ierror = 1
-	return
-	end
 #endif /* FP_IS_DOUBLE */
 #endif	/* HAVE_ROW_SPARSE_JACOBIAN__ILUCG */
+	return
+	end
diff --git a/src/elliptic/lapack_wrapper.F77 b/src/elliptic/lapack_wrapper.F77
index 14f5cbf..95c3375 100644
--- a/src/elliptic/lapack_wrapper.F77
+++ b/src/elliptic/lapack_wrapper.F77
@@ -12,10 +12,14 @@ c norm_int = (in) 0 ==> infinity-norm
 c                 1 ==> 1-norm
 c
 
+c
+c Note that the Compaq f90 compiler complains about empty (or all-comment)
+c files, so we still have to define an empty subroutine even if it's never
+c called.  :( Oh well, memory's cheap...
+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)
@@ -26,6 +30,8 @@ c
 	real WORK(*)
 	integer iwork(*)
 	integer info
+#ifdef HAVE_DENSE_JACOBIAN__LAPACK
+#ifdef FP_IS_FLOAT
 	if	(norm_int .eq. 0) then
 		call sgecon('I', N,A,LDA, anorm,rcond, WORK,IWORK, info)
 	else if (norm_int .eq. 1) then
@@ -33,13 +39,11 @@ c
 	else
 		info = -1;
 	end if
-	return
-	end
 #endif /* FP_IS_FLOAT */
 #endif	/* HAVE_DENSE_JACOBIAN__LAPACK */
+	return
+	end
 
-#ifdef HAVE_DENSE_JACOBIAN__LAPACK
-#ifdef FP_IS_DOUBLE
 	subroutine dgecon_wrapper(norm_int,
      $	                          N, A, LDA, anorm, rcond,
      $	                          WORK, IWORK, info)
@@ -50,6 +54,8 @@ c
 	double precision WORK(*)
 	integer iwork(*)
 	integer info
+#ifdef HAVE_DENSE_JACOBIAN__LAPACK
+#ifdef FP_IS_DOUBLE
 	if	(norm_int .eq. 0) then
 		call dgecon('I', N,A,LDA, anorm,rcond, WORK,IWORK, info)
 	else if (norm_int .eq. 1) then
@@ -57,7 +63,7 @@ c
 	else
 		info = -1;
 	end if
-	return
-	end
 #endif /* FP_IS_DOUBLE */
 #endif	/* HAVE_DENSE_JACOBIAN__LAPACK */
+	return
+	end