---

ElementGatherSeq.hh

#ifndef SEQELEMENTGATHER
#define SEQELEMENTGATHER

#include "ElementGather.hh"
#include "SimpleIntData.hh"

class IndexSetCollector;

struct Range {
  int start, extent;
  Range( int s, int e ) : start(s), extent(e) {}
};

class ElementGatherSeq : public ElementGather {
public:
  ElementGatherSeq( int nx, int ny, int capacity = 0 );
  ElementGatherSeq( int domain, IndexSetCollector & w );
  ElementGatherSeq( int nx, const Range & r );
  static ElementGatherSeq * Identity( int nx );

  virtual IotrMask * asMask();
  virtual ~ElementGatherSeq();
  bool operator==( const ElementGatherSeq & e ) const;
  bool operator!=( const ElementGatherSeq & e ) const
  {
    return !(*this == e );
  }

  virtual void copyFrom( const ElementGatherSeq& elt );
  virtual void copyInto( ElementGatherSeq& elt ) const;

  virtual ElementGather * copy() const;
  virtual void gather( double alpha, const IotrVector & x,
                       double beta,  IotrVector & y ) const;

  virtual void gatherScaled( double alpha, const IotrVector & diag,
                             const IotrVector & x,
                             double beta,
                             IotrVector & y ) const;
  virtual void gatherThenScale( double alpha, const IotrVector & x, 
                                double beta,  IotrVector & y,
                                const IotrVector & d ) const;
  virtual ElementGather * composeWith( const ElementGather & b ) const;

  virtual void unionWith2( const ElementGather & b, const ElementGather & c,
                           ElementGather *& u,
                           ElementGather *& AinU, ElementGather *&BinU,
                           ElementGather *& CinU ) const;


  virtual IotrVector * rangeVector() const;
  virtual IotrVector * rangeVector( const IotrVector & v ) const;
  virtual IotrVector * domainVector() const;
  
  //  specializes ElementGather::scatter()
  virtual void scatter( double alpha, const IotrVector & x, 
                        double beta,  IotrVector & y ) const;
  virtual void scatterInverse( double alpha, const IotrVector & x, 
                               double beta,  IotrVector & y ) const;
  virtual void scatterScaled(  double alpha, const IotrVector & diag,
                               const IotrVector & x,
                               double beta,  IotrVector & y ) const;
  virtual void scatterScaledInverse(  double alpha,
                                      const IotrVector & diag, 
                                      const IotrVector & x, double beta,  
                                      IotrVector & y ) const;
  virtual bool isIdentity() const;
  // specialize ElementGather::scale()
  //virtual void scale( double alpha, IotrVector & x);

  const int * indices() const { return mIndexArray->array(); }
  // Kill the non-const version soon!
  //protected:
  int * indices() { return mIndexArray->array(); }
public:
  int nindices() const { return mIndexArray->size(); }

  virtual ElementTransfer * transferTo( ElementGather * e );
  virtual void unionWith( const ElementGather & b,
                          ElementGather *& abunion,
                          ElementGather *& union2a, 
                          ElementGather *& union2b ) const;

  virtual void copyFromIndices( const int ind[], int & err );

  virtual void inspect( const char * name = 0,
                        const char * filename = 0 ) const;
  virtual void inspectAsMatlab( const char * name = 0,
                                const char * filename = 0 ) const;
  virtual int isValid() const;
protected:
  SimpleIntDataHandle mIndexArray;
};


typedef Handle<ElementGatherSeq> ElementGatherSeqHandle;

class ElementTransferSeq : public ElementTransfer {
protected:
  ElementGatherSeqHandle mFrom, mTo;
public:
  ElementTransferSeq( ElementGatherSeq * from, ElementGatherSeq * to );
  virtual void transferForward( double alpha, IotrVector & x, 
                                double beta,  IotrVector & y );
                                
  virtual void transferBackward( double alpha, IotrVector & x, 
                                 double beta,  IotrVector & y );
                                 
  virtual void inspect( const char * name = 0,
                        const char * filename = 0 );

  virtual void inspectAsMatlab( const char * name = 0,
                                const char * filename = 0 );
};

typedef Handle< ElementTransferSeq> ElementTransferSeqHandle;
#endif

---