Hi Alex, I think we are going to start with hypre, using OMP. This is also new and so it will take some working out, probably. So take your time and when you get it into master, tell us and we can take it from there. Thanks, Mark On Thu, Sep 14, 2017 at 4:45 PM, Kalinkin, Alexander A < [email protected]> wrote:

HI Mark,

Looks like we can implement this functionality (if everything goes correctly) in MKL 2018u2 which will be published during following January

On our side we will test Both MKL triple implementation and wrappers from PETSc to MKL much early so if any petsc user will show interest to try this functionality we can commit PETSc wrapper early and provide engineering build with correspond function.

Thanks,

Alex

*From:* Mark Adams [mailto:[email protected]] *Sent:* Thursday, September 14, 2017 1:18 PM *To:* Kalinkin, Alexander A <[email protected]> *Cc:* Richard Tran Mills <[email protected]>; For users of the development version of PETSc <[email protected]>; Sokolova, Irina < [email protected]> *Subject:* Re: [petsc-dev] hypre + OMP and MKL-AIJ

On Thu, Sep 14, 2017 at 3:25 PM, Kalinkin, Alexander A < [email protected]> wrote:

Hi All,

Of course we can if you need it but I see the problem with testing – we need to implement functionality and test it simultaneously with MatMatMult() wrappers in petsc branch. Any thoughts how it can be realized?

I've warned my users that this new.

If you can get snes/examples/tutorials/ex56 to work then you are probably ready for users. My users will scale up and can work with you to debug.

Thanks,

Thanks,

Alex

*From:* Richard Tran Mills [mailto:[email protected]] *Sent:* Thursday, September 14, 2017 11:45 AM *To:* Mark Adams <[email protected]> *Cc:* For users of the development version of PETSc <[email protected]>; Kalinkin, Alexander A <[email protected]>; Sokolova, Irina < [email protected]> *Subject:* Re: [petsc-dev] hypre + OMP and MKL-AIJ

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?

--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