PETSc зависает на MatView()

Я пытаюсь создать и распечатать MATDENSE матрица в PETSc. Мой MWE выглядит так:

static char help[] = "";

#include <petscmat.h>
#include <petscviewer.h>

int main(int argc,char **args)
{
  Mat            A;
  PetscInt       Istart,Iend,m = 2,n = 2;
  PetscErrorCode ierr;
  PetscMPIInt    rank,size;
  PetscViewer    viewer;

  ierr = PetscInitialize(&argc,&args,(char*)0,help);if (ierr) return ierr;
  ierr = MPI_Comm_rank(PETSC_COMM_WORLD,&rank);CHKERRQ(ierr);
  ierr = MPI_Comm_size(PETSC_COMM_WORLD,&size);CHKERRQ(ierr);
  ierr = PetscOptionsGetInt(NULL,NULL,"-m",&m,NULL);CHKERRQ(ierr);
  ierr = PetscOptionsGetInt(NULL,NULL,"-n",&n,NULL);CHKERRQ(ierr);

  ierr = MatCreateDense(PETSC_COMM_WORLD, PETSC_DECIDE, PETSC_DECIDE, m, n, NULL, &A); CHKERRQ(ierr);
  ierr = MatSetUp(A);CHKERRQ(ierr);

  ierr = MatGetOwnershipRange(A,&Istart,&Iend);CHKERRQ(ierr);

  MatZeroEntries(A);

  printf("rank: %d, istart: %d, iend: %d\n", rank, Istart, Iend);

  ierr = MatAssemblyBegin(A,MAT_FINAL_ASSEMBLY);CHKERRQ(ierr);
  ierr = MatAssemblyEnd(A,MAT_FINAL_ASSEMBLY);CHKERRQ(ierr);

  ierr = MatSetOption(A,MAT_SYMMETRIC,PETSC_TRUE);CHKERRQ(ierr);

  ierr = PetscViewerCreate(PETSC_COMM_WORLD,&viewer);CHKERRQ(ierr);
  ierr = PetscViewerASCIIOpen(PETSC_COMM_WORLD, "out.txt", &viewer);
  ierr = MatView(A, viewer); CHKERRQ(ierr);

  ierr = MatDestroy(&A);CHKERRQ(ierr);
  ierr = PetscViewerDestroy(&viewer);CHKERRQ(ierr);
  ierr = PetscFinalize();
  return ierr;
}

Проблема в том, что при запуске более чем одного процесса программа внезапно зависает на MatView() вместо записи в файл out.txt. Что я здесь делаю неправильно?


c petsc
person Eenoku    schedule 23.10.2019    source источник

Ответы (1)


arrow_upward
0
arrow_downward

Это было вызвано, как ни странно, какой-то ошибкой в ​​MPI-пакетах Lubuntu 18.10. Как только я обновился до Lubuntu 19.04, проблема исчезла.

СОВЕТ. В таком случае всегда пытайтесь запускать make check. Набор тестов будет построен и запущен автоматически, и вы сможете наблюдать за поведением PETSc без длительных ручных проб и ошибок.

person Eenoku    schedule 02.12.2019