/dports/science/siesta/siesta-4.1.5/Src/ |
H A D | printmatrix.F | 30 subroutine printPermMatrix( nuo, ia, ja, perm, invp, argument 34 integer :: nuo, ia(*), ja(0:nuo), perm(nuo), invp(nuo) 40 nnz = ja(nuo)-1 43 call re_alloc( xadj, 0, nuo, name='xadj', 45 call re_alloc( auxx, 1, nuo, name='auxx', 49 do io= 1, nuo 59 xadj(nuo) = ii 60 call printMatrix( nuo, nuo, adj, xadj, filename, title )
|
H A D | diagg.F | 111 & listh(maxnh), numh(nuo), listhptr(nuo), 112 & listd(maxnd), numd(nuo), listdptr(nuo) 120 real(dp) Haux(nuotot,nuo), Saux(nuotot,nuo) 154 do io = 1,nuo 180 do io = 1,nuo 207 call writew(nuotot,nuo,1,k,ispin, 248 if (nuo.gt.0) then 250 nd = listdptr(nuo) + numd(nuo) 263 nuo_LOC(0) = nuo 277 do io= 1, nuo [all …]
|
H A D | diagpol.f | 8 subroutine diagpol( ispin, nspin, nuo, no, nuotot, argument 58 . maxnh, nuotot, no, nspin, nuo, indxuo(no), listh(maxnh), 59 . listhptr(nuo), numh(nuo), ng 63 . xij(3,*), psi(ng,nuotot,nuo), Haux(ng,nuotot,nuo), 64 . Saux(ng,nuotot,nuo) 76 do iuo = 1,nuo 98 call cdiag( Haux, Saux, nuotot, nuo, nuotot, eo, psi, 101 call rdiag( Haux, Saux, nuotot, nuo, nuotot, eo, psi, 111 do iuo = 1,nuo 133 call cdiag( Haux, Saux, nuotot, nuo, nuotot, eo, psi, [all …]
|
H A D | pdosg.F | 8 subroutine pdosg( nspin, nuo, no, maxspn, maxnh, argument 66 . nspin, nuo, no, maxspn, maxnh, 70 . numh(nuo), listhptr(nuo), listh(maxnh), 75 . haux(nuotot,nuo), saux(nuotot,nuo), psi(nuotot,nuo), 101 do iuo = 1, nuo 108 do iuo = 1, nuo 129 do iuo = 1, nuo 135 do iuo = 1, nuo 149 do iuo = 1, nuo 187 do 170 iband = 1, nuo [all …]
|
H A D | diagk.F | 122 . listd(maxnd), numd(nuo), listhptr(nuo), 131 real(dp) Dk(2,nuotot,nuo), Ek(2,nuotot,nuo), 132 . Haux(2,nuotot,nuo), Saux(2,nuotot,nuo) 173 do iuo = 1,nuo 220 nd = listdptr(nuo) + numd(nuo) 242 if (nuo.gt.0) then 280 do iuo = 1,nuo 311 do iuo = 1,nuo 358 do iuo = 1,nuo 390 do iuo = 1,nuo [all …]
|
H A D | diagk_file.F | 12 subroutine diagk_file( nspin, nuo, no, maxspn, argument 119 integer indxuo(no), listh(maxnh), numh(nuo), 120 . listd(maxnd), numd(nuo), listhptr(nuo), 121 . listdptr(nuo) 128 real(dp) Dk(2,nuotot,nuo), Ek(2,nuotot,nuo), 129 . Haux(2,nuotot,nuo), Saux(2,nuotot,nuo), 130 . psi(2,nuotot,nuo) 183 call write_wfs_netcdf(nuotot,nuo,ik,ispin, psi, 218 nd = listdptr(nuo) + numd(nuo) 267 do iuo = 1,nuo [all …]
|
H A D | pdosk.F | 71 . nspin, nuo, no, maxspn, maxnh, NK, 75 . numh(nuo), listhptr(nuo), listh(maxnh), 81 . Haux(2,nuotot,nuo), Saux(2,nuotot,nuo), psi(2,nuotot,nuo), 113 do iuo = 1,nuo 122 do iuo = 1, nuo 152 do iuo = 1,nuo 160 do iuo = 1, nuo 182 do iuo = 1, nuo 189 do iuo = 1,nuo 234 do 170 iband = 1, nuo [all …]
|
H A D | pdos3g.F | 8 subroutine pdos3g( nuo, no, maxuo, maxnh, argument 66 . nuo, no, maxuo, maxnh, 70 . numh(nuo), listhptr(nuo), listh(maxnh), indxuo(no) 75 complex(dp), target :: psi(2,nuotot,2*nuo) 76 complex(dp) Haux(2,nuotot,2,nuo), Saux(2,nuotot,2,nuo) 99 call re_alloc(Spr, 1, nuotot, 1, nuo, name='Spr', 102 do io = 1,nuo 124 do io = 1,nuo 145 do io = 1, nuo 179 do iband = 1, nuo*2 [all …]
|
H A D | pdos2g.F | 8 subroutine pdos2g( nuo, no, maxuo, maxnh, argument 66 . nuo, no, maxuo, maxnh, 70 . numh(nuo), listhptr(nuo), listh(maxnh), indxuo(no) 75 complex(dp), target :: psi(2,nuotot,2*nuo) 77 complex(dp) Haux(2,nuotot,2,nuo), Saux(2,nuotot,2,nuo) 99 call re_alloc(Spr, 1, nuotot, 1, nuo, name='Spr', 103 do io = 1,nuo 126 do io = 1,nuo 147 do io = 1, nuo 181 do iband = 1, nuo*2 [all …]
|
H A D | diagon.F | 160 call GetNodeOrbs(nuotot,Node,Nodes,nuo) 162 nuo = nuotot 172 nhs = nuotot * nuo 176 nhs = 2 * nuotot * nuo 177 npsi = 2 * nuotot * nuo 187 nhs = 2 * (2*nuotot) * (2*nuo) 188 npsi = 2 * (2*nuotot) * (2*nuo) 225 call diagkp( spin%H, nuo, no, nspin, maxnh, 248 call diagk_file(spin%H, nuo, no, nspin, 258 call diagk( spin%H, nuo, no, nspin, [all …]
|
H A D | diag2g.F | 8 subroutine diag2g( nuo, no, maxnh, maxnd, maxo, argument 96 integer no, nuo, nuotot, iscf, neigwanted 98 integer listh(maxnh), numh(nuo), listhptr(nuo) 99 integer listd(maxnd), numd(nuo), listdptr(nuo) 105 complex(dp), dimension(2,nuotot,2*nuo), target :: psi 106 complex(dp), dimension(2,nuotot,2,nuo) :: Haux, Saux 165 do io = 1,nuo 182 call cdiag(Haux,Saux,2*nuotot,2*nuo,2*nuotot,eo,psi, 191 do io = 1,nuo 213 call writew(nuotot,nuo,1,k,1, [all …]
|
H A D | diag3g.F | 8 subroutine diag3g( nuo, no, maxnh, maxnd, maxo, argument 91 integer no, nuo, nuotot, iscf, neigwanted 93 integer listh(maxnh), numh(nuo), listhptr(nuo) 94 integer listd(maxnd), numd(nuo), listdptr(nuo) 100 complex(dp), dimension(2,nuotot,2*nuo), target :: psi 102 complex(dp), dimension(2,nuotot,2,nuo) :: Haux, Saux 154 do io = 1,nuo 171 call cdiag(Haux,Saux,2*nuotot,2*nuo,2*nuotot,eo,psi, 179 do io = 1,nuo 207 call writew(nuotot,nuo,1,k,1, [all …]
|
H A D | detover.F | 10 . no, nuo, xij, maxnh, nuotot, nocc, 62 . nuo, maxnh, nuotot, no, nocc, 63 . listh(maxnh), listhptr(nuo), numh(nuo), indxuo(no), 68 . psiprev(2,nuo,nuotot), dk(3),detr, deti, 70 . psi(2,nuotot,nuo), kpoint(3) 91 call re_alloc( Aux2, 1, nuotot, 1, nuo, name='Aux2', 96 do iuo = 1,nuo 140 do iuo=1,nuo 164 do iuo=1,nuo
|
H A D | transition_rate.F | 12 . indxuo, no, nuo, nuotot, xij, 77 integer nuo, nuotot, maxnh, no, 78 . listh(maxnh), numh(nuo), listhptr(nuo), 83 . psi(ng,nuotot,nuo), kpoint(3), Aux(ng,nuotot,nuo), 84 . Aux2(ng,nuotot,nuo), ek(nuotot), temp, efermi, 85 . intraband(ng,nuo), smooth, wmin, wmax, ekloc(nuotot), 118 do iuo = 1,nuo 131 do iuo = 1,nuo 193 do iuo = 1,nuo 265 do iuo = 1,nuo
|
H A D | savepsi.F | 8 subroutine savepsi(psiprev, psi, nuo, nuotot, nocc) argument 34 integer nuo, nuotot, nocc 37 . psiprev(2,nuo,nuotot), 38 . psi(2,nuotot,nuo) 74 do juo = 1,nuo
|
H A D | diag2k.F | 8 subroutine diag2k( nuo, no, maxnh, maxnd, maxo, argument 108 integer maxuo, maxnd, maxnh, maxo, nk, no, nuo, nuotot, iscf 110 integer indxuo(no), numh(nuo), numd(nuo) 121 complex(dp), dimension(2,nuotot,2*nuo), target :: psi 194 do iuo = 1,nuo 272 do iuo = 1,nuo 296 call cdiag(Haux,Saux,2*nuotot,2*nuo,2*nuotot,caux,psi, 306 do iuo = 1,nuo 326 call cdiag(Haux,Saux,2*nuotot,2*nuo,2*nuotot,caux,psi, 406 do iuo = 1,nuo [all …]
|
H A D | diag3k.F | 8 subroutine diag3k( nuo, no, maxnh, maxnd, maxo, argument 106 integer maxuo, maxnd, maxnh, maxo, nk, no, nuo, nuotot, iscf 108 integer indxuo(no), numh(nuo), numd(nuo) 119 complex(dp), dimension(2,nuotot,2*nuo), target :: psi 192 do iuo = 1,nuo 271 do iuo = 1,nuo 295 call cdiag(Haux,Saux,2*nuotot,2*nuo,2*nuotot,caux,psi, 305 do iuo = 1,nuo 327 call cdiag(Haux,Saux,2*nuotot,2*nuo,2*nuotot,caux,psi, 393 do iuo = 1,nuo [all …]
|
H A D | pdos2k.F | 79 complex(dp), target :: psi(2,nuotot,2*nuo) 81 complex(dp) Haux(2,nuotot,2,nuo), Saux(2,nuotot,2,nuo) 112 call re_alloc(Spr, 1, nuotot, 1, nuo, name='Spr', 121 do iuo = 1,nuo 144 call cdiag( Haux, Saux, 2*nuotot, 2*nuo, 2*nuotot, 153 do iuo = 1,nuo 175 call cdiag( Haux, Saux, 2*nuotot, 2*nuo, 2*nuotot, 182 do iuo = 1,nuo 200 call MPI_AllReduce(nuo,maxnuo,1,MPI_integer,MPI_max, 225 do iband = 1, nuo*2 [all …]
|
H A D | pdos3k.F | 79 complex(dp), target :: psi(2,nuotot,2*nuo) 81 complex(dp) Haux(2,nuotot,2,nuo), Saux(2,nuotot,2,nuo) 112 call re_alloc(Spr, 1, nuotot, 1, nuo, name='Spr', 121 do iuo = 1,nuo 146 call cdiag( Haux, Saux, 2*nuotot, 2*nuo, 2*nuotot, 155 do iuo = 1,nuo 179 call cdiag( Haux, Saux, 2*nuotot, 2*nuo, 2*nuotot, 186 do iuo = 1,nuo 204 call MPI_AllReduce(nuo,maxnuo,1,MPI_integer,MPI_max, 229 do iband = 1, nuo*2 [all …]
|
H A D | optical.F | 8 subroutine optical( nua, na, xa, scell, ucell, nuotot, nuo, no, argument 83 integer maxnh, nuotot, nuo, no, nspin 84 integer indxuo(no), listh(maxnh), numh(nuo) 85 integer listhptr(nuo) 402 maxhs = nuotot*nuo 403 maxpsi = nuotot*nuo 405 maxhs = 2*nuotot*nuo 406 maxpsi = 2*nuotot*nuo 431 call re_alloc( intraband, 1, ngamma, 1, nuo, 483 call diagpol( ispin, nspin, nuo, no, nuotot, [all …]
|
H A D | ksv.f | 12 . nuo, no, nspin, qspin, maxnh, 85 integer maxkpol, maxnh, nuo, nuotot, no, nspin, na 86 integer indxuo(no), listh(maxnh), numh(nuo), nkpol 205 nhs = 2*nuotot*nuo 206 npsi = 2*nuotot*nuo 238 do iuo = 1,nuo 336 call phirphi(nua, na, nuo, no, scell, xa, rmaxo, 355 call diagpol( ispin, nspin, nuo, no, 369 do io = 1,nuo 384 call savepsi(psiprev,psi,nuo,nuotot,nocc(ispin)) [all …]
|
H A D | diag2kspiral.F90 | 105 integer :: maxnd, maxnh, maxo, nk, no, nuo, nuotot, iscf, neigwanted local 107 integer :: indxuo(no), listh(maxnh), numh(nuo), listd(maxnd), numd(nuo), & 108 listhptr(nuo), listdptr(nuo) 115 complex(dp), dimension(2,nuotot,2,nuo) :: Haux, Saux, Dk, Ek 116 complex(dp), dimension(2,nuotot,2*nuo) :: psi 138 nuo2 = nuo * 2 232 do iuo = 1,nuo 261 do iuo = 1,nuo 326 do iuo = 1,nuo
|
H A D | meshdscf.F | 287 subroutine matrixOtoM( maxnd, numd, listdptr, maxndl, nuo, argument 325 integer :: maxnd, maxndl, nspin, nuo, numd(nuo), 326 & listdptr(nuo) 405 do io = 1,nuo 419 subroutine matrixMtoO( maxnvl, maxnv, numVs, listVsptr, nuo, argument 463 & maxnv, maxnvl, nspin, nuo, numVs(nuo), 464 & listVsptr(nuo) 548 do i = 1,nuo 605 & maxnv, maxnvl, nuo, numVs(nuo), 606 & listVsptr(nuo) [all …]
|
/dports/science/siesta/siesta-4.1.5/Util/Denchar/Src/ |
H A D | readwavesx.f | 124 allocate(psi(1,nuo)) 126 allocate(psi(2,nuo)) 149 read(unit) (psi(1:,j), j=1,nuo) 151 RPSI(1:nuo,IK,IWF,ISPIN)=dble(psi(1,1:nuo)) 152 IPSI(1:nuo,IK,IWF,ISPIN)=0.d0 154 RPSI(1:nuo,IK,IWF,ISPIN)=dble(psi(1,1:nuo)) 155 IPSI(1:nuo,IK,IWF,ISPIN)=dble(psi(2,1:nuo))
|
/dports/science/siesta/siesta-4.1.5/Util/Vibra/Src/ |
H A D | outbands.f | 8 subroutine outbands(iunits, nspin, maxuo, nuo, maxk, nk, argument 42 integer iunits, maxk, maxuo, nk, nlines, nuo, nspin 90 do io = 1, nuo 102 write(iu,*) nuo, nspin, nk 110 . path, ((ek(io,ik,ispin)/eV,io=1,nuo),ispin=1,nspin)
|