corr rW (use 0.5 prefactor) in isobrittle and related test upd

This commit is contained in:
Yi Hu 2023-01-23 15:27:44 +01:00
parent 6ee31e2ab3
commit dd74fa2d1d
2 changed files with 2 additions and 2 deletions

@ -1 +1 @@
Subproject commit 547bfe56d67cba358a7bb9582f2f7c0e344befd8
Subproject commit ae2ecb3623b5061d2880e8d9fed8c6fc342e265e

View File

@ -114,7 +114,7 @@ module subroutine isobrittle_deltaState(C, Fe, ph,en)
associate(prm => param(ph), stt => state(ph), dlt => deltaState(ph))
r_W = (2.0_pReal*dot_product(epsilon,matmul(C,epsilon)))/prm%W_crit
r_W = (0.5_pReal*dot_product(epsilon,matmul(C,epsilon)))/prm%W_crit
dlt%r_W(en) = merge(r_W - stt%r_W(en), 0.0_pReal, r_W > stt%r_W(en))
end associate