Форум программистов, компьютерный форум, киберфорум
Наши страницы
C для начинающих
Войти
Регистрация
Восстановить пароль
 
oksanaBM
1 / 1 / 0
Регистрация: 27.01.2011
Сообщений: 91
1

Разложение LDLt

01.12.2012, 21:19. Просмотров 1121. Ответов 0
Метки нет (Все метки)

Имеется программа для разложения матрицы в виде LDLT в ленточном формате. Что нужно исправить в формировании D, чтобы диагональ была равна http://www.cyberforum.ru/cgi-bin/latex.cgi?{d}_{ii}=\pm 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
#include <stdio.h>
#include <conio.h>
 
typedef float type;
 
#define FMT_STR " %f"
#define FMT_STR2 " %f\n"
 
int n, ndiag, nelem;
type **al, *di, *vector;
 
 
void ReadData ()
{
    FILE *F;
    F = fopen ("n8.txt", "r");
    fscanf(F, "%i", &n);
    fscanf(F, "%i", &ndiag);
    nelem=n*ndiag;
    di = new type [n];
    for (int i=0; i<n; i++)
        fscanf(F,FMT_STR,&di[i]);
    
    al = new type* [n];
    for (int i=0; i<n; i++) 
    {
        al[i] = new type[ndiag];
        for (int j=0; j<ndiag; j++)
            fscanf (F, FMT_STR, &al[i][j]);
    }
    vector = new type [n];
    for (int i=0; i<n; i++)
        fscanf (F, FMT_STR, &vector[i]);
    fclose(F);
}
 
void CalcLDLt()
{
    FILE *fl;
    fl = fopen ("matrix2.txt", "w");
    for (int i=0; i<n; i++)
    {
        int j0 = i - ndiag;
        int j = j0;
        double sd = 0;
        for (int jl=0; jl<ndiag; jl++,j++)
        {
            double sl = 0;
            if(j>=0)
            {
                int kk = j0;
                int jlj = i-j;
 
                for (int jlk=0; jlk<jl; jlk++,jlj++,kk++) 
                {
                        sl += al[i][jlk]*al[j][jlj]*di[kk];
                }
                al[i][jl] = (al[i][jl]-sl)/di[j];
                sd += al[i][jl]*al[i][jl]*di[j];
            }
        }
        di[i]=di[i]-sd;
    }
    fprintf(fl, "al:\n");
    for (int i=0; i<n; i++) 
    {
        for (int j=0; j<ndiag; j++)
            fprintf (fl, FMT_STR2, al[i][j]);
    }
    fprintf(fl, "\ndi:\n");
    for(int i=0; i<n; i++)
        fprintf(fl, FMT_STR2, di[i]);
}
 
void CalcYZX()
{
    type *x = vector;
    type *y = vector;
    type *z = vector;
    type *f = vector;
    FILE *fl;
    fl = fopen ("matrix1.txt", "w");
 
    for (int i=0; i<n; i++)
    {
        double s=0;
        int j = i - ndiag;
        for (int jl=0; jl<ndiag; jl++,j++) 
        {
            if(j>=0)
                s += al[i][jl]*y[j];
        }
        y[i] = f[i]-s;
    }
    
        fprintf(fl, FMT_STR, z[i]);
    fprintf(fl, "\n");*/
 
    
    for (int i=n-1; i>=0; i--)
    {
        int j = i-1;
        type xi = z[i] / 1;
        for (int jl=ndiag-1; jl>=0; jl--,j--)
        {
            if(j>=0)
                z[j] -= al[i][jl]*xi;
        }
        x[i] = xi;
 
    }
    for(int i=0; i<n; i++)
        fprintf(fl, FMT_STR2, x[i]);
    fclose(fl);
    
    FILE *fli;
    fli = fopen ("output.txt", "w");
    for(int i=1; i<=n; i++)
        fprintf(fli, FMT_STR2, i-x[i-1]);
    fclose(fli);
}
 
void main ()
{
    ReadData();
    CalcLDLt();
    CalcYZX();
}
0
Надоела реклама? Зарегистрируйтесь и она исчезнет полностью.
Similar
Эксперт
41792 / 34177 / 6122
Регистрация: 12.04.2006
Сообщений: 57,940
01.12.2012, 21:19
Ответы с готовыми решениями:

Разложение числа на слагаемые
Не могу найти найти не рекурсивный алгоритм разложения числа на заданное кол-во...

Разложение числа на множители
Описать функцию Factors(A,N,F), находящую разложение натурального числа А на...

Разложение числа на слагаемые
Здравствуйте, товарищи форумчане! Нужна помощь в доработке программы. ...

Разложение матрицы по строке
Доброго времени суток! Недавно вот делал задачу и понял,что запутался,нужна...

Разложение в ряд Тейлора
Пожалуйста, помогите исправить программу. Условие:...

0
MoreAnswers
Эксперт
37091 / 29110 / 5898
Регистрация: 17.06.2006
Сообщений: 43,301
01.12.2012, 21:19

Разложение на простые множители
Разложение на простые множители

Разложение на простые множители
Требуется вывести представление целого числа N в виде произведения простых...

Разложение функции в ряд фурье
Здравствуйте! Подскажите пожалуйста как можно разложить в ряд Фурье функцию...


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

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

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