Dear Yuyun, Here is the simple example that I wrote to test the function: https://petsc.org/release/src/mat/tests/ex248.c.html You can stick to MatCreate(), but this will only work if the type of your Mat is indeed MatSeqAIJ. If you need this for other types, let us know.
Thanks, Pierre > On 5 Dec 2022, at 9:12 AM, Yuyun Yang <[email protected]> wrote: > > Dear PETSc team, > > Is there an example for using MatSeqAIJKron? I’m using MatCreate for all > matrices in the code, and wonder if I need to switch to MatCreateSeqAIJ in > order to use this function? Just want to compute simple Kronecker products of > a sparse matrix with an identity matrix. > > Thank you, > Yuyun
