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_scf_clscf_h
00029 #define _chemistry_qc_scf_clscf_h
00030
00031 #ifdef __GNUC__
00032 #pragma interface
00033 #endif
00034
00035 #include <chemistry/qc/scf/scf.h>
00036
00037 namespace sc {
00038
00039
00040
00043 class CLSCF: public SCF {
00044 protected:
00045 Ref<PointGroup> most_recent_pg_;
00046 int user_occupations_;
00047 int tndocc_;
00048 int nirrep_;
00049 int *initial_ndocc_;
00050 int *ndocc_;
00051
00052 ResultRefSymmSCMatrix cl_fock_;
00053
00054 public:
00055 CLSCF(StateIn&);
00070 CLSCF(const Ref<KeyVal>&);
00071 ~CLSCF();
00072
00073 void save_data_state(StateOut&);
00074
00075 void print(std::ostream&o=ExEnv::out0()) const;
00076
00077 double occupation(int irrep, int vectornum);
00078
00079 int n_fock_matrices() const;
00084 RefSymmSCMatrix fock(int i);
00086 RefSymmSCMatrix effective_fock();
00087
00088 RefSymmSCMatrix density();
00089
00090 void symmetry_changed();
00091
00092
00093 int spin_polarized();
00094
00095 protected:
00096
00097 RefSymmSCMatrix cl_dens_;
00098 RefSymmSCMatrix cl_dens_diff_;
00099 RefSymmSCMatrix cl_gmat_;
00100
00101 void set_occupations(const RefDiagSCMatrix& evals);
00102
00103
00104 void init_vector();
00105 void done_vector();
00106 void reset_density();
00107 double new_density();
00108 double scf_energy();
00109
00110 Ref<SCExtrapData> extrap_data();
00111
00112
00113 void init_gradient();
00114 void done_gradient();
00115
00116 RefSymmSCMatrix lagrangian();
00117 RefSymmSCMatrix gradient_density();
00118
00119
00120 void init_hessian();
00121 void done_hessian();
00122
00123
00124 void two_body_deriv_hf(double*grad,double exchange_fraction);
00125 };
00126
00127 }
00128
00129 #endif
00130
00131
00132
00133
00134