Ciencia computacional

9
¿Hay algún paquete FEM "ligero"?

Básicamente, FEM parece ser un problema que está prácticamente "resuelto". Existen numerosos marcos potentes, como Trilinos, PETSc, FEniCS, Libmesh o MOOSE. Una cosa que tienen en común: son extremadamente "pesados". Primero, la instalación normalmente es muy dolorosa. En segundo lugar, su...

9
¿Cuál es la razón por la que LAPACK usa

La rutina QR de LAPACK almacena Q como reflectores de cabeza de familia. Escala el vector de reflexión vvv con 1/v11/v11/v_1 , por lo que el primer elemento del resultado se convierte en 111 , por lo que no tiene que almacenarse. Y almacena un vector ττ\tau separado , que contiene los factores de...

8
Prolongar el trabajo de PBS

Es bastante doloroso descubrir que un trabajo de varios días se eliminará prematuramente debido a un error al establecer el límite de tiempo de pared. ¿Hay alguna manera de cambiarlo para un trabajo PBS en ejecución

8
¿Cómo puedo calcular el complemento Schur en PETSc?

¿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...