example_funcconvexhull_isinside.cpp

Programa para el testeo de las funciones.

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

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

Para executar o programa:

./example_funcconvexhull_isinside

Retornando por consola:

P:
-0.94587764374254  -0.074276470613795 
-0.81005665325097  0.32971316311961   
-0.46853830966565  0.7237385412323    
-0.22499284065561  0.91415284802865   
0.13486687659047   0.98928737034522   
0.32727364931594   0.84430528657711   
0.88845428213871   0.40241389833503   
0.98407706012208   0.081744992212274  
0.79740584492563   -0.55141916896748  
0.28988710711239   -0.85939151740605  
-0.095184385820843 -0.97667028055371  
-0.18451674477408  -0.95439511302598  
-0.25752690353316  -0.93611814032128  
-0.6136719848093   -0.63209217071165  
-0.94247473447699  -0.27121230693125  

example_funcconvexhull_isinside_data.png

example_funcconvexhull_isinside_points.png

Código example_funcconvexhull_isinside.cpp:

 
#include <Pds/Ra>
    
int main(void)
{
    // Set a random seed to rand function.
    Pds::Ra::Randomize();
    
    unsigned int L=100;
    
    // Creating the random points
    Pds::Circle C1(0,0,1);
    Pds::Matrix X = C1.RandSamples(L);
 
    // Plotting X
    Pds::Octave::Plot::PointsX2D(X,"testando_data.m","example_funcconvexhull_isinside_data.png");
    
    //Convex hull
    Pds::Matrix P=Pds::ConvexHull::IncrementalAlgorithm(X);
    P.Print("P:\n");
    
    Pds::Circle C2(0,0,2);
    X = C2.RandSamples(10*L);
    Pds::Matrix R=Pds::ConvexHull::IsInside(P,X);
 
    // Plotting X and P
    Pds::Octave::Plot::PointsX2DY(X,R,"testando.m","example_funcconvexhull_isinside_points.png");
    
    return 0;
}