/home/uryu/Cocal_svn_checkinout/Cocal/code/Subroutine/error_surface.f90

Go to the documentation of this file.

subroutine error_surface(pot,pot_bak,error,flag)
  use phys_constant,  only : long
  use grid_parameter, only : nrg, ntg, npg, eps
  implicit none
  real(long), pointer :: pot(:,:), pot_bak(:,:)
  real(long), intent(out) :: error
  integer,    intent(out) :: flag
  real(long) :: error_pot = 0.0d0, small = 1.0d-14
  integer    :: itg, ipg
!
  error = 0.0d0
  flag = 0
  do itg = 0, ntg
    do ipg = 0, npg
      error_pot = 2.0d0*dabs(pot(itg,ipg) -     pot_bak(itg,ipg)) &
    &                 /(dabs(pot(itg,ipg))+dabs(pot_bak(itg,ipg))+small)
      if (error_pot > eps) flag = 1
      error = dmax1(error,error_pot)
    end do
  end do
end subroutine error_surface

---