résolution système avec MPI et ScalaPack

Bonjour
Je travaille sur la parallélisation d'un code Fortran avec MPI
à un moment le code résout un système tri-diagonale A*X=B

donc dans ma décomposition de domaine chaque processus aura accès à un bloc (une partie ) de la matrice global A

j'ai pensé à utiliser la librairie ScalaPack pour paralléliser cette résolution

ScalaPack utilise l'outil BLACS pour créer une grille de processus p*q avec la commande :CALL BLACS_GRIDINIT(icontxt,'c',p,q)

c'est là ou je comprends pas parce que j'ai déjà crée avant une topologie MPI (une grille de processus ) avec

CALL MPI_DIMS_CREATE et CALL MPI_CART_CREATE

et je pense qu'il faut utiliser cette dernière au lieu de créer une autre . comment on fait techniquement ?

Merci