Hi Lawrence,
So you build UU, UP, PU, and PP separately?
Yes, but as you say, I should only build them once, not multiple times, which is incredibly wasteful.
Given they never change, I would assemble them up front and then just use them as and where. If you always do (whereever you need it):
UU = assemble(...)
Then you will pay the cost multiple times.
If you build everything up front then calling mat.M a second time doesn't reassemble (since now the promise to assemble is fulfilled) and the operator is ready.
Ok, good, so that sounds like what I’m going to do will work and I only pay the cost of assembling the mixed operator once. Thanks, Eike
Lawrence -----BEGIN PGP SIGNATURE----- Version: GnuPG v1
iQEcBAEBAgAGBQJVNitHAAoJECOc1kQ8PEYvrB8H/0SW4RYh7IEIsv2UjU9O6vyW DEVlLdiHMJTmbgJu9kJ7N9/6LabdPqBOVr15AB/IpRbm3GD0BahOevCaFbI3Ck+r 7hyhgZyBAsZXzggI+cjXKPOTJtvmXrvCzM4uoK0EHbK496gKwoa4o8xc6JrUKt+5 7J6594g6PH6nudIzd51Gm6YRpF0ldPU6HvJ5Hg9EUuSV1R5mtlPvhM8BU9cxpB4N Cq+BR0EnSvruL238nTvSHMyRitrqPe01eHFPby426xAzPFrcdy3WERLOO6GQXN1k IVuGl1cFV5L8ylsocxOlv/6odNEJX7+06dvPqLizu+haSoqbMtb8t3k2Uqq1LOw= =g9Oo -----END PGP SIGNATURE-----
_______________________________________________ firedrake mailing list firedrake@imperial.ac.uk https://mailman.ic.ac.uk/mailman/listinfo/firedrake