Skip to content
Snippets Groups Projects
Commit 6e3f38e4 authored by Oliver Sander's avatar Oliver Sander Committed by sander
Browse files

Don't abort when no Neumann function is given 

[[Imported from SVN: r9389]]
parent 3c4597f8
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
− 3
View file @ 6e3f38e4
...@@ -551,7 +551,8 @@ energy(const Entity& element, ...@@ -551,7 +551,8 @@ energy(const Entity& element,
// Assemble boundary contributions // Assemble boundary contributions
////////////////////////////////////////////////////////////////////////////// //////////////////////////////////////////////////////////////////////////////
assert(neumannFunction_); if (not neumannFunction_)
return energy;
for (typename Entity::LeafIntersectionIterator it = element.ileafbegin(); it != element.ileafend(); ++it) { for (typename Entity::LeafIntersectionIterator it = element.ileafbegin(); it != element.ileafend(); ++it) {
...@@ -992,8 +993,6 @@ assembleGradient(const Entity& element, ...@@ -992,8 +993,6 @@ assembleGradient(const Entity& element,
// Assemble boundary contributions // Assemble boundary contributions
////////////////////////////////////////////////////////////////////////////// //////////////////////////////////////////////////////////////////////////////
assert(neumannFunction_);
for (typename Entity::LeafIntersectionIterator it = element.ileafbegin(); it != element.ileafend(); ++it) { for (typename Entity::LeafIntersectionIterator it = element.ileafbegin(); it != element.ileafend(); ++it) {
if (not neumannBoundary_ or not neumannBoundary_->contains(*it)) if (not neumannBoundary_ or not neumannBoundary_->contains(*it))
......
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