Skip to content
Snippets Groups Projects
Commit b43abd8d authored by Jö Fahlke's avatar Jö Fahlke
Browse files

lambda_volume: don't use view for r

parent 6a8f2c46
No related branches found
No related tags found
No related merge requests found
......@@ -399,9 +399,8 @@ namespace PPS {
template<typename EG, typename LFSVC>
void assembleVVolume(const EG & eg, const LFSVC & lfsv_cache)
{
rl_view.setWeight(1.0);
Dune::PDELab::LocalAssemblerCallSwitch<LOP,LOP::doLambdaVolume>::
lambda_volume(lop,eg,lfsv_cache.localFunctionSpace(),rl_view);
lambda_volume(lop,eg,lfsv_cache.localFunctionSpace(),rl);
}
template<typename IG, typename LFSUC, typename LFSVC>
......
......@@ -78,7 +78,7 @@ public:
geo.integrationElement(ip.position());
auto f=param.f(eg.entity(),ip.position());
for (size_t i=0; i<lfsv.size(); i++)
r.accumulate(lfsv,i,-f*phihat[i]*factor);
r(lfsv,i) -= f*phihat[i]*factor;
}
}
......
......@@ -74,7 +74,7 @@ public:
geo.integrationElement(ip.position());
auto f=param.f(eg.entity(),ip.position());
for (size_t i=0; i<lfsv.size(); i++)
r.accumulate(lfsv,i,-f*phihat[i]*factor);
r(lfsv,i) -= f*phihat[i]*factor;
}
}
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment