PlayaLoadableVector.hpp

/* @HEADER@ */
//   
 /* @HEADER@ */

#ifndef PLAYA_LOADABLEVECTOR_HPP
#define PLAYA_LOADABLEVECTOR_HPP

#include "PlayaDefs.hpp"

namespace Playa
{
template <class Scalar>
class LoadableVector 
{
public:
  virtual ~LoadableVector() {;}

  virtual void setElement(int globalIndex, const Scalar& value) = 0 ;

  virtual void addToElement(int globalIndex, const Scalar& value) = 0 ;

  virtual void setElements(int numElems, 
    const int* globalIndices, 
    const Scalar* values) ;

  virtual void addToElements(int numElems, 
    const int* globalIndices, 
    const Scalar* values);

  virtual void finalizeAssembly() {;}
};

/* Default implementation of setElements makes multiple calls to
 * setElement(). If at all possible, this should be overridden
 * with a method specialized to the underlying type.  */
template <class Scalar> 
inline void LoadableVector<Scalar>::setElements(int numElems, 
  const int* globalIndices, 
  const Scalar* values)
{
  for (int i=0; i<numElems; i++)
  {
    setElement(globalIndices[i], values[i]);
  }
}

/* Default implementation of addToElements makes multiple calls to
 * addToElement(). If at all possible, this should be overridden
 * with a method specialized to the underlying type.  */
template <class Scalar> 
inline void LoadableVector<Scalar>::addToElements(int numElems, 
  const int* globalIndices, 
  const Scalar* values)
{
  for (int i=0; i<numElems; i++)
  {
    addToElement(globalIndices[i], values[i]);
  }
}

  
  
}

#endif