GaussJacobiQuadCInterp.h File Reference

#include <stdio.h>
#include <stdlib.h>
#include <string.h>

Include dependency graph for GaussJacobiQuadCInterp.h:

This graph shows which files directly or indirectly include this file:

Go to the source code of this file.

Functions
void	gauss_jacobi_rec_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_algo665_c (int *npts, double *alpha, double *beta, double x[], double wts[])
void	gauss_jacobi_c (int *npts, double *alpha, double *beta, double x[], double wts[], const char *method)

Function Documentation

gauss_jacobi_algo665_c()

void gauss_jacobi_algo665_c	(	int *	npts,
		double *	alpha,
		double *	beta,
		double	x[],
		double	wts[]
	)

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.

                                                                                                          {
    if (*npts <= 0) {
        fprintf(stderr, "Error: npts must positive\n");
        exit(EXIT_FAILURE);
    }
 
    if (*alpha <= -1.0f) {
        fprintf(stderr, "Error: alpha must be greater than -1\n");
        exit(EXIT_FAILURE);
    }
 
    if (*beta <= -1.0f) {
        fprintf(stderr, "Error: beta must be greater than -1\n");
        exit(EXIT_FAILURE);
    }
 
    if (strcmp(method, "rec") == 0) {
        gauss_jacobi_rec_c(npts, alpha, beta, x, wts);
    } else if (strcmp(method, "gw") == 0) {
        gauss_jacobi_gw_c(npts, alpha, beta, x, wts);
    } else if (strcmp(method, "algo665") == 0) {
        gauss_jacobi_algo665_c(npts, alpha, beta, x, wts);
    } else {
        fprintf(stderr, "Error: Unknown method specified: %s\n", method);
        fprintf(stderr, "Supported methods: 'rec', 'gw', 'algo665'\n");
        exit(EXIT_FAILURE);
    }
}

gauss_jacobi_gw_c()

void gauss_jacobi_gw_c	(	int *	npts,
		double *	alpha,
		double *	beta,
		double	x[],
		double	wts[]
	)

gauss_jacobi_rec_c()

void gauss_jacobi_rec_c	(	int *	npts,
		double *	alpha,
		double *	beta,
		double	x[],
		double	wts[]
	)