const double aux_nodes_coeff = static_cast<double>(TNumNodes);

const GeometryType& rMasterGeometry,

const array_1d<double, 3>& rSlaveNormal,

const DecompositionType& rDecompGeom,

const NormalDerivativesComputation ConsiderNormalVariation

)

{

const array_1d<double, 3>& rSlaveNormal,

const array_1d<double, 3>& rMasterNormal,

const DecompositionType& rDecompGeom,

const NormalDerivativesComputation ConsiderNormalVariation,

const bool DualLM

)

rDecompGeom.ShapeFunctionsValues( N_decomp, rLocalPointDecomp.Coordinates() );

if constexpr (TDim == 2) {

for (IndexType i_mortar_node = 0; i_mortar_node < TNumNodes; ++i_mortar_node) {

// Projecting points in opposite geometry, defining mortar nodes

}

// Auxiliary point for local coordinates

// Delta local coordinates

array_1d<double, TNumNodes> DeltaXi_slave, DeltaXi_master;

rVariables.Initialize();

for (IndexType i_geom = 0; i_geom < rConditionsPointsSlave.size(); ++i_geom) {

array_1d<BelongType, TDim> belong_array;

for (IndexType i_node = 0; i_node < TDim; ++i_node) {

rSlaveGeometry.GlobalCoordinates(global_point, rConditionsPointsSlave[i_geom][i_node]);

belong_array[i_node] = rConditionsPointsSlave[i_geom][i_node].GetBelong();

}

// Integrating the mortar operators

array_1d<double,3> gp_normal;

GeometryType::CoordinatesArrayType slave_gp_global, projected_gp_local;

for ( IndexType point_number = 0; point_number < integration_points_slave.size(); ++point_number ) {

// We reset the derivatives

rDerivativeData.ResetDerivatives();

// We compute the local coordinates

// Calculate the kinematic variables

// We compute the current configuration

noalias(gp_normal) = MortarUtilities::GaussPointUnitNormal(rVariables.NSlave, rSlaveGeometry);

rSlaveGeometry.GlobalCoordinates( slave_gp_global, local_point_parent );

rMasterGeometry.PointLocalCoordinates( projected_gp_local, projected_gp_global.Coordinates( ) ) ;

using ConditionArrayListType = typename std::vector<ConditionArrayType>;

/// The line definition

/// The triangle definition

/// The geometry for decomposition (line in 2D and triangle for 3D)

using DecompositionType = typename std::conditional<TDim == 2, LineType, TriangleType >::type;

const GeometryType& rMasterGeometry,

const array_1d<double, 3>& rSlaveNormal,

const DecompositionType& rDecompGeom,

const NormalDerivativesComputation ConsiderNormalVariation = NormalDerivativesComputation::NO_DERIVATIVES_COMPUTATION

);

const array_1d<double, 3>& rSlaveNormal,

const array_1d<double, 3>& rMasterNormal,

const DecompositionType& rDecompGeom,

const NormalDerivativesComputation ConsiderNormalVariation = NormalDerivativesComputation::NO_DERIVATIVES_COMPUTATION,

const bool DualLM = false

);