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

Решение СЛАУ большой размерности методом сопряженных градиентов - C++

Восстановить пароль Регистрация
Другие темы раздела
C++ разложить число на цифры и вернуть в массив http://www.cyberforum.ru/cpp-beginners/thread738290.html
//--------------------------------------------------------------------------- #include <iostream> #include <string> using namespace std; int *read(int x) { string s; char *c=new char; itoa(x,c,10); s.append(c);
C++ Иерархия типов помогите сделать программу, хотябы дайте пример. Задача. Создать иерархию типов - устройство ввода, устройство вывода, сенсорный экран. Реализовать во всех классах конструкторы по умолчанию, конструкторы копирования, виртуальные деструкторы, перегруженные функции ввода в поток и вывода из потока. Реализовать не менее 5 функций в каждом классе, 3 из них перегрузить в наследнике http://www.cyberforum.ru/cpp-beginners/thread738282.html
C++ Функции пользователя. Где ошибка?
Вот код программы: #include "stdafx.h" #include <iostream> #include <math.h> #include <conio.h> using namespace std; // Определение функции F float F(float x) {
C++ Двусвязный список
Есть код.Двусвязный список.Помогите допилить функцию findforday(); это удаление элемента .h #include"stdafx.h" #include <stdio.h> #include <iomanip> #include <conio.h> #include <string.h> #include <stdlib.h> #include <locale.h>
C++ Написать программу,которая заполняет линейный массив,размер которого задаются с клавиатуры,случайными натуральными числами в заданном диапазоне http://www.cyberforum.ru/cpp-beginners/thread738273.html
Написать программу,которая заполняет линейный массив,размер которого задаются с клавиатуры,случайными натуральными числами в заданном диапазоне,выводит полученный массив на экран и записывает в файл те элементы, которые являются простыми числами. Очень нужно!! #include "stdafx.h" #include <iostream> #include <iomanip> #include <conio.h> #include <stdlib.h> #include <fstream> using...
C++ Multiset Какой параметр надо задать контейнеру, чтобы он сортировал по убыванию? Или нужен какой-то алгоритм? подробнее

Показать сообщение отдельно
nullpointer
 Аватар для nullpointer
45 / 45 / 5
Регистрация: 30.03.2009
Сообщений: 518
19.12.2012, 17:28     Решение СЛАУ большой размерности методом сопряженных градиентов
Всем првиет! Возникла проблемка с методом сопряженных градиентов. Если задавать самому значения матрицы и правой части, то все решается отлично. Сделал rand'омное заполнение. В итоге получаю, что программа работает только с матрицами размера 2*2 и то не со всеми, даже 3*3 она уже не считает правильно, в результате выводит ответ в виде -1.#IND. А мне нужно, чтобы она работала с матрицами большого порядка, мне надо узнать время выполнения алгоритма. В чем можнт быть проблема? Я даже double заменил на long double и все равно не помогает.
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
#include<stdio.h> 
#include <time.h>
#include <conio.h>
#include <iostream>
using namespace std;
 
#define M 2 //размерность системы
 
#define E 0.1 //точность вычислений
 
 
long double Xk[M], Zk[M];
long double Rk[M], Sz[M], alpha, beta, mf;
long double Spr, Spr1, Spz;
long double A[M][M];//={2.0, -1.0, 0.0, -1.0, 5.0, 3.0, 0.0, 3.0, 1.0};
long double F[M] ;//= {1.0, 7.0, 4.0};
 
int main(void){
    int  i, j, kl = 1;
 
    srand(time(0));
    for (i = 0; i < M; i++)
    {
        for (j = 0; j < M; j++){
     A[i][j] = rand() % 2 - rand() % 1; // заполняем массив случайными значениями в диапазоне от -1 до 1 включительно
     cout << A[i][j] << " "; // печать элементов одномерного массива
        }
        cout << endl;
    }
 
    for (j = 0; j < M; j++){
     F[j] = rand() % 2 - rand() % 2; // заполняем массив случайными значениями в диапазоне от -1 до 1 включительно
     cout << F[j] << " "; // печать элементов одномерного массива 
        }
 
    clock_t time;
    time = clock();//начальное время
 
/* Вычисляем сумму квадратов элементов вектора F*/  
    for(mf = 0,i = 0; i < M; i++){
        mf += F[i] * F[i];
      }
 
 
/* Задаем начальное приближение корней. В Хk хранятся значения корней
 * к-й итерации. */
    for(i = 0; i < M; i++){
        Xk[i] = 0.2;
    } 
 
/* Задаем начальное значение r0 и z0. */
    for(i = 0; i < M; i++){
        for(Sz[i]=0,j = 0; j < M; j++)
            Sz[i] += A[i][j] * Xk[j];
      Rk[i] = F[i] - Sz[i];
      Zk[i] = Rk[i];  
      }
    
    do{ 
  /* Вычисляем числитель и знаменатель для коэффициента
   * alpha = (rk-1,rk-1)/(Azk-1,zk-1) */
        Spz = 0;
        Spr = 0;
        for(i = 0; i < M; i++){ 
            for(Sz[i] = 0, j = 0; j < M; j++){
                Sz[i] += A[i][j] * Zk[j];
            }
          Spz += Sz[i] * Zk[i];
          Spr += Rk[i] * Rk[i];
          }
        alpha = Spr / Spz;             /*  alpha    */
    
 
/* Вычисляем вектор решения: xk = xk-1+ alpha * zk-1, 
    вектор невязки: rk = rk-1 - alpha * A * zk-1 и числитель для betaa равный (rk,rk) */
        Spr1 = 0;
        for(i = 0; i < M; i++){
            Xk[i] += alpha * Zk[i];
            Rk[i] -= alpha * Sz[i];
            Spr1 += Rk[i] * Rk[i];
            cout << "Iter #" << kl;
            cout << " " << "X[" << i << "] = " << Xk[i] << endl;
        }
        cout << endl;
        kl++;
 
/* Вычисляем  beta  */
        beta = Spr1 / Spr; 
 
/* Вычисляем вектор спуска: zk = rk+ beta * zk-1 */
        for(i = 0; i < M; i++)
            Zk[i] = Rk[i] + beta * Zk[i];       
      }
/* Проверяем условие выхода из итерационного цикла  */
    while(Spr1 / mf > E * E);
 
    time = clock() - time;//время выполнения
    cout << "Time result = " << (double)time << "ms" << endl;
    system("pause");
    return 0;
  }
Добавлено через 59 минут
UPD.Понял почему не считает! Матрица же симметричная и положительно определенная должна быть!!! Теперь возникает другой вопрос: как рандомно заполнить симметричную матрицу?

Добавлено через 1 час 1 минуту
UPD.Все. Во всем разобрался
После регистрации реклама в сообщениях будет скрыта и будут доступны все возможности форума.
 
Текущее время: 03:15. Часовой пояс GMT +3.
КиберФорум - форум программистов, компьютерный форум, программирование
Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2016, vBulletin Solutions, Inc.
Рейтинг@Mail.ru