Непонятная ошибка сегментации - C (СИ) - Ответ 1070950

@dasalam · Регистрация: 14.12.2008

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

Собственно код

C

#include <stdio.h>
#include <math.h>
#include <unistd.h>
#include <stdlib.h>
#include <fcntl.h>
#include <sys/types.h>
#include <sys/stat.h>
 
#define MAX_ITERR   100
#define N   3
#define EPS 0.01
 
double SCALAR(double *x, double *y)
{
    int i;
    double res=0.0;
    
    for(i=0; i<=N-1; i++)
        res = res + x[i]*y[i];
        
    return res;
}
 
void matrvec(double **A, double *U, double *AU)
{
    int i, j;
    double res;
    
    for(i=0; i<=N-1; i++) {
        res = 0.0;
        for(j=0; j<=N-1; j++) {
            res = res + A[i][j]*U[j];
            AU[i] = res;
        }
    }
}
 
void discrepancy(double **A, double *U, double *F, double *r)
{
    int i;
    double AUk[N];
    
    matrvec(A, U, AUk);
    for(i=0; i<=N-1; i++) {
        r[i] = AUk[i] - F[i];
    }
}
 
int fvread(double *F, int descriptor)
{
    int nread;
    int i;
    
    for(i=0; i<=N-1; i++) {
        nread=read(descriptor, &F[i], sizeof(double));
        if(nread == -1)
            return -1;
    }
    
    return 0;
}
 
int Amread(double **A, int descriptor)
{
    int nread;
    int i, j=0;
    
    for(i=0; i<=N-1; i++)
        for(j=0; j<=N-1; j++) {
            nread=read(descriptor, &A[i][j], sizeof(double));
            if(nread == -1)
                return -1;
        }
    
    return 0;
}
 
int main()
{
    double *F = (double *)malloc(N * sizeof(double));
    double t[MAX_ITERR];
    double Ark[MAX_ITERR], AUk[MAX_ITERR];
    int i=0, j;
    double err[N];
    double **A, **r, **U;
    
    A = (double **)malloc(N * sizeof(double *));
    if(A == NULL) return -1;
        for(i = 0; i < N; i++) {
            A[i] = (double *)malloc(N * sizeof(double));
            if (A[i] == NULL) return -1;
        }
            
    r = (double **)malloc(N * sizeof(double *));
    if(r == NULL) return -1;
        for(i = 0; i < N; i++) {
            r[i] = (double *)malloc(N * sizeof(double));
            if (r[i] == NULL) return -1;
        }
    
    U = (double **)malloc(N * sizeof(double *));
    if(U == NULL) return -1;
        for(i = 0; i < N; i++) {
            U[i] = (double *)malloc(N * sizeof(double));
            if (U[i] == NULL) return -1;
        }
    
    int d = open("data.bin", O_RDONLY);
    printf("d=%d\n", d);
    
    if(fvread(F, d) == -1) return -1;
    for(i=0; i<=N-1; i++)
        printf("F[%d] = %.2f; ", i, F[i]);
    
    if(Amread(A, d) == -1) return -1;
    for(i=0; i<=N-1; i++) {
        printf("\n");
        printf("i = %d\n", i);;
        for(j=0; j<=N-1; j++)
            printf("A[%d][%d] = %.2f; ", i, j, A[i][j]);
    }
 
    for(i=0; i<=N-1; i++)
        U[0][i] = 1;
        
    matrvec(A, U[0], AUk);
    
    do {
        printf("U[%d] = (", i);
        for(j=0; j<=N-1; j++)
            printf("%.4f, ", U[i][j]);
        printf(")^T\n");
        
        discrepancy(A, U[i], F, r[i]);
        matrvec(A, r[i], Ark);
        t[i] = SCALAR(r[i], r[i]) / SCALAR(Ark, r[i]);
        matrvec(A, U[i], AUk);
        for(j=0; j<=N-1; j++)
            U[i+1][j] = U[i][j] - t[i]*(AUk[j]-F[j]);
        for(j=0; j<=N-1; j++)
            err[j] = U[i+1][j]-U[i][j];
        i++;
        printf("||U_i+1-U_i|| = %f\n", SCALAR(err, err));
    } while((i!=MAX_ITERR-2) || sqrt(SCALAR(err, err))>EPS);
    
    return 0;
}

Вывод программы:

Bash

d=3
F[0] = 1.00; F[1] = 1.00; F[2] = 0.00;
 
i = 0
A[0][0] = 1.00; A[0][1] = 0.00; A[0][2] = 0.00; 
i = 1
A[1][0] = 0.00; A[1][1] = 1.00; A[1][2] = 0.00; 
i = 2

Был бы очень благодарен за скорую помощь.

Вернуться к обсуждению:
Непонятная ошибка сегментации C (СИ)

Следующий ответ

Другие темы раздела
C (СИ) Обьеденить программу с указателями и без в СИ Здравствуйте! Нужна Ваша помощь! Проблема такая: существует две одинаковые программы, одна написана без указателей, а вторая такая же только написана при помощи указателей. Программы выполняют одинаковые действия! Нужно эти две программы обьеденить в одну, чтобы один раз вводились данные, но выводилось сразу два одинаковых результата! 1-ый результат программы без указателей, а второй программы с... https://www.cyberforum.ru/ c-beginners/ thread183984.html	C (СИ) Ошибка в сортировке ПОМОГИТЕ,ПОЖАЛУЙСТА,НАЙТИ ОШИБКУ: Нужно: отсортировать по алфавиту список студенток набравших 10 баллов. #include<stdio.h> #include <string.h> #define kstud 8 typedef struct { char name ; char gender; int sumball;
C (СИ) Преобразование отдельного символа строки в число https://www.cyberforum.ru/ c-beginners/ thread183647.html Доброго времени суток!!! Задана строка, содержащая целые числа и слова, разделенные пробелами (одним или несколькими). Определить сумму целых чисел в строке. Т. е. вводим строку "fhgfh123" получаем сумму 6 #include <stdio.h> int main(int argc, char *argv) {	C (СИ) Организовать ввод пароля звездочками Здравствуйте. Пишу задание по программированию ("C", консоль), возникла необходимость и острое желание реализовать пункт меню с "админскими" правами т.е. под паролем. Подскажите пожалуйста можно ли реализовать cin.getline таким образом чтобы вместо вводимого текста на экране отображались символы? например: " ******** " вместо "65983210" или в крайнем случае не отображалось вообще... https://www.cyberforum.ru/ c-beginners/ thread183588.html
C (СИ) Перевод арабских чисел в римские Программа преобразовывает натуральные арабские числа в римские: Не пойму каким образом функция fprint_rim это реализует? #include <stdio.h> #include <stdlib.h> #include <string.h> #define MAX_RIM 3999 void fprint_rim(FILE , char , int);	C (СИ) Неправильно считает интеграл https://www.cyberforum.ru/ c-beginners/ thread183547.html написал два метода Симпсона и Монте-Карло, но 2-ой не правильно считает, помогите пожалуйстаю Значение у Симпсона получается: 1.3352, у Монте-Карло: 0.05338. Вот метод Монте-Карло: double CalcMonte(int n)// метод Монте-Карло { double part = 0.0; double part1 = 3.14; double integral = 0.0; for (int i = 0; i <= n; i++) { double x = (3.14 / double(n)) * rand();
C (СИ) Структура - найти ошибки в коде вот написал код программы,но выходят ошибки,кто нибудь подскажите пожалуйста что нужно исправить #include <stdio.h> int n; struct flower { int cena; int id; https://www.cyberforum.ru/ c-beginners/ thread183440.html	Сортировка односвязного списка со строками C (СИ) Весь исходный код программы, осталось последняя проблема, отсортировать список по длине и алфавиту, функция ContainerSort варианты которые были откоментил, буду благодарен если укажите на ошибки в них, на данный момент остановился на проблеме: сортировка затирает элемент в списке и отображает либо первое значение либо второе, откуда она их берет ума не приложу, т.к. менял всё что только...
C (СИ) Написать программу. Может ли быть вводимое число быть восьмеричным. Собственно, использовать можно , лишь, цикл while. Массивы и функции преобразования строки нельзя. помогите найти ошибку #include <stdio.h> #include <conio.h> int main (void) { int n=0,i=0; unsigned char c=57; https://www.cyberforum.ru/ c-beginners/ thread183276.html	C (СИ) Получить новую матрицу из элементов исходной https://www.cyberforum.ru/ c-beginners/ thread183237.html Даны действительные числа a_1,…,a_n действительная квадратная матрица порядка n (n≥6). Получить действительную матрицу n x (n+1) вставив в исходную матрицу между пятым и шестым столбцами новый столбец с элементами a_1,…,a_n #include<iostream.h> #include<string.h> #define N 6 void main(void) { float B, A; int i; printf("Введите массив из 6 чисел:\n", N); for (i=0;i<N;i++)
C (СИ) Прошу найти ошибку в коде. #include <stdio.h> #include <conio.h> void f1(float a,int n,float c, int t) {int i; float mmm=0; c=(float)malloc(nsizeof(float)); if (a==NULL) { printf("\n ERROR: Thou Art Out Of Memoreth, Good Sir! \n"); exit(1);	C (СИ) Даны действительные числа https://www.cyberforum.ru/ c-beginners/ thread183193.html Даны действительные числа a_1,…,a_n действительная квадратная матрица порядка n (n≥6). Получить действительную матрицу n x (n+1) вставив в исходную матрицу между пятым и шестым столбцами новый столбец с элементами a_1,…,a_n Добавлено через 50 секунд #include<iostream.h> #include<string.h> #define N 6 void main(void) { float B, A; int i; printf("Введите массив из 6 чисел:\n", N);

	30.10.2010, 19:39