ViewVC Help
View File | Revision Log | Show Annotations | View Changeset | Root Listing
root/OpenMD/branches/development/src/nonbonded/EAM.hpp
(Generate patch)

Comparing branches/development/src/nonbonded/EAM.hpp (file contents):
Revision 1479 by gezelter, Mon Jul 26 19:00:48 2010 UTC vs.
Revision 1505 by gezelter, Sun Oct 3 22:18:59 2010 UTC

# Line 42 | Line 42
42   #ifndef NONBONDED_EAM_HPP
43   #define NONBONDED_EAM_HPP
44  
45 + #include "nonbonded/NonBondedInteraction.hpp"
46   #include "types/AtomType.hpp"
47   #include "UseTheForce/ForceField.hpp"
48   #include "math/Vector3.hpp"
# Line 67 | Line 68 | namespace OpenMD {
68      eamUnknown
69    };
70    
71 <  class EAM {
71 >  class EAM : public MetallicInteraction {
72      
73    public:    
74 <    static EAM* Instance();
75 <    static void setForceField(ForceField *ff) {forceField_ = ff;};
76 <    static void initialize();
77 <    static void addType(AtomType* atomType);
78 <    static void addExplicitInteraction(AtomType* atype1, AtomType* atype2, RealType dr, int nr, std::vector<RealType> phiAB);
79 <    
80 <    static void calcDensity(AtomType* at1, AtomType* at2, const RealType rij, RealType &rho_i_at_j, RealType &rho_j_at_i);
74 >    EAM();
75 >    void setForceField(ForceField *ff) {forceField_ = ff;};
76 >    void addType(AtomType* atomType);
77 >    void addExplicitInteraction(AtomType* atype1, AtomType* atype2, RealType dr, int nr, std::vector<RealType> phiAB);
78 >    void calcDensity(DensityData ddat);
79 >    void calcFunctional(FunctionalData fdat);
80 >    void calcForce(InteractionData idat);
81 >    virtual string getName() { return name_; }
82 >    virtual RealType getSuggestedCutoffRadius(AtomType* at1, AtomType* at2);
83  
81 –    static void calcFunctional(AtomType* at1, const RealType rho, RealType &frho, RealType &dfrhodrho);
82 –
83 –    static void calcForce(AtomType* at1, AtomType* at2, const Vector3d d, const RealType rij, const RealType r2, const RealType sw, RealType &vpair, RealType &pot, Vector3d &f1, const RealType rho1, const RealType rho2, const RealType dfrho1, const RealType dfrho2, RealType &fshift1, RealType &fshift2);
84 –                          
85 –    // Fortran support routines;
86 –    static void calc_eam_prepair_rho(int *atid1, int *atid2, RealType *rij,  RealType* rho_i_at_j, RealType* rho_j_at_i);
87 –    static void calc_eam_preforce_Frho(int *atid1, RealType* rho, RealType* frho, RealType* dfrhodrho);
88 –    static void do_eam_pair(int *atid1, int *atid2, RealType *d, RealType *rij, RealType *r2,  RealType *sw, RealType *vpair, RealType *pot, RealType *f1, RealType* rho1, RealType* rho2, RealType* dfrho1, RealType* dfrho2, RealType* fshift1, RealType* fshift2);
89 –    static void setCutoffEAM(RealType *thisRcut);
90 –    static RealType getEAMcut(int *atid1);
91 –
84    private:
85 <    virtual ~EAM() { }
86 <    // singleton pattern, prevent reconstruction
87 <    EAM() { }
88 <    EAM(EAM const &) {};
89 <    EAM& operator=(EAM const&) {};
90 <    static EAM* _instance;
91 <  
100 <    static EAMParam getEAMParam(AtomType* atomType);
101 <    static CubicSpline* getZ(AtomType* atomType);
102 <    static CubicSpline* getRho(AtomType* atomType);
103 <    static CubicSpline* getF(AtomType* atomType);
104 <    static RealType getRcut(AtomType* atomType);
105 <    static CubicSpline* getPhi(AtomType* atomType1, AtomType* atomType2);
85 >    void initialize();  
86 >    EAMParam getEAMParam(AtomType* atomType);
87 >    CubicSpline* getZ(AtomType* atomType);
88 >    CubicSpline* getRho(AtomType* atomType);
89 >    CubicSpline* getF(AtomType* atomType);
90 >    RealType getRcut(AtomType* atomType);
91 >    CubicSpline* getPhi(AtomType* atomType1, AtomType* atomType2);
92      
93 <    static bool initialized_;
94 <    static std::map<int, AtomType*> EAMlist;
95 <    static std::map<AtomType*, EAMAtomData> EAMMap;
96 <    static std::map<std::pair<AtomType*, AtomType*>, EAMInteractionData> MixingMap;
97 <    static ForceField* forceField_;
98 <    static RealType eamRcut_;
99 <    static EAMMixingMethod mixMeth_;    
93 >    bool initialized_;
94 >    std::map<int, AtomType*> EAMlist;
95 >    std::map<AtomType*, EAMAtomData> EAMMap;
96 >    std::map<std::pair<AtomType*, AtomType*>, EAMInteractionData> MixingMap;
97 >    ForceField* forceField_;
98 >    RealType eamRcut_;
99 >    EAMMixingMethod mixMeth_;    
100 >    string name_;
101      
102    };
103   }

Diff Legend

– Removed lines
+ Added lines
< Changed lines
> Changed lines