example_line2d_create.cpp

Programa para el testeo de las funciones.

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

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

Para executar o programa:

./example_line2d_create

Retornando por consola:

Line1: 0.785398 1
a: 0.707107     b: 0.707107     c: -1

Line2: 0.785398 1

Line3: 0.785398 1

example_line2d_create_randsamples.png

example_line2d_create_getline.png

Código example_line2d_create.cpp:

 
#include <cmath>
#include <Pds/Ra>
    
int main(void)
{
    double a;    double b;    double c;
    Pds::Octave::DAspect=true;
    Pds::Matrix X;
    
    double angle=Pds::Ra::Pi/4;
    Pds::Line2D Line1(angle,1);
    Line1.Print("Line1: ");
    Line1.GetLinearForm(a,b,c);
    std::cout<<"a: "<<a<<"\tb: "<<b<<"\tc: "<<c<<"\n\n";
    
    
    Pds::Line2D Line2(a,b,c);
    Line2.Print("Line2: ");
    
    X=Line1.RandSamplesAngle(32,angle-Pds::Ra::Pi/4,angle+Pds::Ra::Pi/4);
    Pds::Octave::Plot::PointsX2D(X,"testando.m","example_line2d_create_randsamples.png");
    
    std::cout<<std::endl;
    
    Pds::Line2D Line3=Line1.GetLineCovEigMethod(X);
    Pds::Matrix P=Line3.RandSamplesAngle(32,angle-Pds::Ra::Pi/4,angle+Pds::Ra::Pi/4);
    Line3.Print("Line3: ");
    
    Pds::Octave::Plot::PointsX2DY(  Pds::MergeVer({P,X}),
                                    Pds::MergeVer({Pds::Ones(P.Nlin(),1),Pds::Zeros(X.Nlin(),1)}),
                                    "testando.m","example_line2d_create_getline.png");
}