inline boost::shared_ptr< Base > Construct( std::vector< t_OriginalVectorType >& vectors, std::vector< double >& labels, unsigned int const trainingSize, std::vector< std::string > const& kernelParameters, unsigned int const cacheSize = 0 ) { uint64_t sparseSize = 0; if ( t_Traits::SPARSE_ALLOWED ) { typename std::vector< SparseVector< t_Type > >::const_iterator ii = vectors.begin(); typename std::vector< SparseVector< t_Type > >::const_iterator iiEnd = vectors.end(); for ( ; ii != iiEnd; ++ii ) sparseSize += ii->Size(); } uint64_t spanSize = 0; if ( t_Traits::SPAN_ALLOWED ) { typename std::vector< SparseVector< t_Type > >::const_iterator ii = vectors.begin(); typename std::vector< SparseVector< t_Type > >::const_iterator iiEnd = vectors.end(); for ( ; ii != iiEnd; ++ii ) spanSize += SpanVector< t_Type >( *ii ).Size(); } uint64_t denseSize = 0; { typename std::vector< SparseVector< t_Type > >::const_iterator ii = vectors.begin(); typename std::vector< SparseVector< t_Type > >::const_iterator iiEnd = vectors.end(); for ( ; ii != iiEnd; ++ii ) denseSize += DenseVector< t_Type >( *ii ).Size(); } boost::shared_ptr< Base > result; if ( t_Traits::SPARSE_ALLOWED && ( ( ! t_Traits::SPAN_ALLOWED ) || ( sparseSize <= spanSize ) ) && ( sparseSize <= denseSize ) ) { result = boost::shared_ptr< Base >( new t_KernelType< SparseVector< t_Type >, t_Traits >( vectors, labels, trainingSize, kernelParameters, cacheSize ) ); } else if ( t_Traits::SPAN_ALLOWED && ( spanSize <= denseSize ) ) { result = boost::shared_ptr< Base >( new t_KernelType< SpanVector< t_Type >, t_Traits >( vectors, labels, trainingSize, kernelParameters, cacheSize ) ); } else { result = boost::shared_ptr< Base >( new t_KernelType< DenseVector< t_Type >, t_Traits >( vectors, labels, trainingSize, kernelParameters, cacheSize ) ); } return result; }