example_funckernel.cpp

 
#include <iostream>
#include <Pds/Ra>
#include <Pds/Ml>
    
int main(void)
{
    Pds::Matrix X="1 2\n-2 1\n2 -1";
    X.Print("X:\n");
 
    Pds::Matrix Kernel;
    Pds::Matrix ID;
    
    // Polinomial
    Kernel=Pds::Kernel::Polynomial(X,3);
    Kernel.Print("\nKernel Poly\n");
    
    ID=Pds::Kernel::PolynomialId(X.Ncol(),3);
    ID.T().Print("\nID.T()\n");    
    
    // Fourier
    Kernel=Pds::Kernel::FourierSeries(X,Pds::Vector("3 3"),1);
    Kernel.Print("\nKernel Fseries\n");
    
    return 0;
}