Estaba tratando de ejecutar casos de prueba para CG y necesito generar: matrices definidas positivas simétricas de tamaño> 10,000 DENSIDAD COMPLETA Usando solo índices de matriz y, si es necesario, 1 vector (como )A ( i , j )= x ( i ) - x ( j )( i + j )A(i,j)=x(i)−x(j)(i+j)A(i,j) =...
Necesito calcular solo un pequeño número de componentes de Fourier de baja frecuencia de una matriz compleja de 2 dimensiones. Calcularé los mismos componentes de Fourier una y otra vez a medida que cambie la matriz de entrada. Claramente, en el límite donde solo quiero un componente de Fourier,...
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 tratando de hacer MD en vapor de agua. Como sé, existen algunos modelos de agua para agua líquida, como SPC, SPC / E, TIP3P, pero ¿se aplicarán también al estado de vapor del agua? ¿Y cuál es la diferencia de las simulaciones de agua líquida y
Parece que Matlab eigsme está dando malas aproximaciones de los vectores propios más pequeños de una matriz. Supongo que puedo usar algunos métodos más lentos que también serían más precisos ... Estoy buscando encontrar el segundo vector propio más pequeño de una matriz lapalciana (conocido como...
Recibo la advertencia en el asunto de la publicación cuando intento optimizar una función en Python con la función scipy.optimize.fmin_bfgs . La salida completa: Advertencia: el error deseado no se logra necesariamente debido a la pérdida de precisión Current function value: nan Iterations:...
Según tengo entendido, el orden en que se reciben los mensajes MPI punto a punto sin bloqueo (Isend e Irecv) es coherente con el orden en que se envían. ¿Hay alguna técnica para dar prioridad a ciertos mensajes sobre otros? Por ejemplo, tenemos un algoritmo multinivel en el que las soluciones de...
Tengo una matriz hermitiana que depende de dos parámetros, digamos x e y . Cuando lo diagonalizo en dos puntos cercanos ( x 1 , y 1 ) y ( x 2 , y 2 ) obtengo dos valores propios cercanos ( ε 1 y ε 2 ) y dos espacios propios correspondientes ( S 1 y S 2 ) de la misma dimensión .HH\mathbf{H}Xxxyyy(...
Cuando pienso en resolver ecuaciones no lineales, generalmente pienso en linealizar primero, luego aplicar un preacondicionador a la matriz lineal. Se me ocurrió pensar que podría ser posible condicionar las ecuaciones no lineales ANTES de que se linealicen. ¿Existen estrategias generales para...
Cualquiera conoce referencias que resumen la historia y las ideas detrás de la supercomputación, incluida la mención de desarrollos en lenguajes de programación paralelos, aplicaciones, empresas de inicio (algunas se establecieron pero no sobrevivieron), necesidades de la industria,...
Mi desafío es resolver el siguiente sistema de ecuaciones, que describe la combustión de gas en medios porosos: 1) continuidad ε∂ρg∂t+∂∂x(ρgux)=0ε∂ρg∂t+∂∂x(ρgux)=0\varepsilon \frac{\partial \rho_g}{\partial t} +\frac{\partial}{\partial x} \left(\rho_g u_x\right)=0 2) Ley de Darcy...
cuáles serían buenos métodos (y / o paquetes de software) para intentar resolver un problema minimizando una función cuadrática , st 0 ≤ x i ≤ 1 , y hay más restricciones, algunas de las cuales son no lineales (y no diferenciables), por ejemplo, ∑ i x i 1 x i > a < b ?F( x ) = ∑nortei = 1(...
Aprendí que para el solucionador BiCGStab para sistemas lineales dispersos es casi siempre necesario usar un preacondicionador. Ahora me di cuenta de que elegir una buena depende del problema. Navegando por la web descubrí que hay muchos preacondicionadores basados en ILU (ILUT, MILU, etc.) Aquí...
Los diferentes métodos de álgebra lineal numérica y de optimización numérica tienen diferentes regímenes de tamaño donde son una 'buena idea', además de sus propias propiedades. Por ejemplo, para problemas de optimización muy grandes, se utilizan los métodos de gradiente, gradiente estocástico y...
Quiero resolver la ecuación isotérmica de Lane-Emden [PDF, eq. 15.2.9] d2ψdξ2+2ξdψdξ=e−ψd2ψdξ2+2ξdψdξ=e−ψ\frac{d^2 \!\psi}{d \xi^2} + \frac{2}{\xi} \frac{d \psi}{d \xi} = e^{-\psi} con las condiciones iniciales ψ(ξ=0)=0dψdξ∣∣∣ξ=0=0ψ(ξ=0)=0dψdξ|ξ=0=0\psi(\xi = 0) = 0 \quad \left.\frac{d\psi}{d...
Soy nuevo en el campo de CFD. ¿Cuándo se debe optar por una cuadrícula estructurada y cuándo se debe optar por una no estructurada? (Sí, depende mucho de la geometría del problema) Más específicamente, quiero saber la diferencia en la potencia de cálculo requerida, la precisión lograda y los...
Más o menos la pregunta. Dada una matriz general dispersa, no simétrica (tanto numérica como estructuralmente), ¿qué importancia tiene el patrón de dispersión (es decir, permutación de fila / columna de matriz / vector) para los solucionadores iterativos? Puedo ver que se vuelve importante para los...
Me enfrento con el problema de encontrar el elipsoide ( es una matriz definida positiva simétrica) de volumen máximo dentro de un conjunto convexo dado como un conjunto de desigualdades lineales . Comprendí cómo se formaliza como un problema de optimización convexa como se da en "Convex...
Soy nuevo en programación lineal y he formulado un programa lineal (LP) con un orden de variables y restricciones, aunque la matriz de restricciones es extremadamente escasa. 10 131013101310^{13}1013101310^{13} ¿Quería saber si un LP de esta escala es manejable o
Estoy interesado en AMR local en mallas no estructuradas. Actualmente, estoy trabajando con la biblioteca OpenFOAM: admite AMR local completamente no estructurado: los criterios de refinamiento de celda determinan una lista de celdas que se cortan las celdas seleccionadas se refinan: se...