Skip to content
Snippets Groups Projects
Commit f7cd0f9f authored by lisa_julia.nebel_at_tu-dresden.de's avatar lisa_julia.nebel_at_tu-dresden.de
Browse files

Correct mistake in mooneyrivlindensity

parent 5147054d
No related branches found
No related tags found
1 merge request!54Correct mistake in mooneyrivlindensity
...@@ -130,7 +130,7 @@ public: ...@@ -130,7 +130,7 @@ public:
{ {
//mooneyrivlin_energy is "square" //mooneyrivlin_energy is "square"
field_type detFMinus1 = detF - 1; field_type detFMinus1 = detF - 1;
return strainEnergy + mooneyrivlin_k* detFMinus1 * detFMinus1; return strainEnergy + 0.5 * mooneyrivlin_k* detFMinus1 * detFMinus1;
} }
} }
} }
......
...@@ -149,10 +149,10 @@ int main (int argc, char *argv[]) ...@@ -149,10 +149,10 @@ int main (int argc, char *argv[])
parameters["mooneyrivlin_k"] = "75e+6"; parameters["mooneyrivlin_k"] = "75e+6";
parameters["mooneyrivlin_energy"] = "square"; parameters["mooneyrivlin_energy"] = "square";
double expectedEnergy = 177758377; double expectedEnergy = 167636683;
double expectedGradientTwoNorm = 2.13791422e+09; double expectedGradientTwoNorm = 2.10685704e+09;
double expectedGradientInfinityNorm = 599188506; double expectedGradientInfinityNorm = 589695526;
double expectedMatrixFrobeniusNorm = 1.62529884e+11; double expectedMatrixFrobeniusNorm = 1.61941167e+11;
int testSquare = assembleAndCompare(basis, parameters, x, expectedEnergy, expectedGradientTwoNorm, expectedGradientInfinityNorm, expectedMatrixFrobeniusNorm); int testSquare = assembleAndCompare(basis, parameters, x, expectedEnergy, expectedGradientTwoNorm, expectedGradientInfinityNorm, expectedMatrixFrobeniusNorm);
parameters["mooneyrivlin_energy"] = "log"; parameters["mooneyrivlin_energy"] = "log";
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment