с++

[Александе еть я]

подскажите пожалуйста как можно проверить в с++ целое число или дробное... в паскале есть функция транк...если такая в си..заранее спасибо..

[ds.Dante]

Теоретически можно так:

float x;
if (x == (int)x)
    ...

Но есть один нюанс:сопроцессор всегда считает с некоторой погрешностью, поэтому после
x = 0.5+0.5;
в результате может получиться 0.99999999 или 1.00000001. Лучше делать как-нибудь так:

#define epsilon 0.00001

...

float x = 0.5+0.5;
if (abs(x-(int)(x+epsilon)) < 2*epsilon)
    ...

Если сделать просто (int)x, то 0.99999999 округлится до нуля, поэтому тут нужно x+epsilon. Если такая проверка проводится много раз, её можно выделить в отдельную функцию bool IsIntegral (double x).

[Vago]

Сообщение от ds.Dante

#define epsilon 0.00001

В хидере float.h для этого существуют предопределённые константы FLT_EPSILON и DBL_EPSILON .

[ds.Dante]

FLT_EPSILON и DBL_EPSILON - это погрешность при значениях, близких к 1.

То есть, если у нас выражение "x = 1000.5 + 0.5", то эти константы окажутся слишком маленькими. В качестве эпсилон надо взять максимальное число, которое не вносит существенной разницы в результат (которое практически можно считать целым)

[Vago]

Сообщение от ds.Dante

FLT_EPSILON и DBL_EPSILON - это погрешность при значениях, близких к 1.

Совершенно верно. А что, выражение

fabs( x - (double)(int)(x) )

Может стать больше единицы? Тем более, на 10 порядков больше?...

#include <stdio.h>
#include <math.h>
#include <float.h>

bool isInt( double x ) {

   double de = fabs( x - (double)(int)(x) );
   return ( de >= 0. ? de < DBL_EPSILON : (1. - de) < DBL_EPSILON );

}


int main() {

   const char* czResult;

   double a = 
//1      -100.9999999999999995
//2      100.9999999995
//3      -10000.0005
//4      -100000000.0000000001
      100000000.5 + 0.5
      ;

   czResult = isInt( a ) ? "almost integer." : "clearly float.";
   printf( "%.13g is %s\n", a, czResult );

   return 0;

}

Другое дело, что принудительное преобразование в int сужает диапазон чисел, которые можно по такому алгоритму проверить. Но на это мы пока не обращаем внимания...

Added 11:44 CET
А кто, собственно, нас, вообще, заставляет пользоваться int?!...

bool isInt( double x ) {

   double de = fabs( x - floor( x ) );
   return ( de >= 0. ? de < DBL_EPSILON : (1. - de) < DBL_EPSILON );

}

[ds.Dante]

Сообщение от Vago

А что, выражение

fabs( x - (double)(int)(x) )

Может стать больше единицы?

Нет, но для ошибки достаточно, чтобы было больше DBL_EPSILON.


А теперь практическая проверка:

#include <iostream>
#include <cmath>
#include <float.h>
using namespace std;

#define EPSILON 0.00001

bool IsIntegral_byDante (double x)
{
	if (abs(x-(int)(x+EPSILON)) < 2*EPSILON)
		return true;
	else
		return false;
}

bool IsIntegral_byVago (double x)
{
   double de = fabs( x - (double)(int)(x) );
   return ( de >= 0. ? de < DBL_EPSILON : (1. - de) < DBL_EPSILON );
}

void main()
{
	double x=10;
	for (int i=0; i<100; i++)
		x -= 0.1;

	cout << x << endl;
	cout << IsIntegral_byDante (x) << endl;
	cout << IsIntegral_byVago (x) << endl;
	cin.get();
}

Вывод на экран:

1.87905e-014
1
0

Если взять второй твой вариант (с floor() вместо приведения), ничего не изменится. Если в моей функции подставить DBL_EPSILON, она тоже выдаст неверный результат.


UPD: я сделал наиболее точную и "правильную" функцию без грязных хаков. Этот вариант я и предлагаю автору топика (нужно только подобрать подходящий EPSILON).

bool IsIntegral (double x)
{
	double de;

	if (x-(int)x < 0.5)
		de = x - (int)x;
	else
		de = x - (int)(x+1);

	if (de < EPSILON)
		return true;
	else
		return false;
}

[Vago]

Сообщение от ds.Dante

нужно только подобрать подходящий EPSILON

Золотые слова!! ))) Только в реальной инженерной задаче сделать это ой как трудно! Что значит "неверный результат"? )) "Результат с требуемой точностью" - так здесь вернее будет сказать. Ты, просто, успокаиваешься при погрешности чуть меньше 10^-5, а я - 10^-21 Подставим, как компромисс, FLT_EPSILON (10^-7), и обе п/п выдадут "Ура, целый ноль!"