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 #ifdef __GNUG__
00029 #pragma interface
00030 #endif
00031
00032 #ifndef _chemistry_qc_mbptr12_twobodygrid_h
00033 #define _chemistry_qc_mbptr12_twobodygrid_h
00034
00035 #include <vector>
00036 #include <util/ref/ref.h>
00037 #include <util/state/statein.h>
00038 #include <util/state/stateout.h>
00039 #include <math/scmat/vector3.h>
00040 #include <chemistry/molecule/molecule.h>
00041
00042 using namespace std;
00043
00044 namespace sc {
00045
00049 class TwoBodyGrid : virtual public SavableState {
00050
00051 std::string name_;
00052 vector<SCVector3> r1_;
00053 vector<SCVector3> r2_;
00054 SCVector3 O_;
00055
00056 public:
00057 TwoBodyGrid(StateIn&);
00060 TwoBodyGrid(const Ref<KeyVal>&);
00061 ~TwoBodyGrid();
00062
00063 void save_data_state(StateOut&);
00064
00065 const std::string& name() const;
00066 int nelem() const;
00067 const SCVector3& origin() const;
00068
00070 SCVector3 xyz1(int i, const SCVector3& O=SCVector3(0.0)) const;
00072 SCVector3 xyz2(int i, const SCVector3& O=SCVector3(0.0)) const;
00074 SCVector3 rtp1(int i, const SCVector3& O=SCVector3(0.0)) const;
00076 SCVector3 rtp2(int i, const SCVector3& O=SCVector3(0.0)) const;
00078 void print(std::ostream&o=ExEnv::out0()) const;
00079
00080 };
00081
00082
00083 }
00084
00085 #endif
00086
00087
00088
00089
00090
00091
00092