Mi sistema es un problema de FE simétrico con multiplicadores lagrange (por ejemplo, flujo incompresible de Stokes): (ABBTC)(ABTBC)\begin{pmatrix}A & B^T \\ B & C\end{pmatrix} donde es el caso típico (incluso me he asegurado de que las ecuaciones estén numeradas para que los multiplicadores de...
¿Cuánto más lento es petsc4py vs c / c ++ / fortran? Me doy cuenta de que dependerá significativamente del código que se ejecute, pero ¿qué pasa con algo simple como un producto de matriz de
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...
El cómputo paralelo se modela frecuentemente utilizando una tasa local determinista de cómputo, sobrecarga de latencia y ancho de banda de red. En realidad, estos son espacialmente variables y no deterministas. Estudios como Skinner y Kramer (2005) observan distribuciones multimodales, pero el...
Aquí hay un uso típico de los punteros de función en C. Me gustaría hacer algo similar en Fortran. Tengo algunas ideas, pero me gustaría saber si hay alguna forma canónica de hacerlo. Los punteros de función y los contextos pasados por el usuario se almacenan y luego se llaman más...
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á...
Que se calcula más rápido, o log a c o b √unasiaba^bIniciar sesiónunaClogac\log_a c ? un,bycson números reales positivos conb>1.C√sicb\sqrt[b]{c}unaaasibbCccb > 1b>1b>1 ¿Qué tipo de algoritmos usarás en la comparación? ¿Cuáles son sus complejidades? Por ejemplo, cuando o c ≈ a bc ≡...
Soy nuevo en el paquete PETSc. Tengo un ~ 4000x4000 matriz A en formato de mercado de matriz y quiero obtener PETSc para resolver esto usando múltiples procesadores. Sé cómo resolver el sistema en un único procesador, pero no sé cómo distribuir la matriz y los vectores entre diferentes...
Tengo un modelo considerable (~ 5000 líneas) escrito en C. Es un programa en serie, sin generación de números aleatorios en ninguna parte. Utiliza la biblioteca FFTW para funciones que usan FFT: no conozco los detalles de la implementación de FFTW, pero supongo que las funciones allí también son...
Estoy buscando referencias creíbles que indiquen la cantidad de recursos que las supercomputadoras gastan en coordinarse en comparación con el trabajo real relacionado con las tareas. Los recursos podrían estar disponibles en potencia de procesamiento, pero incluso los vatios parecen una unidad...
En los últimos años, han aparecido varias bibliotecas / proyectos de software que ofrecen una forma u otra de paralelismo de memoria compartida de uso general basado en datos. La idea principal es que, en lugar de escribir un código explícitamente enhebrado, los programadores implementan sus...
Soy muy nuevo en la programación GPGPU, así que perdónenme si la pregunta no es particularmente apropiada. Por lo que entiendo, la programación de GPU es un trabajo de ingeniería muy complejo en comparación con la programación habitual de CPU. Hay que tener mucho cuidado con los problemas de...
Existen muchos métodos numéricos bien conocidos para resolver ecuaciones del tipo por ejemplo, método de bisección, método de Newton, etc.f( x ) =0,x ∈ Rnorte,f(X)=0 0,X∈Rnorte, f(x) = 0, \quad x \in \mathbb{R}^n, En mi aplicación, se calcula con un método estocástico (el resultado es un...
¿Compilar PETSc con una biblioteca externa BLAS / LAPACK afecta significativamente el rendimiento en matrices dispersas, o solo usa esas bibliotecas para matemáticas de matriz
Como lo ejemplifica la respuesta de Jed Brown a los Costos de las búsquedas versus los cálculos , el uso de operaciones de punto flotante vectorizado versus no vectorizado resulta en un código mucho más rápido. Muchos compiladores modernos afirman que pueden realizar una vectorización automática....
Necesito resolver Ax = b, pero me doy cuenta de que incluso si es escaso, almacenar los coeficientes matriciales de mi problema requerirá demasiada memoria. Así que ahora estoy considerando usar un método sin matriz, porque los mismos coeficientes aparecen mucho tiempo en la matriz, por lo que...
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...
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...
Me estoy sumergiendo en el fascinante mundo del análisis de elementos finitos y me gustaría resolver un gran problema termomecánico (solo térmico mecánico, sin retroalimentación).→→\rightarrow Para el problema mecánico, ya entendí por la respuesta de Geoff , que necesitaré usar un solucionador...
¿Cuál es la forma más eficiente de obtener una matriz dispersa compleja de mi código Fortran a PETSc? Entiendo que esto depende del problema, por lo que traté de dar tantos detalles relevantes como sea posible a continuación. He estado jugando con el solucionador de valores propios FEAST [1] para...