6
голосов
4ответов
5798 просмотров

std::pow gives a wrong approximation for fractional exponents

Here is what I mean trying to do double x=1.1402 double pow=1/3; std::pow(x,pow) -1; result is 0 but I expect 0.4465 the equation is (1 + x) ^3= 1.1402, find x.

23
голосов
8ответов
53652 просмотров

Implementing the derivative in C/C++

How is the derivative of a f(x) typically calculated programmatically to ensure maximum accuracy? I am implementing the Newton-Raphson method, and it requires taking of the derivative of a function.

6
голосов
4ответов
3626 просмотров

How to make topographic map from sparse sampling data?

I need to make a topographic map of a terrain for which I have only fairly sparse samples of (x, y, altitude) data. Obviously I can't make a completely accurate map, but I would like one that is in some sense "smooth". I need to quantify "smoothness" (probably the reciprocal the average of the ...

3
голосов
8ответов
3602 просмотров

Algorithm for multidimensional optimization / root-finding / something

I have five values, A, B, C, D and E. Given the constraint A + B + C + D + E = 1, and five functions F(A), F(B), F(C), F(D), F(E), I need to solve for A through E such that F(A) = F(B) = F(C) = F(D) = F(E). What's the best algorithm/approach to use for this? I don't care if I have to write it ...

107
голосов
13ответов
173651 просмотров

Как я могу использовать numpy.correlate для автокорреляции?

Мне нужно сделать автокорреляцию набора чисел, что, как я понимаю, просто корреляция набора с самим собой. Я пробовал использовать функцию корреляции numpy, но не верю результату, так как она почти всегда дает вектор, в котором первое число не является самым большим, как должно быть. . И...

7
голосов
6ответов
3248 просмотров

Как выполнить численное интегрирование с волновой функцией квантового гармонического осциллятора?

Как выполнить численное интегрирование (какой численный метод и какие приемы использовать) для одномерного интегрирования в бесконечном диапазоне, когда одна или несколько функций в подынтегральном выражении являются 1d квантовый гармонический осциллятор волновые функции. Среди прочего я хоч...

5
голосов
6ответов
13986 просмотров

Может кто-нибудь объяснить, почему scipy.integrate.quad дает разные результаты для одинаково больших расстояний при интеграции sin (X)?

Я пытаюсь численно интегрировать произвольную (известную, когда я кодирую) функцию в мою программу. с использованием методов численного интегрирования. Я использую Python 2.5.2 вместе с пакетом числовой интеграции SciPy. Чтобы почувствовать это, я решил попробовать интегрировать sin (x) и наблюд...