EigenTools.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 <Eigen/Dense>
 
#include "common/core/CepsTypes.hpp"
 
namespace ceps
{
  // copy from FEAssembler
  // using DofMatrix = CepsMathDynamic2D;
  // using DofVector = CepsMathDynamic1D;
  // using Vertex = Eigen::Matrix<CepsReal, 3, 1>;
  // using DofScalars = CepsMathDynamic1D;
 
  namespace math
  {
 
    void
    computeInverse (const CepsMathDynamic2D &m, CepsMathDynamic2D &inv);
 
    void
    computeInverse (const Eigen::Matrix<CepsReal, 3, 2> &m, Eigen::Matrix<CepsReal, 2, 3> &inv);
 
    void
    computeInverse (const Eigen::Matrix<CepsReal, 3, 1> &m, Eigen::Matrix<CepsReal, 1, 3> &inv);
 
    void
    computeInverse (const Eigen::Matrix<CepsReal, 2, 1> &m, Eigen::Matrix<CepsReal, 1, 2> &inv);
 
    CepsMathDynamic2D
    compose (CepsUInt M, CepsVector<CepsMathDynamic2D *> blocks);
 
    CepsMathDynamic1D
    compose (CepsUInt M, CepsVector<CepsMathDynamic1D *> blocks);
 
    void
    putRowsSumInDiagonal (CepsMathDynamic2D &mat);
 
    CepsMathDynamic2D
    multiplyByBlocks(
      const CepsMathDynamic2D&   A,
      const CepsMathDynamic2D&   B,
      const CepsArray2<CepsUInt> a,
      const CepsArray2<CepsUInt> b
    );
 
// see topLeftCorner implemented method in Eigen library (old version ?)
    // /// @brief Block extraction (C/C from Eigen wiki)
    // template <typename Derived>
    // Eigen::Block<Derived>
    // topLeftCorner (Eigen::MatrixBase<Derived> &m, int rows, int cols)
    // {
    //   return Eigen::Block<Derived> (m.derived (), 0, 0, rows, cols);
    // }
 
    // /// @brief Block extraction (C/C from Eigen wiki)
    // template <typename Derived>
    // const Eigen::Block<const Derived> &
    // topLeftCorner (const Eigen::MatrixBase<Derived> &m, int rows, int cols)
    // {
    //   return Eigen::Block<const Derived> (m.derived (), 0, 0, rows, cols);
    // }
 
  }  // namespace math
}  // namespace ceps