FiniteElement<dim> *
FE_DGQ<dim>::clone() const
{
- return new FE_DGQ<dim>(this->degree);
+ // TODO[Prill] : There must be a better way
+ // to extract 1D quadrature points from the
+ // tensor product FE.
+
+ // Construct a dummy quadrature formula
+ // containing the FE's nodes:
+ std::vector<Point<1> > qpoints(this->degree+1);
+ for (unsigned int i=0; i<=this->degree; ++i)
+ qpoints[i] = Point<1>(this->unit_support_points[i][0]);
+ Quadrature<1> pquadrature(qpoints);
+
+ return new FE_DGQ<dim>(pquadrature);
}
<ol>
+ <li> <p> Fixed: Corrected <code class="member">clone</code> method
+ of <code class="class">FE_DGQ</code> class for non-equidistant
+ support points.
+ <br>
+ (F. Prill 2006/10/31)
+ </p>
+
<li> <p> Improved: The lookup mechanism in <code
class="class">FETools</code>::<code
class="member">get_fe_from_name</code> has been changed, so
In the beginning the Universe was created. This has made a lot of
people very angry and has been widely regarded as a bad move.
Douglas Adams