Estoy tratando de familiarizarme con los SDE y he estado leyendo algunos artículos de revisión sobre el tema. Dejan la impresión de que se ha dedicado una gran cantidad de trabajo a los solucionadores sin derivados. A mi entender, esto significa que para un DDE como las derivadas de f y g no son...
En este momento me quedé con un problema. Parece ser realmente trivial, pero aún así es difícil para mí encontrar una solución adecuada. El problema es: uno tiene dos intervalos y debe encontrar la intersección de ellos. Por ejemplo: La intersección de [0, 3] y [2, 4] es [2, 3] La intersección...
Es bastante doloroso descubrir que un trabajo de varios días se eliminará prematuramente debido a un error al establecer el límite de tiempo de pared. ¿Hay alguna manera de cambiarlo para un trabajo PBS en ejecución
¿Cómo puedo calcular el complemento Schur? S= Kb b- Kb aK- 1a aKa bS=Ksisi-KsiunaKunauna-1Kunasi S = K_{bb} - K_{ba} K_{aa}^{-1} K_{ab} dónde K= ( Ka aKb aKa bKb b)K=(KunaunaKunasiKsiunaKsisi) K=\begin{pmatrix} K_{aa} & K_{ab} \\ K_{ba} & K_{bb} \end{pmatrix} (en algún orden) es una matriz...
Utilizo un ejemplo de la teoría de elementos finitos, pero cualquiera que mantenga una gran estructura de datos y la extienda sucesivamente encontrará algo similar. Suponga que tengo una malla no estructurada de puntos y triángulos, donde los puntos están dados por coordenadas (digamos e ) y los...
Se están escribiendo muchos artículos sobre cómo el recientemente lanzado Intel Xeon Phi robará la cuota de mercado HPC \ Super Computer de los competidores. Intel Knights está equipado con 72 núcleos y 4 zócalos, lo que lo convierte en un sistema de 288 núcleos. Mientras que un solo Gtx980 tiene...
Estoy perfilando el rendimiento de los solucionadores lineales de PETSc. Según lo entiendo, speedup=Sequential TimeParallel Time.speedup=Sequential TimeParallel Time.\text{speedup}=\frac{\text{Sequential Time}}{\text{Parallel Time}}. Sé que ejecutar el código paralelo en un procesador se puede...
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é...
¿Cuáles son los "pros y contras" de estos dos métodos de refinamiento de malla? Ambos parecen ser los métodos predominantes. Naturalmente, puedo imaginar que el refinamiento rojo global es comparativamente fácil de implementar y atractivo para la investigación en métodos de cuadrículas múltiples,...
Cualquiera conoce referencias que resumen la historia y las ideas detrás de la supercomputación, incluida la mención de desarrollos en lenguajes de programación paralelos, aplicaciones, empresas de inicio (algunas se establecieron pero no sobrevivieron), necesidades de la industria,...
Más o menos la pregunta. Dada una matriz general dispersa, no simétrica (tanto numérica como estructuralmente), ¿qué importancia tiene el patrón de dispersión (es decir, permutación de fila / columna de matriz / vector) para los solucionadores iterativos? Puedo ver que se vuelve importante para los...
¿Existe una implementación de código abierto y segura para subprocesos de solucionadores de optimización convexa en C / C ++? Algunas bibliotecas como NLopt , Ipopt , OPT ++ no cumplen mis requisitos. OPT ++ e Ipopt no son seguros para subprocesos, y NLopt no parece tener una forma específica /...
Tengo un problema de matriz generalizado: de un método espectral en un problema de análisis de estabilidad lineal. Mi matriz B es diagonal y positiva semi-definida. A es no ermitaño y complejo.A x = λ B xUNAX=λsiXA x = \lambda B x Mi problema es esencialmente que cuando uso el solucionador de...
¿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...
Implementé un Navier Stokes basado en la discretización FEM y PETSc para resolver el sistema lineal de ecuaciones. Para crear un procedimiento de solución eficiente, sigo el documento "Precondicionamiento eficiente de las ecuaciones de Navier-Stokes linealizadas para flujo incompresible" (Silvester...