Als «sparse-matrix» getaggte Fragen

Fragen zu Speicherung, Montage, Operationen und anderen Aspekten des Umgangs mit spärlichen Matrizen, für die nur Elemente ungleich Null gespeichert sind. Fragen, die sich nicht direkt auf spärliche Matrizen beziehen, sondern andere Mittel zur Verwendung von Sparsity, sollten mit [sparse-operator] gekennzeichnet werden.

5
Beste Wahl des Lösers für ein großes, spärlich symmetrisches (aber nicht positiv definiertes) System
Ich arbeite derzeit an der Lösung sehr großer symmetrischer (aber nicht positiv definierter) Systeme, die von bestimmten Algorithmen generiert werden. Diese Matrizen haben eine schöne Blocksparsity, die zum parallelen Lösen verwendet werden kann. Ich kann mich jedoch nicht entscheiden, ob ich einen direkten Ansatz (wie Multi-Frontal) oder einen iterativen Ansatz …
2
Wie man Dirichlet-Randbedingungen effizient in globalen Finite-Elemente-Steifheitsmatrizen mit geringer Dichte implementiert
Ich frage mich, wie Dirichlet-Randbedingungen in globalen Finite-Elemente-Matrizen mit geringer Dichte tatsächlich effizient implementiert werden. Nehmen wir zum Beispiel an, unsere globale Finite-Elemente-Matrix war: K.= ⎡⎣⎢⎢⎢⎢⎢⎢520- 102410001632- 1037000203⎤⎦⎥⎥⎥⎥⎥⎥und rechter Vektorb = ⎡⎣⎢⎢⎢⎢⎢⎢b 1b 2b 3b 4b 5⎤⎦⎥⎥⎥⎥⎥⎥K.=[520- -102410001632- -1037000203]]und rechter Vektorb=[b1b2b3b4b5]]K = \begin{bmatrix} 5 & 2 & 0 & -1 …
2
Welche neuartigen Datenstrukturen werden in der adaptiven FEM verwendet?
Viele adaptive FEM-Bibliotheken verwenden erweiterte Netzdatenstrukturen, um das Hinzufügen / Entfernen von Knoten, Kanten, Dreiecken, Tetraedern usw. zu handhaben. Beispielsweise verwendet die p4est- Bibliothek Octree-Datenstrukturen zur adaptiven Netzverfeinerung . Sie würden nicht oft Oktrees finden, die für Berechnungen in einem statischen Netz verwendet werden. Was ändert sich auf der Seite …
3
Faustregel für spärliche oder dichte Matrixspeicherung
Angenommen, ich kenne die erwartete Sparsamkeit einer Matrix (dh die Anzahl der Nicht-Nullen / die insgesamt mögliche Anzahl der Nicht-Nullen). Gibt es eine Faustregel (möglicherweise eine ungefähre Angabe) für die Entscheidung, ob ein spärlicher Matrixspeicher (insbesondere ein komprimierter Zeilenspeicher) verwendet oder als dichte Matrix gespeichert werden soll? Geschwindigkeit ist in …
2
Berechnen von Eigenvektorkomponenten eines gegebenen Vektors
Ich habe einen Vektor VVV der in den Eigenraum des hermitischen Sparse-Operators zerlegt werden kann :MMM V=∑ivim^iV=∑ivim^iV = \sum_i v_i \hat{m}_i Gibt es eine Möglichkeit, den (den Eigenvektor selbst) zu finden, der dem größten (in der Größe) entspricht?m^im^i\hat{m}_iviviv_i Ich möchte im Wesentlichen die größten paar Terme der Summe, einschließlich der …
4
Sparse Matrix Implementierung des Kalman Filters?
Ich habe einen Kalman-Filter-basierten Modellierungscode, den ich für eine regionale ionosphärische Kartierungsanwendung in Echtzeit entwickelt habe. Der Code assimiliert Daten von verschiedenen Sensoren mithilfe eines Kalman-Filters in eine Karte (beschrieben durch eine Reihe von Basisfunktionen). Ich versuche, dies auf einen größeren Bereich und mehr Sensoren zu skalieren, aber der Matrixalgebra-Teil …
2
Taxonomie von ILU-Vorkonditionierern
Ich habe gelernt, dass es für den BiCGStab-Solver für spärliche lineare Systeme so gut wie immer erforderlich ist, einen Vorkonditionierer zu verwenden. Inzwischen wurde mir klar, dass die Auswahl eines guten problemabhängig ist. Beim Surfen im Internet habe ich herausgefunden, dass es viele ILU-basierte Vorkonditionierer gibt (ILUT, MILU usw.). Hier …
Durch die Nutzung unserer Website bestätigen Sie, dass Sie unsere Cookie-Richtlinie und Datenschutzrichtlinie gelesen und verstanden haben.
Licensed under cc by-sa 3.0 with attribution required.