// tensor.H
#ifndef tensor_H
#define tensor_H
#include "vtr.h"
#include "../dgtd-performance.hpp"

class tensor {
private: 
  fp_t val[3][3];

public:
  tensor(fp_t = 0.0, fp_t = 0.0, fp_t = 0.0,
	 fp_t = 0.0, fp_t = 0.0, fp_t = 0.0,
	 fp_t = 0.0, fp_t = 0.0, fp_t = 0.0);
  void setEntry(int, int, fp_t );
  void settensor(vtr, vtr, vtr);

  // operators
  tensor operator+(const tensor &) const;
  tensor operator-(const tensor &) const;
  tensor operator*(const fp_t &) const;
  tensor operator*(const tensor &) const;
  tensor operator/(const fp_t &) const;
  tensor operator/(const tensor &) const;
  vtr operator*(vtr &);
  tensor& operator=(const tensor &);
  tensor inverse() const;
  tensor transpose() const;

  // get functions
  fp_t getEntry(int, int);
  
  // print function
  void print() const;
};

#endif