From 74791a6686bef9db2daf35ba5ab9ef11df465342 Mon Sep 17 00:00:00 2001 From: Martin Diehl Date: Fri, 2 Aug 2013 11:48:24 +0000 Subject: [PATCH] corrected invalid error for mismatch_p --- code/DAMASK_spectral_solverAL.f90 | 2 +- code/DAMASK_spectral_solverPolarisation.f90 | 5 +---- 2 files changed, 2 insertions(+), 5 deletions(-) diff --git a/code/DAMASK_spectral_solverAL.f90 b/code/DAMASK_spectral_solverAL.f90 index 708039c3d..30bd7b101 100644 --- a/code/DAMASK_spectral_solverAL.f90 +++ b/code/DAMASK_spectral_solverAL.f90 @@ -590,10 +590,10 @@ subroutine AL_formResidual(in,x_scal,f_scal,dummy,ierr) ! constructing residual err_p = 0.0_pReal do k = 1_pInt, grid(3); do j = 1_pInt, grid(2); do i = 1_pInt, grid(1) - err_p = err_p + sum((math_mul3333xx33(C_scale,F_lambda(1:3,1:3,i,j,k) - math_I3) - residual_F(1:3,1:3,i,j,k))**2.0_pReal) residual_F(1:3,1:3,i,j,k) = math_I3 + math_mul3333xx33(S_scale,residual_F(1:3,1:3,i,j,k)) - & F_lambda(1:3,1:3,i,j,k) & + residual_F_lambda(1:3,1:3,i,j,k) + err_p = err_p + sum((math_mul3333xx33(C_scale,residual_F(1:3,1:3,i,j,k) - residual_F_lambda(1:3,1:3,i,j,k)))**2.0_pReal) enddo; enddo; enddo !-------------------------------------------------------------------------------------------------- diff --git a/code/DAMASK_spectral_solverPolarisation.f90 b/code/DAMASK_spectral_solverPolarisation.f90 index 0eb4ceea5..70017901a 100644 --- a/code/DAMASK_spectral_solverPolarisation.f90 +++ b/code/DAMASK_spectral_solverPolarisation.f90 @@ -593,14 +593,11 @@ subroutine Polarisation_formResidual(in,x_scal,f_scal,dummy,ierr) err_p = 0.0_pReal do k = 1_pInt, grid(3); do j = 1_pInt, grid(2); do i = 1_pInt, grid(1) e = e + 1_pInt - err_p = err_p + sum((math_mul3333xx33(C_scale,F_tau(1:3,1:3,i,j,k) -& - F(1:3,1:3,i,j,k) - residual_F_tau(1:3,1:3,i,j,k)/polarBeta -& - math_I3) - & - residual_F(1:3,1:3,i,j,k))**2.0_pReal) residual_F(1:3,1:3,i,j,k) = math_mul3333xx33(math_invSym3333(materialpoint_dPdF(:,:,:,:,1,e) + C_scale), & residual_F(1:3,1:3,i,j,k) - & math_mul3333xx33(C_scale,F_tau(1:3,1:3,i,j,k) - F(1:3,1:3,i,j,k) - math_I3)) & + residual_F_tau(1:3,1:3,i,j,k) + err_p = err_p + sum((math_mul3333xx33(C_scale,residual_F(1:3,1:3,i,j,k) - residual_F_tau(1:3,1:3,i,j,k)))**2.0_pReal) enddo; enddo; enddo !--------------------------------------------------------------------------------------------------