Dear all,
As I see in most of the tutorials when there's a Laplacian for a variable
like U, what is done in assembling the system_matrix is like:
cell_matrix
<http://www.dealii.org/developer/doxygen/deal.II/namespaceLocalIntegrators_1_1Divergence.html#a1a6d2c938fd71efc441cef321ac3baca>(i,j)
+= (fe_values.shape_grad
<http://www.dealii.org/developer/doxygen/deal.II/classFEValuesBase.html#a0a172cb1d5a0e51364b0284288f727a0>
(i, q_point) *
fe_values.shape_grad
<http://www.dealii.org/developer/doxygen/deal.II/classFEValuesBase.html#a0a172cb1d5a0e51364b0284288f727a0>
(j, q_point) *
fe_values.JxW<http://www.dealii.org/developer/doxygen/deal.II/classFEValuesBase.html#ad097580a2f71878695096cc73b271b9d>(q_point));
Now I'm wondering how it defers if we have the Laplacian term multiplied by
the same variable we have, ex. U*Laplacian(U).
Thank you in advance.
Mohsen
_______________________________________________
dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
