Вызвано исключение по адресу 0x013C1891 в Jacobi.exe: 0xC0000005: нарушение прав доступа при чтении по адресу

@Penelent · 18.01.2016, 10:39. **Показов** 23760. **Ответов** 32

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

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
#include <math.h>
const double eps = 0.001; ///< желаемая точность 
 
/// N - размерность матрицы; A[N][N] - матрица коэффициентов, F[N] - столбец свободных членов,
/// X[N] - начальное приближение, ответ записывается также в X[N];
void Jacobi(int N, double** A, double* F, double* X)
{
    double* TempX = new double[N];
    double norm; // норма, определяемая как наибольшая разность компонент столбца иксов соседних итераций.
 
    do {
        for (int i = 0; i < N; i++) {
            TempX[i] = F[i];
            for (int g = 0; g < N; g++) {
                if (i != g)
                    TempX[i] -= A[i][g] * X[g];
            }
            TempX[i] /= A[i][i];
        }
        norm = fabs(X[0] - TempX[0]);
        for (int h = 0; h < N; h++) {
            if (fabs(X[h] - TempX[h]) > norm)
                norm = fabs(X[h] - TempX[h]);
            X[h] = TempX[h];
        }
    } while (norm > eps);
    delete[] TempX;
}
 
int main()
{
    const int N = 3;
    double A[N][N] = { { 5, -4, -1 }, { 4, 1, -2 }, { 3, 1, -5 } };
    double F[N] = { -2, 8, 10 };
    double X[N] = { 0 };
    Jacobi(N, (double**)A, F, X);
}

при выполнении кода выбивает ошибку на 16 строке

Добавлено через 32 секунды
з.ы. код взят отсюда

@_Valera_ · 18.01.2016, 11:04

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

ошибку на 16 строке

какую ошибку. увидел, вопрос должен быть в теме!!!!

@Penelent · 18.01.2016, 11:05 **[ТС]**

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

какую ошибку

в заголовке же "Вызвано исключение по адресу 0x013C1891 в Jacobi.exe: 0xC0000005: нарушение прав доступа при чтении по адресу"

@_Valera_ · 18.01.2016, 11:09

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

double** A,

скорей всего так массив передавать нельзя.
Да, так и есть.

@Penelent · 18.01.2016, 11:10 **[ТС]**

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

скорей всего так массив передавать нельзя

в смысле? и как же по вашему можно?

@_Valera_ · 18.01.2016, 11:14

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

в смысле? и как же по вашему можно?

не помню уже, на тут почитай:
Как в функцию передать параметр - двухмерный массив

@Penelent · 18.01.2016, 11:18 **[ТС]**

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

не помню уже, на тут почитай:
Как в функцию передать параметр - двухмерный массив

но там нет такого варианта со статическим массивом, че мне теперь чтобы матрицу объявить обязательно динамический объявлять? неудобно же, кода больше намного

@_Valera_ · 18.01.2016, 11:22

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

но там нет такого варианта со статическим массивом, че мне теперь чтобы матрицу объявить обязательно динамический объявлять? неудобно же, кода больше намного

ну поменяй параметры в функции.

@Penelent · 18.01.2016, 11:26 **[ТС]**

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

ну поменяй параметры в функции.

да нафига ж мне их менять, мне надо чтоб так работало, что ты мне мозги паришь, не верю что статический массив туда нельзя передать

zer0mail · 18.01.2016, 11:35

Сделай так:

C++
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
int main()
{
    const int N = 3;
    double B[N][N] = { { 5, -4, -1 }, { 4, 1, -2 }, { 3, 1, -5 } };
    double F[N] = { -2, 8, 10 };
    double X[N] = { 0 };
    double **A=new double*[N];
    for (int i=0;i<N;++i) {
        A[i]=new double[N];
        for (int j=0;j<N;++j)
            A[i][j]=B[i][j];
    }
    Jacobi(N, A, F, X); //B так передавать нельзя, а A - можно 
    return 0;
}

Добавлено через 58 секунд
Статические и динамические 2-мерные массивы устроены по-разному

@Penelent · 18.01.2016, 11:38 **[ТС]**

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

Статические и динамические 2-мерные массивы устроены по-разному

я уже догадался, статические как в C# двумерные [,] т.е. последовательно, а динамические как в C# ступенчатые [][], т.е. память выделена в разных местах, так?

zer0mail · 18.01.2016, 11:44

Еще вариант, изменить описание функции:

C++
1
void Jacobi(int N, double A[][3], double* F, double* X)

Добавлено через 1 минуту
Есть и третий - работать в функции с А, как с одномерным массивом. Приводить неохота.

@Penelent · 18.01.2016, 11:44 **[ТС]**

zer0mail, а можно просто вспомнить что это все таки плюсы и заюзать std::vector<std::vector<double>>

zer0mail · 18.01.2016, 11:46

Penelent, вы уж определитесь:

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

zer0mail, а можно просто вспомнить что это все таки плюсы и заюзать std::vector<std::vector<double>>

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

да нафига ж мне их менять, мне надо чтоб так работало, что ты мне мозги паришь, не верю что статический массив туда нельзя передать

@Penelent · 18.01.2016, 12:11 **[ТС]**

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

Penelent, вы уж определитесь

был не прав, извини _Valera_, после C#, трудно вспоминается C++

Добавлено через 18 минут
а еще вопрос теперь, я переделал под плюсы, почему у меня в X все время нули остаются? или туда нужно передавать заполненный X? что за начальное приближение? кто нибудь разбирался в этом методе?

Добавлено через 14 секунд

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
#include <math.h>
#include <vector>
 
void Jacobi(std::vector<std::vector<double>> A, std::vector<double> F, std::vector<double> X, double eps)
{
    size_t N = A.size();
    std::vector<double> TempX(N, 0);
    double norm(0.0);
 
    do {
        for (int i = 0; i < N; i++) {
            TempX[i] = F[i];
            for (int g = 0; g < N; g++) {
                if (i != g)
                    TempX[i] -= A[i][g] * X[g];
            }
            TempX[i] /= A[i][i];
        }
        norm = fabs(X[0] - TempX[0]);
        for (int h = 0; h < N; h++) {
            if (fabs(X[h] - TempX[h]) > norm)
                norm = fabs(X[h] - TempX[h]);
            X[h] = TempX[h];
        }
    } while (norm > eps);
}
 
int main()
{
    std::vector<std::vector<double>> A{ { 5, -4, -1 },{ 4, 1, -2 },{ 3, 1, -5 } };
    std::vector<double> F{ -2, 8, 10 };
    std::vector<double> X(3, 0);
    Jacobi(A, F, X, 0.001);
}

Добавлено через 5 минут
а все понял, опять путаю все с шарпом, тут же надо явно указывать ссылку -_-

@Penelent · 18.01.2016, 12:14 **[ТС]**

а сейчас?

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
#include <math.h>
#include <vector>
 
void Jacobi(const std::vector<std::vector<double>> &A, const std::vector<double> &F, std::vector<double> &X, double eps)
{
    size_t N = A.size();
    std::vector<double> TempX(N, 0);
    double norm(0.0);
 
    do {
        for (int i = 0; i < N; i++) {
            TempX[i] = F[i];
            for (int g = 0; g < N; g++) {
                if (i != g)
                    TempX[i] -= A[i][g] * X[g];
            }
            TempX[i] /= A[i][i];
        }
        norm = fabs(X[0] - TempX[0]);
        for (int h = 0; h < N; h++) {
            if (fabs(X[h] - TempX[h]) > norm)
                norm = fabs(X[h] - TempX[h]);
            X[h] = TempX[h];
        }
    } while (norm > eps);
}
 
int main()
{
    std::vector<std::vector<double>> A{ { 5, -4, -1 },{ 4, 1, -2 },{ 3, 1, -5 } };
    std::vector<double> F{ -2, 8, 10 };
    std::vector<double> X(3, 0);
    Jacobi(A, F, X, 0.001);
}

почему у меня в X какая то вообще хрень выходит?

nord_v · 18.01.2016, 12:42

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

почему у меня в X какая то вообще хрень выходит?

У тебя, в do while(), norm растёт

@Penelent · 18.01.2016, 12:54 **[ТС]**

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

У тебя, в do while(), norm растёт

в каком месте? там везде только присвоения, роста нет, да и вообще алгоритм с википедии, врятли там не правильно сделали

nord_v · 18.01.2016, 13:52

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

в каком месте?

Перед while (norm > eps) выведи на консоль значение norm и посмотри, что с ним происходит.

nord_v · 18.01.2016, 13:56

У меня, вот что выводит.

Новые блоги и статьи Все статьи Все блоги /
SDL3 для Web (WebAssembly): Реализация движения на Box2D v3 - трение и коллизии с повёрнутыми стенами 8Observer8 20.02.2026 Содержание блога Box2D позволяет легко создать главного героя, который не проходит сквозь стены и перемещается с заданным трением о препятствия, которые можно располагать под углом, как верхнее. . .	Конвертировать закладки radiotray-ng в m3u-плейлист damix 19.02.2026 Это можно сделать скриптом для PowerShell. Использование . \СonvertRadiotrayToM3U. ps1 <path_to_bookmarks. json> Рядом с файлом bookmarks. json появится файл bookmarks. m3u с результатом. # Check if. . .	Семь CDC на одном интерфейсе: 5 U[S]ARTов, 1 CAN и 1 SSI Eddy_Em 18.02.2026 Постепенно допиливаю свою "многоинтерфейсную плату". Выглядит вот так: https:/ / www. cyberforum. ru/ blog_attachment. php?attachmentid=11617&stc=1&d=1771445347 Основана на STM32F303RBT6. На борту пять. . .	Камера Toupcam IUA500KMA Eddy_Em 12.02.2026 Т. к. у всяких "хикроботов" слишком уж мелкий пиксель, для подсмотра в ESPriF они вообще плохо годятся: уже 14 величину можно рассмотреть еле-еле лишь на экспозициях под 3 секунды (а то и больше),. . .
И ясному Солнцу zbw 12.02.2026 И ясному Солнцу, и светлой Луне. В мире покоя нет и люди не могут жить в тишине. А жить им немного лет.	«Знание-Сила» zbw 12.02.2026 «Знание-Сила» «Время-Деньги» «Деньги -Пуля»	SDL3 для Web (WebAssembly): Подключение Box2D v3, физика и отрисовка коллайдеров 8Observer8 12.02.2026 Содержание блога Box2D - это библиотека для 2D физики для анимаций и игр. С её помощью можно определять были ли коллизии между конкретными объектами и вызывать обработчики событий столкновения. . . .	SDL3 для Web (WebAssembly): Загрузка PNG с прозрачным фоном с помощью SDL_LoadPNG (без SDL3_image) 8Observer8 11.02.2026 Содержание блога Библиотека SDL3 содержит встроенные инструменты для базовой работы с изображениями - без использования библиотеки SDL3_image. Пошагово создадим проект для загрузки изображения. . .

@Penelent Заблокирован
	18.01.2016, 11:44 [ТС]
	zer0mail, а можно просто вспомнить что это все таки плюсы и заюзать `std::vector<std::vector<double>>` 0

nord_v 337 / 185 / 80 Регистрация: 22.08.2013 Сообщений: 724
	18.01.2016, 13:56
	У меня, вот что выводит. Миниатюры 0

Вызвано исключение по адресу 0x013C1891 в Jacobi.exe: 0xC0000005: нарушение прав доступа при чтении по адресу

Решение