example_triangle_evalsamples.cpp

 
#include <cmath>
#include <Pds/Ra>
    
int main(void)
{
    Pds::Octave::DAspect=true;
    unsigned int L=10000;
    Pds::Matrix Xc(L,2);
    
    Pds::Point2D p1(1.0,1.0);
    Pds::Point2D p2(1.5,2.0);
    Pds::Point2D p3(2.0,1.0);
    
    Pds::Triangle T(p1,p2,p3);
    T.Print("T: ");
    
    Xc.FillRandU();
    Xc=Xc*3;
    
    Pds::Vector F=T.EvalSamples(Xc);
    std::cout<<"F.Min()"<<F.Min()<<"\n";
    std::cout<<"F.Max()"<<F.Max()<<"\n";
    
    Pds::Octave::Plot::ScatterX2DY(Xc,F,"testando.m","example_triangle_evalsamples_scatterx2dy.png");
    
    std::cout<<"T.Evaluate(p1):"<<T.Evaluate(p1)<<"\n";
    std::cout<<"T.Evaluate(p2):"<<T.Evaluate(p2)<<"\n";
    std::cout<<"T.Evaluate(p3):"<<T.Evaluate(p3)<<"\n";
    std::cout<<"T.Evaluate((p1+p2)/2):"<<T.Evaluate((p1+p2)/2)<<"\n";
    std::cout<<"T.Evaluate((p2+p3)/2):"<<T.Evaluate((p2+p3)/2)<<"\n";
    std::cout<<"T.Evaluate((p3+p1)/2):"<<T.Evaluate((p3+p1)/2)<<"\n";
    std::cout<<"T.Evaluate((p1+p2+p3)/3):"<<T.Evaluate((p1+p2+p3)/3)<<"\n";
    std::cout<<"T.Evaluate({0,0}):"<<T.Evaluate({0,0})<<"\n";
    std::cout<<"T.Evaluate({2,2}):"<<T.Evaluate({2,2})<<"\n";
}