example_ellipse_create.cpp

 
#include <cmath>
#include <Pds/Ra>
    
int main(void)
{
    Pds::Octave::DAspect=true;
    Pds::Matrix Xc;
    
    Pds::Ellipse Ellipse(3,2,1,2,-3*Pds::Ra::Pi/4);
    Ellipse.Print("Ellipse:\n");
    
    std::cout<<"  a: "<<Ellipse.GetMajorAxis()<<"\n";
    std::cout<<"  b: "<<Ellipse.GetMinorAxis()<<"\n";
    std::cout<<"ang: "<<Ellipse.GetAngle()*180/Pds::Ra::Pi<<" grad\n";
    std::cout<<"  a: "<<Ellipse.GetMajorVector().T();
    std::cout<<"  b: "<<Ellipse.GetMinorVector().T();
    std::cout<<" Xc: "<<Ellipse.GetCentroid().T();
    std::cout<<"Xc1: "<<Ellipse.GetCentroid1().T();
    std::cout<<"Xc2: "<<Ellipse.GetCentroid2().T();
    
    
    Xc=Ellipse.RandSamples(1000);
    Pds::Octave::Plot::PointsX2D(Xc,"testando.m","example_ellipse_create_randsamples.png");
    
    Xc=Ellipse.RandPerimeter(64);
    Pds::Octave::Plot::PointsX2D(Xc,"testando.m","example_ellipse_create_randperimeter.png");
    
    std::cout<<"\n";
    
    double A;double B;double C;
    double D;double E;double F;
    Ellipse.GetQuadraticForm(A,B,C,D,E,F);
    
    Pds::Ellipse Ellipse2(A,B,C,D,E,F);
    std::cout<<"  a: "<<Ellipse2.GetMajorAxis()<<"\n";
    std::cout<<"  b: "<<Ellipse2.GetMinorAxis()<<"\n";
    std::cout<<"ang: "<<Ellipse2.GetAngle()*180/Pds::Ra::Pi<<" grad\n";
    std::cout<<"  a: "<<Ellipse2.GetMajorVector().T();
    std::cout<<"  b: "<<Ellipse2.GetMinorVector().T();
    std::cout<<" Xc: "<<Ellipse2.GetCentroid().T();
    std::cout<<"Xc1: "<<Ellipse2.GetCentroid1().T();
    std::cout<<"Xc2: "<<Ellipse2.GetCentroid2().T();
}