← Back to branch summary

~jose-soler/siesta/unfolding

« back to all changes in this revision

Viewing changes to Src/m_ts_weight.F90

  • Committer: Jose M Soler
  • Date: 2018-10-29 18:09:36 UTC
  • mfrom: (718.1.26 trunk)
  • Revision ID: jose.soler@uam.es-20181029180936-7p8z2sbde0mkvxoz
Merged with trunk-744

modified:
  Docs/Contributors.txt
  Docs/siesta.bib
  Docs/siesta.tex
  Docs/tbtrans.tex
  Docs/tex/fdf.tex
  Docs/tex/setup.tex
  RELEASE_NOTES
  Src/final_H_f_stress.F
  Src/grdsam.F
  Src/m_ts_options.F90
  Src/m_ts_tri_init.F90
  Src/m_ts_weight.F90
  Src/save_density_matrix.F90
  Src/state_init.F
  Util/TS/TBtrans/Makefile
  Util/TS/TBtrans/m_tbt_tri_init.F90

Show diffs side-by-side

added added

removed removed

Lines of Context:
Src/m_ts_weight.F90
589
589
            end do
590
590
 
591
591
            ! Calculate error from estimated density
592
 
            e_f  = e_f  - DM(ind,1)
 
592
            e_f = e_f - DM(ind,1)
593
593
            Ee_f = Ee_f - EDM(ind,1)
594
594
            do mu_i = 2 , N_mu
595
595
              tmp = DM(ind,mu_i) - DM(ind,1)
636
636
            end do
637
637
 
638
638
            ! Calculate error from estimated density
639
 
            e_f  = e_f  - DM(ind,1)
 
639
            e_f = e_f - DM(ind,1)
640
640
            do mu_i = 2 , N_mu
641
641
              tmp = DM(ind,mu_i) - DM(ind,1)
642
642
              if ( abs(tmp) > abs(e_f) ) e_f = tmp