Hi,
Sorry to keep asking this question but I tried changing the code using both
VectorTools::interpolate_boundary_values and
DOFTools::make_zero_boundary_constraints and I get the same runtime error. Can
anyone tell me what this error means?
An error occurred in line <1506> of file
</users/stirupat/Research/Software/deal.II/deal.II/include/numerics/vectors.templates.h>
in function
static void dealii::VectorTools::interpolate_boundary_values(const
dealii::Mapping< DH::dimension, DH::space_dimension>&, const DH&, const
typename dealii::FunctionMap< DH::space_dimension>::type&, std::map<unsigned
int, double, std::less<unsigned int>, std::allocator<std::pair<const unsigned
int, double> > >&,const std::vector<bool, std::allocator<bool> >&) [with DH =
dealii::DoFHandler<2, 2>]
The violated condition was:
cell->get_fe().is_primitive (i)
The name and call sequence of the exception was:
ExcMessage ("This function can only deal with requested boundary " "values
that correspond to primitive (scalar) base " "elements")
Additional Information:
This function can only deal with requested boundary values that correspond to
primitive (scalar) base elements
An error occurred in line <5359> of file <source/dofs/dof_tools.cc> in function
static void dealii::DoFTools::make_zero_boundary_constraints(const DH<dim,
spacedim>&, dealii::ConstraintMatrix&, const std::vector<bool,
std::allocator<bool> >&) [with int dim = 2, int spacedim = 2, DH =
dealii::DoFHandler]
The violated condition was:
cell->get_fe().is_primitive (i)
The name and call sequence of the exception was:
ExcMessage ("This function can only deal with requested boundary " "values
that correspond to primitive (scalar) base " "elements")
Additional Information:
This function can only deal with requested boundary values that correspond to
primitive (scalar) base elements
Thanks,
Seshu
-----Original Message-----
From: Wolfgang Bangerth [mailto:[email protected]]
Sent: Tue 3/2/2010 9:17 AM
To: [email protected]
Cc: Tirupathi, Seshu
Subject: Re: [deal.II] Step 21 Boundary Conditions
> Can anyone help me with this error, my approach or a better way to impose
> neumann bcs.
The VectorTools::compute_no_normal_flux function can't deal with the
Raviart-Thomas element you are using.
That said, since the only degrees of freedom on the faces of cells are the
normal fluxes, you could use
VectorTools::interpolate_boundary_values
with a ZeroFunction object, or alternatively
DoFToolsc
This will have the same effect.
Best
W.
--
-------------------------------------------------------------------------
Wolfgang Bangerth email: [email protected]
www: http://www.math.tamu.edu/~bangerth/
_______________________________________________
dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
