Me estoy sumergiendo en el fascinante mundo del análisis de elementos finitos y me gustaría resolver un gran problema termomecánico (solo térmico mecánico, sin retroalimentación).→→\rightarrow Para el problema mecánico, ya entendí por la respuesta de Geoff , que necesitaré usar un solucionador...
Implementé un Navier Stokes basado en la discretización FEM y PETSc para resolver el sistema lineal de ecuaciones. Para crear un procedimiento de solución eficiente, sigo el documento "Precondicionamiento eficiente de las ecuaciones de Navier-Stokes linealizadas para flujo incompresible" (Silvester...
¿Cuál es la forma más eficiente de obtener una matriz dispersa compleja de mi código Fortran a PETSc? Entiendo que esto depende del problema, por lo que traté de dar tantos detalles relevantes como sea posible a continuación. He estado jugando con el solucionador de valores propios FEAST [1] para...
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...
Tengo un problema de matriz generalizado: de un método espectral en un problema de análisis de estabilidad lineal. Mi matriz B es diagonal y positiva semi-definida. A es no ermitaño y complejo.A x = λ B xUNAX=λsiXA x = \lambda B x Mi problema es esencialmente que cuando uso el solucionador de...
Estoy perfilando el rendimiento de los solucionadores lineales de PETSc. Según lo entiendo, speedup=Sequential TimeParallel Time.speedup=Sequential TimeParallel Time.\text{speedup}=\frac{\text{Sequential Time}}{\text{Parallel Time}}. Sé que ejecutar el código paralelo en un procesador se puede...
¿Cómo puedo calcular el complemento Schur? S= Kb b- Kb aK- 1a aKa bS=Ksisi-KsiunaKunauna-1Kunasi S = K_{bb} - K_{ba} K_{aa}^{-1} K_{ab} dónde K= ( Ka aKb aKa bKb b)K=(KunaunaKunasiKsiunaKsisi) K=\begin{pmatrix} K_{aa} & K_{ab} \\ K_{ba} & K_{bb} \end{pmatrix} (en algún orden) es una matriz...