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

Кубические сплайны - C++

Войти
Регистрация
Восстановить пароль
Другие темы раздела
C++ Параллельные запросы к MySQL с использованием OpenMP. http://www.cyberforum.ru/cpp/thread269487.html
Здравствуйте. Помогите пожалуйста разобраться. Программа посылает запросы к базе данных (субд MySQL), нужно их распараллелить. Сначала идут некоторые вычисления, подключение к MySQL, собственно кусок, который нужно распараллелить, чтобы одновременно выполнялись несколько запросов: for (int je = 0; je < j; je++) { for (int iks = 8; iks < 18; iks++) { #pragma omp parallel for ...
C++ enums Kак можно осуществить данную операцию? enum number{null=0x0,one=0x1,two=0x2,three=0x3}; void Main(){ number a=one; number b=two; b|=a;//b=b|a; //ругается здесь //error: this operation on an enumerated type requires an applicable user-defined operator function } Добавлено через 1 час 39 минут http://www.cyberforum.ru/cpp/thread269381.html
Выбор платформы C++
Привет!! Подскажите какая платформа с++ под windows наиболее востребована работодателями?
C++ Изменение готовой программы
Здравствуйте, сразу хочу сказать - С++ для меня неизведанная область, но вопрос, всё же задам. Итак! У нас в университете есть компьютер что-то IBM 486 SX, к нему подключена специальная установка (AC-магнитометр), на компьютере есть специальное приложение для управления установкой и снятия показания. Вопрос! Можно ли внести изменения в программу, она написана на C++ ? Дело в том, что...
C++ символы С++ http://www.cyberforum.ru/cpp/thread269042.html
мне надо в С++ вывести на экран некий символ, как мне это сделать? мне надо █ вставить, простой такой блок. вот я написал программу, которая должна была дать мне код этого символа: #include <iostream> #include <conio.h> using namespace std; int main(){ char f; cin >> f; int i = f;
C++ Прицепить g++ к программе Сабж. Люди, кто-нибудь таким занимался? Пишу, в общем, задачник (олимпиадный) по информатике, хотелось бы чтобы код на С++, который лежит в plain\text в базе данных, можно было скомпилировать, слинковать да запустить. Как это лучше всего организовать? подробнее

Показать сообщение отдельно
maxteam
2 / 2 / 0
Регистрация: 04.07.2010
Сообщений: 66

Кубические сплайны - C++

03.04.2011, 23:09. Просмотров 2977. Ответов 0
Метки (Все метки)

Здраствуйте..

У меня задано таблично {xi,yi} и мне нада с помощъю кубических сплайнов найти Y0 в точке Х0...

x 1 2 3 4 5 6
y 3 5 4 2 1 8

x0=2.4

нужна програмная реализация(желательно на С/С++ или Matcad). Я нашол какой то исходник на вики,но там все как то сложно и запутано...нада просто...

помогите пожалуйста!

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
#include <cstdlib>
#include <cmath>
#include <limits>
 
class cubic_spline
{
private:
        // Структура, описывающая сплайн на каждом сегменте сетки
        struct spline_tuple
        {
                double a, b, c, d, x;
        };
 
        spline_tuple *splines; // Сплайн
        std::size_t n; // Количество узлов сетки
 
        void free_mem(); // Освобождение памяти
 
public:
        cubic_spline(); //конструктор
        ~cubic_spline(); //деструктор
 
        // Построение сплайна
        // x - узлы сетки, должны быть упорядочены по возрастанию, кратные узлы запрещены
        // y - значения функции в узлах сетки
        // n - количество узлов сетки
        void build_spline(const double *x, const double *y, std::size_t n);
 
        // Вычисление значения интерполированной функции в произвольной точке
        double f(double x) const;
};
 
cubic_spline::cubic_spline() : splines(NULL)
{
 
}
 
cubic_spline::~cubic_spline()
{
        free_mem();
}
 
void cubic_spline::build_spline(const double *x, const double *y, std::size_t n)
{
        free_mem();
 
        this->n = n;
 
        // Инициализация массива сплайнов
        splines = new spline_tuple[n];
        for (std::size_t i = 0; i < n; ++i)
        {
                splines[i].x = x[i];
                splines[i].a = y[i];
        }
        splines[0].c = splines[n - 1].c = 0.;
 
        // Решение СЛАУ относительно коэффициентов сплайнов c[i] методом прогонки для трехдиагональных матриц
        // Вычисление прогоночных коэффициентов - прямой ход метода прогонки
        double *alpha = new double[n - 1];
        double *beta = new double[n - 1];
        alpha[0] = beta[0] = 0.;
        for (std::size_t i = 1; i < n - 1; ++i)
        {
                double h_i = x[i] - x[i - 1], h_i1 = x[i + 1] - x[i];
                double A = h_i;
                double C = 2. * (h_i + h_i1);
                double B = h_i1;
                double F = 6. * ((y[i + 1] - y[i]) / h_i1 - (y[i] - y[i - 1]) / h_i);
                double z = (A * alpha[i - 1] + C);
                alpha[i] = -B / z;
                beta[i] = (F - A * beta[i - 1]) / z;
        }
 
        // Нахождение решения - обратный ход метода прогонки
        for (std::size_t i = n - 2; i > 0; --i)
                splines[i].c = alpha[i] * splines[i + 1].c + beta[i];
 
        // Освобождение памяти, занимаемой прогоночными коэффициентами
        delete[] beta;
        delete[] alpha;
 
        // По известным коэффициентам c[i] находим значения b[i] и d[i]
        for (std::size_t i = n - 1; i > 0; --i)
        {
                double h_i = x[i] - x[i - 1];
                splines[i].d = (splines[i].c - splines[i - 1].c) / h_i;
                splines[i].b = h_i * (2. * splines[i].c + splines[i - 1].c) / 6. + (y[i] - y[i - 1]) / h_i;
        }
}
 
double cubic_spline::f(double x) const
{
        if (!splines)
                return std::numeric_limits<double>::quiet_NaN(); // Если сплайны ещё не построены - возвращаем NaN
 
        spline_tuple *s;
        if (x <= splines[0].x) // Если x меньше точки сетки x[0] - пользуемся первым эл-тов массива
                s = splines + 1;
        else if (x >= splines[n - 1].x) // Если x больше точки сетки x[n - 1] - пользуемся последним эл-том массива
                s = splines + n - 1;
        else // Иначе x лежит между граничными точками сетки - производим бинарный поиск нужного эл-та массива
        {
                std::size_t i = 0, j = n - 1;
                while (i + 1 < j)
                {
                        std::size_t k = i + (j - i) / 2;
                        if (x <= splines[k].x)
                                j = k;
                        else
                                i = k;
                }
                s = splines + j;
        }
 
        double dx = (x - s->x);
        return s->a + (s->b + (s->c / 2. + s->d * dx / 6.) * dx) * dx; // Вычисляем значение сплайна в заданной точке по схеме Горнера (в принципе, "умный" компилятор применил бы схему Горнера сам, но ведь не все так умны, как кажутся)
}
 
void cubic_spline::free_mem()
{
        if (splines)
        {
                delete[] splines;
                splines = NULL;
        }
}
Добавлено через 4 часа 26 минут
у меня уже дурка от этого бреда....
После регистрации реклама в сообщениях будет скрыта и будут доступны все возможности форума.
 
КиберФорум - форум программистов, компьютерный форум, программирование
Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2017, vBulletin Solutions, Inc.
Рейтинг@Mail.ru