Supongamos que tengo dos matrices Nx2, Mx2 que representan N, M 2d vectores respectivamente. ¿Hay una manera simple y buena de calcular distancias entre cada par de vectores (n, m)? La manera fácil pero ineficiente es, por supuesto: d = zeros(N, M); for i = 1:N, for j = 1:M, d(i,j) =...
Estoy tratando de averiguar si hay una forma más rápida de calcular todos los valores propios y vectores propios de una matriz de adyacencia muy grande y dispersa que usando scipy.sparse.linalg.eigsh Hasta donde sé, este método solo usa la dispersión y atributos de simetría de la matriz. Una matriz...
Tanto en los métodos de descomposición de dominio (DD) como en multigrid (MG), se puede componer la aplicación de las actualizaciones de bloque o las correcciones generales como aditivas o multiplicativas . Para los solucionadores puntuales, esta es la diferencia entre las iteraciones de Jacobi y...
Considere dos computadoras con diferentes configuraciones de hardware y software. Cuando se ejecuta exactamente el mismo código de Navier-Stokes en serie en cada plataforma, lleva tiempo xey ejecutar una iteración para la computadora 1 y 2, respectivamente. En este caso, , es la diferencia de...
El cómputo paralelo se modela frecuentemente utilizando una tasa local determinista de cómputo, sobrecarga de latencia y ancho de banda de red. En realidad, estos son espacialmente variables y no deterministas. Estudios como Skinner y Kramer (2005) observan distribuciones multimodales, pero el...
Este es un simple problema de enumeración de simetría. Doy todos los antecedentes aquí, pero no se necesita conocimiento de química cuántica. La integral de dos partículas es: Y tiene las siguientes 4 simetrías: Tengo una función que calcula la integral y la almacena en una matriz 1D ,...
Estoy buscando la literatura más importante y las referencias de diapositivas para modelar los costos de energía actuales y futuros de las operaciones de punto flotante y las transferencias de datos a través de la CPU, la memoria, la red y el almacenamiento. He marcado esta pregunta como un wiki...
¿La multiplicación de matrices (Mat * Mat y Mat * Vec) se escala con el número de ceros o con el tamaño de la matriz? O alguna combinación de los dos. ¿Qué pasa con la forma? Por ejemplo, tengo una matriz de 100 x 100 con 100 valores, o una matriz de 1000 x 1000 con 100 valores. Al cuadrar estas...
Tengo que hacer una transformación de coordenadas entre dos sistemas de referencia (ejes). Para eso, se deben multiplicar tres matrices ( ) debido a que se utilizan algunos ejes intermedios. He pensado en dos enfoques para resolver esto:3 × 33×33\times3 Método # 1 : Hacer la multiplicación...
C ++ 11 introduce la semántica de movimiento que puede, por ejemplo, mejorar el rendimiento del código en situaciones donde C ++ 03 necesitaría realizar una construcción de copia o asignación de copia. Este artículo informa que el siguiente código experimenta una velocidad 5 veces mayor cuando se...
Estoy tratando de entender la diferencia entre estas dos tarjetas gráficas para la informática académica, específicamente para el componente DGEMM. Si miramos las estadísticas en bruto, ambas tienen el mismo chip GK110, tienen estadísticas comparables en prácticamente todas las categorías y,...
Tengo un problema en el que necesito encontrar todos los pares propios positivos (ya que el valor propio es positivo) de una matriz no simétrica pequeña (generalmente más pequeña que 60x60). Puedo dejar de calcular cuando el valor propio es menor que un cierto umbral. Sé que los valores propios son...
El título dice la mayor parte. Estoy buscando una biblioteca liviana y fácil de usar que pueda usar para proyectos de Android (NDK). Para cosas densas me gusta usar Eigen pero no he encontrado muchas bibliotecas completas (¡y documentadas!) Para cosas dispersas que "simplemente funcionan" en un...
double trap(double func(double), double b, double a, double N) { double j; double s; double h = (b-a)/(N-1.0); //Width of trapezia double func1 = func(a); double func2; for (s=0,j=a;j<b;j+=h){ func2 = func(j+h); s = s + 0.5*(func1+func2)*h; func1 = func2; } return s; } Lo anterior...
¿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...
He visto muchas publicaciones en revistas de física computacional que utilizan diferentes métricas para el desempeño de su código. Especialmente para el código GPGPU, parece que hay una gran variedad de resultados de tiempo que la gente publica. En particular, he visto: Comparaciones de...
Es común decir que Linux es más rápido y por buenas razones. Pero como se indica en el título, ¿cuánto importa la elección del sistema operativo para el rendimiento del código de computación científica? Para algo en lo que el sistema operativo está más directamente involucrado (como jugar juegos o...
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( z ); (and more) ¿Tiene...
Estoy perfilando el rendimiento de los solucionadores lineales de PETSc. Según lo entiendo, speedup=Sequential TimeParallel Time.speedup=Sequential TimeParallel Time.\text{speedup}=\frac{\text{Sequential Time}}{\text{Parallel Time}}. Sé que ejecutar el código paralelo en un procesador se puede...
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...