Ciencia computacional

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...

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

¿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

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 f (x) , incluso si tiene una...

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...

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 los clústeres de...

Estoy buscando libros o artículos, o publicaciones de blog, o cualquier material publicado en general, que aborde específicamente los usos de las características modernas de C ++ (semántica de movimiento, STL, iteradores, evaluación diferida, etc.) en informática científica. ¿Puedes sugerir...

Estoy usando el paquete de solución no lineal de PETSc SNES para resolver un sistema de ecuaciones no lineales obtenidas al discretizar una ecuación diferencial parcial. ¿Cómo puedo determinar por qué el solucionador no converge y qué puedo hacer para resolver con éxito mis...

Se da una secuencia finita creciente de números positivos . ¿Cuál de los siguientes dos algoritmos es mejor para calcular la suma de los números?z1, z2, . . . . . znortez1,z2,.....znortez_{1} ,z_{2},.....z_{n} s=0; for \ i=1:n s=s + z_{i} ; end O: s=0; for \ i=1:n s=s + z_{n-i+1} ;...

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

Cuál sería una buena discretización de diferencia finita para la siguiente ecuación: ∂ρ∂t+∇⋅(ρu)=0∂ρ∂t+∇⋅(ρu)=0\frac{\partial \rho}{\partial t} + \nabla \cdot \left(\rho u\right)=0 ? Podemos tomar el caso 1D: ∂ρ∂t+ddx(ρu)=0∂ρ∂t+ddx(ρu)=0\frac{\partial \rho}{\partial t} + \frac{d}{dx}\left(\rho...

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 una sola 1 entrada con el resto siendo 0 's, de tal manera que G T G = I . A es enorme,...

¿Cuál es el enfoque preferido y eficiente para interpolar datos multidimensionales? Cosas que me preocupan: rendimiento y memoria para la construcción, evaluación individual / por lotes manejo de dimensiones de 1 a 6 lineal o de orden superior capacidad de obtener gradientes (si no...

Voy a enseñar a los estudiantes de pregrado un curso titulado Introducción a la programación de computadoras. Estoy un poco confundido En Física Computacional, los científicos usan C / C ++ o Python o Fortran, CUDA, etc. ... este es el momento de construir su base. ¿Qué debo usar? Sé que puedes...

El paradigma de programación OpenCL promete ser un estándar abierto libre de regalías para la computación heterogénea. ¿Deberíamos invertir nuestro tiempo en desarrollar software basado en OpenCL? ¿Pros

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...

Tengo una función bidimensional cuyos valores me gustaría muestrear. La función es muy costosa de calcular y tiene una forma compleja, por lo que necesito encontrar una manera de obtener la mayor información sobre su forma utilizando el menor número de puntos de muestra.F( x , y)f(x,y)f(x,y) ¿Qué...

Estoy buscando algoritmos para dibujar gráficos 2D estándar para funciones que pueden tener o no singularidades. El propósito es escribir un "Mini-CAS", por lo que no tengo conocimiento a priori de los tipos de funciones que los usuarios quieren graficar. Este problema es muy antiguo, así que...

Comenzaré con mi experiencia personal en nuestro laboratorio. En los últimos 9 y 10 días, solíamos ser bastante agresivos con las optimizaciones, compilando con -O3 e indicadores específicos del procesador (-xW -xSSE4.2 por ejemplo). Pero comenzando con ifort 11, comenzamos a notar: 1. algunas...

Estoy buscando resolver un problema de optimización restringida donde conozco los límites de algunas de las variables (específicamente una restricción en recuadro). argmintuF( u , x )arg⁡minuf(u,x) \arg \min_u f(u,x) sujeto a a ≤ d ( u , x ) ≤ bc ( u , x ) = 0c(u,x)=0 c(u,x) = 0 a ≤ d( u , x )...