#include <TC2.h>
Public Member Functions | |
TC2 (Tmatrix &F, Tmatrix &DM, const int size, const int noc, const Treal trunc=0, const int maxmm=100) | |
Constructor. | |
Treal | fermi_level (Treal tol=1e-15) const |
Returns the Fermi level. | |
Treal | homo (Treal tol=1e-15) const |
Returns upper bound of the HOMO eigenvalue. | |
Treal | lumo (Treal tol=1e-15) const |
Returns lower bound of the LUMO eigenvalue. | |
int | n_multiplies () const |
Returns the number of used matrix matrix multiplications. | |
void | print_data (int const start, int const stop) const |
virtual | ~TC2 () |
Destructor. | |
Protected Member Functions | |
void | purify () |
Runs purification. | |
Protected Attributes | |
Tmatrix & | X |
Fock / Kohn-Sham matrix at initialization. | |
Tmatrix & | D |
Density matrix after purification. | |
const int | n |
System size. | |
const int | nocc |
Number of occupied orbitals. | |
const Treal | frob_trunc |
Threshold for the truncation. | |
const int | maxmul |
Number of tolerated matrix multiplications. | |
Treal | lmin |
Lower bound for eigenvalue spectrum. | |
Treal | lmax |
Upper bound for eigenvalue spectrum. | |
int | nmul |
Number of used matrix multiplications. | |
int | nmul_firstpart |
Number of used matrix multiplications in the first part of the purification. | |
Treal * | idemerror |
Upper bound of euclidean norm ||D-D^2||_2 before each step. | |
Treal * | tracediff |
The difference between the trace of the matrix and the number of occupied orbitals before each step. | |
int * | polys |
Choices of polynomials 0 for x^2 and 1 for 2x-x^2 Length: nmul. | |
Classes | |
class | Fun |
Help class for bisection root finding calls. More... |
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.
< 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. |
virtual mat::TC2< Treal, Tmatrix >::~TC2 | ( | ) | [inline, virtual] |
Destructor.
Treal mat::TC2< Treal, Tmatrix >::fermi_level | ( | Treal | tol = 1e-15 |
) | const |
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. |
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. |
int mat::TC2< Treal, Tmatrix >::n_multiplies | ( | ) | const [inline] |
Returns the number of used matrix matrix multiplications.
void mat::TC2< Treal, Tmatrix >::print_data | ( | int const | start, | |
int const | stop | |||
) | const |
void mat::TC2< Treal, Tmatrix >::purify | ( | ) | [protected] |
Runs purification.
Run by constructor.
Density matrix after purification.
const Treal mat::TC2< Treal, Tmatrix >::frob_trunc [protected] |
Threshold for the truncation.
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
Upper bound for eigenvalue spectrum.
Lower bound for eigenvalue spectrum.
Number of tolerated matrix multiplications.
System size.
Number of used matrix multiplications.
int mat::TC2< Treal, Tmatrix >::nmul_firstpart [protected] |
Number of used matrix multiplications in the first part of the purification.
Number of occupied orbitals.
Choices of polynomials 0 for x^2 and 1 for 2x-x^2 Length: nmul.
The difference between the trace of the matrix and the number of occupied orbitals before each step.
Length: nmul + 1