anchor.diff:  !$omp end do OMP_SIMD
   !$omp end parallel
     end subroutine
  +
     !====== tetragonal(8) structure =========================!
     subroutine sym8_tetragonal(zutmp,zu_NB,zrho_l,zrhotmp1,zrhotmp2)
       use global_variables
  
positive.code: 
      zfac=1.0d0/4d0
  
  !$omp parallel private(tid)
  !$  tid=omp_get_thread_num()
      call reduce(tid,zfac,zutmp,zu_NB)
  
  ! 1.T_3
  !$omp do OMP_SIMD
  
negative.code:     ! wk(ik)=8.0,(ikx==iky >. wk(ik)=4.0)
      zfac=1.0d0/16d0
  
  #ifdef _OPENACC
      call reduce_acc(zfac,zutmp,zu_NB,zrhotmp(:,0))
  #endif
  
  !$omp parallel private(tid)
  !$  tid=omp_get_thread_num()
  
  #ifndef _OPENACC
      call reduce(tid,zfac,zutmp,zu_NB)
  #endif
  
  ! 1.T_3
  !$omp do OMP_SIMD
  