nannumel.c
No OneTemporary
Actions

Subscribers

None

File Metadata

	#include "nanaccum.h"

	/* What is C a horrible language. To hoops we have to get through to make this
	* work for different data types... Since overloading does not work, we
	* generate type-specific functions: */
	#define fname(name, suffix) name ## _ ## suffix
	#define nanstat_template(TYPE)\
	double fname(nanstat, TYPE)(int n, TYPE *x0, mwSize stride) \
	{\
	int i; int count = 0;\
	for (i = 0; i < n; ++i) {\
	if (!isnan(x0[i * stride]))\
	count += 1;\
	}\
	return count;\
	}

	nanstat_template(float);
	nanstat_template(double);
	nanstat_template(int8_T); nanstat_template(uint8_T);
	nanstat_template(int16_T); nanstat_template(uint16_T);
	nanstat_template(int32_T); nanstat_template(uint32_T);
	nanstat_template(int64_T); nanstat_template(uint64_T);

	#include "nanaccum.c" /* takes care of generic nan-handling of arrays :D. */