example_triangle_evalsamples.cpp

Programa para el testeo de las funciones.

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

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

Para executar o programa:

./example_triangle_evalsamples

Retornando por consola:

T: 1    1       1.5     2       2       1
F.Min()-0.321616
F.Max()3.94703
T.Evaluate(p1):-0
T.Evaluate(p2):-0
T.Evaluate(p3):-0
T.Evaluate((p1+p2)/2):-0
T.Evaluate((p2+p3)/2):-0
T.Evaluate((p3+p1)/2):-0
T.Evaluate((p1+p2+p3)/3):-0.333333
T.Evaluate({0,0}):3
T.Evaluate({2,2}):1

example_triangle_evalsamples_scatterx2dy.png

Código 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";
}