quad_supp_func< K > Class Template Reference

#include <qsc_approximation_fcts.hpp>

List of all members.

Public Member Functions

• quad_supp_func ()
• quad_supp_func (support_function< K > h, normal_vector< K > n)
• Seq< K > Cq ()

---

Detailed Description

template<class K>
class mmx::shape::quad_supp_func< K >

Definition at line 1046 of file qsc_approximation_fcts.hpp.

---

Constructor & Destructor Documentation

quad_supp_func

(

)

[inline]

Definition at line 1049 of file qsc_approximation_fcts.hpp.

{} ;

quad_supp_func	(	support_function< K >	h,
		normal_vector< K >	n
	)			[inline]

Definition at line 1072 of file qsc_approximation_fcts.hpp.

References support_function< K >::eval(), support_function< K >::eval_diff(), normal_vector< K >::eval_unit(), mmx::linear::LUsolve(), normal_vector< K >::middle(), mmx::shape::N, mmx::pow(), and VECT.

{
  // mmout << "b\n";
  VECT N[3];
  N[0]=n.eval_unit(0);
  N[1]=n.eval_unit(1);
  N[2]=n.middle();
  // mmout << "c\n";
  K Vh[5]; // the values of h and its derivative at the control points
  K Mq[25]; // the values of q and its derivative at the control points (expressed in the base of the unknown we want to find
  for (int i=0; i<3; i++) 
    {
      VECT nn=N[i]; // 3 first lines : values at n0, n1, n1/2 
      Vh[i]=h.eval(nn);
      Mq[5*i]=std::pow(nn[0],2); // x^2
      Mq[5*i+1]=std::pow(nn[1],2); // y^2
      Mq[5*i+2]=nn[0]*nn[1]; // x*y CHANGE *2 ??
      Mq[5*i+3]=nn[0]; // x
      Mq[5*i+4]=nn[1]; // y
    }
  // mmout << "d\n";
  for (int i=3; i<5; i++) 
    {
      VECT nn=N[i-3]; // 2 last lines : derivative at n0, n1
      Vh[i]=h.eval_diff(nn);
      Mq[5*i+0]=-2*nn[0]*nn[1]; // -2*x*y
      Mq[5*i+1]=2*nn[0]*nn[1]; // 2*x*y
      Mq[5*i+2]=(std::pow(nn[0],2)-std::pow(nn[1],2)); // x^2-y^2 CHANGE *2 ?
      Mq[5*i+3]=-nn[1]; // -y
      Mq[5*i+4]=nn[0]; // x
    }

  linear::LUsolve(q,Mq,Vh,5); // Mq*q=Vh ; size 5
  // output : q[0]*x^2+q[1]*y^2+q[2]*x*y+q[3]*x+q[4]*y

  std::cout << "Valeur de Vh et Mq" << "\n";
  for (int i=0; i<5; i++){
    std::cout << Vh[i] << " = ";
    for (int j=0; j<5; j++)
      std::cout << Mq[5*i+j] << " ";
    std::cout << "\n";}
    std::cout << "\n";

}

---

Member Function Documentation

Seq< K > Cq

(

)

[inline]

Definition at line 1059 of file qsc_approximation_fcts.hpp.

{
  // mmout << "z\n";
  Seq<K> Q;
  // mmout << Q << "\n";
  for (int i=0; i<5; i++) 
    {
      Q << q[i];
    }
  return Q;
}

---

The documentation for this class was generated from the following file:

• include/shape/qsc_approximation_fcts.hpp