src/IplPaws/DLRepresentation.cpp

Go to the documentation of this file.

// -*- C++ -*-
/***************************************************************************
 *
 * The IPPL Framework
 * 
 * This program was prepared by PSI. 
 * All rights in the program are reserved by PSI.
 * Neither PSI nor the author(s)
 * makes any warranty, express or implied, or assumes any liability or
 * responsibility for the use of this software
 *
 * Visit http://www.acl.lanl.gov/POOMS for more details
 *
 ***************************************************************************/

// -*- C++ -*-
/***************************************************************************
 *
 * The IPPL Framework
 * 
 *
 * Visit http://www.acl.lanl.gov/POOMS for more details
 *
 ***************************************************************************/

#include <iostream.h>
#include "IpplPaws/DLRepresentation.h"


//
// Default DLRepresentation constructor.
//

template<class D, unsigned Dim>
DLRepresentation<D,Dim>::DLRepresentation()
{
  static char* id = "DLRepresentation::DLRepresentation()";
  PDebug(PAWSD_APP, ("%s entering\n", id));
  PDebug(PAWSD_APP, ("%s exiting\n", id));
}


//
// Construct a DLRepresentation given a FieldLayout
// Divide into Global and LocalBlocks
//

template<class D, unsigned Dim>
DLRepresentation<D,Dim>::DLRepresentation(FieldLayout<Dim> &layout)
{
  static char* id = "DLRepresentation::DLRepresentation(FieldLayout)";
  PDebug(PAWSD_APP, ("%s entering\n", id));

  update(layout);

  PDebug(PAWSD_APP, ("%s exiting\n", id));
}


//
// Update the DLRepresentation with a new FieldLayout and num_nodes
// This routine is called when a PawsRemoteData is created with an
// empty, dynamic array.  We want it to match the local domain, but
// to have the global blocks calculated using the number of processors
// on the remote machine
//
template<class D, unsigned Dim>
void
DLRepresentation<D,Dim>::update(FieldLayout<Dim> &layout)
{
  static char* id = "DLRepresentation::update(FieldLayout)";
  PDebug(PAWSD_APP, ("%s entering\n", id));

  // Clear out the current representation and add the new
  clearLocal();
  clearGlobal();

  setStorageOrder(PAWS_COLUMN);
  int total_nodes = Ippl::getNodes();
  int node_rank = Ippl::myNode();
  int i;

  // Number the blocks when building the representation
  // Used to assign the block number of the LField within the structure
  //
  int *block_num = new int[total_nodes];
  for (i = 0; i < total_nodes; i++)
    block_num[i] = 0;

  // Get the whole domain
  //
  int first[Dim];
  int last[Dim];
  int stride[Dim];
  const D &whole_domain = layout.getDomain();
  for (i = 0; i < Dim; i++) {
    first[i] = whole_domain[i].first();
    last[i] = whole_domain[i].last();
    stride[i] = whole_domain[i].stride();
  }
  initialize(Dim, first, last, stride, node_rank);

  // Iterate through the FieldLayout to copy out the local part
  // Iterator returns vmap with second element being a vnode having a domain
  //
  int bnum = 0;
  FieldLayout<Dim>::iterator_iv loc;
  for (loc = layout.begin_iv(); loc != layout.end_iv(); ++loc) {

    // Store the local block within the DLRepresentation
    const D &local_domain = (*loc).second->getDomain();
    for (i = 0; i < Dim; i++) {
      first[i] = local_domain[i].first();
      last[i] = local_domain[i].last();
      stride[i] = local_domain[i].stride();
    }
    insertLocal(Dim, first, last, stride, bnum++, node_rank);
  }

  // Iterate through the FieldLayout to copy out the remote part
  // Iterator returns actual vnodes
  //
  FieldLayout<Dim>::iterator_dv rem;
  for (rem = layout.begin_rdv(); rem != layout.end_rdv(); ++rem) {

    // Store the global block within the DLRepresentation
    int node = (*rem).second->getNode();
    const D &global_domain = (*rem).first;
    for (i = 0; i < Dim; i++) {
      first[i] = global_domain[i].first();
      last[i] = global_domain[i].last();
      stride[i] = global_domain[i].stride();
    }
    insertGlobal(Dim, first, last, stride, block_num[node]++, node);
  }

  delete [] block_num;
  PDebug(PAWSD_APP, ("%s exiting\n", id));
}


/***************************************************************************
 * $RCSfile: DLRepresentation.cpp,v $   $Author: adelmann $
 * $Revision: 1.1.1.1 $   $Date: 2003/01/23 07:40:32 $
 * IPPL_VERSION_ID: $Id: DLRepresentation.cpp,v 1.1.1.1 2003/01/23 07:40:32 adelmann Exp $ 
 ***************************************************************************/

---