Can you come up with a small program that triggers this error? Take a
look at the existing test cases for this function in
http://dealii.org/svn/dealii/trunk/tests/deal.II/no_flux_01.cc
through
http://dealii.org/svn/dealii/trunk/tests/deal.II/no_flux_09.cc
to see how such a test should ideally look like.
Sure, here's what I came up with. From what I can see, it looks like
the difference between this program and test case 8 is the use of
HalfHyperShellBoundary.
-Jennifer
#include <grid/tria.h>
#include <grid/grid_generator.h>
#include <grid/tria_boundary_lib.h>
#include <dofs/dof_handler.h>
#include <fe/fe_system.h>
#include <fe/mapping_q1.h>
#include <numerics/vectors.h>
using namespace dealii;
template <int dim>
void run()
{
Triangulation<dim> triangulation;
FESystem<dim> fe(FE_Q<dim>(1), dim);
DoFHandler<dim> dof_handler (triangulation);
ConstraintMatrix constraints;
MappingQ<dim> mapping(1);
GridGenerator::half_hyper_shell (triangulation,
Point<dim>(),
0.5,
1);
static HalfHyperShellBoundary<dim> boundary;
triangulation.set_boundary (0, boundary);
dof_handler.distribute_dofs (fe);
std::set<unsigned char> no_normal_flux_boundaries;
no_normal_flux_boundaries.insert (0);
VectorTools::compute_no_normal_flux_constraints
(dof_handler, 0,
no_normal_flux_boundaries,
constraints, mapping);
}
int main (int argc, char *argv[])
{
deallog.depth_console (0);
run<2>();
}
_______________________________________________
dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
