Mein System ist ein symmetrisches FE-Problem mit Lagrange-Multiplikatoren (z. B. inkompressibler Stokes-Fluss): ( ABBTC)(EINBTBC)\begin{pmatrix}A & B^T \\ B & C\end{pmatrix} Dabei ist der typische Fall (ich habe sogar darauf geachtet, dass die Gleichungen nummeriert sind, damit die...
Parallele Berechnungen werden häufig unter Verwendung einer deterministischen lokalen Berechnungsrate, des Latenzaufwands und der Netzwerkbandbreite modelliert. In Wirklichkeit sind diese räumlich variabel und nicht deterministisch. Studien wie Skinner und Kramer (2005) beobachten multimodale...
Wie viel langsamer ist Haustierec4py vs c / c ++ / fortran? Mir ist klar, dass dies erheblich vom ausgeführten Code abhängt, aber was ist mit etwas Einfachem wie einem
Ich bekomme eine Matrix , die symmetrisch, invertierbar, positiv definitiv und dicht ist. Ich muss testen, ob wobei J die All-One-Matrix ist.Q det ( Q ) = det ( 12 I - Q - J )12×1212×1212 \times 12QQQdet(Q)=det(12I−Q−J)(1)det(Q)=det(12I−Q−J)(1)\det(Q) = \det(12I-Q-J) \; \; (1)JJJ Ich mache das...
Hier ist eine typische Verwendung von Funktionszeigern in C. Ich möchte etwas Ähnliches in Fortran tun. Ich habe einige Ideen, aber ich würde gerne wissen, ob es einen kanonischen Weg gibt, dies zu tun. Die vom Benutzer übergebenen Funktionszeiger und Kontexte werden gespeichert und später...
Ich möchte die Runge-Kutta-Methode 8. Ordnung (89) in einer in C ++ geschriebenen Anwendung für Himmelsmechanik / Astrodynamik unter Verwendung eines Windows-Computers verwenden. Daher frage ich mich, ob jemand eine gute Bibliothek / Implementierung kennt, die dokumentiert und kostenlos zu...
Ich bin neu im PETSc-Paket. Ich habe eine ~ 4000x4000 Matrix A im Matrix-Market-Format und möchte PETSc dazu bringen, dies mit mehreren Prozessoren zu lösen. Ich weiß, wie man das System auf einem einzelnen Prozessor löst, aber ich weiß nicht, wie man die Matrix und die Vektoren auf...
Was wird schneller berechnet, oder oder ? , und sind positive Realzahlen mit .log a c b √einbaba^bLogeinclogac\log_a c abcb>1c√bcb\sqrt[b]{c}einaabbbcccb > 1b>1b>1 Welche Arten von Algorithmen werden Sie im Vergleich verwenden? Was sind ihre Komplexitäten? Zum Beispiel, wenn oder c ≈ a bc...
Ich habe ein großes Modell (~ 5000 Zeilen) in C geschrieben. Es ist ein serielles Programm, bei dem nirgendwo Zufallszahlen generiert werden. Es nutzt die FFTW-Bibliothek für Funktionen, die FFT verwenden. Ich kenne die Details der FFTW-Implementierung nicht, gehe jedoch davon aus, dass die darin...
In den letzten Jahren sind mehrere Bibliotheken / Softwareprojekte erschienen, die die eine oder andere Form der datengesteuerten Parallelität für gemeinsam genutzten Speicher für allgemeine Zwecke bieten. Die Hauptidee ist, dass Programmierer ihre Algorithmen nicht als explizit getunten Code...
Ich bin sehr neu in der GPGPU-Programmierung. Bitte verzeihen Sie mir, wenn die Frage nicht besonders geeignet ist. Soweit ich weiß, ist die GPU-Programmierung im Vergleich zur üblichen CPU-Programmierung eine sehr komplizierte technische Arbeit. Man muss sehr vorsichtig sein, wenn es um...
Es gibt viele bekannte numerische Methoden zum Lösen von Gleichungen vom Typ z. B. Halbierungsmethode, Newtonsche Methode usw.f(x)=0,x∈Rn,f(x)=0,x∈Rn, f(x) = 0, \quad x \in \mathbb{R}^n, In meiner Anwendung wird mit einer stochastischen Methode berechnet (das Ergebnis ist ein...
Beeinträchtigt das Kompilieren von PETSc mit einer externen BLAS / LAPACK-Bibliothek die Leistung auf dünn besetzten Matrizen erheblich oder werden diese Bibliotheken nur für die Mathematik mit dichter Matrix
Wie aus Jed Browns Antwort auf die Kosten für Suchvorgänge im Vergleich zu Berechnungen hervorgeht , führt die Verwendung von vektorisierten oder nicht vektorisierten Gleitkommaoperationen zu einem viel schnelleren Code. Viele moderne Compiler behaupten, dass sie eine automatische Vektorisierung...
Ich muss Ax = b lösen, aber mir ist klar, dass das Speichern der Matrixkoeffizienten meines Problems zu viel Speicher benötigt, selbst wenn es spärlich ist. Jetzt überlege ich mir, eine matrixfreie Methode zu verwenden, da die gleichen Koeffizienten viel Zeit in der Matrix erscheinen, sodass ich...
Der größte Teil meiner Programmierung besteht aus einmaligen Forschungscodes in C für meinen eigenen Gebrauch. Ich habe niemals Code an andere als enge Mitarbeiter verteilt. Ich habe einen Algorithmus entwickelt, den ich in einer wissenschaftlichen Zeitschrift veröffentliche. Ich möchte den...
Ich suche nach glaubwürdigen Referenzen, die angeben, wie viel Ressourcen Supercomputer für die Koordination im Vergleich zur eigentlichen aufgabenbezogenen Arbeit ausgeben. Ressourcen könnten verfügbare Rechenleistung sein, aber selbst Watt scheint eine gültige Einheit zu sein. Ich glaube, einer...
Grundsätzlich scheint FEM ein Problem zu sein, das so ziemlich "gelöst" ist. Es gibt zahlreiche leistungsstarke Frameworks wie Trilinos, PETSc, FEniCS, Libmesh oder MOOSE. Eines haben sie gemeinsam: Sie sind extrem "schwer". Erstens ist die Installation normalerweise sehr schmerzhaft. Zweitens ist...
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...
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...