Preguntas etiquetadas con fortran

En este programa de muestra estoy haciendo lo mismo (al menos eso creo) de dos maneras diferentes. Estoy ejecutando esto en mi PC con Linux y monitoreando el uso de memoria con top. Usando gfortran encuentro que en la primera forma (entre "1" y "2") la memoria utilizada es de …

13 performance  fortran 

Los estudios de caso en mi tesis doctoral requieren que tenga extensiones de intervalo de subrutinas Fortran en CHEMKIN-II (disculpas por el enlace; es el mejor que pude encontrar para un paquete que Sandia National Laboratory ya no distribuye). CHEMKIN-II es un paquete químico usado en química de combustión; Como …

13 fortran  c  automatic-differentiation  interval-arithmetic 

La conclusión aquí: ¿Cuánto mejor son realmente los compiladores de Fortran? es que gfortran y gcc son tan rápidos para un código simple. Entonces quería probar algo más complicado. Tomé el ejemplo de tiroteo de la norma espectral. Primero precalculo la matriz 2D A (:, :), y luego calculo la …

13 fortran  c 

Me estoy metiendo en Fortran 95 para algunas simulaciones de mecánica cuántica. Honestamente, Octave me ha echado a perder, así que he dado por sentado la exponenciación de matrices. Dada una (pequeña, ) sesgo -matriz hemitiana de tamaño n × n , ¿cuál es la forma más eficiente de usar …

11 fortran  lapack 

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 

Soy nuevo en computación científica. Estoy buscando una implementación Fortran (preferiblemente f90) de un Octree. Mi problema requiere un Octree que divide mi dominio hasta que no haya más de algunas partículas de N (o fuentes donde conozco el valor de densidad que se puede conectar en un método de …

10 fortran  integral-equations 

Hace poco hice una pregunta en la misma línea para las matrices ermitas sesgadas. Inspirado por el éxito de esa pregunta, y después de golpearme la cabeza contra la pared durante un par de horas, estoy mirando la matriz exponencial de las matrices asimétricas reales. La ruta para encontrar los …

10 linear-algebra  fortran  lapack 

Recientemente, he encontrado un extraño problema con FORTRAN95. Inicialicé las variables X e Y de la siguiente manera: X=1.0 Y=0.1 Luego los agrego e imprimo el resultado: 1.10000000149012 Después de examinar las variables, parece que 0.1 no está representado con doble precisión con total precisión. Hay alguna manera de evitar …

10 fortran  floating-point 

Me preguntaba cuál es la mejor manera de aprender a programar en Fortran. He estado buscando libros, los ejemplos que me parecieron bien fueron "Explicado Fortran 95/2003" y "Fortran para científicos e ingenieros" . No sé qué piensan ustedes de esto. Otra cosa que me pregunto es qué compilador debo …

9 fortran 

Como lo ejemplifica la respuesta de Jed Brown a los Costos de las búsquedas versus los cálculos , el uso de operaciones de punto flotante vectorizado versus no vectorizado resulta en un código mucho más rápido. Muchos compiladores modernos afirman que pueden realizar una vectorización automática. ¿Cómo puedo saber qué …

9 fortran  c  compiling  vectorization 

¿Cuál es la forma más eficiente de obtener una matriz dispersa compleja de mi código Fortran a PETSc? Entiendo que esto depende del problema, por lo que traté de dar tantos detalles relevantes como sea posible a continuación. He estado jugando con el solucionador de valores propios FEAST [1] para …

8 matrix  petsc  performance  fortran  sparse-matrix 

He implementado CG en FORTRAN al vincularlo a Intel MKL. Cuando hay declaraciones como: ( Consulte Wikipedia ) p=r; x=x+alpha*p r=r-alpha*Ap; o similares en QMR (en una cantidad mucho mayor) v_tld = r; y = v_tld; rho = norm( y ); w_tld = r; z = w_tld; xi = norm( …

8 performance  fortran  blas 

Estoy investigando el uso de redes neuronales artificiales (ANN) para predecir las velocidades de reacción en mi fluido en lugar de resolver el sistema completo de EDO rígidas. Algunas personas de mi laboratorio ya han trabajado en eso, así que no empiezo desde cero, pero tengo problemas con mis aplicaciones. …

8 performance  fortran  python  data-management