mirror of
https://github.com/QuantumPackage/qp2.git
synced 2024-11-07 14:03:37 +01:00
commit
faaddd438a
@ -696,8 +696,8 @@ def pyscf2QP2(cell,mf, kpts, kmesh=None, cas_idx=None, int_threshold = 1E-8,
|
||||
# k,mo,ao(,2)
|
||||
mo_coef_f = np.array(mo_k.transpose((0,2,1)),order='c',dtype=np.complex128)
|
||||
#mo_coef_blocked=block_diag(*mo_k)
|
||||
#mo_coef_blocked_f = block_diag(*mo_coef_f)
|
||||
#qph5.create_dataset('mo_basis/mo_coef_complex',data=mo_coef_blocked_f.view(dtype=np.float64).reshape((Nk*nmo,Nk*nao,2)))
|
||||
mo_coef_blocked_f = block_diag(*mo_coef_f)
|
||||
qph5.create_dataset('mo_basis/mo_coef_complex',data=mo_coef_blocked_f.view(dtype=np.float64).reshape((Nk*nmo,Nk*nao,2)))
|
||||
qph5.create_dataset('mo_basis/mo_coef_kpts',data=mo_coef_f.view(dtype=np.float64).reshape((Nk,nmo,nao,2)))
|
||||
|
||||
if print_debug:
|
||||
|
@ -174,6 +174,7 @@ for cell_type in cell_types:
|
||||
print (cell_type)
|
||||
print ('===============================')
|
||||
s = diamond.structure.copy()
|
||||
s.change_units('B') #required for supertwists to be in 1/au instead of 1/Å
|
||||
kmap = s.kmap()
|
||||
|
||||
print ('supercell kpoints/twists')
|
||||
|
Loading…
Reference in New Issue
Block a user