[petsc-dev] hypre + OMP and MKL-AIJ

[Kalinkin, Alexander A]

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:mfadams at lbl.gov]
Sent: Thursday, September 14, 2017 1:18 PM
To: Kalinkin, Alexander A <alexander.a.kalinkin at intel.com>
Cc: Richard Tran Mills <rtmills at anl.gov>; For users of the development version of PETSc <petsc-dev at mcs.anl.gov>; Sokolova, Irina <irina.sokolova at intel.com>
Subject: Re: [petsc-dev] hypre + OMP and MKL-AIJ



On Thu, Sep 14, 2017 at 3:25 PM, Kalinkin, Alexander A <alexander.a.kalinkin at intel.com<mailto:alexander.a.kalinkin at intel.com>> 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:rtmills at anl.gov<mailto:rtmills at anl.gov>]
Sent: Thursday, September 14, 2017 11:45 AM
To: Mark Adams <mfadams at lbl.gov<mailto:mfadams at lbl.gov>>
Cc: For users of the development version of PETSc <petsc-dev at mcs.anl.gov<mailto:petsc-dev at mcs.anl.gov>>; Kalinkin, Alexander A <alexander.a.kalinkin at intel.com<mailto:alexander.a.kalinkin at intel.com>>; Sokolova, Irina <irina.sokolova at intel.com<mailto:irina.sokolova at intel.com>>
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 <mfadams at lbl.gov<mailto:mfadams at lbl.gov>> 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


-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20170914/737bdd64/attachment.html>