@rouson Thank you very much for your insight. You are right, presently, we are looking at PETSc for an efficient linear algebra solver: I will try to promote PSBLAS into our "group". Just for your knowledge, we are trying to design a library similar (hopefully better) to OpenFOAM in Fortran standard, OOD-TDD on the good experience obtained by FOODIE (P.S. I hope to finish FOODIE's test on Fermi this week). At some point, when the preliminary design of this new library will be working I would like to invite you to review/join our project. My best regards.