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

С++ для начинающих

Войти
Регистрация
Восстановить пароль
 
Vyzov
6 / 6 / 1
Регистрация: 19.02.2013
Сообщений: 64
#1

Шаг LU разложения - C++

05.02.2014, 19:03. Просмотров 228. Ответов 0
Метки нет (Все метки)

Пишу программу для LU разложения, вот код одного шага:
C++
1
2
3
4
5
6
7
8
9
10
11
12
void stepLU (double **a,int i, int n)
{
    for (int j = i+1; j<n; j++)
        {
            double u = a[j][i]/a[i][i];//для проверки вычислений
            for(int k = 0; k < n;k++)
            {           
                double t = (a[j][k] - ((a[j][i]/a[i][i]) * a[i][k]));//для проверки вычислений
                a[j][k] = t;//a[j][k] - a[j][i]/a[i][i] * a[i][k];
            }
        }
}
вызывается из в мейне из цикла for(int i = 0; i < n-1; i++) где n - размерность матрицы

Проблема в том что код изменяет только первые елементы строки, т.е.
из матрицы
-1, -3, -4
-1, -3, -5
-2, -1, -5

получается
-1, -3, -4
0, -3, -5
0, -1, -5

вместо

-1, -3, -4
0, 0, -1
0, 5, 13

Пройдя по циклу на листке все получается, а вот при запуске на компе вот такая проблема(см.выше).
При пошаговом дебаге видно, что после 1-й элемент строки высчитывается правильно, а последующие просто присваиваются без изменений, сам пока не дошел почему так происходит, может кто раньше поймет...

З.Ы. Шаг LU разложения - тот же метод гауса

Добавлено через 14 минут
Вопрос закрыт все заработало после изменения кода таким образом:
C++
1
2
3
4
5
6
7
8
9
10
11
12
13
14
void stepLU (double **a,int i, int n)
{
    for (int j = i+1; j<n; j++)
        {
            double u = a[j][i]/a[i][i];
            for(int k = 0; k < n;k++)
            {
                double x = a[j][k];
                double y = a[i][k];
                double t = x - u*y;
                a[j][k] = t;
            }
        }
}
От таки дела =\
0
Надоела реклама? Зарегистрируйтесь и она исчезнет полностью.
Similar
Эксперт
41792 / 34177 / 6122
Регистрация: 12.04.2006
Сообщений: 57,940
05.02.2014, 19:03
Здравствуйте! Я подобрал для вас темы с ответами на вопрос Шаг LU разложения (C++):

Шаг компиляции, шаг компоновки, и шаг запуска - C++
Что происходит на шаге компиляции, шаге компоновки, и шаге запуска, с переменными и функциями. что происходит при встрече в коде...

Шаг младенца, шаг великана Трудоемкость - C++
Помогите реализовать данную функцию с трудоемкостью (√P×logP). Как делать именно поиск равных шагов, не сравнивая оба массива. ...

Шаг в 0,5 - C++
Вычислить значение функции y={x}^{5}-\sqrt{x} При 0 \leq X \leq 5 c шагом в 0,5 Т.е. x 0 0.5 ...

Изменить шаг - C++
#include &lt;iostream&gt; #include &lt;conio.h&gt; int main() { for (int i=100;i&lt;=200;i++) {std::cout&lt;&lt;i&lt;&lt;'\n'; ...

Шаг конем - C++
просьба . помогите . ну вот я должен написать код обход коня по всей доске ,программа находит первую точку коня и его первые возможные...

следующий шаг - C++
нужен совет заканчиваю изучать 1-ю книгу по С++ Дайтела Какие лучше книги дальше изучать, по: STL, MFC, OpenGL, Win 32 API, или что-то...

0
MoreAnswers
Эксперт
37091 / 29110 / 5898
Регистрация: 17.06.2006
Сообщений: 43,301
05.02.2014, 19:03
Привет! Вот еще темы с ответами:

Как изменить шаг на 2? - C++
#include &lt;iostream&gt; #include &lt;ctime&gt; int main() { clock_t t = clock(); { std::size_t size = 1000; ...

первый шаг он трудный самый! - C++
Как наверное понятно из заголовка, я решил изучить С++. Скачал книгу, компилятор Dev-C++. В книге написано создать свою первую программу,...

Герберт Шилдт, С++ Шаг за шагом - C++
Приветствую всех! Я чайник в С++, и заказал себе книгу Герберт Шилдт, С++ Шаг за шагом. Стоило того? Будет ли понятна эта книга, если я...

Вывести каждый шаг сортировки - C++
Имеется вот такой код шейкерной сортировки, необходимо, чтобы каждый шаг сортировки выводился на экран, все мои попытки не увенчались...


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

Или воспользуйтесь поиском по форуму:
1
Ответ Создать тему
Опции темы

КиберФорум - форум программистов, компьютерный форум, программирование
Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2017, vBulletin Solutions, Inc.
Рейтинг@Mail.ru