Thanks Markus !
On Mon, Jun 7, 2010 at 4:52 PM, Markus Bürg <[email protected]> wrote:
> 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
>
>
_______________________________________________
dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
