From ba1f4f889c909aeb755ca828c95da95eae9cca14 Mon Sep 17 00:00:00 2001
From: =?UTF-8?q?B=C3=A5rd=20Skaflestad?= <Bard.Skaflestad@sintef.no>
Date: Mon, 9 Aug 2010 11:05:03 +0000
Subject: [PATCH] Disable UMFPACK solver (and hard-coded source terms).

---
 mex_schur_comp_symm.c | 8 ++++++--
 1 file changed, 6 insertions(+), 2 deletions(-)

diff --git a/mex_schur_comp_symm.c b/mex_schur_comp_symm.c
index 5a73807a..ac405bf6 100644
--- a/mex_schur_comp_symm.c
+++ b/mex_schur_comp_symm.c
@@ -204,8 +204,10 @@ mexFunction(int nlhs,       mxArray *plhs[],
 
         for (i = 0; i < nc; i++)      { src[i]   = 0.0; } /* No sources */
         for (i = 0; i < ncf_tot; i++) { gflux[i] = 0.0; } /* No gravity */
+#if 0
         src[0] = 1;
         src[nc-1]=-1;
+#endif
         Binv = mxGetPr(prhs[0]);
         get_nconn(prhs[1], nconn);
 
@@ -230,13 +232,15 @@ mexFunction(int nlhs,       mxArray *plhs[],
         /* double *b = malloc(nf * sizeof *b); */
         double *b;
         hybsys_assemble(nc, nf, nconn, conn, ptr, sys->r, &A, &b);
+#if 0
         double *x = malloc(A.n * sizeof *x);
         
         callMWUMFPACK(A.n, A.ia, A.ja, A.sa, b, x);
 
         int i;for(i=0; i<nf; ++i)mexPrintf("x[%d]=%f\n", i, x[i]);
-        free(A.ia); free(A.ja); free(A.sa); free(x); free(b);
-        
+#endif
+        free(A.ia); free(A.ja); free(A.sa); /* free(x);  */free(b);
+
         ptr = mxGetPr(plhs[1]);
         memcpy(ptr, sys->r, ncf_tot * sizeof *ptr);