Multithreading with new pastix
Concerns branch feature@newpastix.
-
With the new pastix, to enable multithreading, the number of threads in use have to be set before pastixInit. However, the current SDS factorization and sparse precond factorization call XMPH_SDS_set_multithread (setting multithreading) after XMPH_SDS_set_matrix (calling pastixInit). Some investigation are required to find the appropriate order which will preserve both pastix and mumps functionalities. -
About multithreading again, manual binding has to be investigated too. The subroutine pastixInitWithAffinity is going to be reachable from fortran very soon. -
Sparse preconditioning with PILUT: requires pastix perm, unavailable for this moment.