[petsc-users] Tensor product as matrix free method

[Jed Brown]

Mathis Friesdorf <mathisfriesdorf at gmail.com> writes:

> Alright, so your point is that while this construction is in principal
> possible, it will not be efficient. Thanks for clearing this up!

Yes, TAIJ chooses a specific order for simplicity.  A general tensor
library is a big problem that has already gotten a lot of attention, so
we would prefer to use one.  Wrapping as a MatShell is very simple,
though it may be interesting to include a full interface in the library.
-------------- 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/a0a545d6/attachment.pgp>