Suponga que la siguiente matriz se da con su transpuesta . El producto rinde ,[ 0.500 - 0.333 - 0,167 - 0.500 0,667 - 0,167 - 0.500 - 0.333 0.833 ] A T A T A = G [ 0,750 - 0,334 - 0.417 - 0,334 0,667 - 0,333 - 0.417 - 0.333 0.750
Estoy resolviendo la ecuación diferencial con condiciones iniciales u (0) = u (1) = 0 , u '' (0) = u '' (1) = 0 . Aquí \ sigma (x) \ geqslant \ sigma_ {0}> 0 es el parámetro. En forma de operador podemos reescribir la ecuación diferencial como Au = f , donde el operador A es definitivo...
Enumere el paquete de Python (petsc4py, etc.) y los pocos solucionadores directos que admite. Una respuesta (community-wiki) por paquete, por
(Espero que esta pregunta se ajuste a este sitio; si no, acepte mis disculpas). Ejecuté una cierta simulación y obtuve una serie temporal y (t), t = 0, 1, ... 20. Después de probar algunas funciones, descubrí que: y(t) =~ 1 / (A t + B) Donde A y B son coeficientes que calculé usando regresión...
Supongamos que A es una matriz dispersa general, y quiero calcular los valores propios. No sé cómo detectar la multiplicidad de los valores propios. Hasta donde yo sé, para un caso especial, encontrar las raíces polinómicas mediante el método de matriz complementaria, podemos aplicar RRQR para...
Dado el sistema donde A ∈ R n × n , leí que, en caso de que la iteración de Jacobi se use como solucionador, el método no convergerá si b tiene un componente distinto de cero en el espacio nulo de A . Entonces, ¿cómo podría uno declarar formalmente eso, siempre que b tenga un componente distinto de...
¿Cuánto más lento es petsc4py vs c / c ++ / fortran? Me doy cuenta de que dependerá significativamente del código que se ejecute, pero ¿qué pasa con algo simple como un producto de matriz de
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...
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 ,...
Tengo conjuntos de datos de campo multipunto, cada conjunto de datos de punto se relaciona con una sola celda de una malla no estructurada. El objetivo es interpolar los datos al centro celular, directa o indirectamente, de la manera más precisa. Si uso la interpolación ponderada de distancia...
Estoy intentando resolver una ecuación del tipo: (−∂2∂x2−f(x))ψ(x)=λψ(x)(−∂2∂x2−f(x))ψ(x)=λψ(x) \left( -\tfrac{\partial^2}{\partial x^2} - f\left(x\right) \right) \psi(x) = \lambda \psi(x) Donde f(x)f(x)f(x) tiene un polo simple en 000 , para los valores propios y vectores propios más pequeños...
¿Hay algún algoritmo para la detección comunitaria de gráficos bipartitos (redes de 2 modos) implementados en igraph, networkX, R o Python, etc.? En particular, ¿existe tal implementación en la cual uno podría restringir la detección de comunidades solo en uno de los dos
La siguiente ecuación matricial en Σ - para matrices B y C dadas - aparece en mi trabajo como una caracterización de una matriz de covarianza. He aprendido que esta ecuación es conocida, en particular en la teoría del control de tiempo continuo, como la ecuación de Lyapunov , y que existen...
Tengo experiencia en la codificación de OpenMP para máquinas de memoria compartida (tanto en C como en FORTRAN) para llevar a cabo tareas simples como la suma de matrices, la multiplicación, etc. (Solo para ver cómo compite con LAPACK). Conozco OpenMP lo suficiente como para llevar a cabo tareas...
La revista Association for Computing Machinery Transactions on Mathematical Software (ACM TOMS) publica muchos artículos sobre algoritmos numéricos que incluyen implementaciones de software. De acuerdo con su política editorial , la presentación de un documento de algoritmo incluye el código fuente...
¿Hay alguna esperanza en resolver el siguiente sistema lineal de manera eficiente con un método iterativo? A∈Rn×n,x∈Rn,b∈Rn, with n>106A∈Rn×n,x∈Rn,b∈Rn, with n>106A \in \mathbb{R}^{n \times n}, x \in \mathbb{R}^n, b \in \mathbb{R}^n \text{, with } n > 10^6 Ax=bAx=bAx=b con , donde Δ es...
Supongamos que es una matriz simétrica real y su descomposición en autovalores V Λ V T es dado. Es fácil ver qué sucede con los valores propios de la suma A + c I donde c es una constante escalar (vea esta pregunta ). ¿Podemos sacar alguna conclusión en el caso general A + D donde D es una matriz...
Estoy produciendo imágenes 1d y 2d usando Matlab en este momento para comparar la precisión con un modelo dado. Necesito comparar mis métodos con el modelo gaussiano estándar .wfn y lo haré analizando la densidad de una molécula y la laplaciana. Estoy interesado en producir imágenes en 2D de las...
Una pregunta bastante simple: para hacer una integral multidimensional, dado que uno ha decidido que algún tipo de método de Monte Carlo es apropiado, ¿hay alguna ventaja que tenga una integración de MC regular usando números pseudoaleatorios sobre una integración cuasi-Monte Carlo usando una...
nVidia, por ejemplo, tiene CUBLAS, que promete una aceleración de 7-14x. Ingenuamente, esto no está cerca del rendimiento teórico de ninguna de las tarjetas GPU de nVidia. ¿Cuáles son los desafíos para acelerar el álgebra lineal en las GPU, y hay rutas de álgebra lineal más rápidas ya...