454 generalize matrix over container functor

include/alp/reference/blas3.hpp

@@ -1130,7 +1130,7 @@ namespace alp {
 		class Operator
 	>
 	Matrix< typename Operator::D3, structures::General, Density::Dense,
-		view::Functor< std::function< void( InputType1 &, const size_t &, const size_t & ) > >,
+		view::Functor< std::function< void( InputType1 &, const size_t, const size_t ) > >,
 		imf::Id, imf::Id,
 		reference
 	>
@@ -1151,8 +1151,8 @@ namespace alp {
 			"called with a postfactor vector that does not match the first domain "
 			"of the given multiplication operator" );
 
-		std::function< void( typename Operator::D3 &, const size_t &, const size_t & ) > lambda =
-			[ &x, &y, &mul ]( typename Operator::D3 &result, const size_t &i, const size_t &j ) {
+		std::function< void( typename Operator::D3 &, const size_t, const size_t ) > lambda =
+			[ &x, &y, &mul ]( typename Operator::D3 &result, const size_t i, const size_t j ) {
 				//set( ret, alp::identities::zero );
 				internal::apply( result, x[ i ], y[ j ], mul );
 			};
@@ -1161,14 +1161,14 @@ namespace alp {
 			typename Operator::D3,
 			structures::General,
 			Density::Dense,
-			view::Functor< std::function< void( typename Operator::D3 &, const size_t &, const size_t & ) > >,
+			view::Functor< std::function< void( typename Operator::D3 &, const size_t, const size_t ) > >,
 			imf::Id, imf::Id,
 			reference
 			>(
 				internal::getInitialized( x ) && internal::getInitialized( y ), // Temporary solution, pending proper implemention
-				lambda,
-				imf::Id( getLength( x ) ),
-				imf::Id( getLength( y ) )
+				getLength( x ),
+				getLength( y ),
+				lambda
 			);
 
 	}
@@ -1183,7 +1183,7 @@ namespace alp {
 		class Operator
 	>
 	Matrix< typename Operator::D3, structures::Symmetric, Density::Dense,
-		view::Functor< std::function< void( typename Operator::D3 &, const size_t &, const size_t & ) > >,
+		view::Functor< std::function< void( typename Operator::D3 &, const size_t, const size_t ) > >,
 		imf::Id, imf::Id,
 		reference
 	>
@@ -1202,8 +1202,8 @@ namespace alp {
 			"called with a prefactor vector that does not match the first domain "
 			"of the given multiplication operator" );
 
-		std::function< void( typename Operator::D3 &, const size_t &, const size_t & ) > lambda =
-			[ &x, &mul ]( typename Operator::D3 &result, const size_t &i, const size_t &j ) {
+		std::function< void( typename Operator::D3 &, const size_t, const size_t ) > lambda =
+			[ &x, &mul ]( typename Operator::D3 &result, const size_t i, const size_t j ) {
 				//set( ret, alp::identities::zero );
 				internal::apply( result, x[ i ], x[ j ], mul );
 			};
@@ -1212,13 +1212,13 @@ namespace alp {
 			typename Operator::D3,
 			structures::Symmetric,
 			Density::Dense,
-			view::Functor< std::function< void( typename Operator::D3 &, const size_t &, const size_t & ) > >,
+			view::Functor< std::function< void( typename Operator::D3 &, const size_t, const size_t ) > >,
 			imf::Id, imf::Id,
 			reference
 			>(
 				internal::getInitialized( x ),
-				lambda,
-				imf::Id( getLength( x ) )
+				getLength( x ),
+				lambda
 			);
 
 	}

include/alp/reference/io.hpp

@@ -179,6 +179,32 @@ namespace alp {
 		return PANIC;
 	}
 
+	/**
+	 * @brief \a buildVector version.
+	 *
+	 */
+	template< typename InputType, typename Structure, typename View, typename Imf, typename fwd_iterator >
+	RC buildVector(
+		Vector< InputType, Structure, Density::Dense, View, Imf, reference > &v,
+		const fwd_iterator &start,
+		const fwd_iterator &end
+	) noexcept {
+
+		// Temporarily assuming 1-1 mapping with user container
+		internal::setInitialized(v, true);
+
+		InputType * praw, * p;
+
+		size_t len = internal::getLength( internal::getContainer( v ) );
+		praw = p = internal::getRaw( internal::getContainer( v ) );
+
+		for( fwd_iterator it = start; p < praw + len && it != end; ++it, ++p ) {
+			*p = *it;
+		}
+
+		return PANIC;
+	}
+
 } // end namespace ``alp''
 
 #endif // end ``_H_ALP_REFERENCE_IO''