// $Id$
// Version: $Name$
//
-// Copyright (C) 2005, 2006, 2007 by the deal.II authors
+// Copyright (C) 2005, 2006, 2007, 2008 by the deal.II authors
//
// This file is subject to QPL and may not be distributed
// without copyright and license information. Please refer
* but since the multilevel
* method does not have its own
* right hand side, the function
- * values are igneored.
+ * values are ignored.
*
* @arg <tt>boundary_indices</tt>
- * is a vector which contains
- * for each level all indices of
+ * is a vector which on return
+ * contains all indices of
* boundary constraint degrees of
- * freedom.
+ * freedom for each level. Its
+ * length has to match the number
+ * of levels.
*/
template <int dim>
static void make_boundary_list(
// $Id$
// Version: $Name$
//
-// Copyright (C) 1999, 2000, 2001, 2002, 2003, 2004, 2005, 2006, 2007 by the deal.II authors
+// Copyright (C) 1999, 2000, 2001, 2002, 2003, 2004, 2005, 2006, 2007, 2008 by the deal.II authors
//
// This file is subject to QPL and may not be distributed
// without copyright and license information. Please refer
// immediately
if (function_map.size() == 0)
return;
-
- const unsigned int n_components = DoFTools::n_components(dof);
- const bool fe_is_system = (n_components != 1);
+ const unsigned int n_components = DoFTools::n_components(dof);
+ const unsigned int n_levels = dof.get_tria().n_levels();
+ const bool fe_is_system = (n_components != 1);
+
+ AssertDimension (boundary_indices.size(), n_levels);
+
// for (typename FunctionMap<dim>::type::const_iterator i=function_map.begin();
// i!=function_map.end(); ++i)
// Assert (n_components == i->second->n_components,