fei_MatrixTraits_LinProbMgr.hpp

/*--------------------------------------------------------------------*/
/*    Copyright 2005 Sandia Corporation.                              */
/*    Under the terms of Contract DE-AC04-94AL85000, there is a       */
/*    non-exclusive license for use of this work by or on behalf      */
/*    of the U.S. Government.  Export of this program may require     */
/*    a license from the United States Government.                    */
/*--------------------------------------------------------------------*/

#ifndef _fei_MatrixTraits_LinProbMgr_hpp_
#define _fei_MatrixTraits_LinProbMgr_hpp_

//This file defines matrix traits for LinearProblemManager matrix
//representations.
//

#include <fei_LinearProblemManager.hpp>

namespace fei {

  template<>
  struct MatrixTraits<fei::LinearProblemManager> {

    static const char* typeName()
      { return("fei::LinearProblemManager"); }

    static double* getBeginPointer(fei::LinearProblemManager* /*mat*/)
      {
        return NULL;
      }

    static int getOffset(fei::LinearProblemManager* /*mat*/, int /*row*/, int /*col*/)
      {
        return -1;
      }

    static int setValues(fei::LinearProblemManager* mat, double scalar)
      {
  mat->setMatrixValues(scalar);
        return(0);
      }

    static int getNumLocalRows(fei::LinearProblemManager* mat, int& numRows)
    {
      numRows = mat->getLocalNumRows();
      return(0);
    }

    static int getRowLength(fei::LinearProblemManager* mat, int row, int& length)
      {
  length = mat->getRowLength(row);
        if (length < 0) return(length);
        return(0);
      }

    static int copyOutRow(fei::LinearProblemManager* mat,
          int row, int len, double* coefs, int* indices)
      {
  return( mat->copyOutMatrixRow(row, len, coefs, indices) );
      }

    static int putValuesIn(fei::LinearProblemManager* mat,
                           int numRows, const int* rows,
                           int numCols, const int* cols,
                           const double* const* values,
                           bool sum_into)
      {
        return( mat->insertMatrixValues(numRows, rows,
                                        numCols, cols,
                                        values, sum_into) );
      }

    static int globalAssemble(fei::LinearProblemManager* mat)
    {
      return( mat->globalAssemble() );
    }

    static int matvec(fei::LinearProblemManager* mat,
          fei::Vector* x,
          fei::Vector* y)
    {
      return( -1 );
    }
  };//struct MatrixTraits
}//namespace fei

#endif // _fei_MatrixTraits_LinProbMgr_hpp_