// $Id$
// Version: $Name$
//
-// Copyright (C) 2003, 2006, 2007, 2008, 2009 by the deal.II authors
+// Copyright (C) 2003, 2006, 2007, 2008, 2009, 2010 by the deal.II authors
//
// This file is subject to QPL and may not be distributed
// without copyright and license information. Please refer
namespace internal
{
-
+
namespace hp
{
// -------------------------- FEValuesBase -------------------------
present_fe_values_index = TableIndices<3> (fe_index,
mapping_index,
q_index);
-
+
// first check whether we
// already have an object for
// this particular combination
// of indices
- if (fe_values_table(present_fe_values_index) == 0)
+ if (fe_values_table(present_fe_values_index).get() == 0)
fe_values_table(present_fe_values_index)
=
std_cxx1x::shared_ptr<FEValues>
// now there definitely is one!
return *fe_values_table(present_fe_values_index);
- }
+ }
}
}
namespace hp
{
-
+
// -------------------------- FEValues -------------------------
else
real_q_index = 0;
}
-
+
if (real_mapping_index == numbers::invalid_unsigned_int)
{
if (this->mapping_collection->size() > 1)
else
real_mapping_index = 0;
}
-
+
if (real_fe_index == numbers::invalid_unsigned_int)
real_fe_index = cell->active_fe_index();
ExcIndexRange (real_mapping_index, 0, this->mapping_collection->size()));
Assert (real_fe_index < this->fe_collection->size(),
ExcIndexRange (real_fe_index, 0, this->fe_collection->size()));
-
+
// now finally actually get the
// corresponding object and
// initialize it
if (real_q_index == numbers::invalid_unsigned_int)
real_q_index = 0;
-
+
if (real_mapping_index == numbers::invalid_unsigned_int)
real_mapping_index = 0;
ExcIndexRange (real_mapping_index, 0, this->mapping_collection->size()));
Assert (real_fe_index < this->fe_collection->size(),
ExcIndexRange (real_fe_index, 0, this->fe_collection->size()));
-
+
// now finally actually get the
// corresponding object and
// initialize it
if (real_q_index == numbers::invalid_unsigned_int)
real_q_index = 0;
-
+
if (real_mapping_index == numbers::invalid_unsigned_int)
real_mapping_index = 0;
ExcIndexRange (real_mapping_index, 0, this->mapping_collection->size()));
Assert (real_fe_index < this->fe_collection->size(),
ExcIndexRange (real_fe_index, 0, this->fe_collection->size()));
-
+
// now finally actually get the
// corresponding object and
// initialize it
if (real_q_index == numbers::invalid_unsigned_int)
real_q_index = 0;
-
+
if (real_mapping_index == numbers::invalid_unsigned_int)
real_mapping_index = 0;
ExcIndexRange (real_mapping_index, 0, this->mapping_collection->size()));
Assert (real_fe_index < this->fe_collection->size(),
ExcIndexRange (real_fe_index, 0, this->fe_collection->size()));
-
+
// now finally actually get the
// corresponding object and
// initialize it
real_mapping_index,
real_q_index).reinit (cell);
}
-
+
// -------------------------- FEFaceValues -------------------------
else
real_q_index = 0;
}
-
+
if (real_mapping_index == numbers::invalid_unsigned_int)
{
if (this->mapping_collection->size() > 1)
else
real_mapping_index = 0;
}
-
+
if (real_fe_index == numbers::invalid_unsigned_int)
real_fe_index = cell->active_fe_index();
ExcIndexRange (real_mapping_index, 0, this->mapping_collection->size()));
Assert (real_fe_index < this->fe_collection->size(),
ExcIndexRange (real_fe_index, 0, this->fe_collection->size()));
-
+
// now finally actually get the
// corresponding object and
// initialize it
if (real_q_index == numbers::invalid_unsigned_int)
real_q_index = 0;
-
+
if (real_mapping_index == numbers::invalid_unsigned_int)
real_mapping_index = 0;
ExcIndexRange (real_mapping_index, 0, this->mapping_collection->size()));
Assert (real_fe_index < this->fe_collection->size(),
ExcIndexRange (real_fe_index, 0, this->fe_collection->size()));
-
+
// now finally actually get the
// corresponding object and
// initialize it
if (real_q_index == numbers::invalid_unsigned_int)
real_q_index = 0;
-
+
if (real_mapping_index == numbers::invalid_unsigned_int)
real_mapping_index = 0;
ExcIndexRange (real_mapping_index, 0, this->mapping_collection->size()));
Assert (real_fe_index < this->fe_collection->size(),
ExcIndexRange (real_fe_index, 0, this->fe_collection->size()));
-
+
// now finally actually get the
// corresponding object and
// initialize it
if (real_q_index == numbers::invalid_unsigned_int)
real_q_index = 0;
-
+
if (real_mapping_index == numbers::invalid_unsigned_int)
real_mapping_index = 0;
ExcIndexRange (real_mapping_index, 0, this->mapping_collection->size()));
Assert (real_fe_index < this->fe_collection->size(),
ExcIndexRange (real_fe_index, 0, this->fe_collection->size()));
-
+
// now finally actually get the
// corresponding object and
// initialize it
real_mapping_index,
real_q_index).reinit (cell, face_no);
}
-
+
// -------------------------- FESubfaceValues -------------------------
else
real_q_index = 0;
}
-
+
if (real_mapping_index == numbers::invalid_unsigned_int)
{
if (this->mapping_collection->size() > 1)
else
real_mapping_index = 0;
}
-
+
if (real_fe_index == numbers::invalid_unsigned_int)
real_fe_index = cell->active_fe_index();
ExcIndexRange (real_mapping_index, 0, this->mapping_collection->size()));
Assert (real_fe_index < this->fe_collection->size(),
ExcIndexRange (real_fe_index, 0, this->fe_collection->size()));
-
+
// now finally actually get the
// corresponding object and
// initialize it
if (real_q_index == numbers::invalid_unsigned_int)
real_q_index = 0;
-
+
if (real_mapping_index == numbers::invalid_unsigned_int)
real_mapping_index = 0;
ExcIndexRange (real_mapping_index, 0, this->mapping_collection->size()));
Assert (real_fe_index < this->fe_collection->size(),
ExcIndexRange (real_fe_index, 0, this->fe_collection->size()));
-
+
// now finally actually get the
// corresponding object and
// initialize it
if (real_q_index == numbers::invalid_unsigned_int)
real_q_index = 0;
-
+
if (real_mapping_index == numbers::invalid_unsigned_int)
real_mapping_index = 0;
ExcIndexRange (real_mapping_index, 0, this->mapping_collection->size()));
Assert (real_fe_index < this->fe_collection->size(),
ExcIndexRange (real_fe_index, 0, this->fe_collection->size()));
-
+
// now finally actually get the
// corresponding object and
// initialize it
if (real_q_index == numbers::invalid_unsigned_int)
real_q_index = 0;
-
+
if (real_mapping_index == numbers::invalid_unsigned_int)
real_mapping_index = 0;
ExcIndexRange (real_mapping_index, 0, this->mapping_collection->size()));
Assert (real_fe_index < this->fe_collection->size(),
ExcIndexRange (real_fe_index, 0, this->fe_collection->size()));
-
+
// now finally actually get the
// corresponding object and
// initialize it
{
template class FEValuesBase<deal_II_dimension,deal_II_dimension,
dealii::FEValues<deal_II_dimension,deal_II_dimension+1> >;
-// not yet implemented:
+// not yet implemented:
// #if deal_II_dimension == 2
// template class FEValuesBase<deal_II_dimension,deal_II_dimension-1,
// dealii::FEFaceValues<deal_II_dimension> >;
{
template class FEValues<deal_II_dimension, deal_II_dimension+1>;
-// not yet implemented:
+// not yet implemented:
//#if deal_II_dimension == 2
// template class FEFaceValues<deal_II_dimension, deal_II_dimension+1>;
// template class FESubfaceValues<deal_II_dimension, deal_II_dimension+1>;