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

fix index bug in mixed third derivative 

[[Imported from SVN: r9200]]
parent 0bac8d6c
No related branches found
No related tags found
No related merge requests found
Hide whitespace changes
Inline Side-by-side
dune/gfe/hyperbolichalfspacepoint.hh
+ 1
− 1
View file @ 3f07ec39
...@@ -503,7 +503,7 @@ public: ...@@ -503,7 +503,7 @@ public:
for (size_t j=0; j<N; j++) for (size_t j=0; j<N; j++)
for (size_t k=0; k<N; k++) for (size_t k=0; k<N; k++)
result[i][j][k] = alphaPrimePrimePrime * dFdp[i] * dFdq[j] * dFdq[k] result[i][j][k] = alphaPrimePrimePrime * dFdp[i] * dFdq[j] * dFdq[k]
+ alphaPrimePrime * (dFdpdq[i][j] * dFdq[k] + dFdpdq[i][k] * dFdq[j] + dFdqdq[j][k] * dFdq[i]) + alphaPrimePrime * (dFdpdq[i][j] * dFdq[k] + dFdpdq[i][k] * dFdq[j] + dFdqdq[j][k] * dFdp[i])
+ alphaPrime * dFdpdqdq[i][j][k]; + alphaPrime * dFdpdqdq[i][j][k];
return result; return result;
......
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