NAME

kww - Kohlrausch-Williams-Watts function

SYNOPSIS

#include <kww.h>

double kwwc (const double omega, const double beta );

double kwws (const double omega, const double beta );

double kwwp (const double omega, const double beta );

DESCRIPTION

Laplace-Fourier transform of the stretched exponential function exp(-t^beta).

kwwc returns: integral from 0 to +infinity dt cos(omega*t) exp(-t^beta)

kwws returns: integral from 0 to +infinity dt sin(omega*t) exp(-t^beta)

kwwp returns: primitive of kwwc: integral from 0 to omega dw' kwwc(w',beta)

For sufficiently small or large values of |omega|, series expansions are used; otherwise numeric integration is performed using a double-exponential transform.

Allowed parameter range: 0.1 <= beta <= 2.0. However, kwwc is not fully supported for 1.9 < beta < 2.0: For some omega the numeric integration will not attain full accuracy. In these cases, 0 is returned.

ERRORS

On error, a message is written to stderr, and the program exits with errno EDOM (beta out of range) or ENOMEM (workspace allocation failed) or ENOSYS (no convergence in the parameter range that ought to be fully supported; please report to the author).

RESOURCES

Project web site: https://jugit.fz-juelich.de/mlz/kww

REFERENCE

J Wuttke: Fourier Transform of the Stretched Exponential Function: Analytic Error Bounds, Double Exponential Transform, and Open-Source Implementation libkww. http://arxiv.org/abs/0911.4796

BUG REPORTS

Please report bugs to the author.

Failures of the cosine transform in the range 1.9 < beta < 2.0 are not considered bugs. Any other failed convergence should be reported, indicating the values of beta and omega with 17 digits precision.

AUTHOR

Joachim Wuttke <j.wuttke@fz-juelich.de>

COPYING

Copyright (C) 2009, 2012 Joachim Wuttke.

Software: GNU General Public License (GPL).

This documentation: Creative Commons Attribution Share Alike.