/* signed fixed-point resolution 
Given a signed fixed-point format  */
int fixed_resolution(double epsilon)
/* signed fixed-point resolution 
Given a signed fixed-point format  */
int fixed_resolution(double epsilon)
{
	return (int)(epsilon ? ceil (log2(1 / epsilon)) : 0);
}