Commit 47e4ed3e by Francois Gygi

fixed use of info in inverse_det

git-svn-id: http://qboxcode.org/svn/qb/branches/efield@1590 cba15fb0-1239-40c8-b417-11db7ca47a34
parent 44f8b339
......@@ -1960,7 +1960,6 @@ void DoubleMatrix::inverse(void)
////////////////////////////////////////////////////////////////////////////////
double DoubleMatrix::inverse_det(void)
{
int info;
if ( active() )
{
assert(m_==n_);
......@@ -1984,17 +1983,7 @@ double DoubleMatrix::inverse_det(void)
det *= diag[ii];
inverse_from_lu(ipiv);
if(info!=0)
{
cout << " DoubleMatrix::inverse_det, info(getri)=" << info << endl;
#ifdef USE_MPI
MPI_Abort(MPI_COMM_WORLD, 2);
#else
exit(2);
#endif
}
return det;
return det;
}
}
......@@ -2063,7 +2052,6 @@ void ComplexMatrix::inverse(void)
////////////////////////////////////////////////////////////////////////////////
complex<double> ComplexMatrix::inverse_det(void)
{
int info;
if ( active() )
{
assert(m_==n_);
......@@ -2087,17 +2075,7 @@ complex<double> ComplexMatrix::inverse_det(void)
det *= diag[ii];
inverse_from_lu(ipiv);
if(info!=0)
{
cout << " ComplexMatrix::inverse_det, info(getri)=" << info << endl;
#ifdef USE_MPI
MPI_Abort(MPI_COMM_WORLD, 2);
#else
exit(2);
#endif
}
return det;
return det;
}
}
......
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment