GaussJacobiQuadCInterp.c File Reference
+ Include dependency graph for GaussJacobiQuadCInterp.c:

Go to the source code of this file.

Functions

void gauss_jacobi_c (int *npts, double *alpha, double *beta, double x[], double wts[], const char *method)
 

Function Documentation

◆ gauss_jacobi_c()

void gauss_jacobi_c ( int *  npts,
double *  alpha,
double *  beta,
double  x[],
double  wts[],
const char *  method 
)

Definition at line 21 of file GaussJacobiQuadCInterp.c.

21  {
22  if (*npts <= 0) {
23  fprintf(stderr, "Error: npts must positive\n");
24  exit(EXIT_FAILURE);
25  }
26 
27  if (*alpha <= -1.0f) {
28  fprintf(stderr, "Error: alpha must be greater than -1\n");
29  exit(EXIT_FAILURE);
30  }
31 
32  if (*beta <= -1.0f) {
33  fprintf(stderr, "Error: beta must be greater than -1\n");
34  exit(EXIT_FAILURE);
35  }
36 
37  if (strcmp(method, "rec") == 0) {
38  gauss_jacobi_rec_c(npts, alpha, beta, x, wts);
39  } else if (strcmp(method, "gw") == 0) {
40  gauss_jacobi_gw_c(npts, alpha, beta, x, wts);
41  } else if (strcmp(method, "algo665") == 0) {
42  gauss_jacobi_algo665_c(npts, alpha, beta, x, wts);
43  } else {
44  fprintf(stderr, "Error: Unknown method specified: %s\n", method);
45  fprintf(stderr, "Supported methods: 'rec', 'gw', 'algo665'\n");
46  exit(EXIT_FAILURE);
47  }
48 }
void gauss_jacobi_algo665_c(int *npts, double *alpha, double *beta, double x[], double wts[])
void gauss_jacobi_gw_c(int *npts, double *alpha, double *beta, double x[], double wts[])
void gauss_jacobi_rec_c(int *npts, double *alpha, double *beta, double x[], double wts[])