Merge pull request #184 from imh/lapack-mem

Disabled the LAPACK functions that take too much memory

Author: denniszollo

clapack-3.2.1-CMAKE/CMakeLists.txt

@@ -3,6 +3,9 @@ set(CMAKE_C_FLAGS_COVERAGE "${CMAKE_C_FLAGS_DEBUG}")
 set(CMAKE_SHARED_LINKER_FLAGS_COVERAGE "${CMAKE_SHARED_LINKER_FLAGS_DEBUG}")
 set(CMAKE_EXE_LINKER_FLAGS_COVERAGE "${CMAKE_EXE_LINKER_FLAGS_DEBUG}")
 
+# Disable some of the really awful functions that allocate so much memory
+set(CMAKE_C_FLAGS "-DLAPACK_DISABLE_MEMORY_HOGS ${CMAKE_C_FLAGS}")
+
 # Build libraries with PIC flag as they will be linked into shared objects
 if (NOT CMAKE_CROSSCOMPILING)
   set(CMAKE_POSITION_INDEPENDENT_CODE ON)

clapack-3.2.1-CMAKE/SRC/cgbtrf.c

@@ -12,6 +12,7 @@
 
 #include "f2c.h"
 #include "blaswrap.h"
+#include "assert.h"
 
 /* Table of constant values */
 
@@ -41,8 +42,17 @@ static integer c__65 = 65;
 	    complex *, integer *), ccopy_(integer *, complex *, integer *, 
 	    complex *, integer *), cswap_(integer *, complex *, integer *, 
 	    complex *, integer *);
-    complex work13[4160]	/* was [65][64] */, work31[4160]	/* 
-	    was [65][64] */;
+	  #ifdef LAPACK_DISABLE_MEMORY_HOGS
+		  complex work13[1]	/* was [65][64] */, work31[1]	/* 
+		    was [65][64] */;
+	    /** This function uses too much memory, so we stopped allocating the memory
+	     * above and assert false here. */
+	    assert(0 && "cgbtrf_ was called. This function allocates too much"
+	                " memory and has been disabled.");
+	  #else
+	    complex work13[4160]	/* was [65][64] */, work31[4160]	/* 
+		    was [65][64] */;
+	  #endif
     extern /* Subroutine */ int ctrsm_(char *, char *, char *, char *, 
 	    integer *, integer *, complex *, complex *, integer *, complex *, 
 	    integer *), cgbtf2_(integer *, 
@@ -151,6 +161,7 @@ static integer c__65 = 65;
 /*     .. */
 /*     .. Executable Statements .. */
 
+
 /*     KV is the number of superdiagonals in the factor U, allowing for */
 /*     fill-in */
 

clapack-3.2.1-CMAKE/SRC/cgehrd.c

@@ -12,6 +12,7 @@
 
 #include "f2c.h"
 #include "blaswrap.h"
+#include "assert.h"
 
 /* Table of constant values */
 
@@ -32,7 +33,15 @@ static integer c__65 = 65;
 
     /* Local variables */
     integer i__, j;
-    complex t[4160]	/* was [65][64] */;
+    #ifdef LAPACK_DISABLE_MEMORY_HOGS
+      complex t[1]  /* was [65][64] */;
+      /** This function uses too much memory, so we stopped allocating the memory
+       * above and assert false here. */
+      assert(0 && "cgehrd_ was called. This function allocates too much"
+                  " memory and has been disabled.");
+    #else
+      complex t[4160]	/* was [65][64] */;
+    #endif
     integer ib;
     complex ei;
     integer nb, nh, nx, iws;

clapack-3.2.1-CMAKE/SRC/chseqr.c

@@ -12,6 +12,7 @@
 
 #include "f2c.h"
 #include "blaswrap.h"
+#include "assert.h"
 
 /* Table of constant values */
 
@@ -37,7 +38,15 @@ static integer c__49 = 49;
     /* Subroutine */ int s_cat(char *, char **, integer *, integer *, ftnlen);
 
     /* Local variables */
-    complex hl[2401]	/* was [49][49] */;
+    #ifdef LAPACK_DISABLE_MEMORY_HOGS
+      complex hl[1] /* was [49][49] */;
+      /** This function uses too much memory, so we stopped allocating the memory
+       * above and assert false here. */
+      assert(0 && "chseqr_ was called. This function allocates too much"
+                  " memory and has been disabled.");
+    #else
+      complex hl[2401]	/* was [49][49] */;
+    #endif
     integer kbot, nmin;
     extern logical lsame_(char *, char *);
     extern /* Subroutine */ int ccopy_(integer *, complex *, integer *, 

clapack-3.2.1-CMAKE/SRC/cpbtrf.c

@@ -12,6 +12,7 @@
 
 #include "f2c.h"
 #include "blaswrap.h"
+#include "assert.h"
 
 /* Table of constant values */
 
@@ -31,7 +32,15 @@ static integer c__33 = 33;
 
     /* Local variables */
     integer i__, j, i2, i3, ib, nb, ii, jj;
-    complex work[1056]	/* was [33][32] */;
+    #ifdef LAPACK_DISABLE_MEMORY_HOGS
+	    complex work[1]	/* was [33][32] */;
+	    /** This function uses too much memory, so we stopped allocating the memory
+	     * above and assert false here. */
+	    assert(0 && "cpbtrf_ was called. This function allocates too much"
+	                " memory and has been disabled.");
+    #else
+      complex work[1056]	/* was [33][32] */;
+    #endif
     extern /* Subroutine */ int cgemm_(char *, char *, integer *, integer *, 
 	    integer *, complex *, complex *, integer *, complex *, integer *, 
 	    complex *, complex *, integer *), cherk_(char *, 

clapack-3.2.1-CMAKE/SRC/cunmlq.c

@@ -12,6 +12,7 @@
 
 #include "f2c.h"
 #include "blaswrap.h"
+#include "assert.h"
 
 /* Table of constant values */
 
@@ -35,7 +36,15 @@ static integer c__65 = 65;
 
     /* Local variables */
     integer i__;
-    complex t[4160]	/* was [65][64] */;
+    #ifdef LAPACK_DISABLE_MEMORY_HOGS
+      complex t[1]  /* was [65][64] */;
+      /** This function uses too much memory, so we stopped allocating the memory
+       * above and assert false here. */
+      assert(0 && "cunmlq_ was called. This function allocates too much"
+                  " memory and has been disabled.");
+    #else
+      complex t[4160]	/* was [65][64] */;
+    #endif
     integer i1, i2, i3, ib, ic, jc, nb, mi, ni, nq, nw, iws;
     logical left;
     extern logical lsame_(char *, char *);

clapack-3.2.1-CMAKE/SRC/cunmql.c

@@ -12,6 +12,7 @@
 
 #include "f2c.h"
 #include "blaswrap.h"
+#include "assert.h"
 
 /* Table of constant values */
 
@@ -35,7 +36,15 @@ static integer c__65 = 65;
 
     /* Local variables */
     integer i__;
-    complex t[4160]	/* was [65][64] */;
+    #ifdef LAPACK_DISABLE_MEMORY_HOGS
+      complex t[1]  /* was [65][64] */;
+      /** This function uses too much memory, so we stopped allocating the memory
+       * above and assert false here. */
+      assert(0 && "cunmql_ was called. This function allocates too much"
+                  " memory and has been disabled.");
+    #else
+      complex t[4160]	/* was [65][64] */;
+    #endif
     integer i1, i2, i3, ib, nb, mi, ni, nq, nw, iws;
     logical left;
     extern logical lsame_(char *, char *);

clapack-3.2.1-CMAKE/SRC/cunmqr.c

@@ -12,6 +12,7 @@
 
 #include "f2c.h"
 #include "blaswrap.h"
+#include "assert.h"
 
 /* Table of constant values */
 
@@ -35,7 +36,15 @@ static integer c__65 = 65;
 
     /* Local variables */
     integer i__;
-    complex t[4160]	/* was [65][64] */;
+    #ifdef LAPACK_DISABLE_MEMORY_HOGS
+      complex t[1]  /* was [65][64] */;
+      /** This function uses too much memory, so we stopped allocating the memory
+       * above and assert false here. */
+      assert(0 && "cunmqr_ was called. This function allocates too much"
+                  " memory and has been disabled.");
+    #else
+      complex t[4160]	/* was [65][64] */;
+    #endif
     integer i1, i2, i3, ib, ic, jc, nb, mi, ni, nq, nw, iws;
     logical left;
     extern logical lsame_(char *, char *);

clapack-3.2.1-CMAKE/SRC/cunmrq.c

@@ -12,6 +12,7 @@
 
 #include "f2c.h"
 #include "blaswrap.h"
+#include "assert.h"
 
 /* Table of constant values */
 
@@ -35,7 +36,15 @@ static integer c__65 = 65;
 
     /* Local variables */
     integer i__;
-    complex t[4160]	/* was [65][64] */;
+    #ifdef LAPACK_DISABLE_MEMORY_HOGS
+      complex t[1]  /* was [65][64] */;
+      /** This function uses too much memory, so we stopped allocating the memory
+       * above and assert false here. */
+      assert(0 && "cunmrq_ was called. This function allocates too much"
+                  " memory and has been disabled.");
+    #else
+      complex t[4160] /* was [65][64] */;
+    #endif
     integer i1, i2, i3, ib, nb, mi, ni, nq, nw, iws;
     logical left;
     extern logical lsame_(char *, char *);

clapack-3.2.1-CMAKE/SRC/cunmrz.c

@@ -12,6 +12,7 @@
 
 #include "f2c.h"
 #include "blaswrap.h"
+#include "assert.h"
 
 /* Table of constant values */
 
@@ -36,7 +37,15 @@ static integer c__65 = 65;
 
     /* Local variables */
     integer i__;
-    complex t[4160]	/* was [65][64] */;
+    #ifdef LAPACK_DISABLE_MEMORY_HOGS
+      complex t[1]  /* was [65][64] */;
+      /** This function uses too much memory, so we stopped allocating the memory
+       * above and assert false here. */
+      assert(0 && "cunmrz_ was called. This function allocates too much"
+                  " memory and has been disabled.");
+    #else
+      complex t[4160]	/* was [65][64] */;
+    #endif
     integer i1, i2, i3, ib, ic, ja, jc, nb, mi, ni, nq, nw, iws;
     logical left;
     extern logical lsame_(char *, char *);