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

Я портирую существующий код из MATLAB на C ++ и имею линейную систему для решения (вместо более типичной формы )A x = bх = ЬxA=bxA=bA x = bAx=bAx=b Матрица плотная и общего вида, но не больше 1000x1000. Таким образом, в MATLAB решение находится по функции или через косую чертуAAAmrdivide(b,A)x = …

12 linear-solver  c++  lapack 

В статье «Пересмотренные шаблоны выражений: анализ эффективности современных методологий» в SIAM Journal of Scientific Computing упоминается библиотека линейной алгебры «Blaze». Я не слышал об этом раньше, и не могу найти ссылки в Интернете. (Очевидный поиск в Google возвращает вышеуказанную статью.) Так что же это за библиотека и где я могу …

12 linear-algebra  reference-request  c++ 

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

12 libraries  monte-carlo  c++  mathematica 

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

12 linear-algebra  libraries  c++ 

Как я могу заменить метод Эйлера на 4-й порядок Рунге-Кутты, чтобы определить движение свободного падения с не постоянной гравитационной величиной (например, свободное падение с высоты 10 000 км над землей)? До сих пор я писал простую интеграцию по методу Эйлера: while() { v += getMagnitude(x) * dt; x += v …

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

Мне дают матрицу Q которая является симметричной, обратимой, положительно определенной и плотной. Мне нужно проверить, если det, где J является матрицей всех единиц.12 × 1212×1212 \times 12QQQdet ( Q ) = det ( 12 I- Q - J)( 1 )йе(Q)знак равнойе(12я-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 

Я хотел бы использовать метод 8-го порядка Рунге-Кутты (89) в приложении по небесной механике / астродинамике, написанном на C ++, на машине Windows. Поэтому мне интересно, знает ли кто-нибудь хорошую библиотеку / реализацию, которая документирована и бесплатна для использования? Это нормально, если он написан на C, если нет никаких проблем …

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

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

10 finite-element  performance  c++ 

Сейчас я работаю над методом решения дифференциальных уравнений, который называется базис-сплайн-коллокацией. У меня возникли проблемы с созданием метода построения сплайна произвольного порядка с соотношением с начальным условием B 1 i (x)={ 1Вк + 1я( х ) = х - хяИкск + я- хяВКя+ хk + i + 1- хИксk + …

10 c++  b-spline 

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

9 finite-element  python  c++  petsc  fenics 

У меня есть большая проблема с кубическим собственным значением: (A0+ λA1+λ2A2+λ3A3) х =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 как большой: ⎡⎣⎢-A0000я000я⎤⎦⎥⎡⎣⎢ИксYZ⎤⎦⎥= λ⎡⎣⎢A1я0A20яA300⎤⎦⎥⎡⎣⎢ИксYZ⎤⎦⎥,[−A0000I000I][xyz]=λ[A1A2A3I000I0][xyz],\begin{bmatrix} -\mathbf{A}_0 & 0 & 0 …

9 c++  eigensystem  eigenvalues 

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 

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

9 c++  nonlinear-programming 

Поэтому я спорю о том, стоит ли мне заниматься изучением Python. Говоря с моими профессорами, Matlab, кажется, является общим языком, используемым в прикладной математике / вычислительной науке, насколько это касается научных кругов ; в то время как в промышленности мои профессора (особенно те, кто работал в промышленности) говорили, что изучение …

9 matlab  python  c++