added output of current relative tolerance when reporting relative errors

This commit is contained in:
Philip Eisenlohr 2012-11-28 18:44:51 +00:00
parent 732340a4d7
commit 6bb3a475ce
1 changed files with 9 additions and 6 deletions

View File

@ -475,12 +475,15 @@ subroutine AL_converged(snes_local,it,xnorm,snorm,fnorm,reason,dummy,ierr)
reason = 0 reason = 0
endif endif
write(6,'(a,es14.7)') 'error stress BC = ', & write(6,'(a,f12.7,1x,1a,1x,es9.3)') 'error stress BC = ', &
err_stress/min(maxval(abs(P_av))*err_stress_tolrel,err_stress_tolabs) err_stress/min(maxval(abs(P_av))*err_stress_tolrel,err_stress_tolabs),&
write(6,'(a,es14.7)') 'error F = ',& '@',min(maxval(abs(P_av))*err_stress_tolrel,err_stress_tolabs)
err_f/sqrt(sum((F_aim-math_I3)*(F_aim-math_I3)))/err_f_tol write(6,'(a,f12.7,1x,1a,1x,es9.3)') 'error F = ',&
write(6,'(a,es14.7)') 'error P = ', & err_f/sqrt(sum((F_aim-math_I3)*(F_aim-math_I3)))/err_f_tol,&
err_p/sqrt(sum((F_aim-math_I3)*(F_aim-math_I3)))/err_p_tol '@',err_f_tol
write(6,'(a,f12.7,1x,1a,1x,es9.3)') 'error P = ', &
err_p/sqrt(sum((F_aim-math_I3)*(F_aim-math_I3)))/err_p_tol,&
'@',err_p_tol
write(6,'(/,a)') '==========================================================================' write(6,'(/,a)') '=========================================================================='
end subroutine AL_converged end subroutine AL_converged