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...
CVXOPT: http://abel.ee.ucla.edu/cvxopt/index.html OpenOpt: http://openopt.org/Welcome ¿Cuál es la relación entre ellos? ¿Cuáles son las ventajas / desventajas de ellos, respectivamente? Por cierto, ¿hay alguna otra biblioteca de optimización convexa de propósito general de alta calidad para...
¿Existe una lista publicada de las mejores prácticas para garantizar la longevidad del código, con miras a resultados científicos reproducibles? (por ejemplo, código abierto, prácticas de documentación, selección de dependencias, selección de un idioma, máquinas virtuales, etc.). Conozca cualquier...
La mayoría de las bibliotecas matemáticas tienen varias versiones de funciones de logaritmo. La mayoría de las veces suponemos que son perfectos, pero en realidad muchos de ellos solo ofrecen un cierto número de dígitos de precisión. Para algunas funciones, hay variantes numéricamente más...
Hola comunidad scicomp, He trabajado en el área de algoritmos gráficos usando frameworks como NetworkX (Python), JUNG e YFiles (Java). Ahora estoy entrando en el área de computación paralela y de alto rendimiento. Para un nuevo proyecto, estoy buscando una biblioteca de gráficos C ++ con las...
Entonces, el teorema de descomposición de Cholesky establece que cualquier matriz real simétrica positiva definida tiene una descomposición de Cholesky donde es una matriz triangular inferior.MMMM=LL⊤M=LL⊤M= LL^\topLLL Dado , ya sabemos que existen algoritmos rápidos para calcular su factor...
Estoy resolviendo para una enorme matriz definida positiva escasa A usando el método de gradiente conjugado (CG). ¿Es posible calcular el determinante de A utilizando la información producida durante la resolución?A x =
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...
¿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. Soy completamente nuevo en la noción de...
Digamos que tengo una función que deseo integrar sobre un tetraedro T ⊂ R 3 . Si f fuera arbitrario, la cuadratura de Gauss sería una buena solución, pero sé que f es armónico. ¿Cuánto se puede acelerar la cuadratura de Gauss con esta información?f:R3→Rf:R3→Rf : \mathbf{R}^3 \to \mathbf{R}T⊂R3T⊂R3T...
Teóricamente, ¿cómo escala el tiempo para hacer un cálculo de la teoría funcional de densidad (DFT) con el número de electrones? Estoy interesado en implementaciones DFT "típicas" como VASP, ABINIT, etc., no en códigos O
Tengo un sistema lineal no homogéneo A x = bAx=b Ax=b donde es una matriz real n × n con n ≤ 4 . Se garantiza que el espacio nulo de A es de dimensión cero, por lo que la ecuación tiene un inverso único x = A - 1 b . Dado que el resultado entra en el lado derecho de una EDO, que tengo la...
Dada una función desconocida , podemos evaluar su valor en cualquier punto de su dominio, pero no tenemos su expresión. En otras palabras, f es como una caja negra para nosotros.f:Rd→Rf:Rd→Rf:\mathbb R^d \to \mathbb Rfff ¿Cuál es el nombre del problema de encontrar el minimizador de ? ¿Cuáles son...
¿Cómo MATLAB, por ejemplo, calcula la SVD de una matriz dada? Supongo que la respuesta probablemente implica calcular los vectores propios y los valores propios de A*A'. Si ese es el caso, también me gustaría saber cómo los
Tengo un problema físico regido por la ecuación de Poisson en dos dimensiones Tengo mediciones de los dos componentes de gradiente ∂ u / ∂ x y ∂ u / ∂ y a lo largo de alguna parte del límite, Γ m , por lo que me gustaría imponer ∂ u- ∇2u = f( x , y) ,yo nΩ−∇2u=f(x,y),inΩ -\nabla^2 u = f(x,y),...
En el trabajo, esencialmente trabajo como consultor independiente. Para la administración y los clientes, necesito estimar la cantidad de tiempo que llevará desarrollar software como parte de mi investigación en ciencias computacionales. Sin embargo, mis estimaciones de tiempo generalmente están...
Tengo un conjunto de datos de puntos tridimensionales para los que me gustaría construir una malla, usando python. Todo el software que he visto requiere que proporciones los bordes. ¿Hay un programa en Python que toma como entrada un conjunto de puntos en 3D y genera una malla? Si es posible, me...
Comenzando con la ecuación de advección en forma de conservación. ut=(a(x)u)xut=(a(x)u)x u_t = (a(x)u)_x donde es una velocidad que depende del espacio, yu es una concentración de una especie que se conserva.a(x)a(x)a(x)uuu Discretizar el flujo (donde el flujo , se define en los bordes de las...
Actualmente estoy desarrollando un método de descomposición de dominio para la solución del problema de dispersión. Básicamente estoy resolviendo un sistema de Helmholtz BVP de forma iterativa. Discreto las ecuaciones usando el método de elementos finitos sobre mallas triangulares o tetraédricas....
12×1212×1212 \times 12QQQdet(Q)=det(12I−Q−J)(1)det(Q)=det(12I−Q−J)(1)\det(Q) = \det(12I-Q-J) \; \; (1)JJJ Actualmente estoy haciendo esto con la biblioteca de armadillo , pero resulta ser demasiado lento. La cuestión es que necesito hacer esto por un billón de matrices y resulta que calcular los...