--- branches/development/src/parallel/ForceDecomposition.hpp	2011/05/26 13:55:04	1569
+++ branches/development/src/parallel/ForceDecomposition.hpp	2011/06/08 16:05:07	1576
@@ -45,10 +45,15 @@
 #include "brains/SimInfo.hpp"
 #include "brains/SnapshotManager.hpp"
 #include "nonbonded/NonBondedInteraction.hpp"
+#include "nonbonded/Cutoffs.hpp"
+#include "nonbonded/InteractionManager.hpp"
+#include "utils/Tuple.hpp"
 
 using namespace std;
 namespace OpenMD {
   
+  typedef tuple3<RealType, RealType, RealType> groupCutoffs;
+
   /**
    * @class ForceDecomposition 
    *
@@ -91,16 +96,22 @@ namespace OpenMD {
     
     virtual void distributeInitialData() = 0;
     virtual void distributeData() = 0;
+    virtual void zeroWorkArrays() = 0;
     virtual void collectIntermediateData() = 0;
     virtual void distributeIntermediateData() = 0;
     virtual void collectData() = 0;
+    virtual potVec getLongRangePotential() { return longRangePot_; }
 
     // neighbor list routines
     virtual bool checkNeighborList();
     virtual vector<pair<int, int> >  buildNeighborList() = 0;
 
+    // how to handle cutoffs:
+    void setCutoffPolicy(CutoffPolicy cp) {cutoffPolicy_ = cp;}
+    void setUserCutoff(RealType rcut) {userCutoff_ = rcut; userChoseCutoff_ = true; }
+
     // group bookkeeping
-    virtual pair<int, int> getGroupTypes(int cg1, int cg2) = 0;
+    virtual groupCutoffs getGroupCutoffs(int cg1, int cg2) = 0;
 
     // Group->atom bookkeeping
     virtual vector<int> getAtomsInGroupRow(int cg1) = 0; 
@@ -116,14 +127,15 @@ namespace OpenMD {
     virtual Vector3d getInteratomicVector(int atom1, int atom2) = 0;
        
     // atom bookkeeping
-    virtual vector<int> getAtomList() = 0;
-    virtual vector<int> getSkipsForAtom(int atom1) = 0;
+    virtual int getNAtomsInRow() = 0;
+    virtual vector<int> getSkipsForRowAtom(int atom1) = 0;
     virtual bool skipAtomPair(int atom1, int atom2) = 0;
     virtual void addForceToAtomRow(int atom1, Vector3d fg) = 0;
     virtual void addForceToAtomColumn(int atom2, Vector3d fg) = 0;
 
     // filling interaction blocks with pointers
     virtual InteractionData fillInteractionData(int atom1, int atom2) = 0;
+    virtual void unpackInteractionData(InteractionData idat, int atom1, int atom2) = 0;
     virtual InteractionData fillSkipData(int atom1, int atom2) = 0;
     virtual SelfData fillSelfData(int atom1);
     
@@ -131,8 +143,12 @@ namespace OpenMD {
     SimInfo* info_;   
     SnapshotManager* sman_;    
     Snapshot* snap_;
+    ForceField* ff_;
+    InteractionManager* interactionMan_;
+
     int storageLayout_;
     RealType skinThickness_;   /**< Verlet neighbor list skin thickness */    
+    RealType largestRcut_;
 
     map<pair<int, int>, int> topoDist; //< topoDist gives the
                                        //topological distance between
@@ -156,11 +172,19 @@ namespace OpenMD {
                                        //method to fill this.
 
     vector<vector<int> > groupList_;
+
     vector<Vector3i> cellOffsets_;
     Vector3i nCells_;
     vector<vector<int> > cellList_;
     vector<Vector3d> saved_CG_positions_;
+    potVec longRangePot_;
 
+    bool userChoseCutoff_;
+    RealType userCutoff_;
+    CutoffPolicy cutoffPolicy_;
+
+    map<pair<int, int>, tuple3<RealType, RealType, RealType> > gTypeCutoffMap;
+
   };    
 }
 #endif