[petsc-dev] MATOP_MAT_MULT

[Pierre Jolivet]

Hong,
I also removed a couple of deprecated MATOP there https://gitlab.com/petsc/petsc/-/merge_requests/2786 <https://gitlab.com/petsc/petsc/-/merge_requests/2786>
If you git grep MATOP_TRANSPOSE_MAT_MULT, you’ll see that there are still some remaining in src/tao/constrained/impls/admm/admm.c
Actually, in this case, it does not come from a MatHasOperation but from a MatShellSetOperation.
I don’t think there is any example showing how to add an algorithm inside user-code to compute C=A*B. Maybe it would be worth creating such an example? (I’m actually not sure how to do it myself)
That way, it would be quite convenient to try out new algorithms, instead of having to recompile the library itself.

Thanks,
Pierre

> On 9 May 2020, at 1:36 PM, Stefano Zampini <stefano.zampini at gmail.com> wrote:
> 
> Hong
> 
> you may want to take a look at this big commit here https://gitlab.com/petsc/petsc/-/merge_requests/2712/diffs?commit_id=94a81e695f30a0ff4a14c0eb596e09393f9b329f <https://gitlab.com/petsc/petsc/-/merge_requests/2712/diffs?commit_id=94a81e695f30a0ff4a14c0eb596e09393f9b329f>
> 
> Now that you coded MatProduct, I think it is the case of not playing with the destroy/view routines anymore in the various MatMat operations (you know, I have raised this issue already in the past). I think I got it fixed properly. What I did in short, I have added extra slots to the MatProduct struct to keep track of how to destroy/view the extra information you need for the products.
> Also, I have deprecated MatFreeIntermediateStructures, since we now have MatProductClear
> 
> I have spent quite a bit of time on it, and I personally think it is the proper way to go. I think I took care of all the cases already covered. Tests are clean for me
> I would like to hear your opinion.
> 
> 
> Thanks
> Stefano
> 
> Il giorno sab 9 mag 2020 alle ore 13:45 Stefano Zampini <stefano.zampini at gmail.com <mailto:stefano.zampini at gmail.com>> ha scritto:
> Jose
> 
> I have just pushed a test https://gitlab.com/petsc/petsc/-/blob/d64c2bc63c8d5d1a8c689f1abc762ae2722bba26/src/mat/tests/ex69.c <https://gitlab.com/petsc/petsc/-/blob/d64c2bc63c8d5d1a8c689f1abc762ae2722bba26/src/mat/tests/ex69.c>
> See if it fits your framework, and feel free to modify the test to add more checks
> 
> 
> Il giorno ven 8 mag 2020 alle ore 18:48 Jose E. Roman <jroman at dsic.upv.es <mailto:jroman at dsic.upv.es>> ha scritto:
> Attached. Run with -test 1 or -test 2
> 
> > El 8 may 2020, a las 17:14, Stefano Zampini <stefano.zampini at gmail.com <mailto:stefano.zampini at gmail.com>> escribió:
> > 
> > Jose
> > 
> > Just send me a MWE and I’ll fix the case for you
> > 
> > Thanks
> > Stefano
> 
> 
> -- 
> Stefano
> 
> 
> -- 
> Stefano

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