00001
00002
00003
00004
00005
00006
00007
00008
00009
00010
00011
00012
00013
00014
00015
00016
00017
00018
00019
00020
00021
00022
00023
00024
00025
00026
00027
00028 #ifndef _chemistry_qc_dft_clks_h
00029 #define _chemistry_qc_dft_clks_h
00030
00031 #ifdef __GNUC__
00032 #pragma interface
00033 #endif
00034
00035 #include <chemistry/qc/scf/clscf.h>
00036 #include <chemistry/qc/dft/integrator.h>
00037 #include <chemistry/qc/dft/functional.h>
00038
00039 namespace sc {
00040
00041
00042
00043 class CLKS: public CLSCF {
00044 protected:
00045 Ref<DenIntegrator> integrator_;
00046 Ref<DenFunctional> functional_;
00047 RefSymmSCMatrix vxc_;
00048
00049 public:
00050 CLKS(StateIn&);
00051 CLKS(const Ref<KeyVal>&);
00052 ~CLKS();
00053
00054 void save_data_state(StateOut&);
00055
00056 void print(std::ostream&o=ExEnv::out0()) const;
00057
00058 void two_body_energy(double &ec, double &ex);
00059
00060 int value_implemented() const;
00061 int gradient_implemented() const;
00062
00063 RefSymmSCMatrix density();
00064 protected:
00065 void ao_fock(double accuracy);
00066 double exc_;
00067 double scf_energy();
00068 Ref<SCExtrapData> extrap_data();
00069 RefSymmSCMatrix effective_fock();
00070
00071 void init_vector();
00072 void done_vector();
00073
00074 void two_body_deriv(double*);
00075 };
00076
00077 }
00078
00079 #endif
00080
00081
00082
00083
00084