Как оптимизировать код boost :: mpl :: vector?

Добавьте момент, когда я играю с boost :: mpl :: * и попробуйте расширить пример физической единицы. На данный момент у меня есть следующий код:

Простое физическое базовое измерение, представленное mpl :: vector:

template < int Mass, int Length, int Time, int Temperature, int Angle, int Current >
struct base_dimension
{
typedef typename mpl::vector_c< int, Mass, Length, Time, Temperature, Angle, Current >::type type;
};

Скалярный вектор:

typedef base_dimension< 0, 0, 0, 0, 0, 0 >::type base_dimensionless_helper;

Простой тип, хранящий тип base_dimension:

template< class base_dimension >
struct Dimension
{
typedef Dimension< base_dimension > type;
typedef base_dimension              base_dim_type;
};

Скалярное измерение:

typedef Dimension< base_dimensionless_helper > base_dimensionless;

Какой-то простой помощник:

template < class D1, int fac >
struct mul_base_dim_fac_typeof
{
typedef typename base_dimension< fac, fac, fac, fac, fac, fac >::type fac_vec;
typedef typename Detail::multiply_typeof_helper< typename D1, fac_vec >::type type;
};template < class Dim1, class Dim2 >
struct add_dim_typeof_helper
{
typedef typename add_base_dim_typeof_helper< typename Dim1::base_dim_type, typename Dim2::base_dim_type >::type dim;

typedef Dimension< dim > type;
};

И, наконец, производное измерение. Используемая формула:
(D0 * E0) + (D1 * E1) + … + (D5 * E5) … :: type

template < class D0 = base_dimensionless, int E0 = 0,
class D1 = base_dimensionless, int E1 = 0,
class D2 = base_dimensionless, int E2 = 0,
class D3 = base_dimensionless, int E3 = 0,
class D4 = base_dimensionless, int E4 = 0,
class D5 = base_dimensionless, int E5 = 0 >
struct derived_dimension
{
typedef Dimension< typename mul_base_dim_fac_typeof< typename D0::base_dim_type, E0 >::type > d0_type;
typedef Dimension< typename mul_base_dim_fac_typeof< typename D1::base_dim_type, E1 >::type > d1_type;
typedef Dimension< typename mul_base_dim_fac_typeof< typename D2::base_dim_type, E2 >::type > d2_type;
typedef Dimension< typename mul_base_dim_fac_typeof< typename D3::base_dim_type, E3 >::type > d3_type;
typedef Dimension< typename mul_base_dim_fac_typeof< typename D4::base_dim_type, E4 >::type > d4_type;
typedef Dimension< typename mul_base_dim_fac_typeof< typename D5::base_dim_type, E5 >::type > d5_type;

typedef typename Detail::add_dim_typeof_helper< d0_type, d1_type >::type d0_d1_type;
typedef typename Detail::add_dim_typeof_helper< d0_d1_type, d2_type >::type d0_d1_d2_type;
typedef typename Detail::add_dim_typeof_helper< d0_d1_d2_type, d3_type >::type d0_d1_d2_d3_type;
typedef typename Detail::add_dim_typeof_helper< d0_d1_d2_d3_type, d4_type >::type d0_d1_d2_d3_d4_type;
typedef typename Detail::add_dim_typeof_helper< d0_d1_d2_d3_d4_type, d5_type >::type type;
};

Хорошо, это работает как ожидалось. Но я хочу украсить вычисление производного измерения, потому что все эти промежуточные определения типов ужасны. Моя первая идея состояла в том, чтобы поместить входные векторы в один вектор (-> вектор векторов) и умножить их в одном хорошем цикле for_each, но до сих пор безуспешно. Итак, мой вопрос:

Любые советы, как украсить расчет?