Esistono implementazioni ILU multilivello inverse basate su più livelli?

Sono rimasto molto colpito dalle prestazioni seriali dei precondizionatori ILU multilivello basati su inverse , in particolare per l' eterogenea Helmholtz , ma sono sorpreso di non riuscire a trovare implementazioni open source. In particolare, ILUPACK rende i binari liberamente disponibili per gli accademici, ma non sembra che rilascino il loro codice sorgente.

È davvero vero che nessuno ha aperto la propria implementazione?

Hai esaminato il pacchetto di solutori iterativi di Yousef Saad ITSOL ? Include un'altra libreria chiamata Algebraic Recursive Multilevel Solver pARMS , di cui sono abbastanza sicuro di cosa stai parlando; Saad è l'autore dell'articolo che hai collegato.

C'è anche Euclid , a cui PETSc può interfacciarsi, ma non credo che Euclid funzioni in un modo multilivello come PARMS.