diff liboctave/CmplxHESS.cc @ 7482:29980c6b8604

don't check f77_exception_encountered
author John W. Eaton <jwe@octave.org>
date Thu, 14 Feb 2008 21:57:50 -0500 (2008-02-15)
parents a1dbe9d80eee
children eb63fbe60fab
line wrap: on
line diff
--- a/liboctave/CmplxHESS.cc
+++ b/liboctave/CmplxHESS.cc
@@ -88,56 +88,34 @@
 			     n, h, n, ilo, ihi, pscale, info
 			     F77_CHAR_ARG_LEN (1)));
 
-  if (f77_exception_encountered)
-    (*current_liboctave_error_handler) ("unrecoverable error in zgebal");
-  else
-    {
-      Array<Complex> tau (n-1);
-      Complex *ptau = tau.fortran_vec ();
+  Array<Complex> tau (n-1);
+  Complex *ptau = tau.fortran_vec ();
 
-      Array<Complex> work (lwork);
-      Complex *pwork = work.fortran_vec ();
+  Array<Complex> work (lwork);
+  Complex *pwork = work.fortran_vec ();
 
-      F77_XFCN (zgehrd, ZGEHRD, (n, ilo, ihi, h, n, ptau, pwork, lwork, info));
+  F77_XFCN (zgehrd, ZGEHRD, (n, ilo, ihi, h, n, ptau, pwork, lwork, info));
 
-      if (f77_exception_encountered)
-	(*current_liboctave_error_handler) ("unrecoverable error in zgehrd");
-      else
-	{
-	  unitary_hess_mat = hess_mat;
-	  Complex *z = unitary_hess_mat.fortran_vec ();
+  unitary_hess_mat = hess_mat;
+  Complex *z = unitary_hess_mat.fortran_vec ();
 
-	  F77_XFCN (zunghr, ZUNGHR, (n, ilo, ihi, z, n, ptau, pwork,
-				     lwork, info));
+  F77_XFCN (zunghr, ZUNGHR, (n, ilo, ihi, z, n, ptau, pwork,
+			     lwork, info));
 
-	  if (f77_exception_encountered)
-	    (*current_liboctave_error_handler)
-	      ("unrecoverable error in zunghr");
-	  else
-	    {
-	      F77_XFCN (zgebak, ZGEBAK, (F77_CONST_CHAR_ARG2 (&job, 1),
-					 F77_CONST_CHAR_ARG2 (&side, 1),
-					 n, ilo, ihi, pscale, n, z, n, info
-					 F77_CHAR_ARG_LEN (1)
-					 F77_CHAR_ARG_LEN (1)));
+  F77_XFCN (zgebak, ZGEBAK, (F77_CONST_CHAR_ARG2 (&job, 1),
+			     F77_CONST_CHAR_ARG2 (&side, 1),
+			     n, ilo, ihi, pscale, n, z, n, info
+			     F77_CHAR_ARG_LEN (1)
+			     F77_CHAR_ARG_LEN (1)));
 
-	      if (f77_exception_encountered)
-		(*current_liboctave_error_handler)
-		  ("unrecoverable error in zgebak");
-	      else
-		{
-		  // If someone thinks of a more graceful way of
-		  // doing this (or faster for that matter :-)),
-		  // please let me know!
+  // If someone thinks of a more graceful way of
+  // doing this (or faster for that matter :-)),
+  // please let me know!
 
-		  if (n > 2)
-		    for (octave_idx_type j = 0; j < a_nc; j++)
-		      for (octave_idx_type i = j+2; i < a_nr; i++)
-			hess_mat.elem (i, j) = 0;
-		}
-	    }
-	}
-    }
+  if (n > 2)
+    for (octave_idx_type j = 0; j < a_nc; j++)
+      for (octave_idx_type i = j+2; i < a_nr; i++)
+	hess_mat.elem (i, j) = 0;
 
   return info;
 }