C++で2次元あるいは3次元ベクトル使う時って、boost::numeric::ublas::c_vectorを使いますよね。

・・・使いませんか？

私はuBLASとかViennaCLとか使いたいので、使いたいんですよ！
というか、C++用の線形代数ライブラリっていくつかありますけど、boostないC++なんてこの世に存在してはいけないので、逆に言えばuBLASとかC++ではなにもしなくても使えるもののはずなんです！

と、話が逸れましたが本題。

このboost::numeric::ublas::c_vectorって、普通のboost::numeric::ublas::vectorとは違って、配列(array)ベースのベクトルで、サイズを固定で作れるのでN次元のサイズ固定ベクトルのために重宝してるんです。
ただ、コンストラクタが貧弱なんですよね。

引数なし（値不定）：c_vector()
サイズ指定（Nがあるのになにに使うの？）：c_vector(size_type size)
コピーコンストラクタ：c_vector(const c_vector &v)
vector_expression指定：c_vector(const vector_expression< AE > &ae)

まぁ仕方ないんですけど、でも、そうなると初期化時に

normalConstruct.cpp

boost::numeric::ublas::c_vector<double, 2> vec;
vec(0) = 0;
vec(1) = 0;

とかしなきゃいけなくなります。ということは、そう、constにできないんですよね。

constConstruct_NG.cpp

const boost::numeric::ublas::c_vector<double, 2> vec;
vec(0) = 0; // ←できない
vec(1) = 0; // ←できない

list_ofとか使えればなんとかなるかなぁと試行錯誤したんですが、ダメでした。

一応、stackoverflowにも質問が上がってましたけど、結局constにできなさそうでした。

ので、仕方ないので、こんな力技で解決しました。

constVector_OK.cpp


// 次元数
const int DIM = DIMENSTION;

// N次元固定ベクトル
typedef boost::numeric::ublas::c_vector<double, DIM> Vector;

// 定数ベクトル
class VectorConst : public Vector
{
public:
	VectorConst(const double& value)
	{
		std::fill_n(this->begin(), DIM, value);
	}
};

// ゼロベクトル
const Vector VectorZero = VectorConst(0.0);

もっと良い解決方法ないものか・・・ご存知のかた教えて下さい。