ergo
|
Trace correcting purification. More...
#include <TC2.h>
Classes | |
class | Fun |
Help class for bisection root finding calls. More... | |
Public Member Functions | |
TC2 (Tmatrix &F, Tmatrix &DM, const int size, const int noc, const Treal trunc=0, const int maxmm=100) | |
Constructor Initializes everything. More... | |
Treal | fermi_level (Treal tol=1e-15) const |
Returns the Fermi level. More... | |
Treal | homo (Treal tol=1e-15) const |
Returns upper bound of the HOMO eigenvalue. More... | |
Treal | lumo (Treal tol=1e-15) const |
Returns lower bound of the LUMO eigenvalue. More... | |
int | n_multiplies () const |
Returns the number of used matrix matrix multiplications. More... | |
void | print_data (int const start, int const stop) const |
virtual | ~TC2 () |
Destructor. More... | |
Protected Member Functions | |
void | purify () |
Runs purification. More... | |
Protected Attributes | |
Tmatrix & | X |
Fock / Kohn-Sham matrix at initialization. More... | |
Tmatrix & | D |
Density matrix after purification. More... | |
const int | n |
System size. More... | |
const int | nocc |
Number of occupied orbitals. More... | |
const Treal | frob_trunc |
Threshold for the truncation. More... | |
const int | maxmul |
Number of tolerated matrix multiplications. More... | |
Treal | lmin |
Lower bound for eigenvalue spectrum. More... | |
Treal | lmax |
Upper bound for eigenvalue spectrum. More... | |
int | nmul |
Number of used matrix multiplications. More... | |
int | nmul_firstpart |
Number of used matrix multiplications in the first part of the purification. More... | |
Treal * | idemerror |
Upper bound of euclidean norm ||D-D^2||_2 before each step. More... | |
Treal * | tracediff |
The difference between the trace of the matrix and the number of occupied orbitals before each step. More... | |
int * | polys |
Choices of polynomials 0 for x^2 and 1 for 2x-x^2 Length: nmul. More... | |
Trace correcting purification.
This template instantiates the trace correcting purification algorithm developed by Niklasson [Phys. Rev. B 66, 155115 (2002)] with modifications by Rubensson and Rudberg [unpublished]. The template can be used with any matrix class Tmatrix that has the following member functions:
The matrix class should also support the following syntax:
where A and B are of type Tmatrix and alpha and beta are of type Treal.
mat::TC2< Treal, Tmatrix >::TC2 | ( | Tmatrix & | F, |
Tmatrix & | DM, | ||
const int | size, | ||
const int | noc, | ||
const Treal | trunc = 0 , |
||
const int | maxmm = 100 |
||
) |
Constructor Initializes everything.
Constructor.
< Run purification
F | Fock/Kohn-Sham matrix (input/workspace) |
DM | Density matrix (output) |
size | System size (Number of basis functions) |
noc | Number of occupied orbitals. |
trunc | Threshold for truncation in Frobenius norm. |
maxmm | Maximum aloud number of mm-multiplications. |
References mat::TC2< Treal, Tmatrix >::D, mat::TC2< Treal, Tmatrix >::frob_trunc, mat::TC2< Treal, Tmatrix >::idemerror, mat::TC2< Treal, Tmatrix >::lmax, mat::TC2< Treal, Tmatrix >::lmin, mat::TC2< Treal, Tmatrix >::maxmul, mat::TC2< Treal, Tmatrix >::nocc, mat::TC2< Treal, Tmatrix >::polys, mat::TC2< Treal, Tmatrix >::purify(), mat::TC2< Treal, Tmatrix >::tracediff, and mat::TC2< Treal, Tmatrix >::X.
|
inlinevirtual |
Destructor.
References mat::TC2< Treal, Tmatrix >::idemerror, mat::TC2< Treal, Tmatrix >::polys, and mat::TC2< Treal, Tmatrix >::tracediff.
Treal mat::TC2< Treal, Tmatrix >::fermi_level | ( | Treal | tol = 1e-15 | ) | const |
Returns the Fermi level.
Run after call to purify().
tol | Fault-tolerance for result. |
References mat::bisection().
Treal mat::TC2< Treal, Tmatrix >::homo | ( | Treal | tol = 1e-15 | ) | const |
Returns upper bound of the HOMO eigenvalue.
Run after call to purify().
tol | Fault-tolerance for result. |
References mat::bisection(), and template_blas_sqrt().
Treal mat::TC2< Treal, Tmatrix >::lumo | ( | Treal | tol = 1e-15 | ) | const |
Returns lower bound of the LUMO eigenvalue.
Run after call to purify().
tol | Fault-tolerance for result. |
References mat::bisection(), and template_blas_sqrt().
|
inline |
Returns the number of used matrix matrix multiplications.
References mat::TC2< Treal, Tmatrix >::nmul.
void mat::TC2< Treal, Tmatrix >::print_data | ( | int const | start, |
int const | stop | ||
) | const |
|
protected |
Runs purification.
Run by constructor.
References template_blas_sqrt().
Referenced by mat::TC2< Treal, Tmatrix >::TC2().
|
protected |
Density matrix after purification.
Referenced by mat::TC2< Treal, Tmatrix >::TC2().
|
protected |
Threshold for the truncation.
Referenced by mat::TC2< Treal, Tmatrix >::TC2().
|
protected |
Upper bound of euclidean norm ||D-D^2||_2 before each step.
This means: idemerror[i] = norm(D[i]-D[i]^2) where D[0] is the initial matrix and D[i] is the matrix after i steps in the purification. This value is calculated after the step since D[i]^2 or 2D[i] - D[i]^2 is needed. Length: nmul
Referenced by mat::TC2< Treal, Tmatrix >::TC2(), and mat::TC2< Treal, Tmatrix >::~TC2().
|
protected |
Upper bound for eigenvalue spectrum.
Referenced by mat::TC2< Treal, Tmatrix >::TC2().
|
protected |
Lower bound for eigenvalue spectrum.
Referenced by mat::TC2< Treal, Tmatrix >::TC2().
|
protected |
Number of tolerated matrix multiplications.
Referenced by mat::TC2< Treal, Tmatrix >::TC2().
|
protected |
System size.
|
protected |
Number of used matrix multiplications.
Referenced by mat::TC2< Treal, Tmatrix >::n_multiplies().
|
protected |
Number of used matrix multiplications in the first part of the purification.
|
protected |
Number of occupied orbitals.
Referenced by mat::TC2< Treal, Tmatrix >::TC2().
|
protected |
Choices of polynomials 0 for x^2 and 1 for 2x-x^2 Length: nmul.
Referenced by mat::TC2< Treal, Tmatrix >::TC2(), and mat::TC2< Treal, Tmatrix >::~TC2().
|
protected |
The difference between the trace of the matrix and the number of occupied orbitals before each step.
Length: nmul + 1
Referenced by mat::TC2< Treal, Tmatrix >::TC2(), and mat::TC2< Treal, Tmatrix >::~TC2().
|
protected |
Fock / Kohn-Sham matrix at initialization.
Then used as workspace by purify(). Empty after call to purify().
Referenced by mat::TC2< Treal, Tmatrix >::TC2().