On 6/22/23 04:06, 'Jost Arndt' via deal.II User Group wrote:
N_{i,j} F, with N_{i,j,k} = (\phi_i \phi_k, \phi_j)

?

This would be a third order tensor. It would be sparse, because only those entries ijk are nonzero where phi_i, phi_j, and phi_k overlap. But that would still be quite a lot of entries. A back of the envelope calculation would be that it's 10-100 times more memory than the mass matrix, depending on the space dimension and the polynomial degree you use. That's not likely going to fit into your computer's memory for all but the smallest problems :-)


Or do I simply have to use a Newton solver forĀ  (f g, \phi) for f and g?

It does not actually change anything. The N tensor above is just a convenient way to write (f g, phi), but the problem is still nonlinear. If you need to solve for both f and g, you're stuck with a nonlinear problem for which Newton's method is one way to do it.

Best
 W.


--
------------------------------------------------------------------------
Wolfgang Bangerth          email:                 bange...@colostate.edu
                           www: http://www.math.colostate.edu/~bangerth/


--
The deal.II project is located at http://www.dealii.org/
For mailing list/forum options, see 
https://groups.google.com/d/forum/dealii?hl=en
--- You received this message because you are subscribed to the Google Groups "deal.II User Group" group.
To unsubscribe from this group and stop receiving emails from it, send an email 
to dealii+unsubscr...@googlegroups.com.
To view this discussion on the web visit 
https://groups.google.com/d/msgid/dealii/5b372b05-8053-c713-06f3-8881116e4485%40colostate.edu.

Reply via email to