test_vectorfunc.cpp

 
#include <iostream>
#include <Pds/Ra>
 
 
    
int main(void)
{
    
    Pds::Vector W=Pds::LinSpace(0,10,5);
    std::cout<<"W:\n"<<W;
    
 
    Pds::Matrix X2 = Pds::X2D();
    Pds::Matrix Y2 = Pds::Y2D();
    
    X2.Print("\nX2:\n");
    Y2.Print("\nY2:\n");
    
    Pds::Matrix X = Pds::X3D();
    Pds::Matrix Y = Pds::Y3D();
    Pds::Matrix Z = Pds::Z3D();
    
    X.Print("\nX:\n");
    Y.Print("\nY:\n");
    Z.Print("\nZ:\n");
    
    Pds::Vector L="4 1 2 3";
    Pds::Vector P="4 4 1";
    double d;
    
    d=Pds::DistancePointToPlane(P,L,true);
    
    P.T().Print("P.T()");
    L.T().Print("L.T()");
    std::cout<<"d: "<<d<<std::endl;
    
    L="0 0 0 1";
    P="0 0 -1";
    d=Pds::DistancePointToPlane(P,L,true);
    
    P.T().Print("P.T()");
    L.T().Print("L.T()");
    std::cout<<"d: "<<d<<std::endl;
    
    return 0;
}