Matrix.h

Go to the documentation of this file.

/* *******************************************************************
 * Rocstar Simulation Suite                                          *
 * Copyright@2015, Illinois Rocstar LLC. All rights reserved.        *
 *                                                                   *
 * Illinois Rocstar LLC                                              *
 * Champaign, IL                                                     *
 * www.illinoisrocstar.com                                           *
 * sales@illinoisrocstar.com                                         *
 *                                                                   *
 * License: See LICENSE file in top level of distribution package or *
 * http://opensource.org/licenses/NCSA                               *
 *********************************************************************/
/* *******************************************************************
 * THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND,   *
 * EXPRESS OR IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES   *
 * OF MERCHANTABILITY, FITNESS FOR A PARTICULAR PURPOSE AND          *
 * NONINFRINGEMENT.  IN NO EVENT SHALL THE CONTRIBUTORS OR           *
 * COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER       *
 * LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE,   *
 * ARISING FROM, OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE    *
 * USE OR OTHER DEALINGS WITH THE SOFTWARE.                          *
 *********************************************************************/
// $Id: Matrix.h,v 1.4 2008/12/06 08:45:23 mtcampbe Exp $

#ifndef __MATRIX_H__
#define __MATRIX_H__

#include "mopbasic.h"

MOP_BEGIN_NAMESPACE



class Matrix{
  friend std::ostream &operator <<(std::ostream &, const Matrix &);
public:
  
  Matrix() {}

  ~Matrix() {}

  double&  operator() (int, int);

  double   operator() (int, int) const;

  const Matrix &operator = (const Matrix &);


  const Matrix &operator ^= (const Matrix &);

  void info ();
protected:
  int     rows_;
  int     cols_;
  double  data_[9];
};


class J_Matrix : public Matrix {

public:
  
  // Constructor
  J_Matrix() {}

  // Destructor
  ~J_Matrix() {}

  // Constructor from noal coords and number of columns
  //  
  //  \param pnts[] A Vector_3 array of ordered nodal points.
  //  \param cols   Columns in the Jacobian.  Should be 2 for triangles or quadrilaterals
  //                and 3 for tetrahedrons or hexahedrons.
  J_Matrix(Vector_3<double> pnts[],int cols);


  double det();
};

MOP_END_NAMESPACE

#endif