GlobalMasterServer Class Reference

#include <GlobalMasterServer.h>

Public Member Functions
	GlobalMasterServer (ComputeMgr *m, int theNumDataSenders)
virtual	~GlobalMasterServer ()
void	recvData (ComputeGlobalDataMsg *)
void	addClient (GlobalMaster *newClient)

Detailed Description

Definition at line 22 of file GlobalMasterServer.h.

Constructor & Destructor Documentation

GlobalMasterServer::GlobalMasterServer	(	ComputeMgr *	m,
		int	theNumDataSenders
	)

Definition at line 580 of file GlobalMasterServer.C.

References DebugM, SimParameters::lattice, Node::Object(), and Node::simParameters.

                                                              {
  DebugM(3,"Constructing GlobalMasterServer\n");
  myComputeManager = m;
  numDataSenders = 0; // theNumDataSenders;
  numForceSenders = 0;
  latticeCount = 0;
  lattice = Node::Object()->simParameters->lattice;
  recvCount = 0; /* we haven't gotten any messages yet */
  recvPatchesCount = 0;
  firstTime = 1; /* XXX temporary */
  step = -1;
  totalAtomsRequested = 0;
  totalGroupsRequested = 0;
  forceSendEnabled = 0;
  if ( Node::Object()->simParameters->tclForcesOn ) forceSendEnabled = 1;
  if ( Node::Object()->simParameters->colvarsOn ) forceSendEnabled = 1;
  DebugM(3,"done constructing\n");
}

GlobalMasterServer::~GlobalMasterServer ( )

virtual

Definition at line 600 of file GlobalMasterServer.C.

                                        {
  GlobalMaster *m_i = *clientList.begin();
  GlobalMaster *m_e = *clientList.end();
  
  /* delete each of the masters */
  while(m_i != m_e) {
    delete m_i;
    m_i++;
  }
}

Member Function Documentation

void GlobalMasterServer::addClient

(

GlobalMaster *

newClient

)

Definition at line 31 of file GlobalMasterServer.C.

References DebugM, and GlobalMaster::setLattice().

Referenced by ComputeMgr::createComputes().

                                                          {
  DebugM(3,"Adding client\n");
  clientList.add(newClient);
  newClient->setLattice(&lattice);
  DebugM(2,"Added.\n");
}

void GlobalMasterServer::recvData

(

ComputeGlobalDataMsg *

msg

)

Definition at line 38 of file GlobalMasterServer.C.

References ComputeGlobalDataMsg::aid, ResizeArray< T >::begin(), ComputeGlobalDataMsg::count, DebugM, ResizeArray< T >::end(), endi(), ComputeGlobalDataMsg::fid, ComputeGlobalDataMsg::gcom, ComputeGlobalDataMsg::gmass, ComputeGlobalDataMsg::gridobjindex, ComputeGlobalDataMsg::gridobjvalue, ComputeGlobalDataMsg::gtf, ComputeGlobalDataMsg::lat, NAMD_bug(), PatchMap::numPatches(), numPatches, PatchMap::Object(), ComputeGlobalDataMsg::p, ComputeGlobalDataMsg::patchcount, ResizeArray< T >::size(), ComputeGlobalDataMsg::step, and ComputeGlobalDataMsg::tf.

Referenced by ComputeMgr::recvComputeGlobalData().

                                                           {
  DebugM(3,"Storing data (" << msg->aid.size() << " positions) on master\n");

  if ( msg->step != -1 ) step = msg->step;

  /* get the beginning and end of the lists */
  AtomIDList::iterator a_i = msg->aid.begin();
  AtomIDList::iterator a_e = msg->aid.end();
  PositionList::iterator p_i = msg->p.begin();
  PositionList::iterator g_i = msg->gcom.begin();
  PositionList::iterator g_e = msg->gcom.end();
  BigRealList::iterator gm_i = msg->gmass.begin();

  /* iterate over each member of the atom lists */
  for ( ; a_i != a_e; ++a_i, ++p_i ) {
    receivedAtomIDs.add(*a_i);
    receivedAtomPositions.add(*p_i);
  }
  
  /* iterate over each member of "total force" lists */
  a_e = msg->fid.end();
  ForceList::iterator f_i=msg->tf.begin();
  for (a_i=msg->fid.begin() ; a_i!=a_e; ++a_i,++f_i) {
    receivedForceIDs.add(*a_i);
    receivedTotalForces.add(*f_i);
  }

  /* iterate over each member of the group position list */
  int i=0;
  for ( ; g_i != g_e; ++g_i, ++gm_i ) {
    DebugM(1,"Received center of mass "<<*g_i<<"\n");
    if(i >= totalGroupsRequested) NAMD_bug("Received too many groups.");
    receivedGroupPositions[i] += (*g_i);
    receivedGroupMasses[i] += (*gm_i);
    i++;
  }
  if(i!=totalGroupsRequested) NAMD_bug("Received too few groups.");

  /* iterate over each member of group total force lists */
  int ntf = msg->gtf.size();
  if ( ntf && ntf != receivedGroupTotalForces.size() ) NAMD_bug("Received wrong number of group forces.");
  ForceList::iterator gf_i=msg->gtf.begin();
  ForceList::iterator gf_e=msg->gtf.end();
  for ( i=0 ; gf_i != gf_e; ++gf_i, ++i ) {
    receivedGroupTotalForces[i] += (*gf_i);
  }

  // Get values of the GridForce objects
  int ngov = msg->gridobjvalue.size();
  IntList::iterator goi_i = msg->gridobjindex.begin();
  BigRealList::iterator gov_i = msg->gridobjvalue.begin();
  BigRealList::iterator gov_e = msg->gridobjvalue.end();
  for ( i=0 ; gov_i != gov_e; gov_i++, goi_i++, i++ ) {
    receivedGridObjIndices[i] = (*goi_i);
    receivedGridObjValues[i] += (*gov_i);
  }
  if (ngov && ngov != receivedGridObjValues.size()) {
    NAMD_bug("Received wrong number of grid objects.");
  }

  if ( msg->lat.size() ) {
    if ( latticeCount ) {
      NAMD_bug("GlobalMasterServer::recvData received lattice twice.");
    }
    lattice = msg->lat[0];
    latticeCount = 1;
  }

  recvCount += msg->count;
  recvPatchesCount += msg->patchcount;
  int numPatches = PatchMap::Object()->numPatches();

  /* done with the message, delete it */
  delete msg;

  /* check whether we've gotten all the expected messages */

  DebugM(3, "received " << recvCount << " atom messages and "
         << recvPatchesCount << " patch messages (out of " << numPatches
         << " patches) from the ComputeGlobals\n" << endi);

  int numNeededPatches = totalGridObjsRequested ? numPatches : 0;
  if (firstTime) numNeededPatches = 0;

  if (recvCount > numDataSenders + numForceSenders + 1) {
    NAMD_bug("GlobalMasterServer::recvData recvCount too high.");
  }

  if (recvPatchesCount > numPatches) {
    NAMD_bug("GlobalMasterServer::recvData too many patches received.");
  }

  if ((recvCount == numDataSenders + numForceSenders + 1) &&
      (recvPatchesCount == numNeededPatches)) {

    if ( ! latticeCount ) {
      NAMD_bug("GlobalMasterServer::recvData did not receive lattice.");
    }

    DebugM(3, "Received all needed data\n" << endi);

#ifdef DEBUGM
    for (size_t ig = 0; ig < receivedGroupMasses.size(); ig++) {
      // The groups' positions will be normalized later, in callClients()
      DebugM(3, "Group mass " << ig << " = "
             << receivedGroupMasses[ig] << "\n" << endi);
    }
    for (size_t ig = 0; ig < receivedGridObjValues.size(); ig++) {
      DebugM(3, "Grid object " << receivedGridObjIndices[ig] << " = "
             << receivedGridObjValues[ig] << "\n" << endi);
    }
#endif

    int oldTotalGroupsRequested = totalGroupsRequested;

    int resendCoordinates = callClients();

    if (resendCoordinates) {
      DebugM(3, "resendCoordinates\n" << endi);
    }

    /* now restart */
    step = -1;
    receivedAtomIDs.resize(0);
    receivedAtomPositions.resize(0);
    receivedGroupPositions.resize(totalGroupsRequested);
    receivedGroupPositions.setall(Vector(0,0,0));
    receivedGroupMasses.resize(totalGroupsRequested);
    receivedGroupMasses.setall(0);
    receivedGridObjIndices.resize(totalGridObjsRequested);
    receivedGridObjIndices.setall(-1);
    receivedGridObjValues.resize(totalGridObjsRequested);
    receivedGridObjValues.setall(0);
    latticeCount = 0;
    recvCount = 0;
    recvPatchesCount = 0;
    if ( resendCoordinates ) {
      recvCount += numForceSenders;
    } else {
      receivedGroupTotalForces.resize(totalGroupsRequested);
      receivedGroupTotalForces.setall(0);
      receivedForceIDs.resize(0);
      receivedTotalForces.resize(0);
    }
  }
}

---

The documentation for this class was generated from the following files:

• GlobalMasterServer.h
• GlobalMasterServer.C