Skip to content
Snippets Groups Projects
Commit e42c8408 authored by Oliver Sander's avatar Oliver Sander Committed by sander@FU-BERLIN.DE
Browse files

enable all three parts of the gradient 

[[Imported from SVN: r8570]]
parent 5c649bd2
No related branches found
No related tags found
No related merge requests found
Hide whitespace changes
Inline Side-by-side
dune/gfe/cosseratenergystiffness.hh
+ 2
− 2
View file @ e42c8408
......@@ -816,7 +816,7 @@ assembleGradient(const Entity& element,
if (gridDim==2) {
typename TargetSpace::EmbeddedTangentVector tmp(0);
longQuadraticMembraneEnergyGradient(tmp,R,dR_dv,derivative,derOfGradientWRTCoefficient,U);
//embeddedLocalGradient[i].axpy(weight * thickness_, tmp);
embeddedLocalGradient[i].axpy(weight * thickness_, tmp);
tmp = 0;
curvatureEnergyGradient(tmp,R,DR,dDR_dv);
......@@ -824,7 +824,7 @@ assembleGradient(const Entity& element,
tmp = 0;
bendingEnergyGradient(tmp,R,dR_dv,DR,dDR_dv);
//embeddedLocalGradient[i].axpy(weight * std::pow(thickness_,3) / 12.0, tmp);
embeddedLocalGradient[i].axpy(weight * std::pow(thickness_,3) / 12.0, tmp);
} else if (gridDim==3) {
assert(gridDim==2); // 3d not implemented yet
// energy += weight * quadraticMembraneEnergyGradient(U);
......
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment