¿Cómo puedo reemplazar el método de Euler por el cuarto orden de Runge-Kutta para determinar el movimiento de caída libre en una magnitud gravitacional no constante (por ejemplo, caída libre desde 10 000 km sobre el suelo)? Hasta ahora escribí integración simple por el método de Euler: while() {...
Tengo mi propia subrutina para la integración numérica (cuadratura), que es una adaptación en C ++ de un programa ALGOL publicado por Bulirsch & Stoer en 1967 (Numerische Mathematik, 9, 271-278). Me gustaría actualizar a un algoritmo más moderno (adaptativo) y preguntarme si hay bibliotecas C...
C ++ 11 introduce la semántica de movimiento que puede, por ejemplo, mejorar el rendimiento del código en situaciones donde C ++ 03 necesitaría realizar una construcción de copia o asignación de copia. Este artículo informa que el siguiente código experimenta una velocidad 5 veces mayor cuando se...
En este momento estoy trabajando en un método de resolución de ecuaciones diferenciales llamado colocación de base-spline. Con lo que tengo problemas es con la construcción de un método para construir una spline de orden arbitrario, con la relación con la condición inicial B 1 i (x)={...
Me gustaría usar el método de octavo orden Runge-Kutta (89) en una aplicación de mecánica / astrodinámica celestial, escrita en C ++, usando una máquina Windows. Por lo tanto, me pregunto si alguien conoce una buena biblioteca / implementación documentada y de uso gratuito. Está bien si está...
Tengo un gran problema de valor propio cúbico: (A0+λA1+λ2A2+λ3A3)x=0.(UNA0 0+λUNA1+λ2UNA2+λ3UNA3)X=0.\left(\mathbf{A}_0 + \lambda\mathbf{A}_1 + \lambda^2\mathbf{A}_2 + \lambda^3\mathbf{A}_3\right)\mathbf{x} = 0. Podría resolver esto convirtiendo a un problema de valor propio lineal, pero daría...
Básicamente, FEM parece ser un problema que está prácticamente "resuelto". Existen numerosos marcos potentes, como Trilinos, PETSc, FEniCS, Libmesh o MOOSE. Una cosa que tienen en común: son extremadamente "pesados". Primero, la instalación normalmente es muy dolorosa. En segundo lugar, su...
Estoy usando Armadillo para hacer multiplicaciones matriciales muy intensas con longitudes laterales , donde puede ser de hasta 20 o incluso más. Estoy usando Armadillo con OpenBLAS para la multiplicación de matrices, que parece estar haciendo un muy buen trabajo en núcleos paralelos, excepto que...
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; } Lo anterior...
Para mi trabajo de investigación, he estado desarrollando una pequeña biblioteca de C ++ destinada a facilitar la comunicación entre los códigos computacionales de C ++ y Octave / Matlab (cuando este último se utiliza con fines de procesamiento posterior). Me gustaría lanzar una biblioteca de este...
La mayor parte de mi programación son códigos de investigación únicos en C para mi propio uso. Nunca he distribuido ningún código a otros que no sean colaboradores cercanos. He desarrollado un algoritmo que estoy publicando en una revista científica. Quiero proporcionar el código fuente y quizás el...
¿Hay alguna directiva de preprocesador que pueda usarse para usar la función polylog? ¿O está incluido en cmath? Si es así, ¿lo llamas Li o Polylog? EDITAR: lo que realmente estoy tratando de hacer es dar un valor analítico para la integral indefinida de la función X3miX- 1x3ex−1...
Así que he estado debatiendo si debería molestarme o no en aprender Python. Al hablar con mis profesores, Matlab parece ser el lenguaje común utilizado en matemáticas aplicadas / ciencias computacionales en lo que respecta a la academia ; mientras que en la industria , mis profesores (especialmente...
Actualmente estoy tratando de resolver el problema de minimización restringida no lineal implementado en la función matlab "fmincon". Mis expectativas son, minimizar (fun1, x0, uB, lB, fun2) donde x0 es el estado inicial, fun1 es una función que debe minimizarse, uB son límites superiores, lB son...
Estoy buscando una biblioteca C ++ conveniente y gratuita para gráficos que incluya un solucionador para el problema de subgrafo común máximo (MCS). Soy consciente de Boost Graph Library y LEMON , pero ninguno incluye un solucionador
Al escribir el código MPI 3.0, ¿debería pasar argcy argva la MPI_Initllamada o no, y por qué? EDITAR: Dado que la respuesta de Bill Barth planteó algunas preguntas, quiero hacer algunos comentarios: Pasar argc/ argvno es necesario desde MPI 1.1. La pregunta es específicamente sobre por qué...
Mi proyecto actual es una reprogramación de un modelo de plegamiento de proteínas que involucra la solución de miles de EDO en C ++. He estado deteniendo y comenzando a progresar mientras escribo el solucionador para que se ejecute totalmente en la GPU. Finalmente logré integrarlo, pero cuando...
Soy un estudiante de doctorado en Informática Científica y en los últimos meses, pasé una buena cantidad de tiempo aprendiendo Python y C ++ de la manera correcta. Siento que he aprendido bien C ++ y puedo usar Python para hacer lo que quiero si mantengo un buen libro de referencia. También...