// --------------------------------------------------------------------
static matrix_type _create_b_vec(
const matrix_type & current_values,
const size_t row_padding)
{
assert(current_values.is_vector());
const auto K = current_values.get_length();
matrix_type b_vec (K + K + row_padding, 1);
for (size_t k = 0; k < K; k++)
{
b_vec[k] = current_values[k];
b_vec[k + K] = value_type(1) - current_values[k];
}
for (size_t k = K + K; k < b_vec.get_height(); k++)
b_vec[k] = value_type(0);
return b_vec;
}
