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_scfops_h
00029 #define _chemistry_qc_scf_scfops_h
00030
00031 #ifdef __GNUC__
00032 #pragma interface
00033 #endif
00034
00035 #include <math/scmat/elemop.h>
00036 #include <math/scmat/blocked.h>
00037
00038 #include <chemistry/qc/scf/scf.h>
00039
00040 namespace sc {
00041
00042 class SCFEnergy : public SCElementOp2 {
00043 private:
00044 double eelec;
00045 int deferred_;
00046
00047 public:
00048 SCFEnergy();
00049 ~SCFEnergy();
00050
00051 int has_collect();
00052 void defer_collect(int h);
00053 void collect(const Ref<MessageGrp>&grp);
00054 double result();
00055 void reset();
00056
00057 void process(SCMatrixBlockIter&i, SCMatrixBlockIter&j);
00058 };
00059
00060 class LevelShift : public BlockedSCElementOp {
00061 protected:
00062 SCF *scf_;
00063 double shift;
00064
00065 public:
00066 LevelShift(SCF*);
00067 ~LevelShift();
00068
00069 int has_side_effects();
00070 void set_shift(double);
00071
00072 void process(SCMatrixBlockIter&);
00073 };
00074
00075 class ALevelShift : public LevelShift {
00076 public:
00077 ALevelShift(SCF*);
00078 ~ALevelShift();
00079 void process(SCMatrixBlockIter&);
00080 };
00081
00082 class BLevelShift : public LevelShift {
00083 public:
00084 BLevelShift(SCF*);
00085 ~BLevelShift();
00086 void process(SCMatrixBlockIter&);
00087 };
00088
00089
00090
00091
00092
00093
00094
00095
00096
00097 class MOLagrangian : public BlockedSCElementOp2 {
00098 private:
00099 SCF *scf_;
00100
00101 public:
00102 MOLagrangian(SCF* s);
00103 ~MOLagrangian();
00104
00105 int has_side_effects();
00106
00107 void process(SCMatrixBlockIter& bi1, SCMatrixBlockIter& bi2);
00108 };
00109
00110 }
00111
00112 #endif
00113
00114
00115
00116
00117