Hi
I apologise if the question's too simplistic, but does this give the correct
behaviour:
MatrixCreator::create_mass_matrix (dof_handler, QGauss<dim>(3),
mass_matrix);
MatrixCreator::create_laplace_matrix (dof_handler, QGauss<dim>(3),
laplace_matrix);
system_matrix.copy_from (mass_matrix);
system_matrix.add (theta * time_step * tensor_coefficient,
laplace_matrix);
What I'm worried about is the tensor_coefficient, which is a 3x3 matrix. Is
it properly contracted with the laplace_matrix, i.e, is it equivalent to :
cell_matrix(i,j) += ( tensor_coefficient * fe_values.shape_grad
<http://www.dealii.org/6.3.1/doxygen/deal.II/classFEValuesBase.html#a0a172cb1d5a0e51364b0284288f727a0>
(i, q_point) *
fe_values.shape_grad
<http://www.dealii.org/6.3.1/doxygen/deal.II/classFEValuesBase.html#a0a172cb1d5a0e51364b0284288f727a0>
(j, q_point) *
fe_values.JxW
<http://www.dealii.org/6.3.1/doxygen/deal.II/classFEValuesBase.html#ad097580a2f71878695096cc73b271b9d>
(q_point) );
--
Best regards,
Ted
_______________________________________________
dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
