tensor.h

Go to the documentation of this file.

 
#ifndef TENSOR_H
#define TENSOR_H
 
#include <iostream>
#include <cmath>
#include <fstream>
#include "linalg.h"
using namespace std;
using namespace softsusy;
 
class Tensor;
ostream & operator <<(ostream &left, const Tensor &);
 
class Tensor {
private:
  DoubleMatrix A1, A2, A3; 
  
  void setTensor(const DoubleMatrix &, const DoubleMatrix &, const
                 DoubleMatrix &); 
  
public:
  Tensor(); 
  Tensor(const Tensor &);
  
  void threecheck(const DoubleMatrix &, const DoubleMatrix &, 
                  const DoubleMatrix &); 
  
  double & operator () (int, int, int);
  DoubleMatrix & operator () (int);
  const DoubleMatrix & display(int) const;
  double display(int, int, int) const;
  void checkOut(double) const;
  void set(int i, int j, int k, double f);
 
  DoubleVector trace(int) const;
  Tensor transpose() const;
  DoubleMatrix dotProd(const DoubleVector & v, int i) const;
  Tensor operator*(double) const;
  Tensor operator/(double) const;
  Tensor operator*(const DoubleMatrix &) const;
  Tensor operator+(const Tensor &) const;
  Tensor operator-(const Tensor &) const;
  Tensor product(const DoubleMatrix &) const;
  Tensor swap(int);
  Tensor raise(const DoubleMatrix & M) const;
};
Tensor operator*(double, const Tensor &);
Tensor operator*(const DoubleMatrix &, const Tensor &);
Tensor outerProduct(const DoubleVector &, const DoubleMatrix &, int);
DoubleMatrix matrixify(const Tensor &, const Tensor &);
// \f$ M_{ij} = A^{lik} B^{lkj} \f$
DoubleMatrix sumProd(const Tensor & A, const Tensor & B);
 
#endif