Re: [petsc-users] Kronecker Product

2023-01-30 Thread Jose E. Roman
oKit. You have to supply a matrix into the Slepc MFN > operator to set the problem up as far as I know. > > Tyler > > From: Matthew Knepley > Date: Monday, January 30, 2023 at 12:24 PM > To: Guglielmo, Tyler Hardy > Cc: Barry Smith , petsc-users@mcs.anl.gov > >

Re: [petsc-users] Kronecker Product

2023-01-30 Thread Guglielmo, Tyler Hardy via petsc-users
Cc: Barry Smith , petsc-users@mcs.anl.gov Subject: Re: [petsc-users] Kronecker Product On Mon, Jan 30, 2023 at 3:08 PM Guglielmo, Tyler Hardy via petsc-users mailto:petsc-users@mcs.anl.gov>> wrote: I would need the Kronecker product to be explicitly available to perform matrix expone

Re: [petsc-users] Kronecker Product

2023-01-30 Thread Matthew Knepley
st, > > Tyler > > > > *From: *Barry Smith > *Date: *Monday, January 30, 2023 at 12:01 PM > *To: *Guglielmo, Tyler Hardy > *Cc: *petsc-users@mcs.anl.gov > *Subject: *Re: [petsc-users] Kronecker Product > > > > What is large? If A and B have dimensions of

Re: [petsc-users] Kronecker Product

2023-01-30 Thread Guglielmo, Tyler Hardy via petsc-users
To: Guglielmo, Tyler Hardy Cc: petsc-users@mcs.anl.gov Subject: Re: [petsc-users] Kronecker Product What is large? If A and B have dimensions of 1000, then the Kronecker product is of size 1,000,000. Do you want the Kronecker product to be explicitly formed or just available as matrix vector products

Re: [petsc-users] Kronecker Product

2023-01-30 Thread Barry Smith
gt; > Cc: Barry Smith mailto:bsm...@petsc.dev>>, > petsc-users@mcs.anl.gov <mailto:petsc-users@mcs.anl.gov> > mailto:petsc-users@mcs.anl.gov>> > Subject: Re: [petsc-users] Kronecker Product > > On Mon, Jan 30, 2023 at 2:24 PM Guglielmo, Tyler Hardy via petsc-use

Re: [petsc-users] Kronecker Product

2023-01-30 Thread Guglielmo, Tyler Hardy via petsc-users
-users@mcs.anl.gov Subject: Re: [petsc-users] Kronecker Product On Mon, Jan 30, 2023 at 2:24 PM Guglielmo, Tyler Hardy via petsc-users mailto:petsc-users@mcs.anl.gov>> wrote: Thanks Barry, I saw that function, but wasn’t sure how to apply it since the documentation says that S and T are

Re: [petsc-users] Kronecker Product

2023-01-30 Thread Matthew Knepley
To: *Guglielmo, Tyler Hardy > *Cc: *petsc-users@mcs.anl.gov > *Subject: *Re: [petsc-users] Kronecker Product > > > >Do you need the explicit sparse representation of the Kronecker > product? Or do you want to apply it as an operator or solve systems with > it? If t

Re: [petsc-users] Kronecker Product

2023-01-30 Thread Guglielmo, Tyler Hardy via petsc-users
at 11:12 AM To: Guglielmo, Tyler Hardy Cc: petsc-users@mcs.anl.gov Subject: Re: [petsc-users] Kronecker Product Do you need the explicit sparse representation of the Kronecker product? Or do you want to apply it as an operator or solve systems with it? If the latter you can use https

Re: [petsc-users] Kronecker Product

2023-01-30 Thread Barry Smith
Do you need the explicit sparse representation of the Kronecker product? Or do you want to apply it as an operator or solve systems with it? If the latter you can use https://petsc.org/release/docs/manualpages/Mat/MatCreateKAIJ/#matcreatekaij Barry > On Jan 30, 2023, at 12:53 PM,

[petsc-users] Kronecker Product

2023-01-30 Thread Guglielmo, Tyler Hardy via petsc-users
Hi all, I am wondering if there is any functionality for taking Kronecker products of large sparse matrices that are parallel? MatSeqAIJKron is as close as I have found, but it seems like this does not work for parallel matrices. Any ideas here? An option could be to make A and B

[petsc-users] Kronecker product

2019-02-25 Thread Yuyun Yang via petsc-users
Hello team, I'd like to ask whether PETSc has a function to compute the Kronecker product of a sparse matrix with an identity matrix? A Google search didn't lead me to a manual page (like most of the other PETSc functions), so I'm wondering if this has been implemented yet. Thanks very much!

Re: [petsc-users] Kronecker product

2019-02-25 Thread Jed Brown via petsc-users
MatCreateMAIJ does that (implicitly). https://www.mcs.anl.gov/petsc/petsc-current/docs/manualpages/Mat/MatCreateMAIJ.html If you want a Kronecker product with a non-identity matrix, this PR may be of interest. https://bitbucket.org/petsc/petsc/pull-requests/1334/rmills-mat-kaij/diff Yuyun Yang