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 +=...
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...
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 compilado com 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...
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,...
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...
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...
Estou usando o Tatu para fazer multiplicações de matriz muito intensivas com comprimentos laterais , onde pode ser de até 20 ou mais. Estou usando o Armadillo com o OpenBLAS para multiplicação de matrizes, o que parece estar fazendo um trabalho muito bom em núcleos paralelos, exceto que tenho um...
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...
Para o meu trabalho de pesquisa, desenvolvi uma pequena biblioteca C ++ destinada a facilitar a comunicação entre os códigos computacionais C ++ e o Octave / Matlab (quando este último é usado para fins de pós-processamento). Gostaria de liberar essa biblioteca sob a GPL em alguns sites de...
A maior parte da minha programação são códigos de pesquisa pontuais em C para uso próprio. Nunca distribuí nenhum código para outros colaboradores que não sejam colaboradores próximos. Eu desenvolvi um algoritmo que estou publicando em uma revista científica. Desejo fornecer o código fonte e talvez...
Existe alguma diretiva de pré-processador que possa ser usada para usar a função polylog? Ou está incluído no cmath? Se sim, você o chama por Li ou por polylog? EDIT: O que realmente estou tentando fazer é fornecer um valor analítico para a integral indefinida da função x3ex- 1x3ex−1...
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...
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...
Estou procurando uma biblioteca C ++ gratuita e conveniente para gráficos que incluem um solucionador para o problema do máximo subgrafo comum (MCS). Estou ciente da Boost Graph Library e do LEMON , mas nenhuma delas inclui um solucionador de
Ao escrever código MPI 3.0, que deve passar argce argvà MPI_Initchamada ou não, e por quê? EDIT: Como a resposta de Bill Barth levantou algumas perguntas, quero fazer algumas observações: A aprovação argc/ argvnão é necessária desde o MPI 1.1. A pergunta é especificamente sobre por que deveria /...
Meu projeto atual é uma reprogramação de um modelo de dobramento de proteínas envolvendo a solução de milhares de ODEs em C ++. Tenho parado e começado a progredir enquanto escrevia o solucionador para rodar totalmente na GPU. Finalmente consegui integrar, mas quando tento resolver dC / dt = -C...
Sou aluno de doutorado em Computação Científica e, nos últimos meses, passei bastante tempo aprendendo Python e C ++ da maneira certa. Sinto que aprendi bem C ++ e posso usar o Python para fazer o que quero, se mantiver um bom livro de referência. Eu também conheço o MATLAB o suficiente para...