14.08.2013 Views

Tutorial slides (PDF) - Clemson University

Tutorial slides (PDF) - Clemson University

Tutorial slides (PDF) - Clemson University

SHOW MORE
SHOW LESS

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

GPAW code structure<br />

Not simply a Python wrapper on legacy Fortran/C code<br />

Python for coding the high-level algorithm, lots of OO<br />

C for coding numerical intense operations, built on NumPy<br />

Use BLAS and LAPACK whenever possible<br />

Here is some pseudo code for iterative eigensolver:<br />

for i in xrange(max SCF):<br />

for n in xrange(number of bands):<br />

R_ng = apply(H_gg,Psi_ng) % Compute residuals<br />

rk(1.0, R_ng, 0.0, H_mn) % construct Hamiltonian<br />

KS-DFT algorithms are well-known and computationally<br />

intensive parts are known a priori.

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!