Ciencia computacional

Un problema muy común en Markov Chain Monte Carlo implica calcular probabilidades que son la suma de términos exponenciales grandes, ea1+ea2+...ea1+ea2+... e^{a_1} + e^{a_2} + ... donde los componentes de lata pueden variar de muy pequeños a muy grandes. Mi enfoque ha sido factorizar el término exponencial más grande para …

24 monte-carlo  floating-point  statistics 

Supongamos que estoy trabajando en un código científico en C ++. En una discusión reciente con un colega, se argumentó que las plantillas de expresión podrían ser algo realmente malo, lo que podría hacer que el software sea compilable solo en ciertas versiones de gcc. Supuestamente, este problema ha afectado …

24 performance  c++  automatic-differentiation  interval-arithmetic 

Por lo que puedo decir, los dos grandes marcos genéricos de software de ciencia computacional del Departamento de Energía de EE . UU . Son PETSc y Trilinos . Parecen similares a primera vista, más allá de las diferencias en el lenguaje (C versus C ++). ¿Cuáles son las principales …

24 petsc  software  trilinos 

No estoy ejecutando ningún código paralelo en este momento, pero preveo ejecutar código paralelo en el futuro usando un híbrido de OpenMP y MPI. Los depuradores han sido herramientas invaluables para mí al ejecutar proyectos en serie. ¿Alguien puede recomendar un depurador paralelo (o múltiples depuradores) para usar para depurar …

24 software  parallel-computing 

Descripción del experimento: En la interpolación de Lagrange, la ecuación exacta se muestrea en puntos (orden polinomial ) y se interpola en 101 puntos. Aquí es variado de 2 a 64. Cada vez , y parcelas de error se preparan. Se ve que, cuando la función se muestrea en puntos …

24 finite-difference  interpolation  error-estimation  polynomials 

Según tengo entendido, hay dos categorías principales de métodos iterativos para resolver sistemas lineales de ecuaciones: Métodos estacionarios (Jacobi, Gauss-Seidel, SOR, Multigrid) Métodos del subespacio de Krylov (gradiente conjugado, GMRES, etc.) Entiendo que la mayoría de los métodos estacionarios funcionan relajando iterativamente (suavizando) los modos de Fourier del error. Como …

24 linear-algebra  convergence  krylov-method  conjugate-gradient 

Un clúster de cómputo Amazon EC2 cuesta alrededor de $ 800- $ 1000 (dependiendo del ciclo de trabajo) por núcleo físico de la CPU en el transcurso de 3 años. En nuestra última ronda de adquisición de hardware, mi laboratorio recogió 48 núcleos de hardware muy similares a los de …

23 cloud-computing  hpc 

Tengo una serie de puntos de datos que espero (aproximadamente) sigan una función que asíntota a una línea en grande . Esencialmente, acerca a cero como , y probablemente se puede decir lo mismo de todas las derivadas , , etc. Pero no sé cuál es la forma funcional de …

23 algorithms 

Muchos enfoques numéricos para CFD pueden extenderse a un orden arbitrariamente alto (por ejemplo, métodos discontinuos de Galerkin, métodos WENO, diferenciación espectral, etc.). ¿Cómo debo elegir un orden de precisión apropiado para un problema dado?

23 fluid-dynamics  accuracy 

La reproducibilidad se ha vuelto cada vez más importante en la investigación en ciencias computacionales. (Por ejemplo, vea este artículo de Roger Peng en Science ; también conozco otros artículos y sitios web de este tipo). Sin embargo, no tengo claro cuánta información debo incluir en un artículo de revista …

23 software  publications  documentation  reproducibility 

¿Cuál es el estado del arte en la aproximación de integrales altamente oscilatorias en una dimensión y dimensiones superiores con precisión arbitraria?

23 quadrature  precision  numerical-analysis 

Esta pregunta se migró de Cross Validated porque se puede responder en Computational Science Stack Exchange. Migrado hace 7 años . Quiero aprender la programación de MATLAB para poder realizar algunas investigaciones / análisis por mi cuenta y también, para poder estudiar / modificar algunos scripts de MATLAB que he …

23 software  matlab  languages  octave 

Tal como está actualmente, esta pregunta no es adecuada para nuestro formato de preguntas y respuestas. Esperamos que las respuestas sean respaldadas por hechos, referencias o experiencia, pero esta pregunta probablemente solicitará debate, argumentos, encuestas o discusión extendida. Si cree que esta pregunta se puede mejorar y posiblemente volver a …

22 languages 

Como sabemos, los métodos de transformaciones ortogonales (rotaciones de Givens y reflexiones de Housholder) para sistemas de ecuaciones lineales son más caros que la eliminación gaussiana, pero en teoría tienen mejores propiedades de estabilidad en el sentido de que no cambian el número de condición del sistema. Aunque conozco solo …

22 linear-algebra  reference-request 

Tengo matrices y G . A es escaso y es n × n con n muy grande (puede ser del orden de varios millones). G es una matriz alta de n × m con m bastante pequeño ( 1 < m < 1000 ) y cada columna solo puede tener …

22 linear-algebra  linear-solver  sparse-matrix  krylov-method