! ---

BEGIN_PROVIDER [double precision, htilde_matrix_elmt_bi_ortho, (N_det,N_det)]

  BEGIN_DOC
  !
  ! htilde_matrix_elmt_bi_ortho(j,i) = <J| H^tilde |I> 
  !
  ! WARNING !!!!!!!!! IT IS NOT HERMITIAN !!!!!!!!!
  !
  END_DOC
 
  implicit none
  integer          :: i, j
  double precision :: t1, t2
  double precision :: htot
 
  PROVIDE N_int
  PROVIDE psi_det
  PROVIDE three_e_3_idx_term

  if(noL_standard) then
    PROVIDE noL_0e
    print*, "noL_0e =", noL_0e
    PROVIDE noL_1e
    PROVIDE noL_2e
  endif

  print *, ' PROVIDING htilde_matrix_elmt_bi_ortho ...'
  call wall_time(t1)

  call provide_all_three_ints_bi_ortho()

  i = 1
  j = 1
  call htilde_mu_mat_opt_bi_ortho_tot(psi_det(1,1,j), psi_det(1,1,i), N_int, htot)


  !$OMP PARALLEL            &
  !$OMP DEFAULT(NONE)       &
  !$OMP PRIVATE(i, j, htot) &
  !$OMP SHARED (N_det, psi_det, N_int, htilde_matrix_elmt_bi_ortho)
  !$OMP DO 
  do i = 1, N_det
    do j = 1, N_det
      ! < J |Htilde | I >
      call htilde_mu_mat_opt_bi_ortho_tot(psi_det(1,1,j), psi_det(1,1,i), N_int, htot)

      htilde_matrix_elmt_bi_ortho(j,i) = htot
    enddo
  enddo
  !$OMP END DO
  !$OMP END PARALLEL

  call wall_time(t2)
  print *, ' wall time for htilde_matrix_elmt_bi_ortho (min) =', (t2-t1)/60.d0 

END_PROVIDER 

! ---

BEGIN_PROVIDER [double precision, htilde_matrix_elmt_bi_ortho_tranp, (N_det,N_det)]

  implicit none
  integer ::i,j

  do i = 1, N_det
    do j = 1, N_det
      htilde_matrix_elmt_bi_ortho_tranp(j,i) = htilde_matrix_elmt_bi_ortho(i,j)
    enddo
  enddo
END_PROVIDER 

! ---