From bdc46a88287fec8c0acc02bb95f44bce9f51669a Mon Sep 17 00:00:00 2001
From: Pratheek Shanthraj
Date: Mon, 7 May 2012 08:12:13 +0000
Subject: [PATCH] using better measure (norm) of phase contrast for
preconditioning
---
code/DAMASK_spectral.f90 | 2 +-
1 file changed, 1 insertion(+), 1 deletion(-)
diff --git a/code/DAMASK_spectral.f90 b/code/DAMASK_spectral.f90
index 155fda2b4..95c38e5d8 100644
--- a/code/DAMASK_spectral.f90
+++ b/code/DAMASK_spectral.f90
@@ -760,7 +760,7 @@ C_ref = C * wgt
temperature(i,j,k),timeinc,ielem,1_pInt,sigma,dsde, &
P_real(i,j,k,1:3,1:3),dPdF)
CPFEM_mode = 2_pInt
- phase_cont(i,j,k) = maxval(dPdF)/maxval(C_ref)
+ phase_cont(i,j,k) = sqrt(sum(dPdF*dPdF)/sum(C_ref*C_ref))
C = C + dPdF
enddo; enddo; enddo
call debug_info()