Me gustaría modelar el flujo laminar de agua desde las raíces hasta el tallo de una planta. Al final de las raíces, los tubos varían de milímetro a centímetro de diámetro y longitud. A medida que nos acercamos al tallo, las raíces se hacen más grandes en longitud y diámetro. Quiero crear dominios...
Muchos algoritmos utilizados en computación científica tienen una estructura inherente diferente a los algoritmos comúnmente considerados en formas menos intensivas en matemáticas de ingeniería de software. En particular, los algoritmos matemáticos individuales tienden a ser muy complejos, a menudo...
Estoy buscando una biblioteca que realice operaciones matriciales en grandes matrices dispersas sin sacrificar la estabilidad numérica. Las matrices serán 1000+ por 1000+ y los valores de la matriz estarán entre 0 y 1000. Realizaré el algoritmo de cálculo del índice, por lo que generaré vectores de...
¿Existe algún algoritmo que enumere las gráficas que corresponden a alguna teselación de puntos de Delaunay en 3D? Si es así, ¿hay una parametrización eficiente de las geometrías que corresponden a cualquier "gráfico de Delaunay"? Estoy buscando enumerar sistemáticamente todas las geometrías...
Estoy tratando de invertir numéricamente la siguiente transformación integral: F(y) = ∫∞0 0yExp[ - 12( y2+ x2) ] Yo0 0( x y) f( x )d xF(y)=∫0∞yexp[−12(y2+x2)]I0(xy)f(x)dxF(y) = \int_{0}^{\infty} y\exp{\left[-\frac{1}{2}(y^2 + x^2)\right]} I_0\left(xy\right)f(x)\;\mathrm{d}x Entonces, para una...
¿Existe un algoritmo SVD truncado que calcula los valores singulares uno a la vez? Mi problema: me gustaría calcular los primeros kkk valores singulares (y vectores singulares) de una gran matriz densa MMM , pero no sé cuál sería un valor apropiado de kkk . MMM es grande, por lo que por razones de...
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...
Tengo una pregunta que es similar a esta antes, excepto en 3D, y solo necesito el volumen, no la forma real del casco. Más precisamente, se me da un pequeño conjunto de puntos (digamos, 10-15) en 3D, todos los cuales se sabe que se encuentran en el casco convexo del conjunto de puntos (por lo que...
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...
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 trabajando en un problema que se puede paralelizar usando una sola operación mpi_allgather o una operación mpi_scatter y una mpi_gather. Estas operaciones se llaman dentro de un ciclo while, por lo que se pueden llamar muchas veces. En la implementación con un esquema MPI_allgather, estoy...
Estoy haciendo una diagonalización de Lanczos de una gran matriz dispersa (~ 2 millones de elementos). Casi todos los pasos del algoritmo de Lanzcos se realizan en paralelo en la GPU, excepto la diagonalización de la matriz de Lanczos para verificar la convergencia. Para eso, he estado usando el...
Pregunta ¿Cómo ordenaría una nube de puntos con respecto a una malla no estructurada de células hexaédricas? Cada celda tiene un centro y una etiqueta única para representarlo. Básicamente, hay dos puntos de nube (nube de puntos original y una nube de puntos de los centros celulares), pero la...
Dada una matriz simétrica definida positiva, ¿cuál es el algoritmo más rápido para calcular la matriz inversa y su determinante? Para los problemas que me interesan, la dimensión de la matriz es de 30 o menos. La alta precisión y velocidad es realmente necesaria. (se realizan millones de...
Cerrado. Esta pregunta está fuera de tema . Actualmente no está aceptando respuestas. ¿Quieres mejorar esta pregunta? Actualice la pregunta para que sea sobre el tema de Computational Science Stack Exchange. Cerrado hace 5 años . ¿Existe un grado de...
Tengo una matriz cuadrada simétrica real densa. La dimensión es de aproximadamente 1000x1000. Necesito calcular el primer componente principal y preguntarme cuál sería el mejor algoritmo para hacerlo. Parece que MATLAB usa los algoritmos Arnoldi / Lanczos (para eigs). Pero al leer sobre ellos, no...
¿Qué situaciones de análisis numérico se vuelven más / menos estables, tienen una convergencia más rápida / más lenta o son bastante diferentes cuando se trata de funciones de variable compleja en lugar de funciones de una variable real?
Me gustaría saber si hay un algoritmo que da un conjunto de puntos y un ángulo calcula el casco convexo si el ángulo es y dado un \ alpha> 0 calcula una envoltura que sigue más de cerca el "perímetro ".α=0α=0\alpha = 0α>0α>0\alpha > 0 Y si hay una definición de un perímetro no...
Estoy tratando de implementar el método Broyden-Fletcher-Goldfarb-Shanno para encontrar el mínimo de una función. Necesito dos conjeturas iniciales y y una aproximación inicial de la matriz de Hesse . Los únicos requisitos que encuentro para es que si el Hessian es simétrico positivo definido,...
La escalada parece ser una herramienta muy poderosa para la optimización. Sin embargo, cómo generar los "vecinos" de una solución siempre me desconcierta. Por ejemplo, estoy optimizando una solución . Aquí está dentro del rango , está dentro del rango , está dentro del rango . ¿Cuál es la mejor...