/dports/science/octopus/octopus-10.5/src/sternheimer/ |
H A D | kdotp_calc_inc.F90 | 30 integer :: ik, ist, ist2, idir1, idir2, pdim 58 do ik = sys%st%d%kpt%start, sys%st%d%kpt%end 74 eff_mass_inv(idir1, idir2, ist, ik) = eff_mass_inv(idir2, idir1, ist, ik) 87 … orth_mask(ist2) = .not. (abs(sys%st%eigenval(ist2, ik) - sys%st%eigenval(ist, ik)) < degen_thres) 143 integer :: ip, idim, ist, ist2, ik 169 if (abs(sys%st%eigenval(ist2, ik) - sys%st%eigenval(ist, ik)) < degen_thres) cycle 178 kdotp_lr%X(dl_psi)(ip, idim, ist, ik) = kdotp_lr%X(dl_psi)(ip, idim, ist, ik) + & 179 psi2(ip, idim)*mtxel/(sys%st%eigenval(ist, ik) - sys%st%eigenval(ist2, ik)) 182 kdotp_lr%X(dl_psi)(ip, idim, ist2, ik) = kdotp_lr%X(dl_psi)(ip, idim, ist2, ik) + & 183 psi1(ip, idim)*R_CONJ(-mtxel)/(sys%st%eigenval(ist2, ik) - sys%st%eigenval(ist, ik)) [all …]
|
/dports/science/dakota/dakota-6.13.0-release-public.src-UI/examples/compiled_interfaces/Scilab/Bridge/macros/ |
H A D | parse.sci | 57 ik = wkey(nk) 61 if keyfound(ik) == 1 then 74 if keytype(ik) == 'real' then 77 for ir = 1:1:abskeynumber(ik) 84 if keynumber(ik) < 0 then 95 keyfound(ik) = 1 99 elseif keytype(ik) == 'integer' then 117 if keynumber(ik) < 0 then 127 keyfound(ik) = 1 141 if keynumber(ik) < 0 then [all …]
|
/dports/science/octopus/octopus-10.5/src/states/ |
H A D | linear_response_inc.F90 | 58 if(abs(st%eigenval(ist, ik) - st%eigenval(jst, ik)) < CNST(1e-6)) then 67 if(abs(st%occ(ist, ik) * st%occ(jst, ik)) > M_EPSILON) then 89 xx = (st%eigenval(ist, ik) - st%eigenval(jst, ik))/dsmear 96 alpha_j = lr_alpha_j(st, jst, ik) 97 delta_e = st%eigenval(ist, ik) - st%eigenval(jst, ik) - omega 167 do ik = st%d%kpt%start, st%d%kpt%end 251 integer :: ist, ik 256 call X(lr_orth_vector) (mesh, st, lr%X(dl_psi)(:,:, ist, ik), ist, ik, omega) 271 integer :: ik, idim, ist local 278 do ik = 1, st%d%nspin [all …]
|
/dports/science/quantum-espresso/q-e-qe-6.7.0/Modules/ |
H A D | bspline.f90 | 290 work(ik) = bcoef(leftx+ik-kx) 291 dl(ik) = x - xknot(leftx+ik-kx) 292 dr(ik) = xknot(leftx+ik) - x 374 work(ik) = bcoef(leftx+ik-kx) 376 dr(ik) = xknot(leftx+ik) - x 398 dr(ik) = xknot(leftx+ik) - x 399 dl(ik) = x - xknot(leftx+1-ik) 411 work(ik) = bcoef(leftx+ik-kx) 412 dr(ik) = xknot(leftx+ik) - x 655 work(ik) = bcoef(leftx+ik-kx) [all …]
|
/dports/science/elk/elk-7.2.42/src/ |
H A D | gwspecf.f90 | 28 call getevalsv(filext,ik,vkl(:,ik),evalsv(:,ik)) 29 call getoccsv(filext,ik,vkl(:,ik),occsv(:,ik)) 52 if (mod(ik-1,np_mpi).ne.lp_mpi) cycle 54 write(*,'("Info(gwspecf): ",I6," of ",I6," k-points")') ik,nkpt 57 call getgwsefm(ik,se) 59 call dysonr(ik,wr,se,sf) 62 call writegwsf(ik,sf) 66 sft(:)=sft(:)+wkpt(ik)*sf(:)
|
H A D | writehmlbse.f90 | 32 call getevalsv(filext,ik,vkl(:,ik),evalsv(:,ik)) 33 call getoccsv(filext,ik,vkl(:,ik),occsv(:,ik)) 57 if (mod(ik-1,np_mpi).ne.lp_mpi) cycle 58 jk=ivkik(ivk(1,ik),ivk(2,ik),ivk(3,ik)) 60 ist=istbse(i,ik) 62 jst=jstbse(j,ik) 63 a=ijkbse(i,j,ik)
|
H A D | gradwf2.f90 | 35 call match(ngk(ispn,ik),vgkc(:,:,ispn,ik),gkc(:,ispn,ik), & 50 if (abs(occsv(j,ik)).lt.epsocc) cycle 51 wo=wkpt(ik)*occsv(j,ik) 97 if (abs(occsv(j,ik)).lt.epsocc) cycle 98 wo=wkpt(ik)*occsv(j,ik)/omega 109 do igk=1,ngk(jspn,ik) 110 ifg=igfft(igkig(igk,jspn,ik)) 118 do igk=1,ngk(1,ik) 119 ifg=igfft(igkig(igk,1,ik)) 128 do igk=1,ngk(jspn,ik) [all …]
|
/dports/math/linpack/linpack-1.0_14/ |
H A D | zhpsl.f | 57 ik = (n*(n - 1))/2 59 kk = ik + k 84 ik = ik - k 90 ikm1 = ik - (k - 1) 110 km1k = ik + k - 1 111 kk = ik + k 121 ik = ik - (k + 1) - k 129 ik = 0 150 ik = ik + k 162 ikp1 = ik + k [all …]
|
H A D | zspsl.f | 57 ik = (n*(n - 1))/2 59 kk = ik + k 84 ik = ik - k 90 ikm1 = ik - (k - 1) 110 km1k = ik + k - 1 111 kk = ik + k 121 ik = ik - (k + 1) - k 129 ik = 0 150 ik = ik + k 162 ikp1 = ik + k [all …]
|
H A D | chpsl.f | 57 ik = (n*(n - 1))/2 59 kk = ik + k 84 ik = ik - k 90 ikm1 = ik - (k - 1) 110 km1k = ik + k - 1 111 kk = ik + k 121 ik = ik - (k + 1) - k 129 ik = 0 150 ik = ik + k 162 ikp1 = ik + k [all …]
|
H A D | cspsl.f | 57 ik = (n*(n - 1))/2 59 kk = ik + k 84 ik = ik - k 90 ikm1 = ik - (k - 1) 110 km1k = ik + k - 1 111 kk = ik + k 121 ik = ik - (k + 1) - k 129 ik = 0 150 ik = ik + k 162 ikp1 = ik + k [all …]
|
H A D | dspsl.f | 57 ik = (n*(n - 1))/2 59 kk = ik + k 84 ik = ik - k 90 ikm1 = ik - (k - 1) 110 km1k = ik + k - 1 111 kk = ik + k 121 ik = ik - (k + 1) - k 129 ik = 0 150 ik = ik + k 162 ikp1 = ik + k [all …]
|
H A D | sspsl.f | 57 ik = (n*(n - 1))/2 59 kk = ik + k 84 ik = ik - k 90 ikm1 = ik - (k - 1) 110 km1k = ik + k - 1 111 kk = ik + k 121 ik = ik - (k + 1) - k 129 ik = 0 150 ik = ik + k 162 ikp1 = ik + k [all …]
|
H A D | zspco.f | 135 kk = ik + k 175 ik = ik - k 176 if (ks .eq. 2) ik = ik - (k + 1) 185 ik = 0 201 ik = ik + k 202 if (ks .eq. 2) ik = ik + (k + 1) 252 ik = ik - k 253 if (ks .eq. 2) ik = ik - (k + 1) 263 ik = 0 279 ik = ik + k [all …]
|
H A D | chpco.f | 132 kk = ik + k 172 ik = ik - k 173 if (ks .eq. 2) ik = ik - (k + 1) 182 ik = 0 198 ik = ik + k 199 if (ks .eq. 2) ik = ik + (k + 1) 249 ik = ik - k 250 if (ks .eq. 2) ik = ik - (k + 1) 260 ik = 0 276 ik = ik + k [all …]
|
H A D | cspco.f | 131 kk = ik + k 171 ik = ik - k 172 if (ks .eq. 2) ik = ik - (k + 1) 181 ik = 0 197 ik = ik + k 198 if (ks .eq. 2) ik = ik + (k + 1) 248 ik = ik - k 249 if (ks .eq. 2) ik = ik - (k + 1) 259 ik = 0 275 ik = ik + k [all …]
|
H A D | zhpco.f | 136 kk = ik + k 176 ik = ik - k 177 if (ks .eq. 2) ik = ik - (k + 1) 186 ik = 0 202 ik = ik + k 203 if (ks .eq. 2) ik = ik + (k + 1) 253 ik = ik - k 254 if (ks .eq. 2) ik = ik - (k + 1) 264 ik = 0 280 ik = ik + k [all …]
|
H A D | chpdi.f | 92 ik = 0 94 kk = ik + k 108 ikp1 = ik + k 140 ik = ik + k 148 ik = 0 151 kk = ik + k 152 ikp1 = ik + k 221 ksj = ik + ks 224 jk = ik + j 237 ik = ik + k [all …]
|
/dports/math/sisl/SISL-SISL-4.6.0-44-g9114631/src/ |
H A D | sh1927.c | 59 int ik; 149 if (pcurve->ik != ik) goto err106; 156 if ((sbder=newarray(ik*ik,DOUBLE)) == SISL_NULL) goto err101; 163 if (ilend > ik) ilend = ik; 165 if (irend > ik) irend = ik; 180 s1220(etau,ik,in,&kleft,etau[ik-1],ilend-1,sbder,&kstat); 186 if (km == ik || etau[ik-km-1] < etau[ik-km]) break; 195 if (km < ik) 291 for (krh=0, kstop=ik-km, kr=km+1; kr<ik-1; kr++) 303 krh = ik - km; [all …]
|
/dports/math/libxsmm/libxsmm-1.16.3/src/template/ |
H A D | libxsmm_dnn_rnncell_st_gru_bwdupd_nc_ck_generic.tpl.c | 348 ik = ikb*bk; 385 for (ik = 0, ikb = 0; ikb < KB_BLOCKS; ik += bk, ikb++) { 393 for (ik = 0, ikb = 0; ikb < KB_BLOCKS; ik += bk, ikb++) { 400 for (ik = 0, ikb = 0; ikb < KB_BLOCKS; ik += bk, ikb++) { 414 ik = ikb*bk; 596 for (ik = thr_begin_k; ik < thr_end_k; ik++) { 598 dbi[ik] += LIBXSMM_VLA_ACCESS(2, di, in, ik, K); 599 dbc[ik] += LIBXSMM_VLA_ACCESS(2, dc, in, ik, K); 600 dbf[ik] += LIBXSMM_VLA_ACCESS(2, df, in, ik, K); 613 ik = ikb*bk; [all …]
|
/dports/science/berkeleygw/BGW-2.0.0/BSE/ |
H A D | genwf_kernel.f90 | 85 if (ik .ne. ikold .and. ik .ne. ikold2 .and. ifirst .ne. 0) then 137 if (ik .ne. ikold .and. ik .ne. ikold2) then 147 if (ik .ne. ikold .and. ik .ne. ikold2) then 157 qk(1:3)=kg%f(1:3,ik) 171 if (ik .ne. ikold .and. ik .ne. ikold2) then 242 if (ik .ne. ikold .and. ik .ne. ikold2 .and. ifirst .ne. 0) then 294 if (ik .ne. ikold .and. ik .ne. ikold2) then 304 if (ik .ne. ikold .and. ik .ne. ikold2) then 318 if (ik .ne. ikold .and. ik .ne. ikold2) then 364 if (ik .ne. ikold .and. ik .ne. ikold2) then [all …]
|
/dports/math/scilab/scilab-6.1.1/scilab/modules/overloading/macros/ |
H A D | %s_i_s.sci | 64 if type(ik)==2 |type(ik)==129 then // size implicit subscript $... 65 ik=round(horner(ik,dims(k))) // explicit subscript 68 ik=find(ik) 72 ik=1:dims(k) 95 ik=1 106 ik=round(ik) 111 if size(ik,"*")>1 then 112 ik=ik(:) 114 I=(dims1(k)*I).*.ones(ik)+ones(I).*.(ik-1) 119 I=dims1(k)*I+ik-1 [all …]
|
/dports/science/quantum-espresso/q-e-qe-6.7.0/EPW/src/ |
H A D | grid.f90 | 53 INTEGER :: ik local 149 xkf_(idir, ik) = xkf_(idir, ik) + 1.0d0 299 xkfd(:, ik, 1) = xkf(:, ik) + (/deltaq, 0.d0, 0.d0/) 300 xkfd(:, ik, 2) = xkf(:, ik) - (/deltaq, 0.d0, 0.d0/) 301 xkfd(:, ik, 3) = xkf(:, ik) + (/0.d0, deltaq, 0.d0/) 1734 map_rebal((kpt_in(ik) + 1) / 2) = ik 1735 map_rebal_inv(ik) = (kpt_in(ik) + 1) / 2 1738 map_rebal((kpt_out(ik) + 1) /2) = ikpt + ik 1774 map_rebal(map_rebal_inv(ik)) = ik 1860 nsym(ik) = nsym(ik) + 1 [all …]
|
/dports/science/nwchem-data/nwchem-7.0.2-release/src/ddscf/ |
H A D | fock_2e_lab.F | 460 fik(ik) = fik(ik) + g*dlj(lj) 634 fik(1,ik) = fik(1,ik) + g*dlj(1,lj) 640 fik(2,ik) = fik(2,ik) + g*dlj(2,lj) 815 fik(1,ik) = fik(1,ik) + g*dlj(1,lj) 821 fik(2,ik) = fik(2,ik) + g*dlj(2,lj) 827 fik(3,ik) = fik(3,ik) + g*dlj(3,lj) 1002 fik(1,ik) = fik(1,ik) + g*dlj(1,lj) 1008 fik(2,ik) = fik(2,ik) + g*dlj(2,lj) 1275 fik(ik) = fik(ik) + g*dlj(lj) 1717 fock(ik) = fock(ik) + gk*dens(lj) [all …]
|
/dports/science/nwchem/nwchem-7b21660b82ebd85ef659f6fba7e1e73433b0bd0a/src/ddscf/ |
H A D | fock_2e_lab.F | 460 fik(ik) = fik(ik) + g*dlj(lj) 634 fik(1,ik) = fik(1,ik) + g*dlj(1,lj) 640 fik(2,ik) = fik(2,ik) + g*dlj(2,lj) 815 fik(1,ik) = fik(1,ik) + g*dlj(1,lj) 821 fik(2,ik) = fik(2,ik) + g*dlj(2,lj) 827 fik(3,ik) = fik(3,ik) + g*dlj(3,lj) 1002 fik(1,ik) = fik(1,ik) + g*dlj(1,lj) 1008 fik(2,ik) = fik(2,ik) + g*dlj(2,lj) 1275 fik(ik) = fik(ik) + g*dlj(lj) 1717 fock(ik) = fock(ik) + gk*dens(lj) [all …]
|