diff --git a/examples/dense/testFrontMPI.cpp b/examples/dense/testFrontMPI.cpp index 74ee0a1c..e30ea97e 100644 --- a/examples/dense/testFrontMPI.cpp +++ b/examples/dense/testFrontMPI.cpp @@ -13,8 +13,10 @@ MPI_Init(&argc, &argv); { strumpack::MPIComm c; + strumpack::BLACSGrid grid(c); strumpack::DenseMatrix Aseq; strumpack::DistributedMatrix A; + int m = 0; if (!strumpack::mpi_rank()) { std::string filename; @@ -26,10 +28,12 @@ MPI_Init(&argc, &argv); std::ifstream file(filename, std::ifstream::binary); file.read(reinterpret_cast(&m), sizeof(int)); Aseq = strumpack::DenseMatrix(m, m); - file.read(reinterpret_cast(Aseq.data()), sizeof(double)*m*m); + file.read(reinterpret_cast(Aseq.data()), sizeof(double)*m*m); + std::cout << "# Matrix dimension read from file: " << m << std::endl; } MPI_Bcast(&m, 1, strumpack::mpi_type(), 0, MPI_COMM_WORLD); + A = strumpack::DistributedMatrix(&grid, m, m); A.scatter(Aseq);