Tensor.hpp

Go to the documentation of this file.

#ifndef TENSOR
#define TENSOR
 
#include <mecacell/mecacell.h>
 
struct Tensor {
private:
    double matrix[3][3] = {{0.,0.,0.},{0.,0.,0.},{0.,0.,0.}}; 
public:
    Tensor() {}
 
    Tensor(const MecaCell::Vec &f, const MecaCell::Vec &r){
        for(int i = 0; i < 3; ++i){
            for(int j = 0; j < 3; ++j){
                matrix[i][j] = f.coords[i] * r.coords[j];
            }
        }
    }
 
 
    inline double get(int i, int j) const { return matrix[i][j]; }
 
 
 
 
    Tensor& operator/=(double div) {
        for(int i = 0; i < 3; ++i){
            for(int j = 0; j <3; ++j){
                matrix[i][j] /= div;
            }
        }
        return(*this);
    }
 
 
    Tensor& operator+=(Tensor t){
        for(int i = 0; i < 3; ++i){
            for(int j = 0; j <3; ++j){
                matrix[i][j] += t.get(i,j);
            }
        }
        return(*this);
    }
 
 
    string to_string(){
        ostringstream str;
        for(int i = 0; i < 3; ++i){
            for(int j = 0; j <3; ++j){
               str << matrix[i][j] << ", ";
            }
            str << "\n";
        }
        return str.str();
    }
 
};
 
#endif // TENSOR