Ciencia computacional

Para algunos dominios convexos simples en 2D, tenemos algunos u ( x ) que satisfacen la siguiente ecuación: - d i v ( A ∇ u ) + c u n = f con ciertas condiciones de contorno de Dirichlet y / o Neumann. Que yo sepa, aplicar el método …

11 pde  finite-element 

Me gustaría visualizar una vibración molecular que no es un modo normal. Me gustaría presentar una representación estática, vectorial del movimiento, y me gustaría cierta flexibilidad en el estilo vectorial (tamaño, color, etc.). También estoy interesado en producir un video de la vibración. ¿Cuáles son algunos buenos recursos para mostrar …

11 computational-chemistry 

¿Alguien puede ayudarme a encontrar los libros sobre soluciones numéricas (diferencia finita y métodos de Crank-Nicolson) de Poisson y ecuaciones de difusión que incluyen ejemplos sobre geometría irregular, como un dominio que consiste en el área entre un rectángulo y un círculo (especialmente libros o enlaces en ejemplos de código …

11 pde  finite-difference 

Simultáneamente a mi investigación sobre los métodos del subespacio de Krylov, tengo la opción de explorar las matemáticas detrás de HPC un paso adelante o la teoría de la computación (hardware, sistema operativo, compiladores, etc.). Actualmente, sé tanto suficiente para simplemente sobrevivir. Por ejemplo, sé cómo derivar las ecuaciones para …

11 parallel-computing  education 

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 …

11 petsc  fortran 

Al resolver PDE dependientes del tiempo usando el método de elementos finitos, por ejemplo, digamos la ecuación de calor, si usamos pasos de tiempo explícitos, entonces tenemos que resolver un sistema lineal debido a la matriz de masa. Por ejemplo, si nos atenemos al ejemplo de la ecuación de calor, …

11 finite-element  time-integration  navier-stokes 

¿Existe un algoritmo SVD truncado que calcula los valores singulares uno a la vez? Mi problema: me gustaría calcular los primeros kkk valores singulares (y vectores singulares) de una gran matriz densa MMM , pero no sé cuál sería un valor apropiado de kkk . MMM es grande, por lo …

11 linear-algebra  algorithms  svd 

Leí en alguna parte que Pandas se desarrolló por primera vez para el mundo financiero, al menos no especialmente para las ciencias naturales (física, biología, etc.), entonces, ¿hay algún paquete similar de Python de análisis de datos más "orientado a las ciencias naturales"? Acabo de comenzar a usar Pandas y …

11 python  data-analysis  data-management 

Tengo una pregunta que es similar a esta antes, excepto en 3D, y solo necesito el volumen, no la forma real del casco. Más precisamente, se me da un pequeño conjunto de puntos (digamos, 10-15) en 3D, todos los cuales se sabe que se encuentran en el casco convexo del …

11 algorithms  reference-request  computational-geometry 

Para un proyecto, tengo que implementar estos dos métodos y comparar cómo funcionan en diferentes funciones. Parece que el método de gradiente conjugado está destinado a resolver sistemas de ecuaciones lineales de Ax=bAx=b A\mathbf{x} = \mathbf{b} Donde es una matriz n-por-n que es simétrica, positiva definida y real.AAA Por otro …

11 optimization  conjugate-gradient 

Estoy implementando el documento " Transporte masivo óptimo para el registro y la deformación ", mi objetivo es ponerlo en línea ya que simplemente no puedo encontrar ningún código de transporte masivo euleriano en línea y esto sería interesante al menos para la comunidad de investigadores en el procesamiento de …

11 pde  finite-difference  matlab  fluid-dynamics  image-processing 

Descubrí que gran parte de mi programación de ciencias computacionales tiene requisitos de prueba que no están cubiertos por los marcos de prueba estándar: Prueba de tiempo de computación Para asegurarse de que los algoritmos no sean más lentos. Podría hacer algo como, assureSmallerEqual(RuntimeWrapper(algorithm),53)pero me gustaría que el umbral de …

10 testing 

Cada matriz real se puede reducir a un formato Schur real de usando un ortogonal similiary transformada . Aquí la matriz es una forma cuasi triangular con 1 por 1 o 2 por 2 bloques en la diagonal principal. Cada 1 por 1 bloque corresponde a un valor propio real …

10 eigenvalues  eigensystem  lapack  numerics  scalapack 

Los algoritmos existentes para resolver ODE manejan funciones , dondey∈Rn. Pero en muchos sistemas físicos, la ecuación diferencial es autónoma, entoncesdyreyret= f( y, t )reyret=F(y,t)\frac{dy}{dt} = f(y, t)y∈ Rnortey∈Rnortey \in \mathbb R^n,y∈Rn, con latfuera. Con este supuesto simplificador, ¿qué mejoras se pueden ver en los métodos numéricos existentes? Por ejemplo, …

10 ode  numerics 

Deje ser matrices reales, cuadradas y densas. G y Q son simétricos. DejarA , G , QA,G,QA, G, QsolGGQQQ H= [ A- Q- G- AT]H=[A−G−Q−AT]H = \begin{bmatrix} A & -G \\ -Q &-A^T \end{bmatrix} ser una matriz hamiltoniana. Quiero calcular la exponencial de la matriz . Necesito la matriz exponencial …

10 linear-algebra  matrix  dense-matrix