contact_searching.h

vs.lib is a math library in C++ with a set of linear algebra and integrable / differentiable objects

#ifndef __CONTACT_SEARCHING_H
#define __CONTACT_SEARCHING_H

class Contact_Omega_eh : public Omega_eh {
	int the_associate_element[2];
   double the_projection_parameter[4];
   double the_normal_vector[4];
public:
	Contact_Omega_eh(int* ae, double* pp, double* nv, Omega_eh& oh);
   int* associate_element() { return the_associate_element; }
   double* projection_parameter() { return the_projection_parameter; }
   double* normal_vector() { return the_normal_vector; }
};

class Contact_Omega_h_i : public Omega_h_i {
public:
	Contact_Omega_h_i(int i); // client program suplied constructor
   void contact_searching_algorithm(Contact_Omega_h_i &interface_oh,
   	Global_Discretization_Gamma_h_i &a_gd, Global_Discretization_Gamma_h_i &b_gd,
   	U_h& lambda_h, gh_on_Gamma_h_i& interface_gh, double tolerance_gap = 1.e-3);
   void contact_searching_algorithm(Contact_Omega_h_i &interface_oh,
   	Global_Discretization_Gamma_h_i &a_gd, Global_Discretization_Gamma_h_i &b_gd,
   	U_h& lambda_h, gh_on_Gamma_h_i& interface_gh, Matrix_Representation& mr_K_0,
      Matrix_Representation& mr_K_1, Matrix_Representation_Couple& mrc_Q_0,
      Matrix_Representation_Couple& mrc_Q_1,
      double tolerance_gap = 1.e-3); // defining the "active set" step in the context of constraint optimization
};
#endif