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

Куда уходят поезда ? Вывод значений матрицы - C++

Восстановить пароль Регистрация
 
Qazan
211 / 59 / 9
Регистрация: 30.04.2013
Сообщений: 778
Записей в блоге: 10
26.04.2014, 16:11     Куда уходят поезда ? Вывод значений матрицы #1
C++
1
2
    
                std::cout << "  [ " << (A[i][j]) << "  ]       ";
вот такой кусочек кода для вещественной матрицы

выводит единицу на некотором http://www.cyberforum.ru/cgi-bin/latex.cgi?(i_{0},j_{0})
привожу к int выводит ноль , Wtf(What the fenomen) ?
Similar
Эксперт
41792 / 34177 / 6122
Регистрация: 12.04.2006
Сообщений: 57,940
26.04.2014, 16:11     Куда уходят поезда ? Вывод значений матрицы
Посмотрите здесь:

C++ вывод значений Enum
умножить элемент каждой строки первой матрицы на наибольшее из значений элементов соответствующей строки второй матрицы C++
C++ Вывод значений на экран
Получить новую матрицу умножением элементов первой матрицы на наибольшее из значений элементов соответствующий строки второй матрицы C++
C++ Вычислить значение заданной функци на данном интервале с заданным шагом. Осуществить вывод значений аргумента и результатов вычислений значений функ
После регистрации реклама в сообщениях будет скрыта и будут доступны все возможности форума.
Instigator
 Аватар для Instigator
2 / 2 / 1
Регистрация: 26.04.2014
Сообщений: 12
26.04.2014, 16:18     Куда уходят поезда ? Вывод значений матрицы #2
По-видимому, у тебя где-то есть деление, результат которое округляется до нуля, так как пользуешься целочисленным значением. Будет круто, если ты покажешь остальной код
Qazan
211 / 59 / 9
Регистрация: 30.04.2013
Сообщений: 778
Записей в блоге: 10
26.04.2014, 16:25  [ТС]     Куда уходят поезда ? Вывод значений матрицы #3
Instigator, Вот


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
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
// Numerical.cpp: определяет точку входа для консольного приложения.
//
 
#include "stdafx.h"
 
#include <iostream>
#include <fstream>
 
using namespace std;
 
 
void Create(double ** &A,int n);
void Destroy(double ** &A,int n);
 
void Show(double ** A,int n);
 
void SwapString(double ** &A,int i,int k);
void SwapColumn(double ** &A,int n,int k,int l);
 
double ** &E(double ** &A,int n); // make Identity matrix
 
void JS(double ** A,double ** &A_inv,int n)
{
    double max = A[0][0];
    int imax = 0,jmax = 0;
 
    for(int i = 0;i < n;i++)
        for(int j = 0;j < n;j++)
            if(A[i][j] > max)
            {
                imax = i;
                jmax = j;
                max = A[i][j];
            }
 
//  SwapString(A,imax,0);
//  SwapColumn(A,n,jmax,0);
 
    for(int i = 0;i < n;i++)
    {
 
        if(A[i][i] != 0)
        {
            double t = 1./A[i][i];
            for(int j = i;j < n;j++)
            {
                A[i][j] = A[i][j]*t;
            }
 
            for(int j = 0;j < n;j++)
            {
                A_inv[i][j] = A_inv[i][j]*t;
            }
 
//          system("pause");
            for(int k = 0;k < i;k++)
            {
                double t = A[k][i];
                for(int l = i;l < n;l++)
                {
                    A[k][l] = A[k][l] - A[i][l]*t;
                }
 
                for(int l = 0;l < n;l++)
                {
                    A_inv[k][l] = A_inv[k][l] - A_inv[i][l]*t;
                }
            }
 
            for(int k = i + 1;k < n;k++)
            {
                double t = A[k][i];
                for(int l = i;l < n;l++)
                {
                    A[k][l] = A[k][l] - A[i][l]*t;
                }
 
                for(int l = 0;l < n;l++)
                {
                    A_inv[k][l] = A_inv[k][l] - A_inv[i][l]*t;
                }
            }
        }
    }
    
//  SwapString(A,imax,0);
//  SwapColumn(A,n,jmax,0);
    
//  SwapColumn(A_inv,n,imax,0);
//  SwapString(A_inv,jmax,0);
        
}
 
double ** &Fill(double ** &A,int n) // make Identity matrix
{
    for(int i = 0;i < n;i++)
    {
        for(int j = 0;j < n;j++)
            {
                A[i][j] = rand()%32;// = rand()%20 * (1.0/2.0)
            }
//      std::cout << std::endl;
    }
    return A;
}
 
void MultL(double ** A,double ** B,double ** &C,int n)
{
    for(int i = 0;i < n;i++)
        for(int j = 0;j < n;j++)
        {
            double s = 0;
            for(int k = 0;k < n;k++)
            {
                s = s + A[i][k] * B[k][j];
            }
            C[i][j] = s;
        }
    
}
 
void Read(double ** &A,int n,char * fname)
{
    std::ifstream fin(fname);
    if(fin)
    {
        for(int i = 0;i < n;i++)
            for(int j = 0;j < n;j++)
            {
                int t;
                fin >> t;
                A[i][j] = (double)t;
            }
        fin.close();
    }
}
 
int main()
{
    const int n = 4;
    double ** p, **p_inv,**res; 
    Create(p,n);
    Create(res,n);
 
    Create(p_inv,n);
        E(p_inv,n);
 
/*  p[0][0] = 0;
    p[0][1] = -2;
    p[1][0] = 0.5;
    p[1][1] = 0;*/
//  Fill(p,n);
 
        Read(p,n,"input");
 
        Show(p,n);
 
        JS(p,p_inv,n);
//          E(p,n);
//      JS(p_inv,p,n);
 
        cout << endl;
 
        Read(p,n,"input");
 
        MultL(p,p_inv,res,n);
        Show(res,n);
        cout << endl;
        Show(p_inv,n);
 
    Destroy(p,n);
    system("pause");
    return 0;
}
 
 
void Create(double ** &A,int n)
{
    A = new double*[n];
    for(int i = 0;i < n;i++)
        A[i] = new double[n];
}
 
void Destroy(double ** &A,int n)
{
    for(int i = 0;i < n;i++)
        delete [] A[i];
    delete [] A;
}
 
double ** &E(double ** &A,int n)
{
    for(int i = 0;i < n;i++)
        for(int j = i + 1;j < n;j++)
            A[i][j] = 0;
 
    for(int i = 0;i < n;i++)
        A[i][i] = 1;
 
    for(int i = 0;i < n;i++)
        for(int j = 0;j < i;j++)
            A[i][j] = 0;
    return A;
}
 
void SwapString(double ** &A,int i,int k)
{
    if(i == k) return;
    double * p = A[i];
        A[i] = A[k];
        A[k] = p;
}
 
void SwapColumn(double ** &A,int n,int k,int l)
{
    if(k == l) return;
 
    for(int i = 0;i < n;i++)
    {
        double t = A[i][k];
        A[i][k] = A[i][l];
        A[i][l] = t;
    }
}
 
void Show(double ** A,int n)
{
    for(int i = 0;i < n;i++)
    {
        for(int j = 0;j < n;j++)
            {
                std::cout.width(10);
                std::cout << "  [ " << (A[i][j]) << "  ]       ";
            }
        std::cout << std::endl;
    }
}
Instigator
 Аватар для Instigator
2 / 2 / 1
Регистрация: 26.04.2014
Сообщений: 12
26.04.2014, 16:54     Куда уходят поезда ? Вывод значений матрицы #4
Как я и предполагал. Если значение A[i][i] целое, то при делении на неё единицы получается ноль.

Цитата Сообщение от Qazan Посмотреть сообщение
C++
1
2
3
4
5
double t = 1./A[i][i];
for(int j = i;j < n;j++)
{
    A[i][j] = A[i][j]*t;
}
Qazan
211 / 59 / 9
Регистрация: 30.04.2013
Сообщений: 778
Записей в блоге: 10
26.04.2014, 17:07  [ТС]     Куда уходят поезда ? Вывод значений матрицы #5
Instigator, Тобишь целое на еденицу это ноль ??
Instigator
 Аватар для Instigator
2 / 2 / 1
Регистрация: 26.04.2014
Сообщений: 12
26.04.2014, 17:19     Куда уходят поезда ? Вывод значений матрицы #6
Наоборот, единица на целое — это ноль
Qazan
211 / 59 / 9
Регистрация: 30.04.2013
Сообщений: 778
Записей в блоге: 10
26.04.2014, 17:25  [ТС]     Куда уходят поезда ? Вывод значений матрицы #7
Instigator, т е еденица на еденицу ноль ?
как с таким бороться ?

Добавлено через 15 секунд
или это уже погрешность метода ?
Instigator
 Аватар для Instigator
2 / 2 / 1
Регистрация: 26.04.2014
Сообщений: 12
26.04.2014, 17:41     Куда уходят поезда ? Вывод значений матрицы #8
Qazan, нет, если единицу поделишь на единицу, это будет единственный случай, когда он выдаст не ноль. Но я сомневаюсь, что у тебя единица в знаменателе. Обычный запуск программы показал, что там не всё так корректно работает, например — не присваиваются самые первые значения. Глубже вникать пока не могу
Qazan
211 / 59 / 9
Регистрация: 30.04.2013
Сообщений: 778
Записей в блоге: 10
27.04.2014, 21:08  [ТС]     Куда уходят поезда ? Вывод значений матрицы #9
Instigator, Простите , где не присваиваются значения ?
я что то использую без предварительной инициализаций ?

Добавлено через 3 минуты
Instigator, Возможно вы об этом
C++
1
   Read(p,n,"input");
В папке исполнения программы должден еще лежать файл с перечисленными целыми значениями (Это так для теста ) в количестве n*n штук (достаточным для заполненния квадратной матрицы)
MoreAnswers
Эксперт
37091 / 29110 / 5898
Регистрация: 17.06.2006
Сообщений: 43,301
27.04.2014, 22:30     Куда уходят поезда ? Вывод значений матрицы
Еще ссылки по теме:

Матрицы. Отрицательные элементы матрицы, расположенные выше главной диагонали, заменить на квадраты их значений C++
C++ Написать и протестировать подпрограмму, которая переставляет в обратном порядке элементы главной диагонали квадратной матрицы - вывод всей матрицы
C++ Задана матрица Т(N,M). Сформировать массив R(N,M-2),куда переписать все элементы матрицы T

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

Или воспользуйтесь поиском по форуму:
Instigator
 Аватар для Instigator
2 / 2 / 1
Регистрация: 26.04.2014
Сообщений: 12
27.04.2014, 22:30     Куда уходят поезда ? Вывод значений матрицы #10
Qazan
Хм, действительно. К сожалению, не могу сейчас всё проверить, но как вариант могу предложить пройти пошаговым дебагом, особенно те места, где задаются значения матрицы A и функция JS, содержащаю в себе строчку
C++
1
double t = 1./A[i][i];
Yandex
Объявления
27.04.2014, 22:30     Куда уходят поезда ? Вывод значений матрицы
Ответ Создать тему
Опции темы

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