Skip to content
Snippets Groups Projects
Commit c5d9b7f8 authored by Lars Bilke's avatar Lars Bilke
Browse files

Merge branch 'MCASAnisotropicBoundsAccessFix' into 'master'

[MatL/S] Fix a bug in MCASAnisotropic in 2D case.

See merge request ogs/ogs!3260
parents 33d7a7ae 359663b6
No related branches found
No related tags found
No related merge requests found
......@@ -96,7 +96,11 @@ cs.setEntryName("ShearFactor");
const int step = sig.size() + 1;
Hill[0] = cn;
Hill[step] = Hill[step * 2] = 1. / cn;
Hill[step * 3] = Hill[step * 4] = cs;
Hill[step * 3] = cs;
if (step == 7) // 3D case
{
Hill[step * 5] = cs;
}
// Compute initial elastic strain
const auto S = invert(D);
......@@ -146,7 +150,11 @@ cs.setEntryName("ShearFactor");
const int step = sig.size() + 1;
Hill[0] = cn;
Hill[step] = Hill[step * 2] = 1. / cn;
Hill[step * 3] = Hill[step * 4] = cs;
Hill[step * 3] = cs;
if (step == 7) // 3D case
{
Hill[step * 5] = cs;
}
if (F)
{
const auto sig_scaled = Hill * sig;
......
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