message: fbcsr matrix generator is generalized for non-square matrices
diff: @@ -153,17 +153,16 @@ std::unique_ptr<matrix::Fbcsr<ValueType, IndexType>> generate_fbcsr_from_csr(
  template <typename ValueType, typename IndexType, typename RandEngine>
  -std::unique_ptr<matrix::Fbcsr<ValueType, IndexType>>
  -generate_random_square_fbcsr(std::shared_ptr<ReferenceExecutor> ref,
  - RandEngine engine, const IndexType nbrows,
  - const int mat_blk_sz, const bool diag_dominant,
  - const bool unsort)
  +std::unique_ptr<matrix::Fbcsr<ValueType, IndexType>> generate_random_fbcsr(
  + std::shared_ptr<const ReferenceExecutor> ref, RandEngine engine,
  + const IndexType nbrows, const IndexType nbcols, const int mat_blk_sz,
  + const bool diag_dominant, const bool unsort)
  {
  using real_type = gko::remove_complex<ValueType>;
  std::unique_ptr<matrix::Csr<ValueType, IndexType>> rand_csr_ref =
  generate_random_matrix<matrix::Csr<ValueType, IndexType>>(
  - nbrows, nbrows,
  - std::uniform_int_distribution<IndexType>(0, nbrows - 1),
  + nbrows, nbcols,
  + std::uniform_int_distribution<IndexType>(0, nbcols - 1),
  std::normal_distribution<real_type>(0.0, 1.0), std::move(engine),
  ref);
  gko::kernels::reference::factorization::add_diagonal_elements(
  