dPdF B_left_coef col major

phys-sim-book · May 14, 2024 · 7fbe672 · 7fbe672
1 parent d74d923
commit 7fbe672
Show file tree

Hide file tree

Showing 3 changed files with 12 additions and 12 deletions.
diff --git a/6_inv_free/NeoHookeanEnergy.py b/6_inv_free/NeoHookeanEnergy.py
@@ -72,10 +72,10 @@ def d2Psi_div_dF2(F, mu, lam):
  rs = s * 2 + r
  dP_div_dF[ij, rs] = M[0, 0] * U[i, 0] * VT[0, j] * U[r, 0] * VT[0, s] \
  + M[0, 3] * U[i, 0] * VT[0, j] * U[r, 1] * VT[1, s] \
- + M[1, 1] * U[i, 0] * VT[1, j] * U[r, 0] * VT[1, s] \
- + M[1, 2] * U[i, 0] * VT[1, j] * U[r, 1] * VT[0, s] \
- + M[2, 1] * U[i, 1] * VT[0, j] * U[r, 0] * VT[1, s] \
- + M[2, 2] * U[i, 1] * VT[0, j] * U[r, 1] * VT[0, s] \
+ + M[1, 1] * U[i, 1] * VT[0, j] * U[r, 1] * VT[0, s] \
+ + M[1, 2] * U[i, 1] * VT[0, j] * U[r, 0] * VT[1, s] \
+ + M[2, 1] * U[i, 0] * VT[1, j] * U[r, 1] * VT[0, s] \
+ + M[2, 2] * U[i, 0] * VT[1, j] * U[r, 0] * VT[1, s] \
  + M[3, 0] * U[i, 1] * VT[1, j] * U[r, 0] * VT[0, s] \
  + M[3, 3] * U[i, 1] * VT[1, j] * U[r, 1] * VT[1, s]
  return dP_div_dF

diff --git a/7_self_contact/NeoHookeanEnergy.py b/7_self_contact/NeoHookeanEnergy.py
@@ -71,10 +71,10 @@ def d2Psi_div_dF2(F, mu, lam):
  rs = s * 2 + r
  dP_div_dF[ij, rs] = M[0, 0] * U[i, 0] * VT[0, j] * U[r, 0] * VT[0, s] \
  + M[0, 3] * U[i, 0] * VT[0, j] * U[r, 1] * VT[1, s] \
- + M[1, 1] * U[i, 0] * VT[1, j] * U[r, 0] * VT[1, s] \
- + M[1, 2] * U[i, 0] * VT[1, j] * U[r, 1] * VT[0, s] \
- + M[2, 1] * U[i, 1] * VT[0, j] * U[r, 0] * VT[1, s] \
- + M[2, 2] * U[i, 1] * VT[0, j] * U[r, 1] * VT[0, s] \
+ + M[1, 1] * U[i, 1] * VT[0, j] * U[r, 1] * VT[0, s] \
+ + M[1, 2] * U[i, 1] * VT[0, j] * U[r, 0] * VT[1, s] \
+ + M[2, 1] * U[i, 0] * VT[1, j] * U[r, 1] * VT[0, s] \
+ + M[2, 2] * U[i, 0] * VT[1, j] * U[r, 0] * VT[1, s] \
  + M[3, 0] * U[i, 1] * VT[1, j] * U[r, 0] * VT[0, s] \
  + M[3, 3] * U[i, 1] * VT[1, j] * U[r, 1] * VT[1, s]
  return dP_div_dF

diff --git a/8_self_friction/NeoHookeanEnergy.py b/8_self_friction/NeoHookeanEnergy.py
@@ -71,10 +71,10 @@ def d2Psi_div_dF2(F, mu, lam):
  rs = s * 2 + r
  dP_div_dF[ij, rs] = M[0, 0] * U[i, 0] * VT[0, j] * U[r, 0] * VT[0, s] \
  + M[0, 3] * U[i, 0] * VT[0, j] * U[r, 1] * VT[1, s] \
- + M[1, 1] * U[i, 0] * VT[1, j] * U[r, 0] * VT[1, s] \
- + M[1, 2] * U[i, 0] * VT[1, j] * U[r, 1] * VT[0, s] \
- + M[2, 1] * U[i, 1] * VT[0, j] * U[r, 0] * VT[1, s] \
- + M[2, 2] * U[i, 1] * VT[0, j] * U[r, 1] * VT[0, s] \
+ + M[1, 1] * U[i, 1] * VT[0, j] * U[r, 1] * VT[0, s] \
+ + M[1, 2] * U[i, 1] * VT[0, j] * U[r, 0] * VT[1, s] \
+ + M[2, 1] * U[i, 0] * VT[1, j] * U[r, 1] * VT[0, s] \
+ + M[2, 2] * U[i, 0] * VT[1, j] * U[r, 0] * VT[1, s] \
  + M[3, 0] * U[i, 1] * VT[1, j] * U[r, 0] * VT[0, s] \
  + M[3, 3] * U[i, 1] * VT[1, j] * U[r, 1] * VT[1, s]
  return dP_div_dF