 
 	custom_vector<real> lm(number_of_contacts);
 
-#pragma omp parallel for
+//#pragma omp parallel for
 		for (uint index = 0; index < number_of_contacts; index++) {
 			int2 id_ = ids[index];
 			real3 _mx;


		shurA(mx);

#pragma omp parallel for
		for (uint index = 0; index < number_of_contacts; index++) {
			real3 temp = R3(0);
			int2 id_ = ids[index];

	real theta_k=1.0;
	real theta_k1=theta_k;
	real beta_k1=0.0;
//custom_vector<real> ml = x;
		Project(x);
//custom_vector<real> ml_candidate = ml;

		custom_vector<real> mg(x.size());
		ShurProduct(x,mg);
