Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Tweak bookkeeping in LOBPCG #980

Open
wants to merge 3 commits into
base: master
Choose a base branch
from
Open

Tweak bookkeeping in LOBPCG #980

wants to merge 3 commits into from

Conversation

antoine-levitt
Copy link
Member

@mfherbst does that fix #979?

@antoine-levitt antoine-levitt marked this pull request as ready for review June 5, 2024 07:00
@mfherbst
Copy link
Member

mfherbst commented Jun 6, 2024

No, unfortunately not, here is a reproducer (you'll need the Artifacts.toml from dftk in your active project).

using LazyArtifacts
using DFTK
using AtomsIO
using Random

lattice = [0.0133767   3.95139   3.81366;
           3.70274    -0.113618  3.76523;
           3.78117     3.8285    0.0147007]
atoms = [ElementPsp(:Al; psp=load_psp(artifact"pd_nc_sr_pbe_stringent_0.4.1_upf/Al.upf"))]
positions = [zeros(3)]
model  = model_PBE(lattice, atoms, positions; 
                   smearing=Smearing.Gaussian(), temperature=1e-3)
basis  = PlaneWaveBasis(model; Ecut=10, kgrid=[6,6,6]);

Random.seed!(20)
scfres = self_consistent_field(basis; maxiter=2);

@mfherbst
Copy link
Member

@antoine-levitt ping

@antoine-levitt
Copy link
Member Author

Hrmf it's still super annoying, I'll give it a shot but no promises...

No, unfortunately not, here is a reproducer (you'll need the Artifacts.toml from dftk in your active project).

Sorry I'm dense, but how do I do that? activate DFTK doesn't seem to do

@mfherbst
Copy link
Member

mfherbst commented Aug 2, 2024

If you activate / instantiate the lobpcg_fix branch and have AtomsIO installed in your general environment, it should work, since the DFTK source tree already includes the Artifacts.toml i.e.

using Pkg; Pkg.add(["AtomsIO"])
Pkg.activate("/path/to/your/location/of/lobpcg_fix/branch")

... paste code above ...

@mfherbst
Copy link
Member

Updated reproducer with current master:

using DFTK
using AtomsIO
using Random
using PseudoPotentialData

lattice = [0.0133767   3.95139   3.81366;
           3.70274    -0.113618  3.76523;
           3.78117     3.8285    0.0147007]
atoms = [ElementPsp(:Al, PseudoFamily("dojo.nc.sr.pbe.v0_4_1.oncvpsp3.stringent.upf"))]
positions = [zeros(3)]
model  = model_DFT(lattice, atoms, positions; functionals=LDA(), smearing=Smearing.Gaussian(), temperature=1e-3)
basis  = PlaneWaveBasis(model; Ecut=10, kgrid=[6,6,6])

Random.seed!(20)
scfres = self_consistent_field(basis; maxiter=2)

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

Successfully merging this pull request may close these issues.

Eigensolver not converged
2 participants