Форум программистов, компьютерный форум CyberForum.ru

При наличие корней в системе уравнений выдает ошибку - C++

Восстановить пароль Регистрация
 
 
Рейтинг: Рейтинг темы: голосов - 10, средняя оценка - 4.90
Larri
0 / 0 / 0
Регистрация: 29.09.2011
Сообщений: 16
29.10.2011, 20:07     При наличие корней в системе уравнений выдает ошибку #1
После выполнения программы (она решает систему линейных уравнений) для случая наличия корней появляется ошибка.
Вот код:
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
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
#include <stdio.h>
#include<stdlib.h>
#include<time.h>
#define N 4
typedef struct
{int*a;
bool l;//нулевой ли вектор
}lama;
int* WC(int *w1,int f,int n)//alfa na vector
{
    int i=0;
    int *w;
    w=(int*)malloc(n*sizeof(int));
    for(i=0;i<n;i++)
        w[i]=0;
    for(i=0;i<n;i++)
        w[i]=(w[i]+w1[i]*f)%2;
    return w;
}
int WW(int *w1,int *w2,int n)//skal proizv
{
    int i=0;
    int w=0;
    for(i=0;i<n;i++)
        w=(w+w1[i]*w2[i])%2;
    return w;
}
int* MW(int n, int *a, int *w)
{
    int i=0,j=0;
    int *w1;
    w1=(int*)malloc(n*sizeof(int));
    for (j=0;j<n;j++)
        w1[j]=0;
    for (j=0;j<n;j++)
        for(i=0;i<n;i++)
            w1[j]=(w1[j]+a[j*n+i]*w[i])%2;
    return w1;
}
lama form(int n,int *a,int *w1, int *w0)
{
    int i=0;
    int *w2;
    lama w;
    w.l=false;  
    int q,q1,e=0;
    w.a=(int*)malloc(n*sizeof(int));
    w2=(int*)malloc(n*sizeof(int));
    w2=MW(n,a,w1);
    q=(WW(MW(n,a,w2),w1,n))*(WW(w2,w1,n));  
    if(w0==NULL)
        for(i=0;i<n;i++)
        {
            w.a[i]=(w2[i]+WC(w1,q,n)[i])%2; //в двоичном сложении 1+1=1-1=0, если не учитывать перенос разряда
            e+=w.a[i];
        }
    else
    {
        q1=(WW(MW(n,a,w2),w0,n))*(WW(MW(n,a,w0),w0,n));
        for(i=0;i<n;i++)
        {
            w.a[i]=(w2[i]+WC(w1,q,n)[i]+WC(w0,q1,n)[i])%2;
            e+=w.a[i];          
        }
    }   
    if(e==0) w.l=true;
    free(w2);
    return w;
}
int*XX(int*a,int*w,int k,int n)//nepr
{
    int j,i,q;
    int*x;
    x=(int*)malloc(n*sizeof(int));  
    for(i=0;i<n;i++)
        x[i]=0;
    for(i=0;i<=k;i++)
    {       
        q=WW((w+i*n),w,n)*WW(MW(n,a,(w+i*n)),(w+i*n),n);
        for(j=0;j<n;j++)
            x[j]=(x[j]+q*w[i*n+j])%2;
    }   
    return x;
}
bool prov(int*a,int*x,int*b,int n)
{   
    int*b1,i,j=0;
    b1=(int*)malloc(n*sizeof(int));
    b1=MW(n,a,x);
    for(i=0;i<n;i++)
        if(b1[i]!=b[i]) j++;
    free(b1);
    if(j==0) return true;
    else return false;
}
int main()
{
    int *a,*b,*w,*x,i,j,k=0;
    lama p; 
    a=(int*)malloc(N*N*sizeof(int));
    b=(int*)malloc(N*sizeof(int));
    w=(int*)malloc(N*N*sizeof(int));
    /*srand(time(0));
    for(i=0;i<N;i++)
    {
        for(j=i;j<N;j++)
        {
            a[i*N+j]=rand()%2;
            if(i!=j) a[j*N+i]=a[i*N+j];
        }
        b[i]=rand()%2;
    }
    for (i=0; i<N; i++)
    {
        for (j=0; j<N; j++)
            printf("\t %i", a[i*N+j]);
        printf ("\n");
        }   */
    for(i=0;i<N;i++)//-пользовательское задание матрицы А
        for(j=0;j<N;j++)
        {
            printf("Enter a[%i,%i]\n",i,j);
            scanf("%i",&a[i*N+j]);
        }
    for(i=0;i<N;i++)
    {
        printf("Enter b[%i]\n",i);
        scanf("%i",&b[i]);
    }
    for (i=0; i<N; i++)
    {
        for (j=0; j<N; j++)
            printf("\t %i", a[i*N+j]);
        printf ("\n");
    }
    for(i=0;i<N;i++)    
        printf("b[%i]=%i\n",i,b[i]);
    for(j=0;j<N;j++)
        w[j]=b[j];
    i=1;
    do
    {
        if(i==1) p=form(N,a,b,NULL);
        else p=form(N,a,(w+(i-1)*N),(w+(i-2)*N));
        for(j=0;j<N;j++)//было j<=N
        {
            w[i*N+j]=p.a[j];//запись вектора W[i] в матрицу омег            
        }
        k++;
        i++;
    }
    while(!p.l);
    if(p.l)
    {
        for (i=0; i<=k; i++)
        {
            for (j=0; j<N; j++)
                printf("\t %i", w[i*N+j]);
            printf ("\n");
        }
        x=XX(a,w,k,N);      
        if(!prov(a,x,b,N)) printf("answer does not exist\n");
        else 
        {
            printf("answer is true\n");
            printf("answer: ");
            for(i=0;i<N;i++)
                printf("%i ",x[i]);
            printf("\n");
        }
    }
    else printf("algorithm does not work\n");
    free(a);
    free(b);
    free(w);
    return 0;
}
Лучшие ответы (1)
После регистрации реклама в сообщениях будет скрыта и будут доступны все возможности форума.
-=ЮрА=-
Заблокирован
Автор FAQ
29.10.2011, 20:20     При наличие корней в системе уравнений выдает ошибку #2
Цитата Сообщение от Larri Посмотреть сообщение
x=(int*)malloc(n*sizeof(int));
Цитата Сообщение от Larri Посмотреть сообщение
return x;
- после выхода из функции всё, что выделено внутри маллоком уничтожится. Правидь код не особо хочется, поэтому предлагаю взамен уже проверенный алгоритм Методом Гаусса решить систему n линейных алгебраических уравнений
Deviaphan
Делаю внезапно и красиво
Эксперт C++
 Аватар для Deviaphan
1283 / 1217 / 50
Регистрация: 22.03.2011
Сообщений: 3,744
29.10.2011, 20:33     При наличие корней в системе уравнений выдает ошибку #3
Цитата Сообщение от -=ЮрА=- Посмотреть сообщение
после выхода из функции всё, что выделено внутри маллоком уничтожится.
Юра, ты простыл?

C++
1
2
3
4
 for(i=0;i<n;i++)
                w[i]=0;
        for(i=0;i<n;i++)
                w[i]=(w[i]+w1[i]*f)%2;
Из-за первого цикла, уторой можно удалить за ненадобностью. (или исправить логику).
Дальше даже читать не буду. Исправляй логику.
Larri
0 / 0 / 0
Регистрация: 29.09.2011
Сообщений: 16
29.10.2011, 20:35  [ТС]     При наличие корней в системе уравнений выдает ошибку #4
Дело в том, что мне нужно решить систему конкретным алгоритмом
Deviaphan
Делаю внезапно и красиво
Эксперт C++
 Аватар для Deviaphan
1283 / 1217 / 50
Регистрация: 22.03.2011
Сообщений: 3,744
29.10.2011, 20:35     При наличие корней в системе уравнений выдает ошибку #5
Всё таки глянул немного... Всё, что выделяется в вотерклозете(WC), никогда не удаляется. Велкам, утечки.
Larri
0 / 0 / 0
Регистрация: 29.09.2011
Сообщений: 16
29.10.2011, 20:43  [ТС]     При наличие корней в системе уравнений выдает ошибку #6
Цитата Сообщение от Deviaphan Посмотреть сообщение
Из-за первого цикла, уторой можно удалить за ненадобностью. (или исправить логику).
Дальше даже читать не буду. Исправляй логику.
Там все нормально. Второй цикл - поиск скалярного произведения, а если убрать первый, программа работает неправильно.

Добавлено через 1 минуту
Цитата Сообщение от Deviaphan Посмотреть сообщение
Всё таки глянул немного... Всё, что выделяется в вотерклозете(WC), никогда не удаляется. Велкам, утечки.
А как исправить? Я так думаю, на это и ругается
-=ЮрА=-
Заблокирован
Автор FAQ
29.10.2011, 20:43     При наличие корней в системе уравнений выдает ошибку #7

Не по теме:

Deviaphan, я не простыл, мне влом с тобой холиворить и что то приводить ты не прав вот и всё


Larri, какой именно метод решения нужен?
Deviaphan
Делаю внезапно и красиво
Эксперт C++
 Аватар для Deviaphan
1283 / 1217 / 50
Регистрация: 22.03.2011
Сообщений: 3,744
29.10.2011, 20:46     При наличие корней в системе уравнений выдает ошибку #8
Цитата Сообщение от Larri Посмотреть сообщение
Там все нормально
w[i]=0;
w[i]=( 0 +w1[i]*f)%2;

Я сперва прочитал, вместо w1 тоже w, думал, вообще в ноль всё выражение. А так у тебя просто избыточность. Если w не должно суммироваться, то можешь смело удалить цикл с обнулением.

Вызов WC некорректен: выделенная память не освобождается. При каждом вызове WC формируется неудаляемый массив, но используется только один элемент... В общем, исправляй логику.

Добавлено через 1 минуту
Цитата Сообщение от -=ЮрА=- Посмотреть сообщение
ты не прав вот и всё
Юра, при разрушении указателя не освобождается память, на которую он указывает. И никакого холивара, так работает С++. И Си.
-=ЮрА=-
Заблокирован
Автор FAQ
29.10.2011, 20:50     При наличие корней в системе уравнений выдает ошибку #9
Larri, опишите алгоритм нахождения корней или дайте название метода иначе покидаю топик...
alkagolik
 Аватар для alkagolik
1510 / 616 / 79
Регистрация: 15.07.2011
Сообщений: 3,552
29.10.2011, 20:54     При наличие корней в системе уравнений выдает ошибку #10
Цитата Сообщение от -=ЮрА=- Посмотреть сообщение
- после выхода из функции всё, что выделено внутри маллоком уничтожится.
Цитата Сообщение от Deviaphan Посмотреть сообщение
Юра, ты простыл?
похоже да, и серьезно. Юрец, опровергни
C
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
#include <stdio.h>
#include <stdlib.h>
 
int *f(int n)
{
    int *x = malloc(n * sizeof(int));
    return x;
}
 
int main ()
{
    int *a, n = 10, i = 0;
    a = f(n);
 
    while (i < n)
        printf("%i ", a[i++]);
 
    free(a);
    return 0;
}
Larri
0 / 0 / 0
Регистрация: 29.09.2011
Сообщений: 16
29.10.2011, 21:02  [ТС]     При наличие корней в системе уравнений выдает ошибку #11
Цитата Сообщение от -=ЮрА=- Посмотреть сообщение
Larri, какой именно метод решения нужен?
Метод Ланцоша

Добавлено через 5 минут
Цитата Сообщение от Deviaphan Посмотреть сообщение
w[i]=0;
w[i]=( 0 +w1[i]*f)%2;

Я сперва прочитал, вместо w1 тоже w, думал, вообще в ноль всё выражение. А так у тебя просто избыточность. Если w не должно суммироваться, то можешь смело удалить цикл с обнулением.

Вызов WC некорректен: выделенная память не освобождается. При каждом вызове WC формируется неудаляемый массив, но используется только один элемент... В общем, исправляй логику.
Я немного ошиблась WC - произведение вектора на число. Используется весь массив, созданный в WC. Но спасибо, что обратил мое внимание на это
-=ЮрА=-
29.10.2011, 21:06
  #12

Не по теме:

alkagolik, код ниже тебе, а ведь не меняю указатель, пытаюсь работать далее с массивом, мне сложно пока выложить код который даст сбой и с malloc - но думаю когда стек кончится то и с этим грабли вылезут

код
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
#include <stdio.h>
#include <stdlib.h>
 
int * f(int n)
{
    int * x = (int *)malloc(n*sizeof(int));
    for(int i = 0; i < n; i++)
        x[i] = i;
    return x;
}
 
int * f2(int * x, int m)
{
    x = (int *)realloc((void *)x,m);
    return x;
}
 
 
int main ()
{
    int i, * a = f(10);
    printf("f1 : \n");
    for(i = 0; i < 10; i++)
        printf("%d ",a[i]);
    a = f2(a, 5);
    printf("\nf2 : \n");
    for(i = 0; i < 5; i++)
        printf("%d ",a[i]);
    free(a);
    return 0;
}

-=ЮрА=-
Заблокирован
Автор FAQ
29.10.2011, 21:07     При наличие корней в системе уравнений выдает ошибку #13
Цитата Сообщение от Larri Посмотреть сообщение
Метод Ланцоша
- за свою практику не встречал, поищу где нибудь матоописание, потом отпишусь...
Larri
0 / 0 / 0
Регистрация: 29.09.2011
Сообщений: 16
29.10.2011, 21:11  [ТС]     При наличие корней в системе уравнений выдает ошибку #14
Цитата Сообщение от -=ЮрА=- Посмотреть сообщение
- за свою практику не встречал, поищу где нибудь матоописание, потом отпишусь...
Точнее сказать, у меня его модификация. Система уравнений над кольцом по модулю 2. Ещё раз перечитала ошибку, которую выводит... Я так думаю, что проблема именно с выделением памяти.
-=ЮрА=-
29.10.2011, 21:20
  #15

Не по теме:

Larri, а у вас нет ничего из теории???Это ускорило бы реализацию алгоритма
PS:Мне нужно матоописание не потому что не хочу править Ваш код, а потому что в силу как я понял сложности метода возможно по самому алгоритму что-то упустили, вот Я не хочу обидеть(мол не так написали) просто хочу 100% как логическую так и ошибку реализации исключить

alkagolik
 Аватар для alkagolik
1510 / 616 / 79
Регистрация: 15.07.2011
Сообщений: 3,552
29.10.2011, 21:25     При наличие корней в системе уравнений выдает ошибку #16
Сообщение было отмечено автором темы, экспертом или модератором как ответ
-=ЮрА=-, просто руками писать надо, а не чем-либо еще
тест твоей программы
deadman@kubenta:~/documents/projects/system/bin/Debug$ valgrind ./system
==2486== Memcheck, a memory error detector
==2486== Copyright (C) 2002-2010, and GNU GPL'd, by Julian Seward et al.
==2486== Using Valgrind-3.6.1-Debian and LibVEX; rerun with -h for copyright info
==2486== Command: ./system
==2486==
f1 :
0 1 2 3 4 5 6 7 8 9
f2 :
==2486== Invalid read of size 4
==2486== at 0x8048595: main (main.c:28)
==2486== Address 0x41c5084 is 4 bytes inside a block of size 5 alloc'd
==2486== at 0x402896C: realloc (vg_replace_malloc.c:525)
==2486== by 0x80484FA: f2 (main.c:14)
==2486== by 0x804856F: main (main.c:25)
==2486==
0 1 0 0 0 ==2486==
==2486== HEAP SUMMARY:
==2486== in use at exit: 0 bytes in 0 blocks
==2486== total heap usage: 2 allocs, 2 frees, 45 bytes allocated
==2486==
==2486== All heap blocks were freed -- no leaks are possible
==2486==
==2486== For counts of detected and suppressed errors, rerun with: -v
==2486== ERROR SUMMARY: 4 errors from 1 contexts (suppressed: 11 from 6)

нормальный код
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
#include <stdio.h>
#include <stdlib.h>
 
int * f(int n)
{
    int * x = (int *)malloc(n*sizeof(int));
        for(int i = 0; i < n; i++)
                x[i] = i;
    return x;
}
 
int * f2(int * x, int m)
{
    int *p = x;
    p = realloc(p, m * sizeof(int));
    return p;
}
 
 
int main ()
{
    int i, * a = f(10);
        printf("f1 : \n");
    for(i = 0; i < 10; i++)
                printf("%d ",a[i]);
        a = f2(a, 5);
        printf("\nf2 : \n");
        for(i = 0; i < 5; i++)
                printf("%d ",a[i]);
    free(a);
    return 0;
}
тест моей программы

deadman@kubenta:~/documents/projects/system/bin/Debug$ valgrind ./system
==2623== Memcheck, a memory error detector
==2623== Copyright (C) 2002-2010, and GNU GPL'd, by Julian Seward et al.
==2623== Using Valgrind-3.6.1-Debian and LibVEX; rerun with -h for copyright info
==2623== Command: ./system
==2623==
f1 :
0 1 2 3 4 5 6 7 8 9
f2 :
0 1 2 3 4 ==2623==
==2623== HEAP SUMMARY:
==2623== in use at exit: 0 bytes in 0 blocks
==2623== total heap usage: 2 allocs, 2 frees, 60 bytes allocated
==2623==
==2623== All heap blocks were freed -- no leaks are possible
==2623==
==2623== For counts of detected and suppressed errors, rerun with: -v
==2623== ERROR SUMMARY: 0 errors from 0 contexts (suppressed: 11 from 6)
Larri
0 / 0 / 0
Регистрация: 29.09.2011
Сообщений: 16
29.10.2011, 21:36  [ТС]     При наличие корней в системе уравнений выдает ошибку #17
Цитата Сообщение от -=ЮрА=- Посмотреть сообщение

Не по теме:

Larri, а у вас нет ничего из теории???Это ускорило бы реализацию алгоритма
PS:Мне нужно матоописание не потому что не хочу править Ваш код, а потому что в силу как я понял сложности метода возможно по самому алгоритму что-то упустили, вот Я не хочу обидеть(мол не так написали) просто хочу 100% как логическую так и ошибку реализации исключить

Не по теме:

Вот ссылка на сам алгоритм, но у меня модификация (по ней я вообще ничего не нашла):http://lib.mexmat.ru/books/4633
Вероятность логической ошибки мала, т к 1. программа вычисляет решение правильно, 2. выдаваемая ошибка идет с сообщением Heap Corruption Detected.
А ошибка реализации только если где-то при работе с памятью...

-=ЮрА=-
Заблокирован
Автор FAQ
31.10.2011, 11:55     При наличие корней в системе уравнений выдает ошибку #18
Larri, для начала решил оттестировать ваш алгоритм на счёт простой системы, как видите на скрине в алгоритме скорее всего логическая ошибка, т.к. всегда попадается b1[i] != b[i]
Цитата Сообщение от Larri Посмотреть сообщение
if(b1[i]!=b[i]) j++;
Ваш код с небоьшими изменениями, в плане его сокращения откомпилировал несколько раз в Release, даже на простой системе Ваш алгоритм не работает. Вообщем у меня многое под сомнением, но вот здесь особенно
Цитата Сообщение от Larri Посмотреть сообщение
w1[j]=(w1[j]+a[j*n+i]*w[i])%2;
a[j*n+i] - допустим j = n - 1 i = n тогда вы пытаетесь обратиться к несуществующему элементу
a[(n - 1)*n + (n - 1)] как я понял по программе в а всего n элементов.
Цитата Сообщение от Larri Посмотреть сообщение
http://lib.mexmat.ru/books/4633
Там книжку читать надо, а быстрый поиск по Ланцош ничего не дал, если можно сделайте принт скрин нужных страниц и сюда тогда смогу проверить что к чему

Цитата Сообщение от Larri Посмотреть сообщение
w=(w+w1[i]*w2[i])%2;
- ну и это мне не понятно, вы двоичное сложение хотите сделать или что этим хотели реализовать???Это не двоичное сложение - а остаток от целочисленного деления на два, поясню
2 + 3 = 5 (10 + 11 = 101)
5 % 2 = 1

Код который использовал
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
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
#include <stdio.h>
#include <stdlib.h>
#include <time.h>
#define N 4
 
typedef struct
{
    int*a;
    bool l;//нулевой ли вектор
}
lama;
 
int* WC(int *w1,int f,int n)//alfa na vector
{
    int i=0;
    int *w=(int*)malloc(n*sizeof(int));
    for(i=0;i<n;i++)
        w[i]=0;
    for(i=0;i<n;i++)
        w[i]=(w[i]+w1[i]*f)%2;
    return w;
}
 
int WW(int *w1,int *w2,int n)//skal proizv
{
    int i=0;
    int w=0;
    for(i=0;i<n;i++)
        w=(w+w1[i]*w2[i])%2;
    return w;
}
 
int* MW(int n, int *a, int *w)
{
    int i,j;
    int *w1=(int*)malloc(n*sizeof(int));
    for (j=0;j<n;j++)
        w1[j]=0;
    for (j=0;j<n;j++)
    for (i=0;i<n;i++)
        w1[j]=(w1[j]+a[j*n+i]*w[i])%2;
    return w1;
}
 
lama form(int n,int *a,int *w1, int *w0)
{
        int i=0;
        int *w2;
        lama w;
        w.l=false;      
        int q,q1,e=0;
        w.a=(int*)malloc(n*sizeof(int));
        w2=(int*)malloc(n*sizeof(int));
        w2=MW(n,a,w1);
        q=(WW(MW(n,a,w2),w1,n))*(WW(w2,w1,n));  
        if(w0==NULL)
        {
                for(i=0;i<n;i++)
                {
                        w.a[i]=(w2[i]+WC(w1,q,n)[i])%2; //в двоичном сложении 1+1=1-1=0, если не учитывать перенос разряда
                        e+=w.a[i];
                }
        }
        else
        {
                q1=(WW(MW(n,a,w2),w0,n))*(WW(MW(n,a,w0),w0,n));
                for(i=0;i<n;i++)
                {
                        w.a[i]=(w2[i]+WC(w1,q,n)[i]+WC(w0,q1,n)[i])%2;
                        e+=w.a[i];                      
                }
        }       
        if(e==0) 
            w.l=true;
        free(w2);
        return w;
}
int*XX(int*a,int*w,int k,int n)//nepr
{
        int j,i,q;
        int*x =(int*)malloc(n*sizeof(int));  
        for(i=0;i<n;i++)
            x[i]=0;
        for(i=0;i<=k;i++)
        {               
                q=WW((w+i*n),w,n)*WW(MW(n,a,(w+i*n)),(w+i*n),n);
                for(j=0;j<n;j++)
                    x[j]=(x[j]+q*w[i*n+j])%2;
        }       
        return x;
}
bool prov(int*a,int*x,int*b,int n)
{       
    int i,j = 0;
    int * b1=MW(n,a,x);
    for(i=0;i<n;i++)
    if(b1[i]!=b[i])
        j++;
    return (j == 0 ? true : false);
}
int main()
{
        int *x,i,j,k=0;
        lama p; 
        int * a=(int*)malloc(N*N*sizeof(int));
        int * b=(int*)malloc(N*sizeof(int));
        int * w=(int*)malloc(N*N*sizeof(int));
        for(i=0;i<N;i++)//-пользовательское задание матрицы А
                for(j=0;j<N;j++)
                {
                        printf("Enter a[%i,%i] : ",i,j);
                        scanf("%i",&a[i*N+j]);
                }
        for(i=0;i<N;i++)
        {
                printf("Enter b[%i] : ",i);
                scanf("%i",&b[i]);
        }
        for (i=0; i<N; i++)
        {
                for (j=0; j<N; j++)
                        printf("\t %i", a[i*N+j]);
                printf ("\n");
        }
        for(i=0;i<N;i++)        
                printf("b[%i]=%i\n",i,b[i]);
        for(j=0;j<N;j++)
                w[j]=b[j];
        i=1;
        do
        {
                if(i==1) 
                    p=form(N,a,b,NULL);
                else 
                    p=form(N,a,(w+(i-1)*N),(w+(i-2)*N));
                for(j=0;j<N;j++)//было j<=N
                {
                        w[i*N+j]=p.a[j];//запись вектора W[i] в матрицу омег                    
                }
                k++;
                i++;
        }
        while(!p.l);
        if(p.l)
        {
                for (i=0; i<=k; i++)
                {
                        for (j=0; j<N; j++)
                                printf("\t %i", w[i*N+j]);
                        printf ("\n");
                }
                x=XX(a,w,k,N); 
                prov(a,x,b,N);
               if(!prov(a,x,b,N)) 
                    printf("answer does not exist\n");
                else 
                {
                        printf("answer is true\n");
                        printf("answer: ");
                        for(i=0;i<N;i++)
                                printf("%i ",x[i]);
                        printf("\n");
                }
        }
        else 
            printf("algorithm does not work\n");
 /*       free(a);
        free(b);
        free(w);*/
        return 0;
}
Миниатюры
При наличие корней в системе уравнений выдает ошибку  
Larri
0 / 0 / 0
Регистрация: 29.09.2011
Сообщений: 16
31.10.2011, 20:10  [ТС]     При наличие корней в системе уравнений выдает ошибку #19
Дело в том, что система на писана для системы над кольцом вычетов по модулю 2. Соответственно, там есть только числа 0 и 1, поэтому приведенный тест, естественно, приведет к неправильному результату. Я проверяла для систем над нужным концом, результат либо верный, либо программа падает. В а n*n элементов, т.к. это квадратная матрица
MoreAnswers
Эксперт
37091 / 29110 / 5898
Регистрация: 17.06.2006
Сообщений: 43,301
01.11.2011, 10:22     При наличие корней в системе уравнений выдает ошибку
Еще ссылки по теме:

C++ Программа запускается, но при компиляции выдает ошибку
Выдает ошибку при создании модуля C++
C++ Выдает ошибку при построении

Искать еще темы с ответами

Или воспользуйтесь поиском по форуму:
-=ЮрА=-
Заблокирован
Автор FAQ
01.11.2011, 10:22     При наличие корней в системе уравнений выдает ошибку #20
Цитата Сообщение от Larri Посмотреть сообщение
Дело в том, что система на писана для системы над кольцом вычетов по модулю 2. Соответственно, там есть только числа 0 и 1, поэтому приведенный тест, естественно, приведет к неправильному результату.
- т.е работаете с Системой н-го порядка с 1-цами и нулями, а как же тогда вектор свободных членов, что тоже только 1 и 0 в нём???


Цитата Сообщение от Larri Посмотреть сообщение
Я проверяла для систем над нужным концом, результат либо верный, либо программа падает.
- предлагаю заменить все массивы статическими, вы же для системы 4-го порядка алгоритм затачиваете
Yandex
Объявления
01.11.2011, 10:22     При наличие корней в системе уравнений выдает ошибку
Ответ Создать тему
Опции темы

Текущее время: 04:59. Часовой пояс GMT +3.
КиберФорум - форум программистов, компьютерный форум, программирование
Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2016, vBulletin Solutions, Inc.
Рейтинг@Mail.ru