Ciencia computacional

Me gustaría realizar una PCA en un conjunto de datos compuesto por aproximadamente 40 000 muestras, cada muestra mostrando alrededor de 10 000 características. El uso de la función princomp de Matlab lleva más de media hora, momento en el que elimino el proceso. Me gustaría encontrar una implementación / …

11 high-dimensional  data-analysis 

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 secuencia cuasialeatoria? ? Si es así, …

11 monte-carlo 

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 …

11 linear-algebra 

¿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 una matriz muy escasa con algunas diagonales, …

11 linear-solver  dense-matrix  linear-system 

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 …

11 matlab  computational-chemistry  publications 

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 para una implementación del algoritmo descrito en dicho documento. …

11 software  documentation 

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 simples …

11 parallel-computing  python 

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 …

11 matrix  matrix-equations 

¿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 modos?

11 python  graph-theory 

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 NNN. Las condiciones de contorno son: ψ(0)=0ψ(0)=0\psi(0) = 0 y ψ(R)=0ψ(R)=0\psi(R)=0 , y solo …

11 finite-difference  ode  accuracy 

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 inversa, en …

11 interpolation  unstructured-mesh 

¿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 vectores?

11 petsc  python 

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 detectar la multiplicidad de …

11 linear-algebra 

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 …

11 linear-algebra  optimization  matrix  iterative-method  linear-solver 

(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 …

11 regression  publications