diff liboctave/dbleHESS.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
parents a1dbe9d80eee
children eb63fbe60fab
line wrap: on
line diff
--- a/liboctave/dbleHESS.cc
+++ b/liboctave/dbleHESS.cc
@@ -87,58 +87,36 @@
 			     n, h, n, ilo, ihi, pscale, info
 			     F77_CHAR_ARG_LEN (1)));
 
-  if (f77_exception_encountered)
-    (*current_liboctave_error_handler) ("unrecoverable error in dgebal");
-  else
-    {
-      Array<double> tau (n-1);
-      double *ptau = tau.fortran_vec ();
+  Array<double> tau (n-1);
+  double *ptau = tau.fortran_vec ();
 
-      Array<double> work (lwork);
-      double *pwork = work.fortran_vec ();
+  Array<double> work (lwork);
+  double *pwork = work.fortran_vec ();
 
-      F77_XFCN (dgehrd, DGEHRD, (n, ilo, ihi, h, n, ptau, pwork,
-				 lwork, info));
+  F77_XFCN (dgehrd, DGEHRD, (n, ilo, ihi, h, n, ptau, pwork,
+			     lwork, info));
 
-      if (f77_exception_encountered)
-	(*current_liboctave_error_handler) ("unrecoverable error in dgehrd");
-      else
-	{
-	  unitary_hess_mat = hess_mat;
-	  double *z = unitary_hess_mat.fortran_vec ();
+  unitary_hess_mat = hess_mat;
+  double *z = unitary_hess_mat.fortran_vec ();
 
-	  F77_XFCN (dorghr, DORGHR, (n, ilo, ihi, z, n, ptau, pwork,
-				     lwork, info));
+  F77_XFCN (dorghr, DORGHR, (n, ilo, ihi, z, n, ptau, pwork,
+			     lwork, info));
 
-	  if (f77_exception_encountered)
-	    (*current_liboctave_error_handler)
-	      ("unrecoverable error in dorghr");
-	  else
-	    {
-	      F77_XFCN (dgebak, DGEBAK, (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 (dgebak, DGEBAK, (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 dgebak");
-	      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;
 }