/[escript]/trunk/escript/src/AbstractSystemMatrix.cpp
ViewVC logotype

Diff of /trunk/escript/src/AbstractSystemMatrix.cpp

Parent Directory Parent Directory | Revision Log Revision Log | View Patch Patch

revision 97 by jgs, Tue Dec 14 05:39:33 2004 UTC revision 123 by jgs, Fri Jul 8 04:08:13 2005 UTC
# Line 58  Data AbstractSystemMatrix::vectorMultipl Line 58  Data AbstractSystemMatrix::vectorMultipl
58  {  {
59       if (isEmpty())       if (isEmpty())
60            throw SystemMatrixException("Error - Matrix is empty.");            throw SystemMatrixException("Error - Matrix is empty.");
      if (right.getFunctionSpace()!=getColumnFunctionSpace())  
           throw SystemMatrixException("Error - column function space and function space of input data do not match.");  
61       if (right.getDataPointSize()!=getColumnBlockSize())       if (right.getDataPointSize()!=getColumnBlockSize())
62            throw SystemMatrixException("Error - column block size and input data size do not match.");            throw SystemMatrixException("Error - column block size and input data size do not match.");
63       DataArrayView::ShapeType shape;       DataArrayView::ShapeType shape;
64       if (getRowBlockSize()>1) shape.push_back(getRowBlockSize());       if (getRowBlockSize()>1) shape.push_back(getRowBlockSize());
65    
66       Data out=Data(0.,shape,getRowFunctionSpace(),true);       Data out=Data(0.,shape,getRowFunctionSpace(),true);
67       ypAx(out,right);       ypAx(out,Data(right,getColumnFunctionSpace()));
68       return out;       return out;
69  }  }
70    
# Line 97  void AbstractSystemMatrix::saveMM(const Line 95  void AbstractSystemMatrix::saveMM(const
95  {  {
96      throw SystemMatrixException("Error - Matrix Market interface not available.");      throw SystemMatrixException("Error - Matrix Market interface not available.");
97  }  }
98  void AbstractSystemMatrix:: setValue(const double value) const  void AbstractSystemMatrix::saveHB(const std::string& fileName) const
99  {  {
100        throw SystemMatrixException("Error - Harwell-Boeing interface not available.");
101    }
102    void AbstractSystemMatrix::setValue(const double value) const
103    {
104        throw SystemMatrixException("Error - setValue is not implemented.");
105    }
106    void AbstractSystemMatrix::resetSolver() const
107    {
108        throw SystemMatrixException("Error - resetSolver not implemented.");
109  }  }
110    
111  }  // end of namespace  }  // end of namespace

Legend:
Removed from v.97  
changed lines
  Added in v.123

  ViewVC Help
Powered by ViewVC 1.1.26