kernel_cram_coefficients.cpp File Reference

Coefficients for CRAM depletion solver. More...

#include <stdio.h>
#include "Origen/Solver/cram/src/kernel_cram.hpp"
#include "Standard/Interface/Reporter.h"

Macros
#define	NAME   "kernel_cram_coefficients.cpp"
	identifier for error messages More...

Functions
int	kernel_cram_coefficients (doublecomplex *alpha0, doublecomplex *alpha, doublecomplex *theta, int cram_order)

Detailed Description

Coefficients for CRAM depletion solver.

Author

Aarno Isotalo Moved to separate file because this gets long when/if high order coefficients for alternate forms become available.

Macro Definition Documentation

#define NAME   "kernel_cram_coefficients.cpp"

identifier for error messages

Referenced by kernel_cram_coefficients().

Function Documentation

int kernel_cram_coefficients	(	doublecomplex *	alpha0,
		doublecomplex *	alpha,
		doublecomplex *	theta,
		int	cram_order
	)

Set CRAM coefficients for partial fraction decomposition form.

Sets the poles and residues of chebyshev rational approximation on the negative real axis. Available orders are 4, 6, 8, 10, 12, 14, and 16.

Order 14 and 16 coefficients are from: Maria Pusa, "Correction to partial fraction decomposition coefficients for Chebyshev rational approximation on the negative real axis." arXiv:1206.2880v1 [math.NA] (2013). Others have been provided by Maria Pusa directly.

Parameters

[out]	alpha0,alpha,theta	The coefficients. alpha and theta must be preallocated to hold at least order/2 values
[in]	cram_order	The desired order of the CRAM approximation.

Returns

EXIT_SUCCESS or EXIT_FAILURE. On failed exit the results are unspecified.

References NAME.

Referenced by kernel_cram().