29 #ifndef __lebesgueConstant_H
30 #define __lebesgueConstant_H 1
37 namespace ublas = boost::numeric::ublas;
39 template<
class Convex,
41 template<
class, u
int16_type,
class>
class PointSetType,
43 value_type lebesgueConstant()
45 static const uint16_type Dim = Convex::nDim;
47 typedef typename mpl::if_< mpl::bool_< Convex::is_simplex >, Simplex<Dim, 1>, Hypercube<Dim, 1> >::type convex_type;
49 typedef Lagrange<Order,Scalar,PointSetType>::template apply<Dim, value_type, Convex>::type basis_type;
52 PointSetEquiSpaced<convex_type, 100 - 80*( Dim-2 ), value_type> test;
54 return ublas::norm_1( M_basis.evaluate( test.points() ) );