Preguntas etiquetadas con c++

Estoy transfiriendo un código existente de MATLAB a C ++ y tengo un sistema lineal para resolver (en lugar de la forma más típica )A x = bxA=bxA=bxA=bAx=bAx=bAx=b La matriz es densa y de forma general, pero no es mayor que 1000x1000. Entonces, en MATLAB, la solución o la notación …

12 linear-solver  c++  lapack 

El documento "Plantillas de expresión revisadas: un análisis de rendimiento de metodologías actuales" en SIAM Journal of Scientific Computing hace referencia a la biblioteca de álgebra lineal "Blaze". No he oído hablar de él antes, y parece que no puedo encontrar referencias en línea. (Las búsquedas obvias de Google están …

12 linear-algebra  reference-request  c++ 

Tengo un programa de Mathematica que realiza algunas integrales en 3 o 4 dimensiones usando el QuasiMonteCarlométodo. El problema es que lleva mucho tiempo ejecutarlo, hasta el punto de que algunos de estos cálculos no se pueden completar en el tiempo de trabajo máximo disponible en nuestro clúster HPC. Así …

12 libraries  monte-carlo  c++  mathematica 

Estoy resolviendo un problema físico usando un esquema numérico implícito. Esto me lleva a resolver una ecuación lineal con matriz tridiagonal. He codificado este algoritmo de Wikipedia. Me pregunto si hay una biblioteca eficiente que permita resolver este tipo de ecuaciones de manera optimizada. Una nota importante es que la …

12 linear-algebra  libraries  c++ 

¿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() { v …

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 Actualmente estoy haciendo esto con la biblioteca de armadillo , pero resulta ser demasiado lento. La cuestión es que necesito hacer esto por un billón de matrices y resulta que calcular los dos determinantes es el cuello de botella de mi programa. …

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

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 ++ (gratuitas) que lo …

10 c++  quadrature 

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á escrito en …

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

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 compila con …

10 finite-element  performance  c++ 

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)={ 1Bk+1i(x)=x−xixk+i−xiBki+xk+i+1−xxk+i+1−xi+1Bki+1(x)Bik+1(x)=x−xixk+i−xiBik+xk+i+1−xxk+i+1−xi+1Bi+1k(x) B^{k+1}_{i}(x)= \frac{x-x_i}{x_{k+i}-x_i}B^k_i …

10 c++  b-spline 

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 interfaz / API es gruesa y pesada: debe traducir …

9 finite-element  python  c++  petsc  fenics 

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 como resultado un sistema tan grande:32323^2 ⎡⎣⎢- A0 00 00 00 0yo0 00 00 0yo⎤⎦⎥⎡⎣⎢Xyz⎤⎦⎥= λ ⎡⎣⎢UNA1yo0 0UNA20 0yoUNA30 …

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; } Lo anterior es mi código C …

9 c++  performance 

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 límites inferiores y fun2 …

9 c++  nonlinear-programming 

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 aquellos que …

9 matlab  python  c++