From: Wolfgang Bangerth Date: Thu, 23 Feb 2006 00:59:35 +0000 (+0000) Subject: Implement all these reinit() functions also for non-hp iterators. X-Git-Tag: v8.0.0~12221 X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=14cfa7a6d3240a33a6d9c948ce244ecdabb903dd;p=dealii.git Implement all these reinit() functions also for non-hp iterators. git-svn-id: https://svn.dealii.org/trunk@12461 0785d39b-7218-0410-832d-ea1e28bc413d --- diff --git a/deal.II/deal.II/include/fe/hp_fe_values.h b/deal.II/deal.II/include/fe/hp_fe_values.h index b270eb2071..b550f8c96b 100644 --- a/deal.II/deal.II/include/fe/hp_fe_values.h +++ b/deal.II/deal.II/include/fe/hp_fe_values.h @@ -350,6 +350,112 @@ namespace hp */ void reinit (const typename hp::DoFHandler::cell_iterator &cell, + const unsigned int q_index = deal_II_numbers::invalid_unsigned_int, + const unsigned int mapping_index = deal_II_numbers::invalid_unsigned_int, + const unsigned int fe_index = deal_II_numbers::invalid_unsigned_int); + + /** + * Like the previous function, + * but for non-hp + * iterators. The reason this + * (and the other non-hp + * iterator) function exists is + * so that one can use + * hp::FEValues not only for + * hp::DoFhandler objects, but + * for all sorts of DoFHandler + * objects, and triangulations + * not associated with + * DoFHandlers in general. + * + * Since + * cell-@>active_fe_index() + * always returns zero for + * non-hp iterators, this + * function chooses the zero-th + * finite element, mapping, and + * quadrature object from the + * relevant constructions + * passed to the constructor of + * this object. The only + * exception is if you specify + * a value different from the + * default value for any of + * these last three arguments. + */ + void + reinit (const typename ::DoFHandler::cell_iterator &cell, + const unsigned int q_index = deal_II_numbers::invalid_unsigned_int, + const unsigned int mapping_index = deal_II_numbers::invalid_unsigned_int, + const unsigned int fe_index = deal_II_numbers::invalid_unsigned_int); + + /** + * Like the previous function, + * but for non-hp + * iterators. The reason this + * (and the other non-hp + * iterator) function exists is + * so that one can use + * hp::FEValues not only for + * hp::DoFhandler objects, but + * for all sorts of DoFHandler + * objects, and triangulations + * not associated with + * DoFHandlers in general. + * + * Since + * cell-@>active_fe_index() + * always returns zero for + * non-hp iterators, this + * function chooses the zero-th + * finite element, mapping, and + * quadrature object from the + * relevant constructions + * passed to the constructor of + * this object. The only + * exception is if you specify + * a value different from the + * default value for any of + * these last three arguments. + */ + void + reinit (const typename MGDoFHandler::cell_iterator &cell, + const unsigned int q_index = deal_II_numbers::invalid_unsigned_int, + const unsigned int mapping_index = deal_II_numbers::invalid_unsigned_int, + const unsigned int fe_index = deal_II_numbers::invalid_unsigned_int); + + /** + * Like the previous function, + * but for non-hp + * iterators. The reason this + * (and the other non-hp + * iterator) function exists is + * so that one can use + * hp::FEValues not only for + * hp::DoFhandler objects, but + * for all sorts of DoFHandler + * objects, and triangulations + * not associated with + * DoFHandlers in general. + * + * Since + * cell-@>active_fe_index() + * doesn't make sense for + * triangulation iterators, + * this function chooses the + * zero-th finite element, + * mapping, and quadrature + * object from the relevant + * constructions passed to the + * constructor of this + * object. The only exception + * is if you specify a value + * different from the default + * value for any of these last + * three arguments. + */ + void + reinit (const typename Triangulation::cell_iterator &cell, const unsigned int q_index = deal_II_numbers::invalid_unsigned_int, const unsigned int mapping_index = deal_II_numbers::invalid_unsigned_int, const unsigned int fe_index = deal_II_numbers::invalid_unsigned_int); @@ -518,6 +624,115 @@ namespace hp */ void reinit (const typename hp::DoFHandler::cell_iterator &cell, + const unsigned int face_no, + const unsigned int q_index = deal_II_numbers::invalid_unsigned_int, + const unsigned int mapping_index = deal_II_numbers::invalid_unsigned_int, + const unsigned int fe_index = deal_II_numbers::invalid_unsigned_int); + + /** + * Like the previous function, + * but for non-hp + * iterators. The reason this + * (and the other non-hp + * iterator) function exists is + * so that one can use + * hp::FEValues not only for + * hp::DoFhandler objects, but + * for all sorts of DoFHandler + * objects, and triangulations + * not associated with + * DoFHandlers in general. + * + * Since + * cell-@>active_fe_index() + * always returns zero for + * non-hp iterators, this + * function chooses the zero-th + * finite element, mapping, and + * quadrature object from the + * relevant constructions + * passed to the constructor of + * this object. The only + * exception is if you specify + * a value different from the + * default value for any of + * these last three arguments. + */ + void + reinit (const typename ::DoFHandler::cell_iterator &cell, + const unsigned int face_no, + const unsigned int q_index = deal_II_numbers::invalid_unsigned_int, + const unsigned int mapping_index = deal_II_numbers::invalid_unsigned_int, + const unsigned int fe_index = deal_II_numbers::invalid_unsigned_int); + + /** + * Like the previous function, + * but for non-hp + * iterators. The reason this + * (and the other non-hp + * iterator) function exists is + * so that one can use + * hp::FEValues not only for + * hp::DoFhandler objects, but + * for all sorts of DoFHandler + * objects, and triangulations + * not associated with + * DoFHandlers in general. + * + * Since + * cell-@>active_fe_index() + * always returns zero for + * non-hp iterators, this + * function chooses the zero-th + * finite element, mapping, and + * quadrature object from the + * relevant constructions + * passed to the constructor of + * this object. The only + * exception is if you specify + * a value different from the + * default value for any of + * these last three arguments. + */ + void + reinit (const typename MGDoFHandler::cell_iterator &cell, + const unsigned int face_no, + const unsigned int q_index = deal_II_numbers::invalid_unsigned_int, + const unsigned int mapping_index = deal_II_numbers::invalid_unsigned_int, + const unsigned int fe_index = deal_II_numbers::invalid_unsigned_int); + + /** + * Like the previous function, + * but for non-hp + * iterators. The reason this + * (and the other non-hp + * iterator) function exists is + * so that one can use + * hp::FEValues not only for + * hp::DoFhandler objects, but + * for all sorts of DoFHandler + * objects, and triangulations + * not associated with + * DoFHandlers in general. + * + * Since + * cell-@>active_fe_index() + * doesn't make sense for + * triangulation iterators, + * this function chooses the + * zero-th finite element, + * mapping, and quadrature + * object from the relevant + * constructions passed to the + * constructor of this + * object. The only exception + * is if you specify a value + * different from the default + * value for any of these last + * three arguments. + */ + void + reinit (const typename Triangulation::cell_iterator &cell, const unsigned int face_no, const unsigned int q_index = deal_II_numbers::invalid_unsigned_int, const unsigned int mapping_index = deal_II_numbers::invalid_unsigned_int, @@ -666,6 +881,118 @@ namespace hp */ void reinit (const typename hp::DoFHandler::cell_iterator &cell, + const unsigned int face_no, + const unsigned int subface_no, + const unsigned int q_index = deal_II_numbers::invalid_unsigned_int, + const unsigned int mapping_index = deal_II_numbers::invalid_unsigned_int, + const unsigned int fe_index = deal_II_numbers::invalid_unsigned_int); + + /** + * Like the previous function, + * but for non-hp + * iterators. The reason this + * (and the other non-hp + * iterator) function exists is + * so that one can use + * hp::FEValues not only for + * hp::DoFhandler objects, but + * for all sorts of DoFHandler + * objects, and triangulations + * not associated with + * DoFHandlers in general. + * + * Since + * cell-@>active_fe_index() + * always returns zero for + * non-hp iterators, this + * function chooses the zero-th + * finite element, mapping, and + * quadrature object from the + * relevant constructions + * passed to the constructor of + * this object. The only + * exception is if you specify + * a value different from the + * default value for any of + * these last three arguments. + */ + void + reinit (const typename ::DoFHandler::cell_iterator &cell, + const unsigned int face_no, + const unsigned int subface_no, + const unsigned int q_index = deal_II_numbers::invalid_unsigned_int, + const unsigned int mapping_index = deal_II_numbers::invalid_unsigned_int, + const unsigned int fe_index = deal_II_numbers::invalid_unsigned_int); + + /** + * Like the previous function, + * but for non-hp + * iterators. The reason this + * (and the other non-hp + * iterator) function exists is + * so that one can use + * hp::FEValues not only for + * hp::DoFhandler objects, but + * for all sorts of DoFHandler + * objects, and triangulations + * not associated with + * DoFHandlers in general. + * + * Since + * cell-@>active_fe_index() + * always returns zero for + * non-hp iterators, this + * function chooses the zero-th + * finite element, mapping, and + * quadrature object from the + * relevant constructions + * passed to the constructor of + * this object. The only + * exception is if you specify + * a value different from the + * default value for any of + * these last three arguments. + */ + void + reinit (const typename MGDoFHandler::cell_iterator &cell, + const unsigned int face_no, + const unsigned int subface_no, + const unsigned int q_index = deal_II_numbers::invalid_unsigned_int, + const unsigned int mapping_index = deal_II_numbers::invalid_unsigned_int, + const unsigned int fe_index = deal_II_numbers::invalid_unsigned_int); + + /** + * Like the previous function, + * but for non-hp + * iterators. The reason this + * (and the other non-hp + * iterator) function exists is + * so that one can use + * hp::FEValues not only for + * hp::DoFhandler objects, but + * for all sorts of DoFHandler + * objects, and triangulations + * not associated with + * DoFHandlers in general. + * + * Since + * cell-@>active_fe_index() + * doesn't make sense for + * triangulation iterators, + * this function chooses the + * zero-th finite element, + * mapping, and quadrature + * object from the relevant + * constructions passed to the + * constructor of this + * object. The only exception + * is if you specify a value + * different from the default + * value for any of these last + * three arguments. + */ + void + reinit (const typename Triangulation::cell_iterator &cell, const unsigned int face_no, const unsigned int subface_no, const unsigned int q_index = deal_II_numbers::invalid_unsigned_int,