← Back to branch summary

~nickpapior/siesta/tddft-work

« back to all changes in this revision

Viewing changes to Util/TS/TBtrans/m_tbt_tri_init.F90

  • Committer: Rafi Ullah
  • Date: 2017-08-30 14:09:10 UTC
  • mfrom: (611.1.19 trunk)
  • Revision ID: rraffiu@gmail.com-20170830140910-bhu0osuh4d59wn8e
Merged with trunk-630

Show diffs side-by-side

added added

removed removed

Lines of Context:
Util/TS/TBtrans/m_tbt_tri_init.F90
125
125
    use create_Sparsity_Union
126
126
 
127
127
    use m_ts_rgn2trimat
128
 
    use m_ts_tri_common, only : nnzs_tri_i8b, ts_pivot_tri_sort_El
 
128
    use m_ts_tri_common, only : ts_pivot_tri_sort_El
 
129
    use m_ts_tri_common, only : nnzs_tri_i8b, nnzs_tri_dp
129
130
    use m_ts_electype
130
131
#ifdef TRANSIESTA_DEBUG
131
132
    use m_ts_debug
259
260
 
260
261
    if ( IONode ) then
261
262
       ! Print out stuff
262
 
       call rgn_print(DevTri, seq_max = 10 , collapse = .false.)
 
263
       call rgn_print(DevTri, seq_max = 8 , repeat = .true.)
263
264
       ! Print out memory estimate
264
265
       els = nnzs_tri_i8b(DevTri%n,DevTri%r)
 
266
       ! check if there are overflows
 
267
       if ( els < int(nnzs_tri_dp(DevTri%n, DevTri%r), i8b) ) then
 
268
          call die('tbtrans: Memory consumption is too large, try &
 
269
               &another pivoting scheme.')
 
270
       end if
265
271
       write(*,'(a,i0)') 'tbtrans: Matrix elements in BTD: ', els
266
272
 
267
273
       write(*,'(/,a)') 'tbtrans: Electrodes tri-diagonal matrices'
268
274
       do i = 1 , N_Elec
269
 
          call rgn_print(ElTri(i), seq_max = 10 , collapse = .false.)
 
275
          call rgn_print(ElTri(i), seq_max = 8 , repeat = .true.)
270
276
       end do
271
277
    end if
272
278