used dp instead of E to initialize arrays

Author: felipeZ

src/davidson.f90

@@ -513,8 +513,8 @@ end function fun_A_gemv
     real(dp), dimension(dim,1) :: tmp_array
 
     do ii = 1,dim
-      tmp_array = 0E0
-      tmp_array(ii,1) = 1.0
+      tmp_array = 0.0_dp
+      tmp_array(ii,1) = 1.0_dp
       tmp_array = fun_A_gemv(tmp_array)
       out(ii) = tmp_array(ii,1)
     end do