Process returned -107. в среде codeblocks как расшифровать ошибку?

@12listov · Регистрация: 22.07.2014

Студворк — интернет-сервис помощи студентам

Помогите понять какую ошибку выдает? стандарт c99, при уменьшении dz до 0.001 работает но метод расчета предполагает большую точность

C++
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
#include <stdio.h>
 
int main()
{
    // Константы
    const double Vm = 22.4; // Молярный объем газа при нормальных условиях, м3/кмоль
    const double rho = 996.5;// Плотность воды при 25 градусах Цельсия, кг/м3
    const double M = 18;        // Молярная масса воды
 
    // Исходные данные
    double Ha = 20.5;           // Высота аппарата, м
    double Da = 0.6;        // Диаметр аппарата, м
    double Gv = 1480;       // Объемный расход воздуха, м3/час
    double y_in = 0.025;        // Концентрация ацетона в воздухе на входе в абсорбер, мольные доли
    double x_in = 0;            // Концентрация ацетона в воде на входе в абсорбер, мольные доли
    double y_out_lim = 0.02*y_in;// Ограничение на концентрацию ацетона в воздухе на выходе из абсорбера, мольные доли
    double KoyV = 0.65*3600;        // Объемный коэффициент массопередачи, отнесенный к движущей силе в газовой фазе, кмоль/(м3*час)
    double m = 1.68;            // Коэффициент распределения
    double Lv = 1;          // Начальное приближение расхода воды, м3/час
    double x_out1 = 0.004;      // Начальное приближение концентрации ацетона в воде на выходе из абсорбера
    double dLv = 0.1;           // Интервал сканирования по расходу воды, м3/час
    double dz = 0.0001;         // Интервал сканирования по высоте абсрбера, м
    double eps = 0.01;      // Заданная точность расчета x_out
int p;
    // Предварительные расчеты
    double S = 3.1416*Da*Da/4;      // Площадь поперечного сечения аппарата, м2
    double G = Gv/Vm*(273.15+25)/273.15;                // Мольный расход воздуха, кмоль/час
    double L = Lv*rho/M;                // Начальной приближение расхода воды, кмоль/час
    double dL = dLv*rho/M;
    int n = 205000;//Ha / dz;               // Число разбиений по высоте абсорбера
 
    // Переменные для расчета
    double y[n+1];          // Массив значений доли ацетона в воздухе по высоте абсорбера
    double x[n+1];          // Массив значений доли ацетона в воде по высоте абсорбера
    double x_out2 = 0;
    double y_out = 0;
 
    // Расчеты
 
    while(1)        // Считаем до тех пор, пока не убдет достигнута заданная выходная концентрация ацетона в воздухе
    {
        //x_out1 = 0.03f;
        y[0] = y_in;
        while(1)    // Считаем, пока не сойдется выходная концентрация ацетона в воде
        {
            x[0] = x_out1;
            for(int i = 0; i < n; i++)      // Цикл по всей высоте абсорбера
            {
                // Для расчетов используем метод Рунге-Кутта 4 порядка
                double k1 = -KoyV * S / G *(y[i] - m*x[i]);
                double k2 = -KoyV * S / G *(y[i] + dz*k1/2 - m*x[i]);
                double k3 = -KoyV * S / G *(y[i] + dz*k2/2 - m*x[i]);
                double k4 = -KoyV * S / G *(y[i] + dz*k3 - m*x[i]);
                //printf("%e %e %e\n", k1, k2, k3);
                y[i+1] = y[i] + dz*(k1/6 + k2/3 + k3/3 + k4/6); // Концентрация ацетона в воздухе на новом шаге
                
 
                k1 = -KoyV * S / L *(y[i] - m*x[i]);
                k2 = -KoyV * S / L *(y[i] - m*(x[i] + dz*k1/2));
                k3 = -KoyV * S / L *(y[i] - m*(x[i] + dz*k2/2));
                k4 = -KoyV * S / L *(y[i] - m*(x[i] + dz*k3));
                x[i+1] = x[i] + dz*(k1/6 + k2/3 + k3/3 + k4/6); // Концентрация ацетона в воде на новом шаге
                //x[i+1] = x[i] + G/L*(y[i]-y[i+1]);
                //printf("%f\n", x[i+1]);
                //getchar();
                if (i=1) {
                       getch();};
            }
            y_out = y[n];
            x_out2 = x_in + G/L*(y_in - y_out);
            printf("y_in = %f; y_out = %f; x_out2 = %f; x_out1 = %f\n", y_in, y_out, x_out2, x_out1);
            //getchar();
            if(((x_out2 - x_out1 >= 0) && (x_out2 - x_out1 <= eps)) ||  ((x_out2 - x_out1 < 0) && (x_out1 - x_out2 <= eps)))
            {
                break;
            }
            x_out1 = x_out2;
        }
        //getchar();
        if(y_out <= y_out_lim)
        {
            break;
        }
        printf("L + dL = %f + %f = ", L, dL);
        L += dL;
        printf("%f\n", L);
    }
    printf("x_out = %f; y_out = %f, L = %f\n", x_out2, y_out, L*M/rho);
    return 0;
}

lss · 23.03.2015, 19:24

Для начала... Вот это не сравнение, а присваивание:

C++
1
2
if (i=1) {
                       getch();};

@12listov · 23.03.2015, 19:37 **[ТС]**

lss, это можно удалить, как то в попыхах добавил,

C++
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
#include <stdio.h>
 
int main()
{
    // Константы
    const double Vm = 22.4; // Молярный объем газа при нормальных условиях, м3/кмоль
    const double rho = 996.5;// Плотность воды при 25 градусах Цельсия, кг/м3
    const double M = 18;        // Молярная масса воды
 
    // Исходные данные
    double Ha = 20.5;           // Высота аппарата, м
    double Da = 0.6;        // Диаметр аппарата, м
    double Gv = 1480;       // Объемный расход воздуха, м3/час
    double y_in = 0.025;        // Концентрация ацетона в воздухе на входе в абсорбер, мольные доли
    double x_in = 0;            // Концентрация ацетона в воде на входе в абсорбер, мольные доли
    double y_out_lim = 0.02*y_in;// Ограничение на концентрацию ацетона в воздухе на выходе из абсорбера, мольные доли
    double KoyV = 0.65*3600;        // Объемный коэффициент массопередачи, отнесенный к движущей силе в газовой фазе, кмоль/(м3*час)
    double m = 1.68;            // Коэффициент распределения
    double Lv = 1;          // Начальное приближение расхода воды, м3/час
    double x_out1 = 0.004;      // Начальное приближение концентрации ацетона в воде на выходе из абсорбера
    double dLv = 0.1;           // Интервал сканирования по расходу воды, м3/час
    double dz = 0.0001;         // Интервал сканирования по высоте абсрбера, м
    double eps = 0.01;      // Заданная точность расчета x_out
int p;
    // Предварительные расчеты
    double S = 3.1416*Da*Da/4;      // Площадь поперечного сечения аппарата, м2
    double G = Gv/Vm*(273.15+25)/273.15;                // Мольный расход воздуха, кмоль/час
    double L = Lv*rho/M;                // Начальной приближение расхода воды, кмоль/час
    double dL = dLv*rho/M;
    int n = 205000;//Ha / dz;               // Число разбиений по высоте абсорбера
 
    // Переменные для расчета
    double y[n+1];          // Массив значений доли ацетона в воздухе по высоте абсорбера
    double x[n+1];          // Массив значений доли ацетона в воде по высоте абсорбера
    double x_out2 = 0;
    double y_out = 0;
 
    // Расчеты
 
    while(1)        // Считаем до тех пор, пока не убдет достигнута заданная выходная концентрация ацетона в воздухе
    {
        //x_out1 = 0.03f;
        y[0] = y_in;
        while(1)    // Считаем, пока не сойдется выходная концентрация ацетона в воде
        {
            x[0] = x_out1;
            for(int i = 0; i < n; i++)      // Цикл по всей высоте абсорбера
            {
                // Для расчетов используем метод Рунге-Кутта 4 порядка
                double k1 = -KoyV * S / G *(y[i] - m*x[i]);
                double k2 = -KoyV * S / G *(y[i] + dz*k1/2 - m*x[i]);
                double k3 = -KoyV * S / G *(y[i] + dz*k2/2 - m*x[i]);
                double k4 = -KoyV * S / G *(y[i] + dz*k3 - m*x[i]);
                //printf("%e %e %e\n", k1, k2, k3);
                y[i+1] = y[i] + dz*(k1/6 + k2/3 + k3/3 + k4/6); // Концентрация ацетона в воздухе на новом шаге
                
 
                k1 = -KoyV * S / L *(y[i] - m*x[i]);
                k2 = -KoyV * S / L *(y[i] - m*(x[i] + dz*k1/2));
                k3 = -KoyV * S / L *(y[i] - m*(x[i] + dz*k2/2));
                k4 = -KoyV * S / L *(y[i] - m*(x[i] + dz*k3));
                x[i+1] = x[i] + dz*(k1/6 + k2/3 + k3/3 + k4/6); // Концентрация ацетона в воде на новом шаге
                //x[i+1] = x[i] + G/L*(y[i]-y[i+1]);
                //printf("%f\n", x[i+1]);
                //getchar();
            
            }
            y_out = y[n];
            x_out2 = x_in + G/L*(y_in - y_out);
            printf("y_in = %f; y_out = %f; x_out2 = %f; x_out1 = %f\n", y_in, y_out, x_out2, x_out1);
            //getchar();
            if(((x_out2 - x_out1 >= 0) && (x_out2 - x_out1 <= eps)) ||  ((x_out2 - x_out1 < 0) && (x_out1 - x_out2 <= eps)))
            {
                break;
            }
            x_out1 = x_out2;
        }
        //getchar();
        if(y_out <= y_out_lim)
        {
            break;
        }
        printf("L + dL = %f + %f = ", L, dL);
        L += dL;
        printf("%f\n", L);
    }
    printf("x_out = %f; y_out = %f, L = %f\n", x_out2, y_out, L*M/rho);
    return 0;
}

Но та же самая ошибка

lss · 23.03.2015, 19:40

У меня ошибку на этом выдаёт (слишком большие массивы):

C++
1
2
double y[n+1];          // Массив значений доли ацетона в воздухе по высоте абсорбера
double x[n+1];

Попробуй динамически память выделить:

C++
1
2
double* y = new double[n+1];          // Массив значений доли ацетона в воздухе по высоте абсорбера
double* x = new double[n+1];

@12listov · 23.03.2015, 19:51 **[ТС]**

lss, Спасибо прям навёл на мысль, я в c++ нуб, кроме как это нужно в программе что то менять?

lss · 23.03.2015, 20:03

Сообщение от 12listov

кроме как это нужно в программе что то менять?

Без понятия. С расчётами сам разбирайся.

@12listov · 23.03.2015, 22:34 **[ТС]**

lss, Я имел ввиду для выделения памяти, ладно когда то я это делал, вспомню.

lss · 23.03.2015, 22:48

Сообщение от 12listov

Я имел ввиду для выделения памяти,

Только это:

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

Попробуй динамически память выделить:
double* y = new double[n+1]; // Массив значений доли ацетона в воздухе по высоте абсорбера
double* x = new double[n+1];

Новые блоги и статьи Все статьи Все блоги /
Воспроизведение звукового файла с помощью SDL3_mixer при касании экрана Android 8Observer8 26.01.2026 Содержание блога SDL3_mixer - это библиотека я для воспроизведения аудио. В отличие от инструкции по добавлению текста код по проигрыванию звука уже содержится в шаблоне примера. Нужно только. . .	Установка Android SDK, NDK, JDK, CMake и т.д. 8Observer8 25.01.2026 Содержание блога Перейдите по ссылке: https:/ / developer. android. com/ studio и в самом низу страницы кликните по архиву "commandlinetools-win-xxxxxx_latest. zip" Извлеките архив и вы увидите. . .	Вывод текста со шрифтом TTF на Android с помощью библиотеки SDL3_ttf 8Observer8 25.01.2026 Содержание блога Если у вас не установлены Android SDK, NDK, JDK, и т. д. то сделайте это по следующей инструкции: Установка Android SDK, NDK, JDK, CMake и т. д. Сборка примера Скачайте. . .	Использование SDL3-callbacks вместо функции main() на Android, Desktop и WebAssembly 8Observer8 24.01.2026 Содержание блога Если вы откроете примеры для начинающих на официальном репозитории SDL3 в папке: examples, то вы увидите, что все примеры используют следующие четыре обязательные функции, а. . .
моя боль iceja 24.01.2026 Выложила интерполяцию кубическими сплайнами www. iceja. net REST сервисы временно не работают, только через Web. Написала за 56 рабочих часов этот сайт с нуля. При помощи perplexity. ai PRO , при. . .	Модель сукцессии микоризы anaschu 24.01.2026 Решили писать научную статью с неким РОманом	http://iceja.net/ математические сервисы iceja 20.01.2026 Обновила свой сайт http:/ / iceja. net/ , приделала Fast Fourier Transform экстраполяцию сигналов. Однако предсказывает далеко не каждый сигнал (см ограничения http:/ / iceja. net/ fourier/ docs ). Также. . .	http://iceja.net/ сервер решения полиномов iceja 18.01.2026 Выкатила http:/ / iceja. net/ сервер решения полиномов (находит действительные корни полиномов методом Штурма). На сайте документация по API, но скажу прямо VPS слабенький и 200 000 полиномов. . .

@12listov 3 / 3 / 0 Регистрация: 22.07.2014 Сообщений: 209
	23.03.2015, 19:51 [ТС]
	lss, Спасибо прям навёл на мысль, я в c++ нуб, кроме как это нужно в программе что то менять? 0

@12listov 3 / 3 / 0 Регистрация: 22.07.2014 Сообщений: 209
	23.03.2015, 22:34 [ТС]
	lss, Я имел ввиду для выделения памяти, ладно когда то я это делал, вспомню. 0