fei_MatrixTraits_LinSysCore.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_LinSysCore_hpp_
#define _fei_MatrixTraits_LinSysCore_hpp_
 
//This file defines matrix traits for LinearSystemCore matrices
//(well, "matrix-views" to be more precise).
//
 
#include <fei_LinearSystemCore.hpp>
 
namespace fei {

  template<>
  struct MatrixTraits<LinearSystemCore> {

    static const char* typeName()
      { return("LinearSystemCore"); }
 
    static double* getBeginPointer(LinearSystemCore* lsc)
      {
         return lsc->getMatrixBeginPointer();
      }
 
    static int getOffset(LinearSystemCore* lsc, int row, int col)
      {
         return lsc->getMatrixOffset(row,col);
      }

    static int setValues(LinearSystemCore* lsc, double scalar)
      {
  return( lsc->resetMatrix(scalar) );
      }

    static int getNumLocalRows(LinearSystemCore* lsc, int& numRows)
    {
      numRows = -1;
      return(-1);
    }

    static int getRowLength(LinearSystemCore* lsc, int row, int& length)
      {
  return( lsc->getMatrixRowLength(row, length) );
      }

    static int copyOutRow(LinearSystemCore* lsc,
          int row, int len, double* coefs, int* indices)
      {
        int dummy;
  return( lsc->getMatrixRow(row, coefs, indices, len, dummy) );
      }

    static int putValuesIn(LinearSystemCore* lsc,
                           int numRows, const int* rows,
                           int numCols, const int* cols,
                           const double* const* values,
                           bool sum_into)
      {
        if (sum_into) {
          return( lsc->sumIntoSystemMatrix(numRows, rows,
                                           numCols, cols, values) );
        }
        else {
    return( lsc->putIntoSystemMatrix(numRows, rows,
                                           numCols, cols, values) );
        }
      }

    static int globalAssemble(LinearSystemCore* lsc)
    {
      return( lsc->matrixLoadComplete() );
    }

    static int matvec(LinearSystemCore* lsc,
          fei::Vector* x,
          fei::Vector* y)
    {
      return( -1 );
    }
  };//struct MatrixTraits
}//namespace fei
 
#endif // _fei_MatrixTraits_LinSysCore_hpp_