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;
}