legendre_p.cu
#include "cuda_runtime.h"
#include <cuda.h>
#include <stdlib.h>
#include <iostream>
#include <sys/time.h>
#include <chrono>
#define M 1000000 // M个点等分[-1, 1]
#define N 100 // Legendre polynomial的阶数
// https://people.sc.fsu.edu/~jburkardt/cpp_src/legendre_polynomial/legendre_polynomial.cpp
// https://people.sc.fsu.edu/~jburkardt/cpp_src/legendre_polynomial/legendre_polynomial.html
float *p_polynomial_value ( int m, int n, float x[] )
{
int i;
int j;
float *v;
if ( n < 0 ) return NULL;
v = new float[m*(n+1)];
for ( i = 0; i < m; i++ ) v[i+0*m] = 1.0f;
if ( n < 1 ) return v;
for ( i = 0; i < m; i++ ) v[i+1*m] = x[i];
for ( j = 2; j <= n; j++ )
{
for ( i = 0; i < m; i++ )
{
v[i+j*m] = ( ( float ) ( 2 * j - 1 ) * x[i] * v[i+(j-1)*m]