Hi,
I used tutorial no. 36 to solve for eigen value problem in electromagnetics
using Nedelec elements. I would like to impose Dirichlet boundary but
'DoFTools::make_zero_boundary_constraints' do not seem to work for vector
elements. Is there any other way ? Any hints would be helpful.
Thanks !
Arun.
Here is the spinets from the code:
FE_Nedelec<dim> fe;
.
..
for (unsigned int q_point=0; q_point<n_q_points; ++q_point)
for (unsigned int i=0; i<dofs_per_cell; ++i)
for (unsigned int j=0; j<dofs_per_cell; ++j)
{
cell_stiffness_matrix (i, j)
+= (fe_values[E].curl (i, q_point) *
fe_values[E].curl (j, q_point)) * fe_values.JxW
(q_point);
cell_mass_matrix (i, j)
+= (fe_values[E].value (i, q_point) *
fe_values[E].value (j, q_point)) * fe_values.JxW
(q_point);
}
cell->get_dof_indices (local_dof_indices);
constraints
.distribute_local_to_global (cell_stiffness_matrix,
local_dof_indices,
stiffness_matrix);
constraints
.distribute_local_to_global (cell_mass_matrix,
local_dof_indices,
mass_matrix);
}
_______________________________________________
dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
