CepsMathOperations.hpp

Go to the documentation of this file.

/*~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
   This file is part of CEPS.
 
   CEPS is free software: you can redistribute it and/or modify
   it under the terms of the GNU General Public License as published by
   the Free Software Foundation, either version 3 of the License, or
   (at your option) any later version.
 
   CEPS is distributed in the hope that it will be useful,
   but WITHOUT ANY WARRANTY; without even the implied warranty of
   MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
   GNU General Public License for more details.
 
   You should have received a copy of the GNU General Public License
   along with CEPS (see file LICENSE at root of project).
   If not, see <https://www.gnu.org/licenses/>.
 
 
   Copyright 2019-2024 Inria, Universite de Bordeaux
 
   Authors, in alphabetical order:
 
     Pierre-Elliott BECUE, Florian CARO, Yves COUDIERE(*), Andjela DAVIDOVIC, 
     Charlie DOUANLA-LONTSI, Marc FUENTES, Mehdi JUHOOR, Michael LEGUEBE(*), 
     Pauline MIGERDITICHAN, Valentin PANNETIER(*), Nejib ZEMZEMI.
     * : currently active authors
 
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~*/
#pragma once
#include "common/core/CepsTypes.hpp"
 
template <class _Left, class _Right>
struct CepsMathOperation;
 
template <>
struct CepsMathOperation<CepsMathScalar, CepsMathScalar>
{
  static CepsMathScalar
  plus(const CepsMathScalar& l, const CepsMathScalar& r)
  {
    return l + r;
  }
 
  static CepsMathScalar
  minus(const CepsMathScalar& l, const CepsMathScalar& r)
  {
    return l - r;
  }
 
  static CepsMathScalar
  mult(const CepsMathScalar& l, const CepsMathScalar& r)
  {
    return l * r;
  }
 
  static CepsMathScalar
  div(const CepsMathScalar& l, const CepsMathScalar& r)
  {
    return l / r;
  }
};
 
template <>
struct CepsMathOperation<CepsMathScalar, CepsMathVertex>
{
  static CepsMathVertex
  plus(const CepsMathScalar& l, const CepsMathVertex& r)
  {
    return l * CepsMathVertex::Ones() + r;
  }
 
  static CepsMathVertex
  minus(const CepsMathScalar& l, const CepsMathVertex& r)
  {
    return l * CepsMathVertex::Ones() - r;
  }
 
  static CepsMathVertex
  mult(const CepsMathScalar& l, const CepsMathVertex& r)
  {
    return l * r;
  }
 
  static CepsMathVertex
  div(const CepsMathScalar& l, const CepsMathVertex& r)
  {
    CEPS_ABORT("operation scalar/vertex is undefined");
    return r;
  }
};
 
template <>
struct CepsMathOperation<CepsMathScalar, CepsMathTensor>
{
  static CepsMathTensor
  plus(const CepsMathScalar& l, const CepsMathTensor& r)
  {
    return l * CepsMathTensor::Ones() + r;
  }
 
  static CepsMathTensor
  minus(const CepsMathScalar& l, const CepsMathTensor& r)
  {
    return l * CepsMathTensor::Ones() - r;
  }
 
  static CepsMathTensor
  mult(const CepsMathScalar& l, const CepsMathTensor& r)
  {
    return l * r;
  }
 
  static CepsMathTensor
  div(const CepsMathScalar& l, const CepsMathTensor& r)
  {
    CEPS_ABORT("operation scalar/tensor is undefined");
    return r;
  }
};
 
template <>
struct CepsMathOperation<CepsMathVertex, CepsMathScalar>
{
  static CepsMathVertex
  plus(const CepsMathVertex& l, const CepsMathScalar& r)
  {
    return l + r * CepsMathVertex::Ones();
  }
 
  static CepsMathVertex
  minus(const CepsMathVertex& l, const CepsMathScalar& r)
  {
    return l - r * CepsMathVertex::Ones();
  }
 
  static CepsMathVertex
  mult(const CepsMathVertex& l, const CepsMathScalar& r)
  {
    return l * r;
  }
 
  static CepsMathVertex
  div(const CepsMathVertex& l, const CepsMathScalar& r)
  {
    return l / r;
  }
};
 
template <>
struct CepsMathOperation<CepsMathVertex, CepsMathVertex>
{
  static CepsMathVertex
  plus(const CepsMathVertex& l, const CepsMathVertex& r)
  {
    return l + r;
  }
 
  static CepsMathVertex
  minus(const CepsMathVertex& l, const CepsMathVertex& r)
  {
    return l - r;
  }
 
  static CepsMathScalar
  mult(const CepsMathVertex& l, const CepsMathVertex& r)
  {
    return l.dot(r);
  }
 
  static CepsMathVertex
  div(const CepsMathVertex& l, const CepsMathVertex& r)
  {
    CEPS_ABORT("operation vertex/vertex is undefined");
    return r;
  }
};
 
template <>
struct CepsMathOperation<CepsMathVertex, CepsMathTensor>
{
  static CepsMathTensor
  plus(const CepsMathVertex& l, const CepsMathTensor& r)
  {
    CEPS_ABORT("operation vertex+tensor is undefined");
    return r;
  }
 
  static CepsMathTensor
  minus(const CepsMathVertex& l, const CepsMathTensor& r)
  {
    CEPS_ABORT("operation vertex-tensor is undefined");
    return r;
  }
 
  static CepsMathVertex
  mult(const CepsMathVertex& l, const CepsMathTensor& r)
  {
    return (l.transpose() * r).transpose();
  }
 
  static CepsMathTensor
  div(const CepsMathVertex& l, const CepsMathTensor& r)
  {
    CEPS_ABORT("operation vertex/tensor is undefined");
    return r;
  }
};
 
template <>
struct CepsMathOperation<CepsMathTensor, CepsMathScalar>
{
  static CepsMathTensor
  plus(const CepsMathTensor& l, const CepsMathScalar& r)
  {
    return l + r * CepsMathTensor::Ones();
  }
 
  static CepsMathTensor
  minus(const CepsMathTensor& l, const CepsMathScalar& r)
  {
    return l - r * CepsMathTensor::Ones();
  }
 
  static CepsMathTensor
  mult(const CepsMathTensor& l, const CepsMathScalar& r)
  {
    return l * r;
  }
 
  static CepsMathTensor
  div(const CepsMathTensor& l, const CepsMathScalar& r)
  {
    return l / r;
  }
};
 
template <>
struct CepsMathOperation<CepsMathTensor, CepsMathVertex>
{
  static CepsMathTensor
  plus(const CepsMathTensor& l, const CepsMathVertex& r)
  {
    CEPS_ABORT("operation tensor+vertex is undefined");
    return l;
  }
 
  static CepsMathTensor
  minus(const CepsMathTensor& l, const CepsMathVertex& r)
  {
    CEPS_ABORT("operation tensor-vertex is undefined");
    return l;
  }
 
  static CepsMathVertex
  mult(const CepsMathTensor& l, const CepsMathVertex& r)
  {
    return l * r;
  }
 
  static CepsMathTensor
  div(const CepsMathTensor& l, const CepsMathVertex& r)
  {
    CEPS_ABORT("operation tensor/vertex is undefined");
    return l;
  }
};
 
template <>
struct CepsMathOperation<CepsMathTensor, CepsMathTensor>
{
  static CepsMathTensor
  plus(const CepsMathTensor& l, const CepsMathTensor& r)
  {
    return l + r;
  }
 
  static CepsMathTensor
  minus(const CepsMathTensor& l, const CepsMathTensor& r)
  {
    return l - r;
  }
 
  static CepsMathTensor
  mult(const CepsMathTensor& l, const CepsMathTensor& r)
  {
    return l * r;
  }
 
  static CepsMathTensor
  div(const CepsMathTensor& l, const CepsMathTensor& r)
  {
    CEPS_ABORT("operation tensor/tensor is undefined");
    return r;
  }
};