Запитання з тегом «c++»

Я переношу існуючий код з MATLAB на C ++ і маю лінійну систему для вирішення (а не більш типову форму )xA=bxA=bxA=bAx=bAx=bAx=b Матриця щільна і загальної форми, але не більше 1000х1000. Отже, у MATLAB рішення знайдено за функцією або позначенням вперед-косою рисоюAAAmrdivide(b,A)x = b/A; Як я можу це вирішити у своєму …

12 linear-solver  c++  lapack 

У роботі "Переглянуті шаблони виразів: аналіз ефективності сучасних методологій" в журналі SIAM Journal of Scientific Comput посилається на бібліотеку лінійних алгебр "Blaze". Я про це раніше не чув і не можу знайти онлайн-посилання. (Очевидні пошукові запити google повертають вищезгаданий папір.) Отже, що це за бібліотека і де я можу дізнатися …

12 linear-algebra  reference-request  c++ 

У мене є програма Mathematica, яка виконує деякі інтеграли в 3 або 4 вимірах за допомогою QuasiMonteCarloметоду. Проблема полягає в тому, що для запуску потрібен прикро довгий час, до тих пір, коли деякі з цих обчислень не зможуть виконатись за максимальний час роботи, наявний у нашому кластері HPC. Тому я …

12 libraries  monte-carlo  c++  mathematica 

Я вирішую фізичну задачу за допомогою неявної числової схеми. Це призводить мене до розв’язування лінійного рівняння з тридіагональною матрицею. Я зашифрував цей алгоритм із Вікіпедії. Цікаво, чи існує ефективна бібліотека, яка дозволяє вирішити цей тип рівнянь оптимізованим способом. Важливе зауваження полягає в тому, що сама матриця змінюється лише тоді, коли …

12 linear-algebra  libraries  c++ 

Як я можу замінити метод Ейлера на Runge-Kutta 4-го порядку, щоб визначити рух вільного падіння в не постійній гравітаційній величині (наприклад, вільне падіння з 10 000 км над землею)? Поки я писав просту інтеграцію методом Ейлера: while() { v += getMagnitude(x) * dt; x += v * dt; time += …

11 numerical-analysis  c++  runge-kutta  integral-equations 

Мені дано матрицю яка симетрична, обернена, позитивно визначена і щільна. Мені потрібно перевірити, чи де J - матриця всіх.12×1212×1212 \times 12QQQdet(Q)=det(12I−Q−J)(1)det(Q)=det(12I−Q−J)(1)\det(Q) = \det(12I-Q-J) \; \; (1)JJJ Зараз я це роблю з бібліотекою броненосців, але це виявляється занадто повільно. Вся справа в тому, що мені потрібно це зробити на трильйон матриць, …

11 linear-algebra  matrix  c++  graph-theory  c 

У мене є своя маленька підпрограма для числової інтеграції (квадратури), яка є адаптацією C ++ програми ALGOL, опублікованої Bulirsch & Stoer у 1967 році (Numerische Mathematik, 9, 271-278). Мені хотілося б перейти до більш сучасного (адаптивного) алгоритму і поцікавитися, чи існують (безкоштовні) бібліотеки C ++, які надають таке. Я мав …

10 c++  quadrature 

Я хотів би використовувати метод Runge-Kutta 8-го порядку (89) у застосуванні небесної механіки / астродинаміки, написаному на C ++, за допомогою машини Windows. Тому мені цікаво, чи хтось знає хорошу бібліотеку / програму, яка є документально підтвердженою та вільною для використання? Це нормально, якщо це написано на C, якщо не …

10 c++  reference-request  libraries  c  runge-kutta 

C ++ 11 вводить семантику переміщення, яка може, наприклад, покращити продуктивність коду в ситуаціях, коли C ++ 03 потрібно буде виконати побудову копії або призначення копії. Ця стаття повідомляє, що наступний код має 5-кратну швидкість, коли компілюється з C + 11: vector<vector<int> > V; for(int k = 0; k < …

10 finite-element  performance  c++ 

На даний момент я працюю над методом розв'язання диференціальних рівнянь, який називається колокалізація базисно-сплайн. З чим у мене виникають проблеми - це побудова методу побудови сплайну довільного порядку із співвідношенням з початкова умова і у мене виникають проблеми, навіть коли я починаю цю проблему, оскільки рекурсивна робота може починатися з …

10 c++  b-spline 

В основному FEM, здається, є проблемою, яка в значній мірі "вирішена". Існує безліч потужних рамок, таких як Trilinos, PETSc, FEniCS, Libmesh або MOOSE. У них є одне спільне: вони надзвичайно «важкі». По-перше, установка зазвичай дуже болісна. По-друге, їх інтерфейс / API товстий і важкий - вам потрібно перекласти всю свою …

9 finite-element  python  c++  petsc  fenics 

double trap(double func(double), double b, double a, double N) { double j; double s; double h = (b-a)/(N-1.0); //Width of trapezia double func1 = func(a); double func2; for (s=0,j=a;j<b;j+=h){ func2 = func(j+h); s = s + 0.5*(func1+func2)*h; func1 = func2; } return s; } Наведене вище - мій C ++ …

9 c++  performance 

У мене є велика проблема кубічного власного значення: (A0+λA1+λ2A2+λ3A3)x=0.(A0+λA1+λ2A2+λ3A3)x=0.\left(\mathbf{A}_0 + \lambda\mathbf{A}_1 + \lambda^2\mathbf{A}_2 + \lambda^3\mathbf{A}_3\right)\mathbf{x} = 0. Я міг би вирішити це, перетворившись на лінійну задачу власного значення, але це призведе до того, що система є великою:32323^2 ⎡⎣⎢−A0000I000I⎤⎦⎥⎡⎣⎢xyz⎤⎦⎥=λ⎡⎣⎢A1I0A20IA300⎤⎦⎥⎡⎣⎢xyz⎤⎦⎥,[−A0000I000I][xyz]=λ[A1A2A3I000I0][xyz],\begin{bmatrix} -\mathbf{A}_0 & 0 & 0 \\ 0 & \mathbf{I} & 0 \\ …

9 c++  eigensystem  eigenvalues 

На даний момент я намагаюся вирішити нелінійну задачу мінімізації з обмеженими можливостями, що реалізується у функції matlab "fmincon". Мої очікування - мінімізувати (fun1, x0, uB, lB, fun2), де x0 - початковий стан, fun1 - це функція, яку потрібно мінімізувати, uB - це верхня межа, lB - нижня межа, а fun2 …

9 c++  nonlinear-programming 

Тому я обговорював, чи варто мені турбуватися над вивченням Python. Якщо говорити з моїми професорами, Матлаб, здається, є загальною мовою, що використовується у прикладній математиці / обчислювальній науці, що стосується наукових шкіл; тоді як у галузі , мої професори (особливо ті, хто працював у галузі) сказали, що навчання с ++ …

9 matlab  python  c++