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

Умные указатели - C++

Войти
Регистрация
Восстановить пароль
Другие темы раздела
C++ Создать однонаправленный список абстрактных данных http://www.cyberforum.ru/cpp-beginners/thread1546455.html
Вы случаем не поможите еще с этой задачей буду счастлив. В каждом варианте необходимо, во-первых, реализовать определенный абстрактный тип данных, с которым можно работать только через...
C++ Определить меньший угол между часовой и минутной стрелками на циферблате часов Доброго времени суток. Прошу помочь с решением. C++ только начал изучать, поэтому решение здесь несложное. Дано время — два целых числа количество часов и минут. Необходимо определить меньший угол... http://www.cyberforum.ru/cpp-beginners/thread1546454.html
C++ Зачем пишется имя класса перед шаблоном?
#include <iostream> template <typename T> class Grid { public: Grid(int mWidth = kDefaultWidth, int mHeight = kDefaultHeight); ~Grid(); Grid(const Grid<T>& rhs); // ??? 1
C++ Ошибка (access violation) при обращении к элементу динамического массива структур
Добрый день! Все прекрасно работает если в структуре статические массивы. Но так уж получилось, что есть необходимость задавать размерность массивов вручную. Итак, есть некая структура с динамическим...
C++ Динамическое программирование http://www.cyberforum.ru/cpp-beginners/thread1546389.html
Столкнулся с такой задачей. Есть 6 фигурок площадью 3. Нужно узнать, сколькими способами можно полностью замостить ими поле n на m, если известно, что nm делится на 3. Я даже не представляю, как...
C++ Передача одного и того же аргумента 13 раз как сделать чтобы в цикле 13 раз передавался аргумент со значением 3 if(int j=0; j<52; j++) //Цикл 52 раза повторяется из них первые 13 аргументом должен быть 3 { funk(чтобы аргументом 13 раз... подробнее

Показать сообщение отдельно
ExModE
2 / 2 / 2
Регистрация: 04.03.2011
Сообщений: 27

Умные указатели - C++

07.10.2015, 11:12. Просмотров 240. Ответов 5
Метки (Все метки)

Кликните здесь для просмотра всего текста
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
#include <fstream>
#include <algorithm>
#include <vector>
#include <iostream>
#include <memory>
using namespace std;
 
template<typename real, typename ext = real>
class gaussmethod
{
    vector<shared_ptr<vector<real>>> mat;
    int dim;
public:
    gaussmethod<real, ext> (string filename)
    {
        ifstream in (filename);
 
        in >> dim;
        mat = vector<shared_ptr<vector<real>>> (dim);
        for (int i = 0; i < dim; i++)
        {
            mat[i] = make_shared<vector<real>> (dim + 1);
            for (int j = 0; j < dim + 1; j++)
                in >> (*mat[i])[j];
        }
    }
 
    gaussmethod<real, ext> (vector<vector<real>> &vv, vector<real> &v)
    {
        dim = vv.size ();
        mat = vector<shared_ptr<vector<real>>> (dim);
        for (int i = 0; i < dim; i++)
        {
            mat[i] = make_shared<vector<real>> (dim + 1);
            for (int j = 0; j < dim; j++)
                (*mat[i])[j] = vv[i][j];
        }
 
        for (int i = 0; i < dim; i++)
            (*mat[i])[dim] = v[i];
    }
 
    vector<real, ext> solve ()
    {
        vector<real> solution (dim);
        for (int col = 0; col < dim; col++)
        {
            auto max = make_pair (0, mat[0][col]);
            for (int i = 1; i < dim; i++)
            if (abs ((*mat[i])[col]) > abs (max.second))
                max = make_pair (i, (*mat[i])[col]);
            swap (mat[max.first], mat[col]);
            for (int j = col; j < dim + 1; j++)
                (*mat[i])[j] /= max.second;
            for (int i = col + 1; i < dim; i++)
            {
                real delta = (*mat[i])[col];
                for (int j = col; j < dim + 1; j++)
                    (*mat[i])[j] -= mat[max.first][j] * delta;
            }
        }
        for (int col = dim - 1; col >= 0; col--)
        {
            for (int j = col + 1; j < dim; j++)
                (*mat[i])[dim] -= solution[j] * (*mat[i])[j];
            solution[col] = (*mat[i])[dim] / (*mat[i])[col];
        }
        return solution;
    }
};
 
class generator
{
public:
    template<typename real, typename ext = real>
    static gaussmethod<real, ext> HilbertGauss (int k)
    {
        vector<vector<real>> vv (n, vector<real> (n));
        vector<real> v (n);
 
        for (int i = 0; i < n; i++)
        for (int j = 0; j < n; j++)
            v[i] += (vv[i][j] = 1. / (i + j + 1)) * (j + 1);
        return gaussmethod<real, ext> (vv, v);
    }
};
 
int main ()
{
    ifstream in ("input.txt");
    ofstream out ("output.txt");
    out << fixed;
    out.precision (15);
    for (int k = 1; k < 10; k++)
    {
        auto gm = generator::HilbertGauss<double> (k);
        auto v = gm.solve ();
        for (auto a : v)
            out << a << endl;
    }
 
    return 0;
}


Почему не компилируется этот код в Visual Studio 2013 ??
0
Надоела реклама? Зарегистрируйтесь и она исчезнет полностью.
 
КиберФорум - форум программистов, компьютерный форум, программирование
Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2017, vBulletin Solutions, Inc.