PlayaSerialVector.hpp

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

#ifndef PLAYA_SERIAL_VECTOR_HPP
#define PLAYA_SERIAL_VECTOR_HPP

#include "PlayaDefs.hpp"
#include "PlayaSingleChunkVector.hpp"
#include "PlayaLoadableVector.hpp"
#include "PlayaVectorSpaceDecl.hpp"
#include "Teuchos_Array.hpp"

namespace Playa
{
using namespace Teuchos;

template <class Scalar> class Vector;

class SerialVector : public SingleChunkVector<double>,
                     public LoadableVector<double>,
                     public Describable
{
public:

  SerialVector(const VectorSpace<double>& vs);

  RCP<const VectorSpaceBase<double> > space() const {return vecSpace_.ptr();}

  void setElement(int globalIndex, const double& value);

  void addToElement(int globalIndex, const double& value);

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


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

  void finalizeAssembly();

  std::string description() const ;



  void getElements(const int* globalIndices, int numElems,
    Array<double>& elems) const ;

  static const SerialVector* getConcrete(const Vector<double>& x);
  static SerialVector* getConcrete(Vector<double>& x);

      
  virtual const double* dataPtr() const {return &(data_[0]);}
  virtual double* dataPtr() {return &(data_[0]);}

  virtual int chunkSize() const {return dim_;}



  
private:

  VectorSpace<double> vecSpace_;

  Array<double> data_;

  int dim_;
};
  
}


#endif