Go to the documentation of this file.00001
00002
00003
00004
00005
00006
00007
00008
00009
00010
00011
00012
00013
00014
00015 #ifndef _VIENNACL_NORM_1_HPP_
00016 #define _VIENNACL_NORM_1_HPP_
00017
00022 #include <math.h>
00023 #include "viennacl/forwards.h"
00024 #include "tag_of.hpp"
00025
00026 namespace viennacl
00027 {
00028
00029
00030
00031
00032
00033 namespace linalg
00034 {
00035
00036 #ifdef VIENNACL_HAVE_UBLAS
00037
00038
00039
00040 template< typename VectorT >
00041 typename VectorT::value_type
00042 norm_1(VectorT const& vector,
00043 typename viennacl::tools::enable_if< viennacl::is_ublas< typename viennacl::traits::tag_of< VectorT >::type >::value
00044 >::type* dummy = 0)
00045 {
00046
00047 return boost::numeric::ublas::norm_1(vector);
00048 }
00049 #endif
00050
00051
00052
00053
00054
00055 template< typename VectorT>
00056 typename VectorT::value_type
00057 norm_1(VectorT const& v1,
00058 typename viennacl::tools::enable_if< viennacl::is_stl< typename viennacl::traits::tag_of< VectorT >::type >::value
00059 >::type* dummy = 0)
00060 {
00061
00062 typename VectorT::value_type result = 0;
00063 for (typename VectorT::size_type i=0; i<v1.size(); ++i)
00064 result += fabs(v1[i]);
00065
00066 return result;
00067 }
00068
00069
00070
00071
00072 template< typename ScalarType, unsigned int alignment >
00073 viennacl::scalar_expression< const viennacl::vector<ScalarType, alignment>,
00074 const viennacl::vector<ScalarType, alignment>,
00075 viennacl::op_norm_1 >
00076 norm_1(viennacl::vector<ScalarType, alignment> const & vector,
00077 typename viennacl::tools::enable_if< viennacl::is_viennacl< typename viennacl::traits::tag_of< viennacl::vector<ScalarType, alignment> >::type >::value
00078 >::type* dummy = 0)
00079 {
00080 return viennacl::scalar_expression< const viennacl::vector<ScalarType, alignment>,
00081 const viennacl::vector<ScalarType, alignment>,
00082 viennacl::op_norm_1 >(vector, vector);
00083 }
00084
00085 }
00086 }
00087 #endif
00088
00089
00090
00091
00092