La investigación reproducible en computación tiene como objetivo hacer que el código necesario para generar los resultados en un documento computacional esté disponible para otros investigadores para que puedan ejecutar este código para reproducir los resultados en ese documento. Me gustaría hacer...
Tengo un poco de dificultad para tratar de entender un artículo. El documento utiliza un método espectral para resolver un valor propio que proviene de un sistema de EDO acopladas. Escribiré solo una ecuación ahora, porque es suficiente para llegar al quid de mi (s) pregunta (s). La ecuación es...
n × n G n n GUNAA es una matriz dispersa simétrica positiva definida (SPD). es una matriz diagonal escasa. es grande ( > 10000) y el número de nonzeros en suele ser de 100 ~ 1000.n × nn×nn \times nsolGGnortennnortennsolGG L D L TUNAA ha sido factorizar en forma de Cholesky como .L D...
Para la solución de grandes sistemas lineales utilizando métodos iterativos, a menudo es interesante introducir el preacondicionamiento, por ejemplo, resolver M - 1 ( A x = b ) , donde M se usa aquí para el preacondicionamiento a la izquierda del sistema. Por lo general, deberíamos tener ese M - 1...
Supongamos que tengo una ecuación matricial definida recursivamente como A[n] = inverse([1 - b[n]A[n+1]]) * a[n] Entonces la ecuación para A [1] se parece a una fracción continua, para la cual existen algunos métodos altamente eficientes que evitan el tedioso recálculo (ver "Recetas numéricas"...
Las cantidades físicas como presión, densidad, energía, temperatura y concentración siempre deben ser positivas, pero los métodos numéricos a veces calculan valores negativos durante el proceso de solución. Esto no está bien porque las ecuaciones computarán valores complejos o infinitos...
Mientras leía algo de literatura sobre solucionadores de PDE, me encontré con el término pseudo paso del tiempo hoy. Parece ser un término común, sin embargo, no pude encontrar una buena definición o un artículo introductorio. Por lo tanto: ¿Qué es el seudo-paso de tiempo y cómo se usa...
Finalmente has logrado descubrir cómo están dispuestos espacialmente los átomos en tu entidad molecular recién descubierta. A través de, por ejemplo, medios espectroscópicos, ahora está en posesión de un montón de coordenadas de átomos, tipos de átomos, longitudes de enlaces, tipos de enlaces y...
Estoy tratando de resolver algunos problemas de optimización no lineal sin restricciones en la GPU (CUDA). La función objetivo es una función no lineal suave, y su gradiente es relativamente barato de calcular analíticamente, por lo que no necesito molestarme con la aproximación numérica. Quiero...
En la mayoría de las implementaciones (¿todas?) Del Método rápido multipolar (FMM), se utilizan octrees para descomponer el dominio relevante. Teóricamente, los octrees proporcionan un límite volumétrico simple, que es útil para probar el tiempo de ejecución O (n) de un FMM. Más allá de esta lógica...
Me gustaría aprender cómo se pueden aplicar los métodos wavelet a PDE, pero desafortunadamente no conozco un buen recurso para aprender sobre este tema. Parece que muchas introducciones a las wavelets se centran en la teoría de la interpolación, por ejemplo, ensamblando una señal mediante una...
Quiero desarrollar un software paralelo de computación científica desde cero. Quiero algunas ideas sobre qué idioma comenzar. El programa implica leer / escribir datos en archivos txt y hacer cálculos pesados en paralelo, con muchas factorizaciones LU y el uso de solucionadores lineales...
¿Cuáles son los mejores algoritmos (escalabilidad y eficiencia) para generar mallas cuádruples no estructuradas en 2D? ¿Dónde puedo encontrar un buen generador de malla cuádruple no estructurado? (código abierto
Me gustaría usar f2pycon Fortran moderno. En particular, estoy tratando de que el siguiente ejemplo básico funcione. Este es el ejemplo útil más pequeño que pude generar. ! alloc_test.f90 subroutine f(x, z) implicit none ! Argument Declarations ! real*8, intent(in) :: x(:) real*8, intent(out)...
Estoy tratando de implementar la siguiente función en coma flotante de doble precisión con un error relativo bajo : logsum(x,y)=log(exp(x)+exp(y))logsum(x,y)=log(exp(x)+exp(y))\mathrm{logsum}(x,y) = \log(\exp(x) + \exp(y)) Esto se usa ampliamente en aplicaciones estadísticas para agregar...
¿Quieres mejorar esta publicación? Proporcione respuestas detalladas a esta pregunta, incluidas citas y una explicación de por qué su respuesta es correcta. Las respuestas sin suficiente detalle pueden ser editadas o eliminadas. Refinamiento de malla adaptable (AMR) es...
Estoy buscando un software de generación de mallas que es gratis y de código abierto, proporciona una interfaz de secuencia de comandos sana para la especificación de dominio, funciona para geometrías complejas, puede generar mallas 2D y 3D, ¿Que opciones
Tengo un no-convexa función 2-D acotada, que me gustaría encontrar el mínimo de. La función es bastante suave. Evaluar es costoso. Un error aceptable es de aproximadamente 3% del dominio de la función en cada eje. He intentado ejecutar la implementación del algoritmo DIRECTO en la biblioteca...
Cuando los cálculos limitados de ancho de banda de memoria se llevan a cabo en entornos de memoria compartida (por ejemplo, roscado a través de OpenMP, Pthreads, o TBB), hay un dilema de cómo garantizar que la memoria se distribuye correctamente a través de física de la memoria, de tal manera que...
Supongamos que tengo alguna función y quiero encontrar tal que . Puede ser que utilice el método de Newton-Raphson. Pero esto requiere que sé que el derivado de la función . Una expresión analítica para puede no estar disponible. Por ejemplo, puede definirse mediante un código informático...