Escribí un programa / biblioteca que usé para obtener resultados en un artículo. (Aquí está , pero mi pregunta es general). Tengo pruebas que ejecuto regularmente ctest(me toma unos minutos ejecutarlas). Para reproducir algunas tablas o figuras en el artículo, tengo que construir un script o un...
Tengo una noción muy básica sobre la aritmética de intervalos (IA), pero parece ser una rama muy interesante de la ciencia computacional, tanto teórica como prácticamente. Está claro que las aplicaciones obvias son problemas informáticos verificados y mal planteados, pero esto es demasiado...
Estaba leyendo esto en Math SE. La pregunta básica es: Suponga que alguien desea estudiar algo avanzado; Una forma de hacerlo sería comenzar desde lo básico y construir. Pero la "imagen más grande" puede perderse en este proceso. Un método más (que prefiero llamar Wiki recursivo) es recoger un...
Estoy muy impresionado con el rendimiento en serie de los preacondicionadores de ILU basados en inversos multinivel , particularmente para Helmholtz heterogéneo , pero me sorprende no poder encontrar ninguna implementación de código abierto. En particular, ILUPACK pone los binarios a disposición...
Tengo un sistema lineal de ecuaciones de tamaño mxm, donde m es grande. Sin embargo, las variables que me interesan son solo las primeras n variables (n es pequeña en comparación con m). ¿Hay alguna manera de aproximar la solución para los primeros valores de m sin tener que resolver todo el...
Me gustaría determinar el número teórico de FLOP (operaciones de punto flotante) que puede hacer mi computadora. ¿Puede alguien ayudarme con esto? (Me gustaría comparar mi computadora con algunas supercomputadoras solo para tener una idea de la diferencia entre
Pregunta: Suponga que tiene dos preacondicionadores diferentes (factorizados) para una matriz simétrica positiva definida : y donde las inversas de los factores son Fácil de aplicar.UNUNAA ≈ BTsiUN≈siTsiA \approx B^TBA ≈ CTC,UN≈CTC,A \approx C^TC,B , BT, C, CTsi,siT,C,CTB, B^T, C, C^T Cuando...
¿Existe algún algoritmo (eficiente) para seleccionar un subconjunto de puntos de un conjunto de puntos ( ) de modo que "cubran" la mayor parte del área (sobre todos los subconjuntos posibles de tamaño )?MMMNNNM<NM<NM < NMMM Supongo que los puntos están en plano 2D. El ingenuo algoritmo es...
¿Cuáles son los métodos más modernos para la solución numérica de EDO con lado derecho discontinuo? Me interesan principalmente las funciones del lado derecho liso por partes, por ejemplo, signo. Estoy tratando de resolver la ecuación de un tipo siguiente: X˙v˙= v= { ( | FexternoEl | - |...
He tenido una idea del análisis numérico (principalmente, métodos numéricos como la búsqueda de raíces, ecuaciones cuadráticas y otras cosas preliminares) en mi clase de cálculo, pero ahora me encuentro con ganas de más sofisticación en mi trabajo. ¿Existe un buen libro que me ayude a comprender...
Estoy optimizando una función de 10-20 variables. La mala noticia es que cada evaluación de la función es costosa, aproximadamente 30 minutos de cálculo en serie. La buena noticia es que tengo un clúster con unas pocas docenas de nodos computacionales a mi disposición. Por lo tanto, la pregunta:...
¡Estoy enseñando una clase de encuesta de análisis numérico y estoy buscando motivación para el método BFGS para estudiantes con experiencia / intuición limitada en optimización! J k ( → x k - → x k - 1 ) = f ( → x k ) - f ( → x k - 1
He estado usando HDF5 durante años, pero a medida que aumenta el tamaño del conjunto de datos, estoy empezando a experimentar los mismos problemas enumerados aquí. http://cyrille.rossant.net/moving-away-hdf5/ ¿Me puede indicar un formato como HDF5 con - Buen soporte para escrituras paralelas -...
Regularmente compito en los llamados "Concursos de programación", donde resuelves problemas algorítmicos difíciles con tu propio código y habilidades de resolución de problemas durante un período de tiempo limitado. Para obtener ejemplos referenciales de cómo se verían, busque competencias como,...
Dejado x, yser dos números de punto flotante. ¿Cuál es la forma correcta de calcular su media? La forma ingenua (x+y)/2puede dar lugar a desbordamientos cuando xy yson demasiado grandes. Creo que 0.5 * x + 0.5 * ytal vez sea mejor, pero implica dos multiplicaciones (que tal vez sea ineficiente), y...
Después de hacer algunas matemáticas relacionadas con la estabilidad de los elementos en el problema de 3D Stokes, me sorprendió un poco darme cuenta de que PAG2- P1PAG2-PAG1P_2-P_1 no es estable para una malla tetraédrica arbitraria. Más precisamente, en caso de que tenga un elemento donde todos...
A veces, mientras se optimiza el código, es necesario cronometrar ciertas partes del código, he estado usando lo siguiente durante años, pero me preguntaba si hay una manera más simple / mejor de hacerlo. call system_clock(count_rate=clock_rate) !Find the time rate call
Estoy tratando de entender algunos resultados y agradecería algunos comentarios generales sobre cómo abordar problemas no lineales. La ecuación de Fisher (una PDE de reacción-difusión no lineal), tut= dtux x+ βu ( 1 - u ) = F( u )ut=duxx+βu(1−u)=F(u) u_t = du_{xx} + \beta u (1 - u) = F(u) en...
Perdón por la larga publicación, pero quería incluir todo lo que pensé que era relevante al principio. Lo que quiero Estoy implementando una versión paralela de los métodos de Krylov subespacio para matrices densas. Principalmente GMRES, QMR y CG. Me di cuenta (después del perfil) que mi rutina...
Me gustaría desarrollar sistemas de resolución de EDO en GPU, en un entorno 'trivialmente paralelo'. Por ejemplo, haciendo un análisis de sensibilidad con 512 conjuntos de parámetros diferentes. Idealmente, quiero resolver ODE con un solucionador de paso de tiempo inteligente y adaptable como...