Homec4science
Diffusion tamaas (petsc) Browse

tamaas/src/bem/bem_grid_polonski.hh
37620cc9f8b9petsc

Last Change
Raw File

bem_grid_polonski.hh
View Options

/**
*
* @author Lucas Frérot <lucas.frerot@epfl.ch>
*
* @section LICENSE
*
* Copyright (©) 2016 EPFL (Ecole Polytechnique Fédérale de
* Lausanne) Laboratory (LSMS - Laboratoire de Simulation en Mécanique des
* Solides)
*
* Tamaas is free software: you can redistribute it and/or modify it under the
* terms of the GNU Lesser General Public License as published by the Free
* Software Foundation, either version 3 of the License, or (at your option) any
* later version.
*
* Tamaas 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 Lesser General Public License for more
* details.
*
* You should have received a copy of the GNU Lesser General Public License
* along with Tamaas. If not, see <http://www.gnu.org/licenses/>.
*
*/
/* -------------------------------------------------------------------------- */
#ifndef __BEM_GRID_POLONSKI_HH__
#define __BEM_GRID_POLONSKI_HH__
/* -------------------------------------------------------------------------- */
#include "bem_grid_kato.hh"
/* -------------------------------------------------------------------------- */
namespace tamaas {
// typedef for swig
using Grid1dReal = Grid<Real, 1>;
/**
* @brief multi component implementation of Polonski algorithm
* This method's main unknown is the pressure.
*/
class BemGridPolonski : public BemGridKato {
public:
/// Constructor
BemGridPolonski(Surface<Real>& surface);
/// Destructor
virtual ~BemGridPolonski();
public:
/* --------------------------------------------------------------------------
*/
/* Actual computations */
/* --------------------------------------------------------------------------
*/
/// Compute equilibrium
virtual void computeEquilibrium(Real tol, const Grid<Real, 1>& p_target);
/// Compute gap functional
Real computeGradientNorm();
/// Compute optimal step size
Real computeOptimalStep();
/// Center field on contact area
void centerOnContactArea(Grid<Real, 2>& field);
/// Update search direction
void updateSearchDirection(Real delta, Real G, Real G_old);
/// Update tractions
Real updateTractions(Real tau);
/// Set pressure on all points
void setPressure(const Grid<Real, 1>& p_target);
protected:
Grid<Real, 2> search_direction;
Grid<Real, 2> projected_search_direction;
};
} // namespace tamaas
#endif // __BEM_GRID_POLONSKI_HH__
c4science · Help