Skip to content
Snippets Groups Projects
Commit b2bd9d11 authored by Jonathan Youett's avatar Jonathan Youett
Browse files

Use ranged-based for-loop

parent 8ce35f94
No related branches found
No related tags found
No related merge requests found
......@@ -74,29 +74,24 @@ public:
{
typename LocalEnergy::ReturnType aEnergy(0.0);
const auto& gridView = this->basis().getGridView();
for (const auto& e : elements(this->basis().getGridView())) {
auto eIt = gridView.template begin<0>();
auto eEndIt = gridView.template end<0>();
for (; eIt != eEndIt; ++eIt) {
const auto& lfe = this->basis().getLocalFiniteElement(*eIt);
const auto& lfe = this->basis().getLocalFiniteElement(e);
// interpolate by local finite element to get the local coefficients
typename LocalEnergy::CoefficientVectorType localCoeff(lfe.localBasis().size());
LocalWrapper fiLocal(*displace,*eIt,0);
LocalWrapper fiLocal(*displace, e, 0);
std::vector<typename LocalWrapper::RangeType> interpolationValues;
for (int i=0; i<dim; i++) {
for (int i=0; i < dim; i++) {
fiLocal.setIndex(i);
lfe.localInterpolation().interpolate(fiLocal,interpolationValues);
for (size_t j=0; j<lfe.localBasis().size(); j++)
lfe.localInterpolation().interpolate(fiLocal, interpolationValues);
for (size_t j=0; j< lfe.localBasis().size(); j++)
localCoeff[j][i] = interpolationValues[j];
}
aEnergy += localEnergy_->energy(*eIt,lfe, localCoeff);
aEnergy += localEnergy_->energy(e,lfe, localCoeff);
}
ReturnType energy;
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment