Segfault on MPI HSS generation with laplace kernel. about strumpack HOT 3 OPEN

#include <cmath>
#include <iostream>

#include "dense/DistributedMatrix.hpp"
#include "HSS/HSSMatrixMPI.hpp"
#include "kernel/Kernel.hpp"
using namespace strumpack;
using namespace strumpack::HSS;

int main(int argc, char * argv[]) {
  MPI_Init(&argc, &argv);
  int N = atoi(argv[1]);

  HSSOptions<double> hss_opts;
  hss_opts.set_verbose(false);
  hss_opts.set_from_command_line(argc, argv);
  {
    BLACSGrid grid(MPI_COMM_WORLD);
    // DistributedMatrix<double> A;

    DenseMatrix<double> vector(N, 1);  // data for a 1D laplace kernel
    for (int i = 0; i < N; ++i) {
      vector(i, 0) = i;
    }

    auto laplace1d = kernel::LaplaceKernel<double>(vector, N, N);

    HSSMatrixMPI<double> HSS(laplace1d, &grid, hss_opts);
  }
  scalapack::Cblacs_exit(1);
  MPI_Finalize();
}

# ANN search quality = 1 after 0 iterations
# ANN search quality = 1 after 0 iterationsterminate called after throwing an instance of '
std::bad_array_new_length'
  what():  std::bad_array_new_length
# ANN search quality = # ANN search quality = 1 after 1 after 0 iterations0 iterations

{    1,    0}:  On entry to 
DESCINIT parameter number    6 had an illegal value 
 ERROR: Could not create DistributedMatrix descriptor!
terminate called after throwing an instance of 'std::bad_array_new_length'
terminate called after throwing an instance of 'std::bad_array_new_length'
  what():  std::bad_array_new_length
  what():  std::bad_array_new_length

===================================================================================
=   BAD TERMINATION OF ONE OF YOUR APPLICATION PROCESSES
=   RANK 0 PID 131297 RUNNING AT g0082
=   KILLED BY SIGNAL: 6 (Aborted)
===================================================================================

===================================================================================
=   BAD TERMINATION OF ONE OF YOUR APPLICATION PROCESSES
=   RANK 1 PID 131298 RUNNING AT g0082
=   KILLED BY SIGNAL: 6 (Aborted)
===================================================================================

===================================================================================
=   BAD TERMINATION OF ONE OF YOUR APPLICATION PROCESSES
=   RANK 2 PID 131299 RUNNING AT g0082
=   KILLED BY SIGNAL: 6 (Aborted)
===================================================================================

===================================================================================
=   BAD TERMINATION OF ONE OF YOUR APPLICATION PROCESSES
=   RANK 3 PID 131300 RUNNING AT g0082
=   KILLED BY SIGNAL: 6 (Aborted)
===================================================================================

#include <cmath>
#include <iostream>
#include "mpi.h"

#include "dense/DistributedMatrix.hpp"
#include "HSS/HSSMatrixMPI.hpp"
#include "kernel/Kernel.hpp"
using namespace strumpack;
using namespace strumpack::HSS;

int main(int argc, char * argv[]) {
  MPI_Init(&argc, &argv);
  int N = atoi(argv[1]);
  int nprocs;
  int rank;

  MPI_Comm_size(MPI_COMM_WORLD, &nprocs);
  MPI_Comm_rank(MPI_COMM_WORLD, &rank);

  std::cout << "nprocs: " << nprocs << " rank: " << rank << std::endl;

  HSSOptions<double> hss_opts;
  hss_opts.set_verbose(false);
  hss_opts.describe_options();
  hss_opts.set_from_command_line(argc, argv);
  
  int nr, nc;
  {
    BLACSGrid grid(MPI_COMM_WORLD, nprocs);
    DenseMatrix<double> vector(N, 1);  // data for a 1D laplace kernel
    for (int i = 0; i < N; ++i) {
      vector(i, 0) = i;
    }

    auto laplace1d = kernel::LaplaceKernel<double>(vector, N, N);

    HSSMatrixMPI<double> HSS(laplace1d, &grid, hss_opts);
  }

  scalapack::Cblacs_exit(1);
  MPI_Finalize();
}