x_rel.compress();
ConstraintMatrix cm(relevant_set);
- DoFTools::make_hanging_node_constraints (* static_cast<DoFHandler<dim>* >(&dofh), cm);
+ DoFTools::make_hanging_node_constraints (dofh, cm);
std::vector<bool> velocity_mask (dim+1, true);
velocity_mask[dim] = false;
- VectorTools::interpolate_boundary_values (static_cast<const DoFHandler<dim>&>(dofh),
+ VectorTools::interpolate_boundary_values (dofh,
0,
ZeroFunction<dim>(dim+1),
cm,
DoFTools::extract_locally_active_dofs (dofh, dof_set);
ConstraintMatrix cm;
- DoFTools::make_hanging_node_constraints (* static_cast<DoFHandler<dim>* >(&dofh), cm);
+ DoFTools::make_hanging_node_constraints (dofh, cm);
ConstraintMatrix cm2(dof_set);
- DoFTools::make_hanging_node_constraints (* static_cast<DoFHandler<dim>* >(&dofh), cm2);
+ DoFTools::make_hanging_node_constraints (dofh, cm2);
{
std::ofstream file((std::string("p4est_2d_constraintmatrix_01/dat.") + Utilities::int_to_string(myid)).c_str());
DoFTools::extract_locally_active_dofs (dofh, dof_set);
ConstraintMatrix cm;
- DoFTools::make_hanging_node_constraints (* static_cast<DoFHandler<dim>* >(&dofh), cm);
+ DoFTools::make_hanging_node_constraints (dofh, cm);
ConstraintMatrix cm2(dof_set);
- DoFTools::make_hanging_node_constraints (* static_cast<DoFHandler<dim>* >(&dofh), cm2);
+ DoFTools::make_hanging_node_constraints (dofh, cm2);
if (myid==0)
{
x_rel.compress();
ConstraintMatrix cm(relevant_set);
- DoFTools::make_hanging_node_constraints (* static_cast<DoFHandler<dim>* >(&dofh), cm);
+ DoFTools::make_hanging_node_constraints (dofh, cm);
std::vector<bool> velocity_mask (dim+1, true);
velocity_mask[dim] = false;
- VectorTools::interpolate_boundary_values (static_cast<const DoFHandler<dim>&>(dofh),
+ VectorTools::interpolate_boundary_values (dofh,
0,
ZeroFunction<dim>(dim+1),
cm,
no_normal_flux_boundaries.insert (1);
- VectorTools::compute_no_normal_flux_constraints (static_cast<const DoFHandler<dim>&>(dofh), 0,
+ VectorTools::compute_no_normal_flux_constraints (dofh, 0,
no_normal_flux_boundaries,
cm);
TrilinosWrappers::MPI::Vector x;
x.reinit(owned_set, MPI_COMM_WORLD);
- VectorTools::interpolate(* static_cast<DoFHandler<dim>* >(&dofh),
+ VectorTools::interpolate(dofh,
TemperatureInitialValues<dim>(),
x);
x.compress();
x_rel.compress();
ConstraintMatrix cm(relevant_set);
- DoFTools::make_hanging_node_constraints (* static_cast<DoFHandler<dim>* >(&dofh), cm);
+ DoFTools::make_hanging_node_constraints (dofh, cm);
/* std::vector<bool> velocity_mask (dim+1, true);
velocity_mask[dim] = false;
- VectorTools::interpolate_boundary_values (static_cast<const DoFHandler<dim>&>(dofh),
+ VectorTools::interpolate_boundary_values (dofh,
0,
ZeroFunction<dim>(1),
cm,
TrilinosWrappers::MPI::Vector x_ref;
x_ref.reinit(owned_set, MPI_COMM_WORLD);
- VectorTools::interpolate(* static_cast<DoFHandler<dim>* >(&dofh),
+ VectorTools::interpolate(dofh,
TemperatureInitialValues<dim>(),
x_ref);
x_ref.compress();
TrilinosWrappers::MPI::Vector x_ref;
x_ref.reinit(owned_set, MPI_COMM_WORLD);
- VectorTools::interpolate(* static_cast<DoFHandler<dim>* >(&dofh),
+ VectorTools::interpolate(dofh,
ConstantFunction<dim> (1.),
x_ref);
x_ref.compress();
// anything
x1 = x_ref;
ConstraintMatrix cm(relevant_set);
- DoFTools::make_hanging_node_constraints (* static_cast<DoFHandler<dim>* >(&dofh), cm);
+ DoFTools::make_hanging_node_constraints (dofh, cm);
cm.distribute(x1);
x1 -= x_ref;
DoFTools::extract_locally_relevant_dofs (dofh, dof_set);
ConstraintMatrix cm;
- DoFTools::make_hanging_node_constraints (* static_cast<DoFHandler<dim>* >(&dofh), cm);
+ DoFTools::make_hanging_node_constraints (dofh, cm);
ConstraintMatrix cm2(dof_set);
- DoFTools::make_hanging_node_constraints (* static_cast<DoFHandler<dim>* >(&dofh), cm2);
+ DoFTools::make_hanging_node_constraints (dofh, cm2);
{
std::ofstream file((std::string("p4est_3d_constraintmatrix_01/dat.") + Utilities::int_to_string(myid)).c_str());
x_rel.compress();
ConstraintMatrix cm(relevant_set);
- DoFTools::make_hanging_node_constraints (* static_cast<DoFHandler<dim>* >(&dofh), cm);
+ DoFTools::make_hanging_node_constraints (dofh, cm);
//x.print(std::cout);
cm.distribute(x);
TrilinosWrappers::MPI::Vector x;
x.reinit(owned_set, MPI_COMM_WORLD);
- VectorTools::interpolate(* static_cast<DoFHandler<dim>* >(&dofh),
+ VectorTools::interpolate(dofh,
TemperatureInitialValues<dim>(),
x);
x.compress();
x_rel.compress();
ConstraintMatrix cm(relevant_set);
- DoFTools::make_hanging_node_constraints (* static_cast<DoFHandler<dim>* >(&dofh), cm);
+ DoFTools::make_hanging_node_constraints (dofh, cm);
/* std::vector<bool> velocity_mask (dim+1, true);
velocity_mask[dim] = false;
- VectorTools::interpolate_boundary_values (static_cast<const DoFHandler<dim>&>(dofh),
+ VectorTools::interpolate_boundary_values (dofh,
0,
ZeroFunction<dim>(1),
cm,
TrilinosWrappers::MPI::Vector x_ref;
x_ref.reinit(owned_set, MPI_COMM_WORLD);
- VectorTools::interpolate(* static_cast<DoFHandler<dim>* >(&dofh),
+ VectorTools::interpolate(dofh,
TemperatureInitialValues<dim>(),
x_ref);
x_ref.compress();
// $Id$
// Version: $Name$
//
-// Copyright (C) 2009, 2010 by the deal.II authors
+// Copyright (C) 2009, 2010, 2011 by the deal.II authors
//
// This file is subject to QPL and may not be distributed
// without copyright and license information. Please refer
IndexSet locally_relevant (dof.n_dofs());
DoFTools::extract_locally_relevant_dofs (dof, locally_relevant);
ConstraintMatrix constraints (locally_relevant);
- DoFTools::make_hanging_node_constraints (static_cast<const DoFHandler<dim>&>(dof), constraints);
+ DoFTools::make_hanging_node_constraints (dof, constraints);
constraints.close();
// print out constraints for each
const IndexSet set
= DoFTools::dof_indices_with_subdomain_association
- (static_cast<const DoFHandler<dim>&>(dofh),
+ (dofh,
tr.locally_owned_subdomain());
deallog << set.n_elements() << std::endl;