uBLAS版に引き続き、ViennaCL版も書いておくと、こんな感じです。
（というか、これがしたくてuBLAS版頑張った）

viennacl_crs.cpp

// 行列の次元数
const int COUNT = 5;

// 非ゼロ要素数
const int NONZERO_COUNT = 15;

// 要素値
double elements[NONZERO_COUNT] ={
6, 2,    1, 4,
2, 7,       2,
      6, 3,   
1,    3, 9,   
4, 2,       8, 
};
// 列番号
unsigned int columnIndeces[NONZERO_COUNT] ={
0, 1,    3, 4,
0, 1,       4,
      2, 3,   
1,    2, 3,   
0, 1,       4, 
};
// 各行先頭の位置
unsigned int rowOffsets[COUNT+1] ={
0,
4,
7,
9,
12,
15 //＝非ゼロ要素数 ※ここ忘れがち
};


// 行列の生成
viennacl::compressed_matrix<double> mat(count, count);

// データの指定
mat.set(rowOffsets, columnIndeces, elements, COUNT, COUNT, NONZERO_COUNT);

あまり解説することもなく、単にsetを使うだけです
コンストラクタでuBLASのcompressed_matrixとかstd::vector>とかも指定できますが、生データがあるならこれが一番速いです。
これでuBLASとデータ共有できますしね！