Положительная граница машинного нуля

@Sh@dow777 · Регистрация: 10.12.2013

Привет всем.
Такое задание: положительную границу машинного нуля E можно определить, реализуя алгоритм: вычислить такое наименьшее $https://www.cyberforum.ru/cgi-bin/latex.cgi?{E}_{k}$ , для которого выполняется условие $https://www.cyberforum.ru/cgi-bin/latex.cgi?{E}_{k}$ > 0, причем $https://www.cyberforum.ru/cgi-bin/latex.cgi?{E}_{k}$ = $https://www.cyberforum.ru/cgi-bin/latex.cgi?{E}_{k-1}$ / 2, $https://www.cyberforum.ru/cgi-bin/latex.cgi?{E}_{0}$ = 1, k = 1,...m. Найти положительную границу машинного нуля для вещественных типов float и double.

Написал программу:

C++

#include <iostream>
#include <vector>
using namespace std;
 
void null(){
    vector<float> float1(10000000);
    vector<double> double1(10000000);
 
    float1[0] = 1.0;
    double1[0] = 1.0;
 
    float min = float1[0];
    double min1 = double1[0];
 
    int k = 1;
    int f = 1;
 
    for(;;){
        float1[k] = float1[k-1] / 2;
        if(float1[k] < 0)
            break;
        min = float1[k];
        k++;
    }
 
    for(;;){
        double1[f] = double1[f-1] / 2;
        if(double1[f] < 0)
            break;
        min1 = double1[f];
        f++;
    }
 
    cout << "Float: " << min << endl;
    cout << "Double: " << min1 << endl;
}
 
int main(){
    null();
 
    system("pause");
 
    return 0;
}

Вроде все верно...бесконечные циклы...break, если граница найдена. Но выдает ошибку "vector subscript out of range".
Что я делаю не так? Понятно, что выхожу за пределы векторов, но как тогда правильно задать их размер?
Заранее спасибо.

@vua72 · 08.11.2015, 07:37

никак, можете использовать push_back, интересно, а какой размер этих векторов вы все-таки предполагаете, когда начинаете алгоритм с 1?
изучайте свой код и алгоритм

C++

#include <iostream>
#include <iomanip>
#include <vector>
using namespace std;
 
void null(){
    vector<float> float1;
   // vector<double> double1();
 
    float1.push_back(1.0);
    //double1.push_back(1.0);
 
    float min = float1[0];
   // double min1 = double1[0];
 
    int k = 1;
   // int f = 1;
 
    for(;;){
        float1.push_back (float1[k-1] / 2);
        if(float1[k] < 0)
            break;
        min = float1[k];
        k++;
        std::cout << k << ":   "<<scientific<<setprecision(8)<<min<< std::endl;
    }
 
    //for(;;){
        //double1[f].push_back(double1[f-1] / 2);
        //if(double1[f] < 0)
            //break;
        //min1 = double1[f];
        //f++;
    //}
 
    cout << "Float: " << min << endl;
    //cout << "Double: " << min1 << endl;
}
 
int main(){
    null();
 
    system("pause");
 
    return 0;
}

Добавлено через 1 минуту

Не по теме:

только ctrl-C вовремя нажимай :)

@Fallenworld · 08.11.2015, 09:30

а не проще ли

C++

#define _FLOAT_
 
#ifdef _FLOAT_
  float current = 1.0, previous = 1.0;
#else
 double current = 1.0, previous = 1.0;
#endif
do{
previous = current;
current/=2.0;
}while(current>0);
cout<<previous;

@Sh@dow777 · 10.11.2015, 05:00 **[ТС]**

Fallenworld, Спасибо, помог ваш вариант. Только я его переделал чуть. А с векторами что-то бесконечно вычисляло

Потом еще всё таки добью вариант с векторами.

@Sh@dow777 17 / 17 / 6 Регистрация: 10.12.2013 Сообщений: 740
	10.11.2015, 05:00 [ТС]	4
	Fallenworld, Спасибо, помог ваш вариант. Только я его переделал чуть. А с векторами что-то бесконечно вычисляло Потом еще всё таки добью вариант с векторами. 0

Опции темы