← Back to branch summary

~siesta-ts/siesta/trunk_ts_soc

« back to all changes in this revision

Viewing changes to Src/m_ts_trik.F90

  • Committer: Nils Wittemeier
  • Date: 2019-02-14 07:45:07 UTC
  • mfrom: (746.1.15 trunk)
  • Revision ID: nils@4wittemeier.de-20190214074507-1mvzbmj9kw19gllr
Merged trunk 761

Show diffs side-by-side

added added

removed removed

Lines of Context:
Src/m_ts_trik.F90
667
667
 
668
668
    call rgn_delete(pvt)
669
669
 
 
670
    ! Nullify external pointers
 
671
    do iEl = 1, N_Elec
 
672
      nullify(Elecs(iEl)%Sigma)
 
673
    end do
 
674
 
670
675
#ifdef TRANSIESTA_DEBUG
671
676
    call write_debug( 'POS transiesta mem' )
672
677
#endif
867
872
       call init_val(spuDM)
868
873
       cE%exist = .true.
869
874
       cE%fake  = .true.
870
 
       cE%e = dcmplx(0._dp, 0._dp)
 
875
       cE%e = cmplx(0._dp, 0._dp,dp)
871
876
       cE%idx = 1
872
877
       cE%idx(1) = 0 ! Signals a fermi-contour
873
878
       if ( Node == Nodes - 1 ) cE%fake = .false.
891
896
          call invert_TriMat(zwork_tri,GF_tri,calc_parts)
892
897
       end if
893
898
 
894
 
       W = dcmplx(kw,0._dp)
 
899
       W = cmplx(kw,0._dp,dp)
895
900
       call add_DM( spuDM, W, spuEDM, W, &
896
901
            GF_tri, r_pvt, pvt, N_Elec, Elecs, DMidx=1)
897
902
 
941
946
 
942
947
    call rgn_delete(pvt)
943
948
 
 
949
    ! Nullify external pointers
 
950
    do iEl = 1, N_Elec
 
951
      nullify(Elecs(iEl)%Sigma)
 
952
    end do
 
953
 
944
954
  end subroutine ts_trik_Fermi
945
955
  
946
956
  subroutine add_DM(DM, DMfact, EDM, EDMfact, &
1133
1143
 
1134
1144
    ! Initialize
1135
1145
!$OMP workshare
1136
 
    GFinv(:) = dcmplx(0._dp,0._dp)
 
1146
    GFinv(:) = cmplx(0._dp,0._dp,dp)
1137
1147
!$OMP end workshare
1138
1148
 
1139
1149
    ! We will only loop in the central region