#include <boost/mpl/assert.hpp>
#include <boost/numeric/ublas/fwd.hpp>
#include <boost/numeric/ublas/matrix.hpp>
#include <boost/numeric/ublas/traits.hpp>
#include <boost/numeric/ublas/traits/layout_type.hpp> // <-- NEW -|
#include <boost/type_traits/is_same.hpp>
#include <complex>

namespace ublas = boost::numeric::ublas;

int main()
{
	typedef double value_type;

	typedef ublas::matrix<value_type, ublas::column_major> matrix_type1;

	typedef ublas::matrix<std::complex<double>, ublas::layout_type<matrix_type1>::type> matrix_type2;
//	typedef ublas::matrix<std::complex<double>, ublas::row_major> matrix_type2;
//	typedef ublas::matrix<std::complex<double>, ublas::column_major> matrix_type2;

	// To check at compile time, uncomment this
	BOOST_MPL_ASSERT((
		   boost::is_same<
			  ublas::matrix_traits<matrix_type1>::orientation_category,
			  ublas::matrix_traits<matrix_type2>::orientation_category
			>
	));

	// To check at runtime time, uncomment this
	bool flag = boost::is_same<
				  ublas::matrix_traits<matrix_type1>::orientation_category,
				  ublas::matrix_traits<matrix_type2>::orientation_category
				>::value;
	if (!flag)
	{
		std::cerr << "Different layout storage." << std::endl;
	}
}