boundary_quadrature = new QGauss2<dim-1>();
break;
case 1:
- fe = new FELinear<dim>();
+ fe = new FEQ1<dim>();
quadrature = new QGauss3<dim>();
boundary_quadrature = new QGauss2<dim-1>();
break;
case 2:
- fe = new FEQuadraticSub<dim>();
+ fe = new FEQ2<dim>();
quadrature = new QGauss4<dim>();
boundary_quadrature = new QGauss3<dim-1>();
break;
case 3:
- fe = new FECubicSub<dim>();
+ fe = new FEQ3<dim>();
quadrature = new QGauss5<dim>();
boundary_quadrature = new QGauss4<dim-1>();
break;
case 4:
- fe = new FEQuarticSub<dim>();
+ fe = new FEQ4<dim>();
quadrature = new QGauss6<dim>();
boundary_quadrature = new QGauss5<dim-1>();
break;
cout << " Distributing degrees of freedom..." << endl;
- FELinear<dim> fe;
+ FEQ1<dim> fe;
dof->distribute_dofs (fe);
cout << " Renumbering degrees of freedom..." << endl;
solution_function = new Solution<dim>::Kink ();
- FELinear<dim> fe;
+ FEQ1<dim> fe;
QGauss3<dim> quadrature;
PoissonEquation<dim> *equation;
boundary_quadrature = new QGauss2<dim-1>();
break;
case 1:
- fe = new FELinear<dim>();
+ fe = new FEQ1<dim>();
quadrature = new QGauss3<dim>();
boundary_quadrature = new QGauss2<dim-1>();
break;
case 2:
- fe = new FEQuadraticSub<dim>();
+ fe = new FEQ2<dim>();
quadrature = new QGauss4<dim>();
boundary_quadrature = new QGauss3<dim-1>();
break;
case 3:
- fe = new FECubicSub<dim>();
+ fe = new FEQ3<dim>();
quadrature = new QGauss5<dim>();
boundary_quadrature = new QGauss4<dim-1>();
break;
case 4:
- fe = new FEQuarticSub<dim>();
+ fe = new FEQ4<dim>();
quadrature = new QGauss6<dim>();
boundary_quadrature = new QGauss5<dim-1>();
break;
if (!set_boundary_values (prm))
return;
- FELinear<dim> fe;
+ FEQ1<dim> fe;
PoissonEquation<dim> equation (*rhs);
QGauss2<dim> quadrature;
boundary_quadrature = new QGauss2<dim-1>();
break;
case 1:
- fe = new FELinear<dim>();
+ fe = new FEQ1<dim>();
quadrature = new QGauss3<dim>();
boundary_quadrature = new QGauss2<dim-1>();
break;
case 2:
- fe = new FEQuadraticSub<dim>();
+ fe = new FEQ2<dim>();
quadrature = new QGauss4<dim>();
boundary_quadrature = new QGauss3<dim-1>();
break;
case 3:
- fe = new FECubicSub<dim>();
+ fe = new FEQ3<dim>();
quadrature = new QGauss5<dim>();
boundary_quadrature = new QGauss4<dim-1>();
break;
case 4:
- fe = new FEQuarticSub<dim>();
+ fe = new FEQ4<dim>();
quadrature = new QGauss6<dim>();
boundary_quadrature = new QGauss5<dim-1>();
break;
cout << " Distributing degrees of freedom..." << endl;
- FELinear<dim> fe;
+ FEQ1<dim> fe;
dof->distribute_dofs (fe);
cout << " Renumbering degrees of freedom..." << endl;
solution_function = new Solution<dim>::Kink ();
- FELinear<dim> fe;
+ FEQ1<dim> fe;
QGauss3<dim> quadrature;
PoissonEquation<dim> *equation;
boundary_quadrature = new QGauss2<dim-1>();
break;
case 1:
- fe = new FELinear<dim>();
+ fe = new FEQ1<dim>();
quadrature = new QGauss3<dim>();
boundary_quadrature = new QGauss2<dim-1>();
break;
case 2:
- fe = new FEQuadraticSub<dim>();
+ fe = new FEQ2<dim>();
quadrature = new QGauss4<dim>();
boundary_quadrature = new QGauss3<dim-1>();
break;
case 3:
- fe = new FECubicSub<dim>();
+ fe = new FEQ3<dim>();
quadrature = new QGauss5<dim>();
boundary_quadrature = new QGauss4<dim-1>();
break;
case 4:
- fe = new FEQuarticSub<dim>();
+ fe = new FEQ4<dim>();
quadrature = new QGauss6<dim>();
boundary_quadrature = new QGauss5<dim-1>();
break;
if (!set_boundary_values (prm))
return;
- FELinear<dim> fe;
+ FEQ1<dim> fe;
PoissonEquation<dim> equation (*rhs);
QGauss2<dim> quadrature;
--- /dev/null
+Make.depend