// // Copyright (c) 2002--2010 // Toon Knapen, Karl Meerbergen, Kresimir Fresl, // Thomas Klimpel and Rutger ter Borg // // Distributed under the Boost Software License, Version 1.0. // (See accompanying file LICENSE_1_0.txt or copy at // http://www.boost.org/LICENSE_1_0.txt) // // THIS FILE IS AUTOMATICALLY GENERATED // PLEASE DO NOT EDIT! // #ifndef BOOST_NUMERIC_BINDINGS_LAPACK_COMPUTATIONAL_HPTRF_HPP #define BOOST_NUMERIC_BINDINGS_LAPACK_COMPUTATIONAL_HPTRF_HPP #include #include #include #include #include #include #include #include #include #include #include // // The LAPACK-backend for hptrf is the netlib-compatible backend. // #include #include namespace boost { namespace numeric { namespace bindings { namespace lapack { // // The detail namespace contains value-type-overloaded functions that // dispatch to the appropriate back-end LAPACK-routine. // namespace detail { // // Overloaded function for dispatching to // * netlib-compatible LAPACK backend (the default), and // * float value-type. // template< typename UpLo > inline std::ptrdiff_t hptrf( const UpLo, const fortran_int_t n, float* ap, fortran_int_t* ipiv ) { fortran_int_t info(0); LAPACK_SSPTRF( &lapack_option< UpLo >::value, &n, ap, ipiv, &info ); return info; } // // Overloaded function for dispatching to // * netlib-compatible LAPACK backend (the default), and // * double value-type. // template< typename UpLo > inline std::ptrdiff_t hptrf( const UpLo, const fortran_int_t n, double* ap, fortran_int_t* ipiv ) { fortran_int_t info(0); LAPACK_DSPTRF( &lapack_option< UpLo >::value, &n, ap, ipiv, &info ); return info; } // // Overloaded function for dispatching to // * netlib-compatible LAPACK backend (the default), and // * complex value-type. // template< typename UpLo > inline std::ptrdiff_t hptrf( const UpLo, const fortran_int_t n, std::complex* ap, fortran_int_t* ipiv ) { fortran_int_t info(0); LAPACK_CHPTRF( &lapack_option< UpLo >::value, &n, ap, ipiv, &info ); return info; } // // Overloaded function for dispatching to // * netlib-compatible LAPACK backend (the default), and // * complex value-type. // template< typename UpLo > inline std::ptrdiff_t hptrf( const UpLo, const fortran_int_t n, std::complex* ap, fortran_int_t* ipiv ) { fortran_int_t info(0); LAPACK_ZHPTRF( &lapack_option< UpLo >::value, &n, ap, ipiv, &info ); return info; } } // namespace detail // // Value-type based template class. Use this class if you need a type // for dispatching to hptrf. // template< typename Value > struct hptrf_impl { typedef Value value_type; typedef typename remove_imaginary< Value >::type real_type; // // Static member function, that // * Deduces the required arguments for dispatching to LAPACK, and // * Asserts that most arguments make sense. // template< typename MatrixAP, typename VectorIPIV > static std::ptrdiff_t invoke( MatrixAP& ap, VectorIPIV& ipiv ) { namespace bindings = ::boost::numeric::bindings; typedef typename result_of::uplo_tag< MatrixAP >::type uplo; BOOST_STATIC_ASSERT( (bindings::is_mutable< MatrixAP >::value) ); BOOST_STATIC_ASSERT( (bindings::is_mutable< VectorIPIV >::value) ); BOOST_ASSERT( bindings::size_column(ap) >= 0 ); return detail::hptrf( uplo(), bindings::size_column(ap), bindings::begin_value(ap), bindings::begin_value(ipiv) ); } }; // // Functions for direct use. These functions are overloaded for temporaries, // so that wrapped types can still be passed and used for write-access. In // addition, if applicable, they are overloaded for user-defined workspaces. // Calls to these functions are passed to the hptrf_impl classes. In the // documentation, most overloads are collapsed to avoid a large number of // prototypes which are very similar. // // // Overloaded function for hptrf. Its overload differs for // template< typename MatrixAP, typename VectorIPIV > inline std::ptrdiff_t hptrf( MatrixAP& ap, VectorIPIV& ipiv ) { return hptrf_impl< typename bindings::value_type< MatrixAP >::type >::invoke( ap, ipiv ); } } // namespace lapack } // namespace bindings } // namespace numeric } // namespace boost #endif