SF_INTEGRATE

Description

SciFortran module for function integration

Quick access

Routines:

trapz(), trapz2d(), simps(), simps2d(), quad(), kronig(), get_quadrature_weights()

Used modules

  • gauss_quadrature

Subroutines and functions

interface  sf_integrate/trapz(f, a, b, dh, x[, a, b, n])
Parameters:

  • f (various shapes) [real, complex]

  • a [real]

  • b [real]

  • dh [real]

  • x (size(f)) [real]

Options:

n [integer]

interface  sf_integrate/trapz2d(func, xrange, yrange, nx, ny[, n0, iterative, threshold, dhx, dhy, xrange, yrange])
Parameters:

  • func (various shapes) [real, complex, external]

  • xrange (2) [real]

  • yrange (2) [real]

  • nx [integer]

  • ny [integer]

Options:

  • n0 [integer]

  • iterative [logical]

  • threshold [real]

  • dhx [real]

  • dhy [real]

interface  sf_integrate/simps(f, dh, a, b, x[, a, b, n])

SAMPLE:

Parameters:

  • f (various shapes) [real, complex]

  • dh [real]

  • a [real]

  • b [real]

  • x (size(f)) [real]

Options:

n [integer]

interface  sf_integrate/simps2d(func, xrange, yrange, nx, ny[, n0, iterative, threshold, dhx, dhy, xrange, yrange])
Parameters:

  • func (various shapes) [real, complex, external]

  • xrange (2) [real]

  • yrange (2) [real]

  • nx [integer]

  • ny [integer]

Options:

  • n0 [integer]

  • iterative [logical]

  • threshold [real]

  • dhx [real]

  • dhy [real]

interface  sf_integrate/quad(func, a, result, fsample, b[, b, epsabs, epsrel, key, inf, singular_endpoint, singular_points, cpole, alfa, beta, omega, weight_func, verbose, strict, ninterp])
Parameters:

  • func [external]

  • a [real]

  • result [real]

  • fsample (•) [real]

Options:

  • b [real]

  • epsabs [real]

  • epsrel [real]

  • key [integer] – =[1:6]”

  • inf [integer] – =[-1,1,2]”

  • singular_endpoint [logical] – T if singular_endpoint exists (QAGS)

  • singular_points (•) [real] – location of singular points in QAGP

  • cpole [real] – location of the pole in QAWC f(x)/x-cpole

  • alfa [real]

  • beta [real] – parameters for QAWS

  • omega [real] – frequency of the weight funcion in QAWF,QAWO

  • weight_func [integer] – =[1,2]”

  • verbose [logical]

  • strict [logical]

  • ninterp [integer]

function  sf_integrate/kronig(fi, wr, m)
Parameters:

  • fi (m) [real]

  • wr (m) [real]

Options:

m [integer]

Return:

fr (m) [real]

subroutine  sf_integrate/get_quadrature_weights(wt[, nrk])
Parameters:

wt (•) [real]

Options:

nrk [integer] – = 2,4”