std::vector< Tensor<1,dim> > &viq = values[maps[i][q]];
const unsigned int s = vgrads[q].size();
viq.resize(s);
- for (unsigned int l=0;l<s;l++)
+ for (unsigned int l=0; l<s; l++)
viq[l] = vgrads[q][l];
}
}
IDScratchData<dim,spacedim>::IDScratchData (const IDScratchData &data)
:
x_fe_values(data.x_fe_values.get_mapping_collection(),
- data.x_fe_values.get_fe_collection(),
- data.x_fe_values.get_quadrature_collection(),
- data.x_fe_values.get_update_flags())
+ data.x_fe_values.get_fe_collection(),
+ data.x_fe_values.get_quadrature_collection(),
+ data.x_fe_values.get_update_flags())
{}
template <int dim, int spacedim>
&shape_derivatives[shape_function_data[shape_function].row_index][0];
for (unsigned int q_point=0; q_point<n_quadrature_points; ++q_point)
derivatives[q_point] += value *
- typename ProductType<Number,dealii::Tensor<order,spacedim> >::type(*shape_derivative_ptr++);
+ typename ProductType<Number,dealii::Tensor<order,spacedim> >::type(*shape_derivative_ptr++);
}
}
&shape_gradients[snc][0];
for (unsigned int q_point=0; q_point<n_quadrature_points; ++q_point)
symmetric_gradients[q_point] += value *
- typename ProductType<Number,dealii::SymmetricTensor<2,spacedim> >::type (symmetrize_single_row(comp, *shape_gradient_ptr++));
+ typename ProductType<Number,dealii::SymmetricTensor<2,spacedim> >::type (symmetrize_single_row(comp, *shape_gradient_ptr++));
}
else
for (unsigned int q_point=0; q_point<n_quadrature_points; ++q_point)
const UpdateFlags update_flags)
:
FEValuesBase<dim,spacedim> (q.size(),
- fe.dofs_per_cell,
- update_default,
- mapping,
- fe),
+ fe.dofs_per_cell,
+ update_default,
+ mapping,
+ fe),
quadrature (q)
{
initialize (update_flags);
const UpdateFlags update_flags)
:
FEValuesBase<dim,spacedim> (q.size(),
- fe.dofs_per_cell,
- update_default,
- StaticMappingQ1<dim,spacedim>::mapping,
- fe),
+ fe.dofs_per_cell,
+ update_default,
+ StaticMappingQ1<dim,spacedim>::mapping,
+ fe),
quadrature (q)
{
initialize (update_flags);
const Quadrature<dim-1>& quadrature)
:
FEValuesBase<dim,spacedim> (n_q_points,
- dofs_per_cell,
- update_default,
- mapping,
- fe),
+ dofs_per_cell,
+ update_default,
+ mapping,
+ fe),
quadrature(quadrature)
{}
const UpdateFlags update_flags)
:
FEFaceValuesBase<dim,spacedim> (quadrature.size(),
- fe.dofs_per_cell,
- update_flags,
- mapping,
- fe, quadrature)
+ fe.dofs_per_cell,
+ update_flags,
+ mapping,
+ fe, quadrature)
{
initialize (update_flags);
}
const UpdateFlags update_flags)
:
FEFaceValuesBase<dim,spacedim> (quadrature.size(),
- fe.dofs_per_cell,
- update_flags,
- StaticMappingQ1<dim,spacedim>::mapping,
- fe, quadrature)
+ fe.dofs_per_cell,
+ update_flags,
+ StaticMappingQ1<dim,spacedim>::mapping,
+ fe, quadrature)
{
initialize (update_flags);
}
const UpdateFlags update_flags)
:
FEFaceValuesBase<dim,spacedim> (quadrature.size(),
- fe.dofs_per_cell,
- update_flags,
- mapping,
- fe, quadrature)
+ fe.dofs_per_cell,
+ update_flags,
+ mapping,
+ fe, quadrature)
{
initialize (update_flags);
}
const UpdateFlags update_flags)
:
FEFaceValuesBase<dim,spacedim> (quadrature.size(),
- fe.dofs_per_cell,
- update_flags,
- StaticMappingQ1<dim,spacedim>::mapping,
- fe, quadrature)
+ fe.dofs_per_cell,
+ update_flags,
+ StaticMappingQ1<dim,spacedim>::mapping,
+ fe, quadrature)
{
initialize (update_flags);
}
{
namespace Vector
{
- template void copy_vector<int,double> (const dealii::Vector<int>&,
- dealii::Vector<double>&);
+ template void copy_vector<int,double> (const dealii::Vector<int> &,
+ dealii::Vector<double> &);
}
}
template
-void Vector<int>::reinit<double>(const Vector<double>&, const bool);
+void Vector<int>::reinit<double>(const Vector<double> &, const bool);
// do a few functions that currently don't fit the scheme because they have
//FIXME
std::vector<typename VectorType::value_type> tmp (patch_values.size());
fe_patch_values.get_function_values (*vector, tmp);
- for (unsigned int i = 0; i < tmp.size();i++)
+ for (unsigned int i = 0; i < tmp.size(); i++)
patch_values[i] = tmp[i];
}
{
//FIXME
std::vector<std::vector<Tensor<1,DH::space_dimension,typename VectorType::value_type> > > tmp(patch_gradients_system.size());
- for (unsigned int i = 0; i < tmp.size();i++)
+ for (unsigned int i = 0; i < tmp.size(); i++)
tmp[i].resize(patch_gradients_system.size());
fe_patch_values.get_function_gradients (*vector, tmp);
- for (unsigned int i = 0; i < tmp.size();i++)
- for (unsigned int j = 0; j < tmp[i].size();j++)
+ for (unsigned int i = 0; i < tmp.size(); i++)
+ for (unsigned int j = 0; j < tmp[i].size(); j++)
patch_gradients_system[i][j] = tmp[i][j];
}
std::vector<Tensor<1,DH::space_dimension,typename VectorType::value_type> > tmp;
tmp.resize(patch_gradients.size());
fe_patch_values.get_function_gradients (*vector, tmp);
- for (unsigned int i = 0; i < tmp.size();i++)
+ for (unsigned int i = 0; i < tmp.size(); i++)
patch_gradients[i] = tmp[i];
}
else
{
std::vector<std::vector<Tensor<1,dim,typename InputVector::value_type>> > values
- (1, std::vector<Tensor<1,dim,typename InputVector::value_type>>(fe_values.get_fe().n_components()));
+ (1, std::vector<Tensor<1,dim,typename InputVector::value_type>>(fe_values.get_fe().n_components()));
fe_values.get_function_gradients (solution, values);
return ProjectedDerivative(values[0][component]);
};
else
{
std::vector<std::vector<Tensor<2,dim,typename InputVector::value_type>> > values
- (1, std::vector<Tensor<2,dim,typename InputVector::value_type>>(fe_values.get_fe().n_components()));
+ (1, std::vector<Tensor<2,dim,typename InputVector::value_type>>(fe_values.get_fe().n_components()));
fe_values.get_function_hessians (solution, values);
return ProjectedDerivative(values[0][component]);
};
finite_element (fe),
face_quadratures (face_quadratures),
fe_face_values_cell (mapping,
- finite_element,
- face_quadratures,
- update_gradients |
- update_JxW_values |
- (need_quadrature_points ?
- update_quadrature_points :
- UpdateFlags()) |
- update_normal_vectors),
+ finite_element,
+ face_quadratures,
+ update_gradients |
+ update_JxW_values |
+ (need_quadrature_points ?
+ update_quadrature_points :
+ UpdateFlags()) |
+ update_normal_vectors),
fe_face_values_neighbor (mapping,
- finite_element,
- face_quadratures,
- update_gradients),
+ finite_element,
+ face_quadratures,
+ update_gradients),
fe_subface_values (mapping,
- finite_element,
- face_quadratures,
- update_gradients),
+ finite_element,
+ face_quadratures,
+ update_gradients),
phi (n_solution_vectors,
- std::vector<std::vector<number> >
- (face_quadratures.max_n_quadrature_points(),
- std::vector<number> (fe.n_components()))),
+ std::vector<std::vector<number> >
+ (face_quadratures.max_n_quadrature_points(),
+ std::vector<number> (fe.n_components()))),
psi (n_solution_vectors,
- std::vector<std::vector<Tensor<1,spacedim,number> > >
- (face_quadratures.max_n_quadrature_points(),
- std::vector<Tensor<1,spacedim,number> > (fe.n_components()))),
+ std::vector<std::vector<Tensor<1,spacedim,number> > >
+ (face_quadratures.max_n_quadrature_points(),
+ std::vector<Tensor<1,spacedim,number> > (fe.n_components()))),
neighbor_psi (n_solution_vectors,
- std::vector<std::vector<Tensor<1,spacedim,number> > >
- (face_quadratures.max_n_quadrature_points(),
- std::vector<Tensor<1,spacedim,number> > (fe.n_components()))),
+ std::vector<std::vector<Tensor<1,spacedim,number> > >
+ (face_quadratures.max_n_quadrature_points(),
+ std::vector<Tensor<1,spacedim,number> > (fe.n_components()))),
normal_vectors (face_quadratures.max_n_quadrature_points()),
coefficient_values1 (face_quadratures.max_n_quadrature_points()),
coefficient_values (face_quadratures.max_n_quadrature_points(),
- dealii::Vector<double> (fe.n_components())),
+ dealii::Vector<double> (fe.n_components())),
JxW_values (face_quadratures.max_n_quadrature_points()),
subdomain_id (subdomain_id),
material_id (material_id),