clear ();
tria = new Triangulation<dim>();
- dof = new DoFHandler<dim> (tria);
+ dof = new DoFHandler<dim> (*tria);
set_tria_and_dof (tria, dof);
};
tria = new Triangulation<dim>();
GridGenerator::hyper_cube(*tria);
- dof = new DoFHandler<dim> (tria);
+ dof = new DoFHandler<dim> (*tria);
};
clear ();
tria = new Triangulation<dim>();
- dof = new DoFHandler<dim> (tria);
+ dof = new DoFHandler<dim> (*tria);
set_tria_and_dof (tria, dof);
boundary = new HyperBallBoundary<dim> ();
};
clear ();
tria = new Triangulation<dim>();
- dof = new MGDoFHandler<dim> (tria);
+ dof = new MGDoFHandler<dim> (*tria);
};
clear ();
tria = new Triangulation<dim>();
- dof = new DoFHandler<dim> (tria);
+ dof = new DoFHandler<dim> (*tria);
set_tria_and_dof (tria, dof);
};
* Constructor. Take #tria# as the
* triangulation to work on.
*/
- DoFHandler (Triangulation<dim> *tria);
-
- /**
- * Constructor. Does the same
- * thing as the one above, but
- * takes a reference instead of a
- * pointer. The pointer form is
- * the older one, but the one
- * encouraged for use now is this
- * one.
- */
DoFHandler (Triangulation<dim> &tria);
/**
* Constructor. Take #tria# as the
* triangulation to work on.
*/
- MGDoFHandler (Triangulation<dim> *tria);
+ MGDoFHandler (Triangulation<dim> &tria);
/**
* Destructor
* Constructor. Take #tria# as the
* triangulation to work on.
*/
- MGDoFHandler (Triangulation<dim> *tria);
+ MGDoFHandler (Triangulation<dim> &tria);
/**
* Destructor
-template <int dim>
-DoFHandler<dim>::DoFHandler (Triangulation<dim> *tria) :
- tria(tria),
- used_dofs (0)
-{
- Assert (tria != 0, ExcInvalidTriangulation());
-};
-
-
template <int dim>
DoFHandler<dim>::DoFHandler (Triangulation<dim> &tria) :
tria(&tria),
/* ------------------------ MGDoFHandler ------------------------------------- */
template <int dim>
-MGDoFHandler<dim>::MGDoFHandler (Triangulation<dim> *tria) :
+MGDoFHandler<dim>::MGDoFHandler (Triangulation<dim> &tria) :
DoFHandler<dim> (tria)
{};
/* ------------------------ MGDoFHandler ------------------------------------- */
template <int dim>
-MGDoFHandler<dim>::MGDoFHandler (Triangulation<dim> *tria) :
+MGDoFHandler<dim>::MGDoFHandler (Triangulation<dim> &tria) :
DoFHandler<dim> (tria)
{};
...
// then associate a DoFHandler object with this triangulation...
-DoFHandler<dim> dof_handler (&tria);
+DoFHandler<dim> dof_handler (tria);
// ...and distribute degrees of freedom for linear elements
FEQ1<dim> fe;
clear ();
tria = new Triangulation<dim>();
- dof = new DoFHandler<dim> (tria);
+ dof = new DoFHandler<dim> (*tria);
set_tria_and_dof (tria, dof);
};
tria = new Triangulation<dim>();
GridGenerator::hyper_cube(*tria);
- dof = new DoFHandler<dim> (tria);
+ dof = new DoFHandler<dim> (*tria);
};
clear ();
tria = new Triangulation<dim>();
- dof = new DoFHandler<dim> (tria);
+ dof = new DoFHandler<dim> (*tria);
set_tria_and_dof (tria, dof);
boundary = new HyperBallBoundary<dim> ();
};
clear ();
tria = new Triangulation<dim>();
- dof = new MGDoFHandler<dim> (tria);
+ dof = new MGDoFHandler<dim> (*tria);
};
clear ();
tria = new Triangulation<dim>();
- dof = new DoFHandler<dim> (tria);
+ dof = new DoFHandler<dim> (*tria);
set_tria_and_dof (tria, dof);
};
make_tria (tria, step);
GridOut().write_gnuplot (tria, logfile);
- DoFHandler<3> dof (&tria);
+ DoFHandler<3> dof (tria);
dof.distribute_dofs (*fe);
ConstraintMatrix constraints;
tria = new Triangulation<dim>();
GridGenerator::hyper_cube(*tria);
- dof = new DoFHandler<dim> (tria);
+ dof = new DoFHandler<dim> (*tria);
};
{
Triangulation<dim> tr;
GridGenerator::hyper_cube(tr,-1,1);
- DoFHandler<dim> dof(&tr);
+ DoFHandler<dim> dof(tr);
dof.distribute_dofs(fe);
StraightBoundary<dim> boundary;
DoFHandler<dim>::active_cell_iterator cell = dof.begin_active();
tria.begin_active()->vertex(2)(0) = 2;
FEQ1<2> fe;
- DoFHandler<2> dof(&tria);
+ DoFHandler<2> dof(tria);
dof.distribute_dofs(fe);
QTrapez<2> q;
for (unsigned int degree=1;degree<=3;degree++)
{
Triangulation<2> tr;
- MGDoFHandler<2> mgdof(&tr);
+ MGDoFHandler<2> mgdof(tr);
DoFHandler<2>& dof(mgdof);
GridGenerator::hyper_cube(tr,-1.,1.);
for (unsigned int degree=1;degree<=1;degree++)
{
Triangulation<2> tr;
- MGDoFHandler<2> mgdof(&tr);
+ MGDoFHandler<2> mgdof(tr);
DoFHandler<2>& dof(mgdof);
GridGenerator::hyper_cube(tr,-M_PI_2,M_PI_2);
GridGenerator::hyper_cube (tria,0,1);
FEQ1<2> fe;
- DoFHandler<2> dof(&tria);
+ DoFHandler<2> dof(tria);
dof.distribute_dofs(fe);
StraightBoundary<2> b;
sweep_info->get_timers().grid_generation.start();
- dof_handler = new DoFHandler<dim>(tria);
+ dof_handler = new DoFHandler<dim>(*tria);
dof_handler->distribute_dofs (fe);
if (parameters.renumber_dofs)