mirror/eigen
2
0
Fork 0
mirror of https://gitlab.com/libeigen/eigen.git synced 2024-12-21 07:19:46 +08:00
Code Issues Packages Projects Releases Wiki Activity

Digamma cleanup

Browse Source 
* Added permission from cephes author to use his code
* Cleanup in ArrayCwiseUnaryOps
This commit is contained in:
Eugene Brevdo 2016-01-15 16:32:21 -08:00
parent f2471f31e0
commit 6a75e7e0d5
2 changed files with 23 additions and 12 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

32
Eigen/src/Core/SpecialFunctions.h
View File

@ -13,6 +13,29 @@
namespace Eigen {
namespace internal {
// Parts of this code are based on the Cephes Math Library.
//
// Cephes Math Library Release 2.8: June, 2000
// Copyright 1984, 1987, 1992, 2000 by Stephen L. Moshier
//
// Permission has been kindly provided by the original author
// to incorporate the Cephes software into the Eigen codebase:
//
// From: Stephen Moshier
// To: Eugene Brevdo
// Subject: Re: Permission to wrap several cephes functions in Eigen
//
// Hello Eugene,
//
// Thank you for writing.
//
// If your licensing is similar to BSD, the formal way that has been
// handled is simply to add a statement to the effect that you are incorporating
// the Cephes software by permission of the author.
//
// Good luck with your project,
// Steve
namespace cephes {
/* polevl (modified for Eigen)
@ -178,11 +201,6 @@ struct digamma_impl<float> {
* message condition value returned
* psi singularity x integer <=0 INFINITY
*/
/*
Cephes Math Library Release 2.2: June, 1992
Copyright 1984, 1987, 1992 by Stephen L. Moshier
Direct inquiries to 30 Frost Street, Cambridge, MA 02140
*/
EIGEN_DEVICE_FUNC
static float run(float xx) {
float p, q, nz, x, s, w, y, z;
@ -297,10 +315,6 @@ struct digamma_impl<double> {
* psi singularity x integer <=0 INFINITY
*/
/*
* Cephes Math Library Release 2.8: June, 2000
* Copyright 1984, 1987, 1992, 2000 by Stephen L. Moshier
*/
double p, q, nz, s, w, y, z;
bool negative;

3
Eigen/src/plugins/ArrayCwiseUnaryOps.h
View File

@ -320,9 +320,6 @@ lgamma() const
}
/** \returns an expression of the coefficient-wise digamma (psi, derivative of lgamma).
*
* Example: \include Cwise_digamma.cpp
* Output: \verbinclude Cwise_digamma.out
*
* \sa cos(), sin(), tan()
*/