On Thu, Sep 14, 2017 at 2:45 PM, Richard Tran Mills <[email protected]> wrote:

Hi Mark,

Regarding MatMatMult() support, I think it should be fairly trivial for me to add that to AIJMKL. However, MKL doesn't currently have a matrix triple product routine (though the possibility of adding this has come up before in discussions with Alex Kalinkin from the MKL sparse team), so for now we'd need to implement that using MatMatMult().

Perhaps Alex or Irina (cc'd) can say something about possible support for sparse matrix triple products in MKL?

I don't mind if you wrap two Mat-Mat into the P'AP. Both of my apps are linear so this is amortized. Thanks,

--Richard

On Thu, Sep 14, 2017 at 11:32 AM, Mark Adams <[email protected]> wrote:

...

I recall Barry saying that he updated the hypre interface after the last hypre release, which includes OpenMP. But, I am not finding the email. Can someone tell me the status of this?

Note, I have two users that are interested in using threads with AMG. I think we would be interested in testing the new MKL-AIJ matrix but as I understand it we are waiting on the Mat-Mat Mult method(s) to get migrated. I should have at least one user that would be interested in testing this when it is ready.

Thanks, Mark