GaussJacobiQuadCCompat.f90

Go to the documentation of this file.

! BEGIN_HEADER
! -----------------------------------------------------------------------------
! Gauss-Jacobi Quadrature Implementation
! Authors: Rohit Goswami <rgoswami[at]ieee.org>
! Source: GaussJacobiQuad Library
! License: MIT
! GitHub Repository: https://github.com/HaoZeke/GaussJacobiQuad
! Date: 2023-08-28
! Commit: c442f77
! -----------------------------------------------------------------------------
! This code is part of the GaussJacobiQuad library, providing an efficient
! implementation for Gauss-Jacobi quadrature nodes and weights computation.
! -----------------------------------------------------------------------------
! To cite this software:
! Rohit Goswami (2023). HaoZeke/GaussJacobiQuad: v0.1.0.
! Zenodo: https://doi.org/10.5281/ZENODO.8285112
! ---------------------------------------------------------------------
! END_HEADER
!
!! @details This module provides a C-compatible interface to the Gauss-Jacobi quadrature methods
!! implemented in the GaussJacobiQuad module. It's designed to make the Fortran implementations
!! accessible from C/C++ code. It wraps the original Fortran subroutines with ones that are C-compatible,
!! using the `iso_c_binding` module for the necessary type conversions.
module gaussjacobiquadccompat
use iso_c_binding, only: c_double, c_int
use gaussjacobiquad, only: gauss_jacobi_rec, gauss_jacobi_gw, gauss_jacobi_algo665
implicit none
contains
 
subroutine gauss_jacobi_rec_c(npts, alpha, beta, x, wts) bind(C, name="gauss_jacobi_rec_c")
    integer(c_int), intent(in) :: npts
    real(c_double), intent(in) :: alpha, beta
    real(c_double), intent(out) :: x(npts), wts(npts)
    call gauss_jacobi_rec(npts, alpha, beta, x, wts)
end subroutine gauss_jacobi_rec_c
 
subroutine gauss_jacobi_gw_c(npts, alpha, beta, x, wts) bind(C, name="gauss_jacobi_gw_c")
    integer(c_int), intent(in) :: npts
    real(c_double), intent(in) :: alpha, beta
    real(c_double), intent(out) :: x(npts), wts(npts)
    call gauss_jacobi_gw(npts, alpha, beta, x, wts)
end subroutine gauss_jacobi_gw_c
 
subroutine gauss_jacobi_algo665_c(npts, alpha, beta, x, wts) bind(C, name="gauss_jacobi_algo665_c")
    integer(c_int), intent(in) :: npts
    real(c_double), intent(in) :: alpha, beta
    real(c_double), intent(out) :: x(npts), wts(npts)
    call gauss_jacobi_algo665(npts, alpha, beta, x, wts)
end subroutine gauss_jacobi_algo665_c
 
end module gaussjacobiquadccompat