// -------------------------------------------------------------------- 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; }