Fix all ph_base/ tests (kpoints case fixed)

LR_Modules/response_kernels.f90

@@ -170,12 +170,15 @@ SUBROUTINE sternheimer_kernel(first_iter, time_reversed, npert, lrdvpsi, iudvpsi
       !
       IF (nksq > 1 .OR. (noncolin .AND. domag)) THEN
          IF (lgamma) THEN
+            !civn: in this case evq is a pointer to evc
             CALL get_buffer(evc, lrwfc, iuwfc, ikmk)
             !$acc update device(evc)
          ELSE
+            !civn: in this case evq is allocated separately and needs to be updated on device
             CALL get_buffer(evc, lrwfc, iuwfc, ikmk)
             !$acc update device(evc)
             CALL get_buffer(evq, lrwfc, iuwfc, ikmkmq)
+            !$acc update device(evq)
          ENDIF
       ENDIF
       !

PHonon/PH/allocate_phq.f90

@@ -63,6 +63,7 @@ subroutine allocate_phq
      !  q!=0 : evq is allocated and calculated at point k+q
      !
      allocate (evq ( npwx*npol , nbnd))
+     !$acc enter data create(evq)
   endif
   !
   allocate (dvpsi ( npwx*npol , nbnd))

PHonon/PH/deallocate_phq.f90

@@ -54,6 +54,7 @@ subroutine deallocate_phq
   if (lgamma) then
      if(associated(evq)) nullify(evq)
   else
+     !$acc exit data delete(evq)
      if(associated(evq)) deallocate(evq)   !why not if allocated? 
   end if
 

PHonon/PH/solve_linter.f90

@@ -267,6 +267,7 @@ SUBROUTINE solve_linter (irr, imode0, npe, drhoscf)
         !
         IF (nksq > 1 .OR. nsolv == 2) THEN
            CALL get_buffer(evc, lrwfc, iuwfc, ikmk)
+           !$acc update device(evc)
         ENDIF
         !
         DO ipert = 1, npe