--- blas2_overloads.hpp.orig	2009-03-03 03:25:10.000000000 +0000
+++ blas2_overloads.hpp	2009-03-03 03:30:14.000000000 +0000
@@ -5,6 +5,7 @@
 // (See accompanying file LICENSE_1_0.txt or copy at
 // http://www.boost.org/LICENSE_1_0.txt)
 //
+// Modified to include dtrsv by Kian Ming A. Chai (3 Mar 2009)
 
 #ifndef BOOST_NUMERIC_BINDINGS_BLAS_BLAS2_OVERLOADS_HPP
 #define BOOST_NUMERIC_BINDINGS_BLAS_BLAS2_OVERLOADS_HPP
@@ -29,6 +30,13 @@ namespace boost { namespace numeric { na
   void ger( const int& m, const int& n, const float  & alpha, const float  * x_ptr, const int& incx, const float  * y_ptr, const int& incy, float  * a_ptr, const int& lda ) { BLAS_SGER( &m, &n, &alpha, x_ptr, &incx, y_ptr, &incy, a_ptr, &lda ) ; }
   inline
   void ger( const int& m, const int& n, const double & alpha, const double * x_ptr, const int& incx, const double * y_ptr, const int& incy, double * a_ptr, const int& lda ) { BLAS_DGER( &m, &n, &alpha, x_ptr, &incx, y_ptr, &incy, a_ptr, &lda ) ; }
+
+  // Added by Kian Ming A. Chai (3 Mar 2009)
+  inline
+  void trsv( char UPLO, char TRANS, char DIAG, const int& n, const double *a_ptr, const int& lda, double *x_ptr, const int& incx) {
+    BLAS_DTRSV( &UPLO, &TRANS, &DIAG, &n, a_ptr, &lda, x_ptr, &incx );
+  }
+
 /*  
   inline
   void geru( const int& m, const int& n, const complex_f & alpha, const complex_f * x_ptr, const int& incx, complex_f * y_ptr, const int& incy, complex_f * a_ptr, const int& lda ) { BLAS_CGERU( &m, &n, complex_ptr( &alpha ), complex_ptr( x_ptr ), &incx, complex_ptr( y_ptr ), &incy, complex_ptr( a_ptr ), &lda ) ; }