Hello Arun,
unfortunately there is no function, which can handle boundary conditions
for Nédélec's elements, at the moment. (I am writing one, but it will
still take a few days until it is ready to release.) Anyway, if you only
want to impose homogeneous Dirichlet boundary conditions, you can simply
set the degrees of freedom located at the boundary to zero using the
constraint matrix.
Best Regards,
Markus
Am 07.06.10 22:18, schrieb arun jaganathan:
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
_______________________________________________
dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
