Re: [starpu-devel] Mix two solvers on top of StarPU

[Philippe SWARTVAGHER <philippe.swartvagher@inria.fr>]

Hello Marek,

(switching to English for the archive of the mailing-list)

Le 25/10/2022 à 16:23, Marek Felšöci a écrit :
> Nous essayons actuellement de développer une application qui combine 
> deux solveurs de systèmes linéaires à base de StarPU. Il s'agit d'un 
> solveur creux en mémoire partagée et d'un solveur dense parallèle 
> distribué faisant appel à StarPU-MPI (ce qui n'est pas le cas du 
> solveur creux). Nous avons modifié les APIs de ces solveurs afin de 
> pouvoir passer des "handles" StarPU d'un solveur à l'autre. Néanmoins, 
> nous aurions besoin d'aide pour faire cohabiter l'instance de StarPU 
> du solveur creux avec celle du solveur dense. Lors de nos premiers 
> essaies, uniquement le solveur qui initialise StarPU en premier semble 
> fonctionner. Les routines du deuxième solveur restent bloquées.
>
> Quelle serait la bonne façon d'initialiser StarPU dans cette 
> configuration ? Est-ce possible de combiner l'utilisation d'un solveur 
> faisant appel à StarPU-MPI et celle d'un solveur qui ne le fait pas ?

I'm not sure I really got your problem... The problem is to make

starpu_init();

/* sparse solver stuff */

starpu_shutdown();

starpu_mpi_init();

/* dense solver stuff */

starpu_mpi_shutdown();

with the same handles?

If so, I guess you can directly call starpu_mpi_init() instead of 
starpu_init(), and then you don't need to restart StarPU specifically 
for the MPI support.

We can discuss more in details offline.

--
Philippe SWARTVAGHER

PhD Student
TADaaM team, Inria Bordeaux Sud-Ouest