sergiotarxz/libvips
1
0
Fork 0
Code Issues Pull Requests Releases Wiki Activity

make log/log10 zero-avoiding

Browse Source 
stops -inf appearing in images
This commit is contained in:
John Cupitt 2015-08-11 13:59:11 +01:00
parent b3085f9f09
commit 35232e751d
1 changed files with 9 additions and 2 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

11
libvips/arithmetic/math.c
View File

@ -23,6 +23,8 @@
* - im_sintra() adapted to make math.c
* 4/11/11
* - redone as a class
* 11/8/15
* - log/log10 zero-avoid
*/
/*
@ -140,6 +142,11 @@ vips_math_build( VipsObject *object )
*/
#define EXP10( X ) (pow( 10.0, (X) ))
/* Zero-avoiding log, cf. zero-avoiding behaviour of /.
*/
#define LOGZ( X ) ((X) == 0.0 ? 0.0 : log( X ))
#define LOGZ10( X ) ((X) == 0.0 ? 0.0 : log10( X ))
static void
vips_math_buffer( VipsArithmetic *arithmetic,
VipsPel *out, VipsPel **in, int width )
@ -157,8 +164,8 @@ vips_math_buffer( VipsArithmetic *arithmetic,
case VIPS_OPERATION_MATH_ASIN: SWITCH( ADSIN ); break;
case VIPS_OPERATION_MATH_ACOS: SWITCH( ADCOS ); break;
case VIPS_OPERATION_MATH_ATAN: SWITCH( ADTAN ); break;
case VIPS_OPERATION_MATH_LOG: SWITCH( log ); break;
case VIPS_OPERATION_MATH_LOG10: SWITCH( log10 ); break;
case VIPS_OPERATION_MATH_LOG: SWITCH( LOGZ ); break;
case VIPS_OPERATION_MATH_LOG10: SWITCH( LOGZ10 ); break;
case VIPS_OPERATION_MATH_EXP: SWITCH( exp ); break;
case VIPS_OPERATION_MATH_EXP10: SWITCH( EXP10 ); break;