basis.h

//
// Created by milinda on 1/12/17.
//

#ifndef SFCSORTBENCH_JACOBIPOLYNOMIAL_H
#define SFCSORTBENCH_JACOBIPOLYNOMIAL_H


#include "assert.h"
#include <iostream>
#include <math.h>
#include <cmath>
#include "cstring"
#ifdef WITH_BLAS_LAPACK
#include "lapac.h"
#endif

#ifndef WITH_BLAS_LAPACK
template <typename T>
void printArray_1D(T *a, int length)
{
    for (int i = 0; i < length; i++) { std::cout<<a[i]<<" "; }
    std::cout<<std::endl;
}


template <typename T>
void printArray_2D(T *a, int length1,int length2)
{
    for (int i = 0; i < length1; i++) {
        for (int j = 0; j < length2; j++) {
            std::cout << a[i * length2 + j] << " ";
        }
        std::cout<<std::endl;
    }
    std::cout<<std::endl;
}

#endif



namespace basis {

    void jacobip(double alpha, double beta, unsigned int N,
                 double *x, double *p , unsigned int np);

    void gradjacobip(double alpha, double beta, int N,double *x, double *dp, unsigned int np);

    void jacobiglq(double alpha, double beta, int N,
                   double *x, double *w);

    void jacobigq(double alpha, double beta, int N,
                  double *x, double *w);



    void lagrange(const double * x0, int N, int at, const double * x, double* px,  int m); 



}

#endif //SFCSORTBENCH_JACOBIPOLYNOMIAL_H