example_circle_create.cpp

Programa para el testeo de las funciones.

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

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

Para executar o programa:

./example_circle_create

Retornando por consola:

Circ: 1 3       2
 r: 2
Xc: 1   3

Circ.IsInside(1,3): 1
Circ.IsInside(3,3): 1
Circ.IsInside(3,5): 0

F.T(): -2.8370398748975  -2.3460696642042  -3.8452365812432  -0.06878828756262 -3.5315438115844  -1.512507151821   -0.8587791711281  -3.0720849945456  -1.9730318576901  -0.310943392731   

example_circle_create_randsamples.png

example_circle_create_randperimeter.png

Código example_circle_create.cpp:

 
#include <cmath>
#include <Pds/Ra>
    
int main(void)
{
    Pds::Octave::DAspect=true;
    
    Pds::Matrix Xc;
    Pds::Circle Circ(1,3,2.0);
    Circ.Print("Circ: ");
    
    std::cout<<" r: "<<Circ.GetRadius()<<"\n";
    std::cout<<"Xc: "<<Circ.GetCentroid().T()<<"\n";
    std::cout<<"Circ.IsInside(1,3): "<<Circ.IsInside(1,3)<<"\n";
    std::cout<<"Circ.IsInside(3,3): "<<Circ.IsInside(3,3)<<"\n";
    std::cout<<"Circ.IsInside(3,5): "<<Circ.IsInside(3,5)<<"\n";
    std::cout<<"Circ.IsInside(Pds::Point2D(1,3)): "<<Circ.IsInside(Pds::Point2D(1,3))<<"\n";
    std::cout<<"Circ.IsInside(Pds::Point2D(3,3)): "<<Circ.IsInside(Pds::Point2D(3,3))<<"\n";
    std::cout<<"Circ.IsInside(Pds::Point2D(3,5)): "<<Circ.IsInside(Pds::Point2D(3,5))<<"\n";
    std::cout<<"Circ.Evaluate(1,3): "<<Circ.Evaluate(1,3)<<"\n";
    std::cout<<"Circ.Evaluate(3,3): "<<Circ.Evaluate(3,3)<<"\n";
    std::cout<<"Circ.Evaluate(3,5): "<<Circ.Evaluate(3,5)<<"\n";
    std::cout<<"Circ.Evaluate(Pds::Point2D(1,3)): "<<Circ.Evaluate(Pds::Point2D(1,3))<<"\n";
    std::cout<<"Circ.Evaluate(Pds::Point2D(3,3)): "<<Circ.Evaluate(Pds::Point2D(3,3))<<"\n";
    std::cout<<"Circ.Evaluate(Pds::Point2D(3,5)): "<<Circ.Evaluate(Pds::Point2D(3,5))<<"\n";
    
    Xc=Circ.RandSamples(200);
    Pds::Octave::Plot::PointsX2D(Xc,"testando.m","example_circle_create_randsamples.png");
    
    Xc=Circ.RandPerimeter(64);
    Pds::Octave::Plot::PointsX2D(Xc,"testando.m","example_circle_create_randperimeter.png");
    
    std::cout<<std::endl;
    Xc=Circ.RandSamples(10);
    Pds::Vector F=Circ.EvalSamples(Xc);
    F.T().Print("F.T(): ");
}