Add example solving -div(a∇u)=f using DivAgradOperator with general inclusions#1228

Draft

Copilot wants to merge 18 commits into

masterfrom

copilot/implement-simple-div-solve

Commits on Jan 14, 2026

Initial plan
Copilot
committed
Add initial divAgrad_solver example with spherical inclusions

Copilot
and
antoine-levitt
committed
Simplify ene_ops and use custom diagonal preconditioner

Copilot
and
antoine-levitt
committed
Fix DivAGrad factor and projection for G=0 mode

Copilot
and
antoine-levitt
committed
Fix LinearMap posdef flag and occupation array structure

Copilot
and
antoine-levitt
committed
Fix convert_dual to use DFTK namespace

Copilot
and
antoine-levitt
committed
Remove test files that were accidentally committed

Copilot
and
antoine-levitt
committed
Co-authored-by: antoine-levitt <88065+antoine-levitt@users.noreply.github.com>
Copilot
committed
Address code review: parametric types, pseudo-inverse preconditioner, better test case

Copilot
and
antoine-levitt
committed
Fix code review issues: move PseudoInversePreconditioner to module level, correct isposdef flag

Copilot
and
antoine-levitt
committed
Improve preconditioner: zero DC component first in pseudo-inverse logic

Copilot
and
antoine-levitt
committed
Add 1/√a (-Δ)^(-1) 1/√a preconditioner with mathematical proof of I+compact form

Copilot
and
antoine-levitt
committed
Revert "Add 1/√a (-Δ)^(-1) 1/√a preconditioner with mathematical proof of I+compact form"
antoine-levitt
committed
Misc fixes
antoine-levitt
committed
Refactor to support non-spherical inclusions with DivAGradFromReal/Fourier

Copilot
and
antoine-levitt
committed
first attempt at fixing
antoine-levitt
committed
up
antoine-levitt
committed
Refactor from claude
antoine-levitt
committed

Provide feedback

Saved searches

Use saved searches to filter your results more quickly

Add example solving -div(a∇u)=f using DivAgradOperator with general inclusions#1228

Add example solving -div(a∇u)=f using DivAgradOperator with general inclusions#1228
Copilot wants to merge 18 commits into
masterfrom
copilot/implement-simple-div-solve

Commits on Jan 14, 2026

Initial plan

Add initial divAgrad_solver example with spherical inclusions

Simplify ene_ops and use custom diagonal preconditioner

Fix DivAGrad factor and projection for G=0 mode

Fix LinearMap posdef flag and occupation array structure

Fix convert_dual to use DFTK namespace

Remove test files that were accidentally committed

Co-authored-by: antoine-levitt <88065+antoine-levitt@users.noreply.github.com>

Address code review: parametric types, pseudo-inverse preconditioner, better test case

Fix code review issues: move PseudoInversePreconditioner to module level, correct isposdef flag

Improve preconditioner: zero DC component first in pseudo-inverse logic

Add 1/√a (-Δ)^(-1) 1/√a preconditioner with mathematical proof of I+compact form

Revert "Add 1/√a (-Δ)^(-1) 1/√a preconditioner with mathematical proof of I+compact form"

Misc fixes

Refactor to support non-spherical inclusions with DivAGradFromReal/Fourier

first attempt at fixing

up

Refactor from claude