Описание тега numerical-computing

Это взаимосвязанное сочетание информатики и математики.
0 ответов

Ошибка вычисления Javascript (числовая ошибка в простом вычислении)

РЕШИТЬ: Читать: Математика с плавающей точкой нарушена? Как бороться с точностью чисел с плавающей точкой в ​​JavaScript? Просто получите любую среду обработки JavaScript и введите: 150.3 + 150.03 Простой расчет, но, похоже, он не дает правильного р…
15 апр '14 в 07:02
0 ответов

Интеграция с переменными пределами

Мне нужно оценить двойной интеграл, где внутренняя верхняя граница является переменной:integral2 between -5 and 5 ( integral1 between 0 and y f(x)dx )dy, Я застрял в расчете внешнего цикла, который зависит от внутреннего цикла. Мой код работает очен…
2 ответа

Ошибка чисел Python: большое число добавить небольшое число

Это упражнение из курса глубокого обучения Udacity. Может кто-нибудь объяснить, почему окончательный ответ не 1,0? v1 = 1e9 v2 = 1e-6 for i in range(int(1e6)): v1 = v1 + v2 print 'answer is', v1 - 1e9 # answer is 0.953674316406
2 ответа

Решение неквадратной линейной системы с помощью R

Как решить неквадратную линейную систему с R: A X = B? (в случае, если система не имеет решения или бесконечно много решений) Пример: A=matrix(c(0,1,-2,3,5,-3,1,-2,5,-2,-1,1),3,4,T) B=matrix(c(-17,28,11),3,1,T) A [,1] [,2] [,3] [,4] [1,] 0 1 -2 3 [2…
2 ответа

Загрузка разреженной матрицы scipy в существующую плотную матрицу

Скажем, у меня огромная матрица А, занимающая десятки гигабайт. Для выделения этой памяти требуется немалое количество времени. Допустим, у меня также есть коллекция скудных разреженных матриц с теми же размерами, что и у матрицы с пустыми краями. И…
27 янв '12 в 23:03
0 ответов

Почему сочетание обратного БПФ (c2r) и прямого БПФ с реальным вводом (r2c) в один комплексный массив не приводит к созданию исходного комплексного массива?

БПФ реального ввода обычно выражается rfft и обратное irfft во многих языках программирования (например, python или же julia). Как сказали numy docs или julia docs, они делают нормализованный FFT, а не как FFTWНормализация БПФ Следовательно, если у …
08 авг '17 в 15:01
2 ответа

Почему этот код возвращает два разных значения, делающих то же самое?

#include <stdio.h> #include <stdlib.h> #include <math.h> int main() { double a; double b; double q0 = 0.5 * M_PI + 0.5 * -2.1500000405000002; double q1 = 0.5 * M_PI + 0.5 * 0.0000000000000000; double w0 = 0.5 * M_PI + 0.5 * -43000.…
1 ответ

Быстрая и простая установка и использование библиотеки IT++ в Linux/Windows?

IT++ - это популярная библиотека GPL (которую сейчас называют LGPL, но пока без признаков) для научных вычислений, особенно в области телекоммуникаций. Это было рекомендовано мне моими коллегами. Кажется, трудно установить этот кусок библиотеки в м…
10 дек '16 в 16:15
6 ответов

Какова относительная скорость сложения с плавающей запятой и умножения с плавающей запятой?

Десять или два года назад стоило написать числовой код, чтобы не использовать умножения и деления, а вместо этого использовать сложение и вычитание. Хорошим примером является использование прямых разностей для оценки полиномиальной кривой вместо неп…
2 ответа

Лучший способ проверить, равен ли double отрицательная бесконечность в C++

Я нашел это: http://en.cppreference.com/w/cpp/numeric/math/isinf но, похоже, он проверяет положительную или отрицательную бесконечность. Я просто хочу проверить, равно ли значение точно отрицательной бесконечности, или в других словах log(0) Спасибо…
23 фев '15 в 21:22
1 ответ

Как избежать вызовов виртуальных функций в вычислительном графе

Я использую DAG (направленный ациклический граф) для представления и оценки выражений; каждый узел представляет операцию (+, -, /, *, накапливать и т. д.), а оценка всего выражения достигается путем последовательной оценки каждого узла в топологичес…
0 ответов

Как построить трехмерную напряженность электрического поля в MATLAB

В общем, я пытаюсь построить график напряженности электрического поля в данной точке [x,y,z] в MATLAB. Я хочу изобразить напряженность поля в данной точке с точки зрения цвета; высокая напряженность поля в точке [x_1, y_1, z_1] будет, например, пока…
02 авг '18 в 15:04
2 ответа

Библиотека численных вычислений для C

Я ищу хорошую библиотеку численных вычислений для языка программирования Си. Подобно LAPACK с хорошей документацией было бы очень полезно. С уважением, Upul
01 фев '11 в 06:41
5 ответов

C++ std::vector для HPC?

Я перевожу программу, которая выполняет числовое моделирование с FORTRAN на C++. Мне приходится иметь дело с большими матрицами типа double размером 800MB. это double M[100][100][100][100]; выдает ошибку сегментации, потому что стек не такой большой…
17 сен '16 в 12:34
0 ответов

Проверка отмены вычитания (точность с плавающей точкой)

Когда вы читаете о численных методах, то отмена вычитания часто является темой. Простой пример отмены: a = 1; b = 1e-16; a-(a+b) = 0 следовательно, мы теряем всю информацию о б. Тем не менее, я редко читаю о том, как проверить эту проблему. Разве мы…
17 окт '18 в 13:08
2 ответа

Лучший тип коллекции Scala для векторизованных числовых вычислений

Ищете правильный тип данных (например, IndexedSeq[Double]) использовать при разработке предметной библиотеки числовых вычислений. Для этого вопроса я ограничиваю работу с 1-мерными массивами Double, Библиотека будет определять числовые функции, кото…
1 ответ

Является ли вычислительно более сложным вычисление экспоненциальной или бесселевой функции?

Часто желательно получить решение математической задачи в замкнутой форме, то есть в виде выражения, содержащего общепринятые функции, такие как полиномы, рациональные и иррациональные функции, корни, а также экспоненты и логарифмы. Я часто слышу об…
1 ответ

Как вывести 2 или более массивов в функции фортрана?

Я пишу программу, которая вычисляет декомпозицию матрицы LU с частичным поворотом, и я хотел бы, чтобы функция выводила несколько (2 или 3) матриц без запуска программы несколько раз для вывода каждой из них по отдельности, что является пустой трато…
3 ответа

Какой способ является более точным?

Мне нужно разделить числовой диапазон на несколько сегментов, имеющих одинаковую длину. Но я не могу решить, какой путь более точный. Например: double r1 = 100.0, r2 = 1000.0, r = r2 - r1; int n = 30; double[] position = new double[n]; for (int i = …
2 ответа

Небольшая числовая ошибка при расчете среднего веса

Вот часть физического движка.Упрощенная функция centerOfMass вычисляет 1D-центр масс двух твердых тел ( демо ):- #include <iostream> #include <iomanip> float centerOfMass(float pos1,float m1, float pos2,float m2){ return (pos1*m1+pos2*m2…
21 янв '18 в 04:10