Skip to content
Snippets Groups Projects
Commit bc039b90 authored by Christian Engwer's avatar Christian Engwer
Browse files

more paranoia

[[Imported from SVN: r3337]]
parent 2907db6e
No related branches found
No related tags found
No related merge requests found
......@@ -192,13 +192,14 @@ namespace Dune
virtual RT derivativelocal (int comp, const Dune::FieldVector<int,n>& d,
const Entity& e, const Dune::FieldVector<DT,n>& xi) const
{
int dir;
int dir=-1;
int order=0;
for (int i=0; i<n; i++)
{
order += d[i];
if (d[i]>0) dir=i;
}
assert(dir != -1);
if (order!=1) DUNE_THROW(GridError,"can only evaluate one derivative");
RT value=0;
......
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