};
-virtual Point<dim>
+ virtual Point<dim>
get_new_point_on_quad (const typename Triangulation<dim>::quad_iterator &quad) const
{
Point<dim> middle = StraightBoundary<dim>::get_new_point_on_quad(quad);
return middle;
-double x=middle(0),
+ double x=middle(0),
y=middle(1);
if (y<x)
};
+
template <int dim>
Point<dim>
CurvedLine<dim>::get_new_point_on_quad (const typename Triangulation<dim>::quad_iterator &quad) const
};
+
template <int dim>
class TestCases
{
};
+
template <int dim>
TestCases<dim>::TestCases () :
tria(0), dof(0) {};
+
template <int dim>
TestCases<dim>::~TestCases ()
{
};
+
template <int dim>
void TestCases<dim>::create_new ()
{
};
+
template <int dim>
void TestCases<dim>::run (const unsigned int test_case)
{
};
-deallog << " Distributing degrees of freedom..." << endl;
+ deallog << " Distributing degrees of freedom..." << endl;
FEQ1<dim> fe;
dof->distribute_dofs (fe);
dof->max_couplings_between_dofs());
-DoFTools::make_sparsity_pattern (*dof, sparsity);
+ DoFTools::make_sparsity_pattern (*dof, sparsity);
int unconstrained_bandwidth = sparsity.bandwidth();
deallog << " Writing sparsity pattern..." << endl;
sparsity.print_gnuplot (logfile);
-// computing constraints
+ // computing constraints
deallog << " Computing constraints..." << endl;
ConstraintMatrix constraints;
DoFTools::make_hanging_node_constraints (*dof, constraints);
sparsity.print_gnuplot (logfile);
-deallog << endl
- << " Total number of cells = " << tria->n_cells() << endl
- << " Total number of active cells = " << tria->n_active_cells() << endl
- << " Number of DoFs = " << dof->n_dofs() << endl
- << " Number of constraints = " << constraints.n_constraints() << endl
- << " Unconstrained matrix bandwidth= " << unconstrained_bandwidth << endl
- << " Constrained matrix bandwidth = " << sparsity.bandwidth()
- << endl << endl;
+ deallog << endl
+ << " Total number of cells = " << tria->n_cells() << endl
+ << " Total number of active cells = " << tria->n_active_cells() << endl
+ << " Number of DoFs = " << dof->n_dofs() << endl
+ << " Number of constraints = " << constraints.n_constraints() << endl
+ << " Unconstrained matrix bandwidth= " << unconstrained_bandwidth << endl
+ << " Constrained matrix bandwidth = " << sparsity.bandwidth()
+ << endl << endl;
// release the lock that dof has to the
// finite element object
};
+
int main ()
{
deallog.attach(logfile);