Beispiel #1
0
 inline
 int gesv (MatrA& a, MatrB& b) {
   // with 'internal' pivot vector 
   
   // gesv() errors: 
   //   if (info == 0), successful
   //   if (info < 0), the -info argument had an illegal value
   //   -- we will use -101 if allocation fails
   //   if (info > 0), U(i-1,i-1) is exactly zero 
   int info = -101; 
   traits::detail::array<int> ipiv (traits::matrix_size1 (a)); 
   if (ipiv.valid()) 
     info = gesv (a, ipiv, b); 
   return info; 
 }
 inline
 int lu_solve (MatrA& a, MatrB& b) {
   return gesv (a, b); 
 }