SYNOPSIS
#include <density.h>
Inherits sc::Volume.
Public Types
enum { X = 0, Y = 1, Z = 2 }
This gives the elements of gradient arrays.
enum { XX = 0, YX = 1, YY = 2, ZX = 3, ZY = 4, ZZ = 5 }
This gives the elements of hessian arrays.
Public Member Functions
BatchElectronDensity (const Ref< KeyVal > &)
BatchElectronDensity (const Ref< Wavefunction > &, double accuracy=DBL_EPSILON)
BatchElectronDensity (const Ref< BatchElectronDensity > &d, bool reference_parent_data=false)
This will construct copies of this.
virtual void boundingbox (double valuemin, double valuemax, SCVector3 &p1, SCVector3 &p2)
Returns the bounding box.
void clear ()
This will cause all stratch storage to be released.
void compute_density (const SCVector3 &r, double *alpha_density, double *alpha_density_grad, double *alpha_density_hessian, double *beta_density, double *beta_density_grad, double *beta_density_hessian)
This is a alternate to the Volume interface that avoids some of the overhead of that interface.
virtual void init (bool initialize_density_matrices=true)
This is called to finish initialization of the object.
virtual void set_densities (const RefSymmSCMatrix &aden, const RefSymmSCMatrix &bden)
This will fill in the internel copies of the density matrices with new values.
void set_linear_scaling (bool b)
Turn linear scaling algorithm on/off.
void set_accuracy (double a)
Sets the accuracy.
void set_use_dmat_bound (bool b)
Turn use of density matrix bounds on/off.
DFT Support Members.
These return some of the internal data, some of which is only value after a density has been computed.
This data is needed by the density functional theory code.
double * alpha_density_matrix ()
Return the alpha density matrix.
double * beta_density_matrix ()
Return the beta density matrix.
int ncontrib ()
int * contrib ()
int ncontrib_bf ()
int * contrib_bf ()
double * bs_values ()
double * bsg_values ()
double * bsh_values ()
void set_need_basis_gradient (bool b)
To ensure that that the basis functions gradients are computed, use this.
void set_need_basis_hessian (bool b)
Protected Member Functions
virtual void init_common_data (bool initialize_density_matrices)
virtual void init_scratch_data ()
void compute_basis_values (const SCVector3 &r)
void compute_spin_density (const double *dmat, double *rho, double *grad, double *hess)
virtual void compute ()
Recompute at least the results that have compute true and are not already computed.
Protected Attributes
Ref< Wavefunction > wfn_
Ref< GaussianBasisSet > basis_
double * alpha_dmat_
double * beta_dmat_
double * dmat_bound_
ShellExtent * extent_
GaussianBasisSet::ValueData * valdat_
int ncontrib_
int * contrib_
int ncontrib_bf_
int * contrib_bf_
double * bs_values_
double * bsg_values_
double * bsh_values_
int nshell_
int nbasis_
int spin_polarized_
int linear_scaling_
int use_dmat_bound_
bool need_hessian_
bool need_gradient_
bool need_basis_hessian_
bool need_basis_gradient_
bool using_shared_data_
double accuracy_
Detailed Description
This a more highly optimized than ElectronDensity since everything is precomputed.
However, it cannot be used if the density and/or geometry might change between computations of the density or bounding box, unless the obsolete member is called.
Member Enumeration Documentation
anonymous enum
This gives the elements of hessian arrays.
anonymous enum
This gives the elements of gradient arrays.
Constructor & Destructor Documentation
sc::BatchElectronDensity::BatchElectronDensity (const Ref< BatchElectronDensity > & d, bool reference_parent_data = false)
This will construct copies of this. If reference_parent_data is true, then data that do not change, such as the density matrices and shell extent, are referenced rather than copied. In this case, the original object that allocated this items must be valid while copied objects are used to compute densities. Also d must have already been intialized and the resulting copy is already initialized (and cannot be reinitialized).
If reference_parent_data is false, then init must be called on this object before it is used.
Member Function Documentation
double* sc::BatchElectronDensity::alpha_density_matrix () [inline]
Return the alpha density matrix.
double* sc::BatchElectronDensity::beta_density_matrix () [inline]
Return the beta density matrix.
virtual void sc::BatchElectronDensity::boundingbox (double valuemin, double valuemax, SCVector3 & p1, SCVector3 & p2) [virtual]
Returns the bounding box.
Implements sc::Volume.
void sc::BatchElectronDensity::clear ()
This will cause all stratch storage to be released.
virtual void sc::BatchElectronDensity::compute () [protected, virtual]
Recompute at least the results that have compute true and are not already computed. This should only be called by Result's members.
Implements sc::Volume.
void sc::BatchElectronDensity::compute_density (const SCVector3 & r, double * alpha_density, double * alpha_density_grad, double * alpha_density_hessian, double * beta_density, double * beta_density_grad, double * beta_density_hessian)
This is a alternate to the Volume interface that avoids some of the overhead of that interface.
virtual void sc::BatchElectronDensity::init (bool initialize_density_matrices = true) [virtual]
This is called to finish initialization of the object. It must not be called with objects that created in a way that they share parent data, those objects are initialized when they are constructed. This member is usually called automatically, however, if it will be used to initial other objects that share parent data, then it must be initialized first and this return is the way to do that. If initialize_density_matrices is false, then the density matrices will be allocated, but not filled in. They must be later filled in with set_densities.
void sc::BatchElectronDensity::set_accuracy (double a) [inline]
Sets the accuracy.
virtual void sc::BatchElectronDensity::set_densities (const RefSymmSCMatrix & aden, const RefSymmSCMatrix & bden) [virtual]
This will fill in the internel copies of the density matrices with new values. aden is the alpha density matrix and bden is the beta density matrix. bden is ignored if the wavefunction is not spin polarized.
void sc::BatchElectronDensity::set_linear_scaling (bool b) [inline]
Turn linear scaling algorithm on/off. The effect of this will be delayed until the next time init() is called.
void sc::BatchElectronDensity::set_need_basis_gradient (bool b) [inline]
To ensure that that the basis functions gradients are computed, use this.
void sc::BatchElectronDensity::set_use_dmat_bound (bool b) [inline]
Turn use of density matrix bounds on/off.
Author
Generated automatically by Doxygen for MPQC from the source code.