File numerics/src/FrictionContact/gfc3d_Solvers.h

Go to the source code of this file

Subroutines for the resolution of contact problems with friction (3-dimensional case).

Typedefs

typedef void (*ComputeErrorGlobalPtr)(GlobalFrictionContactProblem *, double *, double *, double *, double, SolverOptions *, double, double *)
typedef void (*FreeSolverGlobalPtr)(GlobalFrictionContactProblem *)
typedef void (*PostSolverGlobalPtr)(int, double *)
typedef void (*SolverGlobalPtr)(int, int, double *, int *, double *)

Functions

void gfc3d_ACLMFixedPoint(GlobalFrictionContactProblem *problem, double *reaction, double *velocity, double *globalVelocity, int *info, SolverOptions *options)

Solver based on the fixed-point iteration proposed by Cadoux for friction-contact 3D problem.

Parameters
  • problem: the friction-contact 3D problem to solve
  • velocity: global vector (n), in-out parameter
  • reaction: global vector (n), in-out parameters
  • globalVelocity: global vector (m), in-out parameters
  • info: return 0 if the solution is found
  • options: the solver options : iparam[0] : Maximum iteration number iparam[4] ; local strategy dparam[0] : tolerance dparam[2] : localtolerance dparam[1] : (out) error

int gfc3d_ACLMFixedPoint_setDefaultSolverOptions(SolverOptions *options)
void gfc3d_ADMM(GlobalFrictionContactProblem *problem, double *reaction, double *velocity, double *globalVelocity, int *info, SolverOptions *options)
void gfc3d_ADMM_free(GlobalFrictionContactProblem *problem, SolverOptions *options)
void gfc3d_ADMM_init(GlobalFrictionContactProblem *problem, SolverOptions *options)
int gfc3d_ADMM_setDefaultSolverOptions(SolverOptions *options)
void gfc3d_admm_wr(GlobalFrictionContactProblem *problem, double *reaction, double *velocity, double *globalVelocity, int *info, SolverOptions *options)
int gfc3d_admm_wr_setDefaultSolverOptions(SolverOptions *options)
void gfc3d_AVI_gams_path(GlobalFrictionContactProblem *problem, double *reaction, double *velocity, int *info, SolverOptions *options)

solver using PATH (via GAMS) for friction-contact 3D problem based on an AVI reformulation

Parameters
  • problem: the friction-contact 3D problem to solve
  • velocity: global vector (n), in-out parameter
  • reaction: global vector (n), in-out parameters
  • info: return 0 if the solution is found
  • options: the solver options

void gfc3d_AVI_gams_pathvi(GlobalFrictionContactProblem *problem, double *reaction, double *velocity, int *info, SolverOptions *options)

solver using PATHVI (via GAMS) for friction-contact 3D problem based on an AVI reformulation

Parameters
  • problem: the friction-contact 3D problem to solve
  • velocity: global vector (n), in-out parameter
  • reaction: global vector (n), in-out parameters
  • info: return 0 if the solution is found
  • options: the solver options

int gfc3d_checkTrivialCaseGlobal(int dim, double *q, double *velocity, double *reaction, double *globalVelocity, SolverOptions *options)

Check for trivial solution in the friction-contact 3D problem.

Return
int =0 if a trivial solution has been found, else = -1
Parameters
  • dim: of the problem
  • q: global vector (n)
  • velocity: global vector (n), in-out parameter
  • reaction: global vector (n), in-out parameters
  • globalVelocity: the velocity in global coordinates
  • options: the pointer to the array of options to set

int gfc3d_DeSaxceFixedPoint_setDefaultSolverOptions(SolverOptions *options)
void gfc3d_DeSaxceFixedPoint_wr(GlobalFrictionContactProblem *problem, double *reaction, double *velocity, double *globalVelocity, int *info, SolverOptions *options)

Fixed Point iteration on De Saxe formulation solver with reformulation for friction-contact 3D problem.

Parameters
  • problem: the friction-contact 3D problem to solve
  • velocity: global vector (n), in-out parameter
  • reaction: global vector (n), in-out parameters
  • globalVelocity: global vector (m), in-out parameters
  • info: return 0 if the solution is found
  • options: the solver options : iparam[0] : Maximum iteration number dparam[0] : tolerance dparam[2] : localtolerance dparam[1] : (out) error

void gfc3d_nonsmooth_Newton_AlartCurnier(GlobalFrictionContactProblem *problem, double *reaction, double *velocity, double *globalVelocity, int *info, SolverOptions *options)
void gfc3d_nonsmooth_Newton_AlartCurnier_wr(GlobalFrictionContactProblem *problem, double *reaction, double *velocity, double *globalVelocity, int *info, SolverOptions *options)
int gfc3d_nonsmooth_Newton_AlartCurnier_wr_setDefaultSolverOptions(SolverOptions *options)
void gfc3d_nsgs(GlobalFrictionContactProblem *problem, double *reaction, double *velocity, double *globalVelocity, int *info, SolverOptions *options)

Non-Smooth Gauss Seidel solver for friction-contact 3D problem.

Parameters
  • problem: the friction-contact 3D problem to solve
  • velocity: global vector (n), in-out parameter
  • reaction: global vector (n), in-out parameters
  • globalVelocity: global vector (m), in-out parameters
  • info: return 0 if the solution is found
  • options: the solver options : iparam[0] : Maximum iteration number iparam[4] ; local strategy dparam[0] : tolerance dparam[2] : localtolerance dparam[1] : (out) error

void gfc3d_nsgs_velocity_wr(GlobalFrictionContactProblem *problem, double *reaction, double *velocity, double *globalVelocity, int *info, SolverOptions *options)

Non-Smooth Gauss Seidel solver for friction-contact 3D problem with iteration on velocities.

Parameters
  • problem: the friction-contact 3D problem to solve
  • velocity: global vector (n), in-out parameter
  • reaction: global vector (n), in-out parameters
  • globalVelocity: global vector (m), in-out parameters
  • info: return 0 if the solution is found
  • options: the solver options : iparam[0] : Maximum iteration number dparam[0] : tolerance dparam[2] : localtolerance dparam[1] : (out) error

int gfc3d_nsgs_velocity_wr_setDefaultSolverOptions(SolverOptions *options)
void gfc3d_nsgs_wr(GlobalFrictionContactProblem *problem, double *reaction, double *velocity, double *globalVelocity, int *info, SolverOptions *options)

Non-Smooth Gauss Seidel solver with reformulation for friction-contact 3D problem.

Parameters
  • problem: the friction-contact 3D problem to solve
  • velocity: global vector (n), in-out parameter
  • reaction: global vector (n), in-out parameters
  • globalVelocity: global vector (m), in-out parameters
  • info: return 0 if the solution is found
  • options: the solver options : iparam[0] : Maximum iteration number iparam[4] : localsolver choice 0: projection on Cone, 1: Newton/AlartCurnier, 2: projection on Cone with local iteration, 2: projection on Disk with diagonalization, dparam[0] : tolerance dparam[2] : localtolerance dparam[1] : (out) error

int gfc3d_nsgs_wr_setDefaultSolverOptions(SolverOptions *options)
void gfc3d_proximal_wr(GlobalFrictionContactProblem *problem, double *reaction, double *velocity, double *globalVelocity, int *info, SolverOptions *options)

Proximal point solver with reformulation for friction-contact 3D problem.

Parameters
  • problem: the friction-contact 3D problem to solve
  • velocity: global vector (n), in-out parameter
  • reaction: global vector (n), in-out parameters
  • globalVelocity: global vector (m), in-out parameters
  • info: return 0 if the solution is found
  • options: the solver options : iparam[0] : Maximum iteration number iparam[4] : localsolver choice 0: projection on Cone, 1: Newton/AlartCurnier, 2: projection on Cone with local iteration, 2: projection on Disk with diagonalization, dparam[0] : tolerance dparam[2] : localtolerance dparam[1] : (out) error

int gfc3d_proximal_wr_setDefaultSolverOptions(SolverOptions *options)
void gfc3d_set_internalsolver_tolerance(GlobalFrictionContactProblem *problem, SolverOptions *options, SolverOptions *internalsolver_options, double error)
int gfc3d_setDefaultSolverOptions(SolverOptions *options, int solverId)

set the default solver parameters and perform memory allocation for gfc3d

Parameters
  • options: the pointer to the array of options to set
  • solverId: int identifier of the solver

int gfc3d_TrescaFixedPoint_setDefaultSolverOptions(SolverOptions *options)
void gfc3d_TrescaFixedPoint_wr(GlobalFrictionContactProblem *problem, double *reaction, double *velocity, double *globalVelocity, int *info, SolverOptions *options)

Fied Point iteration on Tresca Friction Cylinder with reformulation for friction-contact 3D problem.

Parameters
  • problem: the friction-contact 3D problem to solve
  • velocity: global vector (n), in-out parameter
  • reaction: global vector (n), in-out parameters
  • globalVelocity: global vector (m), in-out parameters
  • info: return 0 if the solution is found
  • options: the solver options : iparam[0] : Maximum iteration number dparam[0] : tolerance dparam[2] : localtolerance dparam[1] : (out) error

void gfc3d_VI_ExtraGradient(GlobalFrictionContactProblem *problem, double *reaction, double *velocity, double *globalVelocity, int *info, SolverOptions *options)
int gfc3d_VI_ExtraGradient_setDefaultSolverOptions(SolverOptions *options)
void gfc3d_VI_FixedPointProjection(GlobalFrictionContactProblem *problem, double *reaction, double *velocity, double *globalVelocity, int *info, SolverOptions *options)
int gfc3d_VI_FixedPointProjection_setDefaultSolverOptions(SolverOptions *options)