[petsc-users] Tensor product as matrix free method

[Jed Brown]

Mathis Friesdorf <mathisfriesdorf at gmail.com> writes:

> Thanks Jed, this is great. I was trying to implement what was proposed in
> http://lists.mcs.anl.gov/pipermail/petsc-users/2011-September/009991.html,
> by you as I just realised. For this the matrices would be of varying size,
> so for three local systems, I would need:
>
> S \otimes Id \otimes Id ; Id \otimes S \otimes Id ; Id \otimes Id \otimes S.
>
> Realistic systems are of the order of at least 20 local systems, which
> would mean that I need the matrices of the above structure for all 20
> possible positions.

This is something for which you should use CTF or TCE.  When S moves out
of the last position in the tensor product, the operation becomes less
local and various forms of transposition is needed.  This is not
something worth recreating in PETSc unless there is a critical
deficiency in the available libraries.
-------------- next part --------------
A non-text attachment was scrubbed...
Name: not available
Type: application/pgp-signature
Size: 835 bytes
Desc: not available
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20140401/b7073733/attachment.pgp>