Ich tauche in die faszinierende Welt der Finite-Elemente-Analyse ein und möchte ein großes thermomechanisches Problem lösen (nur thermische Pfeilmechanik, kein Feedback).→→\rightarrow Für das mechanische Problem habe ich bereits aus Geoffs Antwort verstanden , dass ich aufgrund der Größe meines...
Ich habe einen Navier Stokes implementiert, der auf FEM-Diskretisierung und PETSc basiert, um das lineare Gleichungssystem zu lösen. Um ein effizientes Lösungsverfahren zu erstellen, folge ich dem Artikel "Effiziente Vorkonditionierung der linearisierten Navier-Stokes-Gleichungen für...
Was ist der effizienteste Weg, um eine komplexe Matrix mit geringer Dichte von meinem Fortran-Code auf PETSc zu übertragen? Ich verstehe, dass dies problemabhängig ist, deshalb habe ich versucht, so viele relevante Details wie möglich unten anzugeben. Ich habe mit dem FEAST-Eigenwertlöser [1] für...
So ziemlich die Frage. Wie wichtig ist bei einer allgemein spärlichen, nicht symmetrischen (sowohl numerisch als auch strukturell) Matrix das Sparsity-Muster (dh die Zeilen- / Spaltenpermutation von Matrix / Vektor) für iterative Löser? Ich kann sehen, dass es für Direktlöser (LU) oder...
Ich habe ein verallgemeinertes Matrixproblem: aus einer Spektralmethode für ein Problem der linearen Stabilitätsanalyse. Meine Matrix B ist diagonal und positiv semidefinit. A ist nicht hermitisch und komplex.Ax=λBxAx=λBxA x = \lambda B x Mein Problem ist im Wesentlichen, dass bei Verwendung des...
Ich stelle die Leistung der linearen Löser von PETSc vor. So wie ich es verstehe, speedup=Sequential TimeParallel Time.speedup=Sequential TimeParallel Time.\text{speedup}=\frac{\text{Sequential Time}}{\text{Parallel Time}}. Ich weiß, dass das Ausführen des parallelen Codes auf einem Prozessor als...
Wie kann ich das Schur-Komplement berechnen: S.= K.b b- K.b aK.- 1a aK.a bS.=K.bb- -K.beinK.einein- -1K.einb S = K_{bb} - K_{ba} K_{aa}^{-1} K_{ab} wo K.= ( K.a aK.b aK.a bK.b b)K.=(K.eineinK.einbK.beinK.bb) K=\begin{pmatrix} K_{aa} & K_{ab} \\ K_{ba} & K_{bb} \end{pmatrix} (in einiger...