Home
last modified time | relevance | path

Searched refs:idir2 (Results 1 – 25 of 31) sorted by relevance

12

/dports/science/octopus/octopus-10.5/src/sternheimer/
H A Dem_resp_inc.F90251 do idir2 = 1, idir
269 call X(inhomog_kB_tot)(sh, sys, idir, magn_dir(idir2, 1), magn_dir(idir2, 2), &
272 k2_lr(max(magn_dir(idir2, 1), idir), min(magn_dir(idir2,1), idir), 1:1),&
273 k2_lr(max(magn_dir(idir2, 2), idir), min(magn_dir(idir2,2), idir), 1:1), inhomog)
341 str_tmp = kdotp_wfs_tag(min(idir, idir2), max(idir, idir2))
392 do idir2 = 1, gr%sb%dim
417 kdotp_lr(magn_dir(idir2, 1), 1:1), kdotp_lr(magn_dir(idir2, 2), 1:1), &
418 k2_lr(max(magn_dir(idir2, 1), idir),min(magn_dir(idir2, 1), idir), 1:1),&
433 … em_vars%lr(idir2, :, ifactor),k2_lr(max(idir2, idir), min(idir2,idir), 1:1), inhomog)
498 em_vars%alpha0(idir, idir2, ifactor) = em_vars%alpha(idir, idir2, ifactor)
[all …]
H A Dem_resp_calc_inc.F901243 do idir2 = 1, ndir
1306 zpol0(idir3, idir2, idir1) = zpol0(idir3, idir2, idir1) - &
1472 do idir2 = 1, ndir
1474 … zpol(idir3, idir1, idir2) = zpol(idir3, idir1, idir2) + weight * zpol_k(idir3, idir1, idir2)
1526 do idir2 = 1, ndir
1528 zpol(idir3, idir1, idir2) = zpol(idir3, idir1, idir2) + zpol0(idir3, idir1, idir2)
1543 do idir2 = 1, ndir
1745 do idir2 = 1, ndir
1772 …lr_kk(max(idir1, idir2), min(idir1, idir2))%X(dl_psi)(:,:, ist, ik), Hdl_kk(:,:, idir1, idir2), is…
1989 magn(idir1, idir2) = magn(idir1, idir2) + M_FOURTH * (&
[all …]
H A Dkdotp.F9094 integer :: idir, idir2, ierr, pdim, ispin local
204 do idir2 = idir, pdim
205 call lr_init(kdotp_vars%lr2(1, idir, idir2))
223 do idir2 = idir, pdim
224 str_tmp = kdotp_wfs_tag(idir, idir2)
290 do idir2 = idir, pdim
295 call pert_setup_dir(pert2, idir2)
300 …kdotp_vars%lr2(1:1, idir, idir2), kdotp_vars%perturbation2, restart_dump, "", kdotp_wfs_tag(idir, …
305 …kdotp_vars%lr2(1:1, idir, idir2), kdotp_vars%perturbation2, restart_dump, "", kdotp_wfs_tag(idir, …
337 do idir2 = idir, pdim
[all …]
H A Dem_resp.F90268 do idir2 = 1, gr%sb%dim
318 do idir2 = 1, gr%sb%dim
321 if(idir2 <= idir) then
418 do idir2 = 1, gr%sb%dim
631 do idir2 = 1, gr%sb%dim
633 if(idir2 <= idir) call lr_dealloc(k2_lr(idir, idir2, 1))
651 do idir2 = 1, gr%sb%dim
1105 integer :: idir, idir2
1125 do idir2 = 1, gr%sb%dim
1126 cross(idir, idir2) = units_from_atomic(units_out%length**2, cross(idir, idir2))
[all …]
H A Dkdotp_calc_inc.F9030 integer :: ik, ist, ist2, idir1, idir2, pdim
70 do idir2 = 1, pdim
73 if (idir1 < idir2) then
74 eff_mass_inv(idir1, idir2, ist, ik) = eff_mass_inv(idir2, idir1, ist, ik)
79 …proj_dl_psi(1:mesh%np, 1:sys%hm%d%dim) = lr(1, idir2)%X(dl_psi)(1:mesh%np, 1:sys%hm%d%dim, ist, ik)
97 eff_mass_inv(idir1, idir2, ist, ik) = M_TWO * R_REAL(term)
99 call pert_setup_dir(perturbation, idir1, idir2)
102 eff_mass_inv(idir1, idir2, ist, ik) = &
103 … eff_mass_inv(idir1, idir2, ist, ik) - R_REAL(X(mf_dotp)(mesh, sys%hm%d%dim, psi, pertpsi2))
H A Dem_resp_calc.F90228 character(len=100) function em_wfs_tag(idir, ifactor, idir2, ipert) result(str)
231 integer, optional, intent(in) :: idir2 local
237 if(present(idir2)) write(str, '(3a)') trim(str), "_", index2axis(idir2)
/dports/math/sisl/SISL-SISL-4.6.0-44-g9114631/src/
H A Dsh1783.c70 int idir1, int idir2, double elast[], double enext[], int *jstat) in sh1783() argument
79 int idir2;
219 tend2 = (idir2 == 1) ? st2[kn2] : st2[kk2 - 1];
237 fevalc2 = (idir2 == 1) ? s1221 : s1227;
260 while (idir1 * tx1 < idir1 * tend1 && idir2 * ty1 < idir2 * tend2)
371 tyincre = idir2 * (st2[kleftc2] - ty1);
381 ty2 = ty1 + idir2 * tyincre;
400 if (idir2 * ty >= idir2 * tend2)
428 ty1 = (idir2 > 0) ? MAX(ty1,tpos) : MIN(ty1,tpos);
445 else if (idir1 * tx1 >= idir1 * tend1 && idir2 * ty1 >= idir2 * tend2)
[all …]
/dports/science/py-abipy/abipy-0.9.0/abipy/dfpt/
H A Dconverters.py387 for idir2 in range(3):
388 born_data[idir1, ipert1, idir2, ipert2] = dm1[idir1, idir2] * 2 * np.pi + 0.0j
389 born_data[idir2, ipert2, idir1, ipert1] = dm1[idir1, idir2] * 2 * np.pi + 0.0j
396 for idir2 in range(3):
397 born_data[idir1, ipert1, idir2, ipert2] = dm1[idir1, idir2] * np.pi * volume_bohr + 0.0j
463 for idir2 in range(3):
464 … q_data[(idir1+1, ipert1+1, idir2+1, ipert2+1)] = dm[ipert1, idir1, ipert2, idir2]
472 … q_data[(idir1+1, ipert1+1, idir2+1, ipert2+1)] = born_data[ipert1, idir1, ipert2, idir2]
473 … q_data[(idir2+1, ipert2+1, idir1+1, ipert1+1)] = born_data[ipert2, idir2, ipert1, idir1]
476 for idir2 in range(3):
[all …]
/dports/science/cp2k-data/cp2k-7.1.0/src/
H A Dqs_linres_epr_ownutils.F124 DO idir2 = 1, 3
125 gsum = gsum + epr_env%g_total(idir1, idir2)
182 DO idir2 = 1, 3
183 g_sym(idir1, idir2) = (epr_env%g_total(idir1, idir2) + &
356 CALL set_vecp(idir1, idir2, idir3)
360 CALL qs_rho_get(jrho1_set(idir2)%rho, rho_r=jrho2_r)
435 CALL set_vecp(idir1, idir2, idir3)
857 DO idir2 = 1, 3
858 IF (idir /= idir2) THEN
864 CALL set_vecp_rev(idir, idir2, idir3)
[all …]
H A Dqs_linres_current.F674 idir2 = 1
2197 DO idir2 = 1, 3
2207 DO idir2 = 1, 3
2254 DO idir2 = 1, 9
2339 IF (idir2 .EQ. ii .OR. idir2 .EQ. iii) THEN
2344 … IF (idir2 .EQ. ind_m2(ii, iiiB) .OR. idir2 .EQ. ind_m2(iii, iiiB) .OR. idir2 .EQ. iiiB) THEN
2349 … IF (idir2 .EQ. ind_m2(ii, iiB) .OR. idir2 .EQ. ind_m2(iii, iiB) .OR. idir2 .EQ. iiB) THEN
2495 DO idir2 = 1, 9
2611 DO idir2 = 1, 3
2621 DO idir2 = 1, 3
[all …]
H A Dqs_epr_hyp.F87 idir2, ig, ikind, ir, iso, jatom, & local
313 DO idir2 = idir1, 3 ! tensor symmetry
316 idir1, idir2)
327 hypaniso(idir1, idir2, iatom) = hypaniso(idir1, idir2, iatom) + esum
349 DO idir2 = idir1, 3
350 hypaniso(idir1, idir2, iatom) = hypaniso(idir1, idir2, iatom)* &
352 IF (idir1 /= idir2) THEN
353 hypaniso(idir2, idir1, iatom) = hypaniso(idir1, idir2, iatom)
H A Dqs_linres_atom_current.F117 INTEGER :: bo(2), handle, iac, iat, iatom, ibc, idir2, ii, iii, ikind, ispin, istat, jatom, & local
176 idir2 = 1
178 CALL set_vecp_rev(idir, iB, idir2)
1210 INTEGER :: ia, idir2, iiB, iiiB, ir, & local
1289 idir2 = 1
1291 CALL set_vecp_rev(idir, iB, idir2)
1292 scale = fac_vecp(idir, iB, idir2)
1311 & + scale*(a(idir2, ia)*r(ir) - g(idir2, ir, ia))*Fr_h(ir, iso)
1319 & + scale*(a(idir2, ia)*r(ir) - g(idir2, ir, ia))*Fr_s(ir, iso)
1333 & + scale*a(idir2, ia)*Fr_h(ir, iso)
[all …]
H A Dqs_linres_nmr_shift.F111 INTEGER :: handle, idir, idir2, idir3, iiB, iiiB, & local
182 DO idir2 = 1, 3
183 IF (idir /= idir2) THEN
186 pw_gspace_work, idir2, 0.0_dp)
189 CALL set_vecp_rev(idir, idir2, idir3)
190 scale_fac = fac_vecp(idir3, idir2, idir)
/dports/science/cp2k/cp2k-2e995eec7fd208c8a72d9544807bd8b8ba8cd1cc/src/
H A Dqs_linres_epr_ownutils.F124 DO idir2 = 1, 3
125 gsum = gsum + epr_env%g_total(idir1, idir2)
182 DO idir2 = 1, 3
183 g_sym(idir1, idir2) = (epr_env%g_total(idir1, idir2) + &
356 CALL set_vecp(idir1, idir2, idir3)
360 CALL qs_rho_get(jrho1_set(idir2)%rho, rho_r=jrho2_r)
435 CALL set_vecp(idir1, idir2, idir3)
857 DO idir2 = 1, 3
858 IF (idir /= idir2) THEN
864 CALL set_vecp_rev(idir, idir2, idir3)
[all …]
H A Dqs_epr_hyp.F87 idir2, ig, ikind, ir, iso, jatom, & local
313 DO idir2 = idir1, 3 ! tensor symmetry
316 idir1, idir2)
327 hypaniso(idir1, idir2, iatom) = hypaniso(idir1, idir2, iatom) + esum
349 DO idir2 = idir1, 3
350 hypaniso(idir1, idir2, iatom) = hypaniso(idir1, idir2, iatom)* &
352 IF (idir1 /= idir2) THEN
353 hypaniso(idir2, idir1, iatom) = hypaniso(idir1, idir2, iatom)
H A Dqs_linres_current.F681 idir2 = 1
2380 DO idir2 = 1, 3
2390 DO idir2 = 1, 3
2437 DO idir2 = 1, 9
2522 IF (idir2 .EQ. ii .OR. idir2 .EQ. iii) THEN
2527 … IF (idir2 .EQ. ind_m2(ii, iiiB) .OR. idir2 .EQ. ind_m2(iii, iiiB) .OR. idir2 .EQ. iiiB) THEN
2532 … IF (idir2 .EQ. ind_m2(ii, iiB) .OR. idir2 .EQ. ind_m2(iii, iiB) .OR. idir2 .EQ. iiB) THEN
2678 DO idir2 = 1, 9
2794 DO idir2 = 1, 3
2804 DO idir2 = 1, 3
[all …]
H A Dqs_linres_atom_current.F117 INTEGER :: bo(2), handle, iac, iat, iatom, ibc, idir2, ii, iii, ikind, ispin, istat, jatom, & local
176 idir2 = 1
178 CALL set_vecp_rev(idir, iB, idir2)
1210 INTEGER :: ia, idir2, iiB, iiiB, ir, & local
1289 idir2 = 1
1291 CALL set_vecp_rev(idir, iB, idir2)
1292 scale = fac_vecp(idir, iB, idir2)
1311 & + scale*(a(idir2, ia)*r(ir) - g(idir2, ir, ia))*Fr_h(ir, iso)
1319 & + scale*(a(idir2, ia)*r(ir) - g(idir2, ir, ia))*Fr_s(ir, iso)
1333 & + scale*a(idir2, ia)*Fr_h(ir, iso)
[all …]
H A Dqs_linres_nmr_shift.F111 INTEGER :: handle, idir, idir2, idir3, iiB, iiiB, & local
182 DO idir2 = 1, 3
183 IF (idir /= idir2) THEN
186 pw_gspace_work, idir2, 0.0_dp)
189 CALL set_vecp_rev(idir, idir2, idir3)
190 scale_fac = fac_vecp(idir3, idir2, idir)
/dports/science/octopus/octopus-10.5/src/system/
H A Doutput_linear_response_inc.F9032 integer :: ik, ist, idim, ierr, is, idir2 local
66 do idir2 = 1, gr%mesh%sb%dim
67 tmp(1:gr%mesh%np) = -gr%mesh%x(1:gr%mesh%np, idir2) * lr%X(dl_rho)(:, is)
69 write(fname, '(4a)') 'alpha_density-', index2axis(idir2), '-', index2axis(idir)
71 … write(fname, '(a,i1,4a)') 'alpha_density-sp', is, '-', index2axis(idir2), '-', index2axis(idir)
84 do idir2 = 1, gr%mesh%sb%dim
86 write(fname, '(4a)') 'lr_current-', index2axis(idir2), '-', index2axis(idir)
88 … write(fname, '(a,i1,4a)') 'lr_current-sp', is, '-', index2axis(idir2), '-', index2axis(idir)
91 gr%mesh, lr%dl_j(:, idir2, is), &
H A Dpert_inc.F90479 integer :: iatom, idir, idir2, ip, idim
521 do idir2 = 1, gr%sb%dim
524 …te_r)(hm%ep%proj(iatom), gr%mesh, gr%der%boundaries, hm%d%dim, idir2, ik, f_in2(:, :, idir2), vrnl)
536 xf(ip, idim) = gr%mesh%x(ip, idir) * f_in2(ip, idim, idir2)
541 …(projector_commute_r)(hm%ep%proj(iatom), gr%mesh, gr%der%boundaries, hm%d%dim, idir2, ik, xf, vrnl)
/dports/science/qbox/qbox-public-rel1_73_3/util/dynmat/
H A Ddynmat.cpp108 for ( int idir2 = 0; idir2 < 3; idir2++ ) in main() local
111 int j = 3*ia2+idir2; in main()
112 double aij = ( f[1][ia1][idir1][ia2][idir2] - in main()
113 f[0][ia1][idir1][ia2][idir2] ) / ( 2.0 * h ); in main()
/dports/science/octopus/octopus-10.5/src/grid/
H A Dsymmetrizer_inc.F90259 integer :: idir1, idir2, idir3, ndir
277 do idir2 = 1, ndir
281 i2 = abs(1 * rot(1,idir2) + 2 * rot(2,idir2) + 3 * rot(3,idir2))
285 tensor(idir1,idir2,idir3) * det
H A Dio_function.F90465 integer :: idir, idir2, iatom, index_
505 mesh%sb%rlattice(idir2, idir)), idir2 = 1, mesh%sb%dim)
H A Dio_function_inc.F901249 integer :: ix, iy, iz, idir, idir2, iatom local
1284 (mesh%spacing(idir)*mesh%sb%rlattice_primitive(idir2, idir), idir2 = 1, 3)
1323 integer :: ix, iy, iz, idir2, ix2, iy2, iz2, my_n(3) local
1363 do idir2 = 1, 3
1364 …lattice_vectors(idir2, idir) = mesh%spacing(idir) * (my_n(idir) - 1) * mesh%sb%rlattice_primitive(…
1382 lattice_vectors(idir2, idir)), idir2 = 1, 3)
/dports/science/octopus/octopus-10.5/src/ions/
H A Dsimul_box.F901007 integer :: idir, idir2, ispec local
1060 …write(message(2+idir),'(9f12.6)') (units_from_atomic(units_out%length, sb%rlattice(idir2, idir)), &
1061 idir2 = 1, sb%dim)
1072 …ge(1+idir),'(3f12.6)') (units_from_atomic(unit_one / units_out%length, sb%klattice(idir2, idir)), &
1073 idir2 = 1, sb%dim)
1093 integer :: idir1, idir2 local
1113 do idir2 = 1, sb%dim
1114 …write(iunit, '(x,f11.6)', advance='no') units_from_atomic(unit_angstrom, sb%rlattice(idir2, idir1))

12