example_ellipse_getellipse.cpp

Programa para el testeo de las funciones.

Programa para el testeo de las funciones.Para compilar o código example_ellipse_getellipse.cpp:

g++ -static -o example_ellipse_getellipse example_ellipse_getellipse.cpp -lpdsramm 

Para executar o programa:

./example_ellipse_getellipse

Retornando por consola:

  a: 3
  b: 1.5
ang: -30 grad
 Xc: 5  3

  a: 3.06393
  b: 1.5494
ang: -30.4754 grad
 Xc: 5.0581924912269    2.9518671743604

example_ellipse_getellipse_elipse1.png

example_ellipse_getellipse_elipse2.png

Código example_ellipse_getellipse.cpp:

 
#include <cmath>
#include <Pds/Ra>
    
    
    
int main(void)
{
    Pds::Ra::Randomize();
    Pds::Octave::DAspect=true;
    Pds::Octave::MarkerSize=1;
    unsigned int L=1000;
    
    Pds::Ellipse E(3,1.5,5,3,-Pds::Ra::Pi/6);
    std::cout<<"  a: "<<E.GetMajorAxis()<<"\n";
    std::cout<<"  b: "<<E.GetMinorAxis()<<"\n";
    std::cout<<"ang: "<<E.GetAngle()*180/Pds::Ra::Pi<<" grad\n";
    std::cout<<" Xc: "<<E.GetCentroid().T();
    
    Pds::Matrix Xeq=E.RandSamples(L);
    
    Pds::Octave::Plot::PointsX2D(Xeq,"testando.m","example_ellipse_getellipse_elipse1.png");
    
    std::cout<<std::endl;
 
    Pds::Ellipse E2=Pds::Ellipse::GetEllipse(Xeq);
    
    std::cout<<"  a: "<<E2.GetMajorAxis()<<"\n";
    std::cout<<"  b: "<<E2.GetMinorAxis()<<"\n";
    std::cout<<"ang: "<<E2.GetAngle()*180/Pds::Ra::Pi<<" grad\n";
    std::cout<<" Xc: "<<E2.GetCentroid().T();
    
    Pds::Matrix P=E2.RandPerimeter(L);
    Pds::Octave::Plot::PointsX2DY(  Pds::MergeVer({P,Xeq}),
                                    Pds::MergeVer({Pds::Ones(P.Nlin(),1),Pds::Zeros(Xeq.Nlin(),1)}),
                                    "testando.m","example_ellipse_getellipse_elipse2.png");
 
    return 0;
}