Вопросы с тегом «precision»

Вопросы, связанные с представлением числовых величин в конечном представлении в заданной базе, отличном от их точного математического значения.

77
Есть ли качественный решатель нелинейного программирования для Python?

У меня есть несколько сложных невыпуклых задач глобальной оптимизации. В настоящее время я использую MATLAB Optimization Toolbox (в частности, fmincon()с алгоритмом = 'sqp'), что довольно эффективно . Тем не менее, большая часть моего кода написана на Python, и я бы тоже хотел провести оптимизацию...

14
Научные вычисления на Python с современными графическими процессорами с двойной точностью

Кто-нибудь здесь использовал научные вычисления двойной точности с графическими процессорами нового поколения (например, K20) через Python? Я знаю, что эта технология быстро развивается, но каков наилучший способ сделать это в настоящее время? GPU не входит в число популярных научных библиотек...

14
Численно устойчивый способ вычисления углов между векторами

При применении классической формулы для угла между двумя векторами: α=arccosv1⋅v2∥v1∥∥v2∥α=arccos⁡v1⋅v2‖v1‖‖v2‖\alpha = \arccos \frac{\mathbf{v_1} \cdot \mathbf{v_2}}{\|\mathbf{v_1}\| \|\mathbf{v_2}\|} обнаруживается, что при очень малых / острых углах наблюдается потеря точности, и результат не...

13
Одинарная или двойная точность с плавающей точкой

Числа с плавающей запятой одинарной точности занимают половину памяти, и на современных машинах (кажется, даже на графических процессорах) операции могут выполняться с ними почти вдвое быстрее, чем с двойной точностью. Многие коды FDTD, которые я обнаружил, используют исключительно арифметику...

13
Почему плохо обусловленные линейные системы могут быть решены точно?

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

11
Числовая производная и конечно-разностные коэффициенты: какое-либо обновление метода Форнберга?

Когда кто-то хочет вычислить числовые производные, метод, представленный Бенгтом Форнбергом здесь (и сообщенный здесь ), очень удобен (точен и прост в реализации). Как оригинальная статья 1988 года, я хотел бы знать, есть ли лучшая альтернатива сегодня (как (или почти) как простая и более...

10
Диагонализация плотных плохо обусловленных матриц

Я пытаюсь диагонализировать некоторые плотные, плохо обусловленные матрицы. В машинной точности результаты являются неточными (возвращая отрицательные собственные значения, собственные векторы не имеют ожидаемой симметрии). Я переключился на функцию Eigensystem [] Mathematica, чтобы использовать...

10
Небольшие, непредсказуемые результаты в прогонах детерминированной модели

У меня есть значительная модель (~ 5000 строк), написанная на C. Это последовательная программа, нигде не генерирующая случайные числа. Она использует библиотеку FFTW для функций, использующих FFT - я не знаю деталей реализации FFTW, но я предполагаю, что функции в ней также являются...

9
Высокоточная арифметика с плавающей точкой в ​​числовом PDE

У меня сложилось впечатление, из очень разных источников и разговоров с исследователями, что существует растущий спрос на высокоточные вычисления в численных уравнениях в частных производных. Здесь высокая точность означает большую точность, чем стандартная 64-битная двойная точность. Интересно о...

9
Представление чисел Эйзенштейна без поплавков

У меня есть проект , где нужно использовать квадратные поля В частности число видаа + б- 3---√a+b−3a + b \sqrt{-3} с a,b∈Qa,b∈Qa,b \in \mathbb{Q}, Например, вот простые числа в числах Эйзенштейна : Я не хочу использовать шалфей. Я хотел бы написать свой собственный тип данных для включения numpy....