matrix.hpp

Go to the documentation of this file.

// Copyright (c) 2010, Lawrence Livermore National Security, LLC. Produced at
// the Lawrence Livermore National Laboratory. LLNL-CODE-443211. All Rights
// reserved. See file COPYRIGHT for details.
//
// This file is part of the MFEM library. For more information and source code
// availability see http://mfem.googlecode.com.
//
// MFEM is free software; you can redistribute it and/or modify it under the
// terms of the GNU Lesser General Public License (as published by the Free
// Software Foundation) version 2.1 dated February 1999.

#ifndef MFEM_MATRIX
#define MFEM_MATRIX

// Abstract data types matrix, inverse matrix

#include "../general/array.hpp"
#include "operator.hpp"

class  MatrixInverse;

class Matrix : public Operator
{
   friend class MatrixInverse;
public:
   explicit Matrix (int s) { size=s; }

   virtual double& Elem (int i, int j) = 0;

   virtual const double& Elem (int i, int j) const = 0;

   virtual MatrixInverse * Inverse() const = 0;

   virtual void Finalize(int) { }

   virtual void Print (ostream & out = cout, int width = 4) const;

   virtual ~Matrix() { }
};


class MatrixInverse : public Operator
{
protected:
   const Matrix *a;

public:
   MatrixInverse (const Matrix &mat) { size = mat.size; a = &mat; }

   virtual ~MatrixInverse() { }
};

#endif