Perguntas com a marcação «c++»

Estou portando um código existente do MATLAB para C ++ e tenho um sistema linear para resolver (em vez da forma mais típica A x = b )x A = bxA=bxA=bA x = bAx=bAx=b A matriz é densa e de forma geral, mas não é maior que 1000x1000. Portanto, no …

12 linear-solver  c++  lapack 

O artigo "Modelos de Expressão Revisitados: Uma Análise de Desempenho das Metodologias Atuais" no SIAM Journal of Scientific Computing faz referência à biblioteca de álgebra linear "Blaze". Eu nunca ouvi falar disso antes e não consigo encontrar referências on-line. (As pesquisas óbvias no Google estão devolvendo o artigo acima.) Então, …

12 linear-algebra  reference-request  c++ 

Eu tenho um programa Mathematica que executa algumas integrais em 3 ou 4 dimensões usando o QuasiMonteCarlométodo O problema é que demora muito tempo para ser executado, a ponto de alguns desses cálculos não serem concluídos no tempo máximo de trabalho disponível em nosso cluster HPC. Então, eu estou pensando …

12 libraries  monte-carlo  c++  mathematica 

Estou resolvendo um problema físico usando esquema numérico implícito. Isso me leva a resolver uma equação linear com matriz tridiagonal. Eu codifiquei esse algoritmo da Wikipedia. Gostaria de saber se existe uma biblioteca eficiente que permita resolver esse tipo de equação de maneira otimizada. Uma observação importante é que a …

12 linear-algebra  libraries  c++ 

Como posso substituir o método de Euler pela 4ª ordem de Runge-Kutta para determinar o movimento de queda livre em magnitude gravimétrica não constante (por exemplo, queda livre de 10.000 km acima do solo)? Até agora escrevi uma integração simples pelo método Euler: while() { v += getMagnitude(x) * dt; …

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

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 Atualmente, estou fazendo isso com a biblioteca do tatu , mas ela acaba sendo muito lenta. O problema é que eu preciso fazer isso por um trilhão de matrizes e acontece que calcular os dois determinantes é o gargalo do meu programa. …

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

Eu tenho minha própria pequena sub-rotina para integração numérica (quadratura), que é uma adaptação em C ++ de um programa ALGOL publicado por Bulirsch & Stoer em 1967 (Numerische Mathematik, 9, 271-278). Gostaria de atualizar para um algoritmo (adaptável) mais moderno e me pergunto se existem bibliotecas C ++ (gratuitas) …

10 c++  quadrature 

Eu gostaria de usar o método de 8ª ordem Runge-Kutta (89) em um aplicativo de mecânica / astrodinâmica celeste, escrito em C ++, usando uma máquina Windows. Portanto, pergunto-me se alguém conhece uma boa biblioteca / implementação documentada e livre para usar? Tudo bem se estiver escrito em C, desde …

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

O C ++ 11 introduz a semântica de movimentação que pode, por exemplo, melhorar o desempenho do código em situações em que o C ++ 03 precisaria executar uma construção de cópia ou atribuição de cópia. Este artigo relata que o código a seguir experimenta uma velocidade de 5x quando …

10 finite-element  performance  c++ 

No momento, estou trabalhando em um método de solução de equações diferenciais chamado colocação de base-spline. Estou com problemas para criar um método para construir um spline de ordem arbitrária, com o relacionamento com o condição inicial e estou tendo problemas até para começar com esse problema, pois é recursivo …

10 c++  b-spline 

Basicamente, o FEM parece ser um problema praticamente "resolvido". Existem inúmeras estruturas poderosas, como Trilinos, PETSc, FEniCS, Libmesh ou MOOSE. Uma coisa que eles têm em comum: são extremamente "pesados". Primeiro, a instalação normalmente é super dolorosa. Segundo, a interface / API deles é espessa e pesada - você precisa …

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; } O código acima é o meu …

9 c++  performance 

Eu tenho um grande problema de autovalor cúbico: (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. Eu poderia resolver isso convertendo para um problema de autovalor linear, mas resultaria em um sistema maior:32323^2 ⎡⎣⎢−A0000I000I⎤⎦⎥⎡⎣⎢xyz⎤⎦⎥=λ⎡⎣⎢A1I0A20IA300⎤⎦⎥⎡⎣⎢xyz⎤⎦⎥,[−A0000I000I][xyz]=λ[A1A2A3I000I0][xyz],\begin{bmatrix} -\mathbf{A}_0 & 0 & 0 \\ 0 & \mathbf{I} & 0 \\ 0 & 0 …

9 c++  eigensystem  eigenvalues 

Atualmente, estou tentando resolver um problema de minimização com restrições não-lineares, conforme implementado na função matlab "fmincon". Minhas expectativas são: minimizar (fun1, x0, uB, lB, fun2) onde x0 é o estado inicial, fun1 é uma função que precisa ser minimizada, uB são limites superiores, lB são limites inferiores e lB …

9 c++  nonlinear-programming 

Então, eu estou debatendo se devo ou não aprender o Python. Ao falar com meus professores, o Matlab parece ser a linguagem comum usada na matemática aplicada / ciência da computação no que diz respeito à academia ; enquanto na indústria , meus professores (especialmente aqueles que trabalharam na indústria) …

9 matlab  python  c++