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_scf_h
00029 #define _chemistry_qc_scf_scf_h
00030
00031 #ifdef __GNUC__
00032 #pragma interface
00033 #endif
00034
00035 #include <util/group/thread.h>
00036
00037 #include <math/optimize/scextrap.h>
00038
00039 #include <chemistry/qc/basis/tbint.h>
00040 #include <chemistry/qc/wfn/accum.h>
00041 #include <chemistry/qc/wfn/obwfn.h>
00042
00043 namespace sc {
00044
00045
00046
00049 class SCF: public OneBodyWavefunction {
00050 protected:
00051 int need_vec_;
00052 int compute_guess_;
00053
00054 int keep_guess_wfn_;
00055 Ref<OneBodyWavefunction> guess_wfn_;
00056
00057 int always_use_guess_wfn_;
00058
00059 Ref<SelfConsistentExtrapolation> extrap_;
00060
00061 Ref<AccumH> accumdih_;
00062 Ref<AccumH> accumddh_;
00063
00064 int maxiter_;
00065 int dens_reset_freq_;
00066 int reset_occ_;
00067 int local_dens_;
00068 size_t storage_;
00069 int print_all_evals_;
00070 int print_occ_evals_;
00071
00072 double level_shift_;
00073
00074 Ref<MessageGrp> scf_grp_;
00075 Ref<ThreadGrp> threadgrp_;
00076 int local_;
00077
00078 Ref<TwoBodyInt>* tbis_;
00079 virtual void init_threads();
00080 virtual void done_threads();
00081
00082
00083 virtual void compute();
00084
00085
00086 virtual double compute_vector(double&);
00087
00088
00089 virtual Ref<SCExtrapError> extrap_error();
00090
00091
00092 virtual void compute_gradient(const RefSCVector&);
00093
00094
00095 virtual void compute_hessian(const RefSymmSCMatrix&);
00096
00097
00098
00099 virtual void savestate_iter(int);
00100
00101
00102 signed char * init_pmax(double *);
00103
00104
00105
00106
00107 enum Access { Read, Write, Accum };
00108 RefSymmSCMatrix get_local_data(const RefSymmSCMatrix&, double*&, Access);
00109
00110
00111
00112
00113
00114 virtual void initial_vector(int needv=1);
00115
00116
00117
00118
00119 void init_mem(int);
00120
00121 void so_density(const RefSymmSCMatrix& d, double occ, int alp=1);
00122
00123
00124
00125 int *read_occ(const Ref<KeyVal> &, const char *name, int nirrep);
00126 public:
00127 SCF(StateIn&);
00180 SCF(const Ref<KeyVal>&);
00181 ~SCF();
00182
00183 void save_data_state(StateOut&);
00184
00185 RefSCMatrix oso_eigenvectors();
00186 RefDiagSCMatrix eigenvalues();
00187
00188 int spin_unrestricted();
00189
00190
00191 virtual int n_fock_matrices() const =0;
00192
00193
00194 virtual RefSymmSCMatrix fock(int) =0;
00195
00196
00197 virtual RefSymmSCMatrix effective_fock() =0;
00198
00199 virtual double one_body_energy();
00200 virtual void two_body_energy(double &ec, double &ex);
00201
00202 void symmetry_changed();
00203
00204 void obsolete();
00205
00206 void print(std::ostream&o=ExEnv::out0()) const;
00207
00208 protected:
00209
00210 RefSCMatrix oso_scf_vector_;
00211 RefSCMatrix oso_scf_vector_beta_;
00212 RefSymmSCMatrix hcore_;
00213
00214
00215
00216
00217
00218 virtual void set_occupations(const RefDiagSCMatrix&) =0;
00219
00220
00221
00222 virtual void init_vector() =0;
00223 virtual void done_vector() =0;
00224
00225
00226 virtual double new_density() =0;
00227
00228
00229 virtual void reset_density() =0;
00230
00231
00232 virtual double scf_energy() =0;
00233
00234
00235 virtual Ref<SCExtrapData> extrap_data() =0;
00236
00237
00238 virtual void ao_fock(double accuracy) =0;
00239
00240
00241
00242 virtual void init_gradient() =0;
00243 virtual void done_gradient() =0;
00244
00245 virtual RefSymmSCMatrix lagrangian() =0;
00246 virtual RefSymmSCMatrix gradient_density() =0;
00247 virtual void two_body_deriv(double*) =0;
00248
00249
00250
00251 virtual void init_hessian() =0;
00252 virtual void done_hessian() =0;
00253 };
00254
00255 }
00256
00257 #endif
00258
00259
00260
00261
00262