From dd74fa2d1d8d32ae7b84f309a2d72b2f899e40ce Mon Sep 17 00:00:00 2001
From: Yi Hu <y.hu@mpie.de>
Date: Mon, 23 Jan 2023 15:27:44 +0100
Subject: [PATCH] corr rW (use 0.5 prefactor) in isobrittle and related test
 upd

---
 PRIVATE                         | 2 +-
 src/phase_damage_isobrittle.f90 | 2 +-
 2 files changed, 2 insertions(+), 2 deletions(-)

diff --git a/PRIVATE b/PRIVATE
index 547bfe56d..ae2ecb362 160000
--- a/PRIVATE
+++ b/PRIVATE
@@ -1 +1 @@
-Subproject commit 547bfe56d67cba358a7bb9582f2f7c0e344befd8
+Subproject commit ae2ecb3623b5061d2880e8d9fed8c6fc342e265e
diff --git a/src/phase_damage_isobrittle.f90 b/src/phase_damage_isobrittle.f90
index cfce3a31c..e9b62a03a 100644
--- a/src/phase_damage_isobrittle.f90
+++ b/src/phase_damage_isobrittle.f90
@@ -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