Als «petsc» getaggte Fragen

PETSc ist eine Open-Source-Suite von Datenstrukturen und -routinen für die skalierbare (parallele) Lösung wissenschaftlicher Anwendungen, die durch partielle Differentialgleichungen modelliert werden.

17
Gibt es einen hochwertigen nichtlinearen Programmierlöser für Python?
Ich habe mehrere herausfordernde nicht konvexe globale Optimierungsprobleme zu lösen. Derzeit verwende ich die Optimization Toolbox von MATLAB (speziell fmincon()mit algorithm = 'sqp'), was sehr effektiv ist . Der größte Teil meines Codes ist jedoch in Python, und ich würde die Optimierung gerne auch in Python durchführen. Gibt es einen …
2
Warum konvergiert mein iterativer linearer Löser nicht?
Was kann schief gehen, wenn mit vorkonditionierten Krylov-Methoden von KSP ( dem linearen Lösungspaket von PETSc ) ein spärliches lineares System gelöst wird, wie es durch Diskretisierung und Linearisierung partieller Differentialgleichungen erhalten wird? Welche Schritte kann ich unternehmen, um festzustellen, was bei meinem Problem schief läuft? Welche Änderungen kann ich …
5
Was sind die Hauptunterschiede zwischen PETSc und Trilinos?
Soweit ich das beurteilen kann, handelt es sich bei den beiden großen generischen IT- Software-Frameworks des US-Energieministeriums um PETSc und Trilinos . Sie scheinen auf den ersten Blick ähnlich zu sein, abgesehen von Sprachunterschieden (C versus C ++). Was sind die Hauptunterschiede zwischen den beiden Frameworks und welche Faktoren sollten …
1
Warum konvergiert Newtons Methode nicht?
Ich verwende das nichtlineare Lösungspaket SNES von PETSc , um ein System nichtlinearer Gleichungen zu lösen, das durch Diskretisieren einer partiellen Differentialgleichung erhalten wird. Wie kann ich feststellen, warum der Löser nicht konvergiert, und was kann ich tun, um meine Gleichungen erfolgreich zu lösen?
2
Welche Bibliotheken unterstützen Multigrid auf hohem Niveau?
Ich plane, mit Multigrid einige Eigenwerte und Vektoren zu berechnen, und mir ist aufgefallen, dass PETSc Multigrid auf hohem Niveau unterstützt. Die PETSc-Dokumentation besagt, dass dieser Teil von PETSc nicht verwendet werden sollte, da er bald ersetzt wird. Welche anderen Bibliotheken unterstützen Multigrid auf hohem Niveau, und wie bald wird …
3
Euklidischer Abstand in Oktave
Ich würde gerne wissen, ob es einen schnellen Weg gibt, den euklidischen Abstand zweier Vektoren in Oktave zu berechnen. Es scheint, dass es dafür keine spezielle Funktion gibt. Soll ich also einfach die Formel mit verwenden sqrt?
2
Wie nützlich ist PETSc für dichte Matrizen?
Wo immer ich es gesehen habe, sagen PETSc-Tutorials / -Dokumente usw., dass es für die lineare Algebra nützlich ist und normalerweise angibt, dass spärliche Systeme davon profitieren. Was ist mit dichten Matrizen? Ich mache mir Sorgen, für dichtes A zu lösen .Ax=bAx=bAx=bAAA Ich habe meinen eigenen Code für CG und …
4
Gibt es Empfehlungen für Unit-Testing-Frameworks, die mit Code / Bibliotheken kompatibel sind, die MPI verwenden?
Normalerweise schreibe ich seriellen Code, und wenn ich das tue, schreibe ich Komponententests mit einem Testframework im xUnit-Stil (MATLAB xUnit, PyUnit / nose oder Googles C ++ - Testframework). Basierend auf einer flüchtigen Google-Suche habe ich nicht viel darüber gesehen, wie Praktiker Unit-Test-Code verwenden, der MPI verwendet. Gibt es dafür …
5
Berechnung der Sparsity-Struktur für Finite-Elemente-Matrizen
Frage: Mit welchen Methoden kann die Sparsity-Struktur einer Finite-Elemente-Matrix genau und effizient berechnet werden? Info: Ich arbeite an einem Poisson-Druckgleichungslöser nach der Methode von Galerkin auf quadratischer Lagrange-Basis, geschrieben in C, und verwende PETSc für die Speicherung von spärlicher Matrix und KSP-Routinen. Um PETSc effizient zu nutzen, muss der globalen …
2
Welche Vorkonditionierer (und Löser) in PETSc für unbestimmte symmetrische Systeme sollte ich verwenden?
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 Lagrange-Multiplikatoren zuletzt erscheinen). Das System ist ziemlich groß (+ 100k Zeilen).C= 0C=0C = …
3
Sparse Linear Solver für viele rechte Seiten
Ich muss das gleiche dünne lineare System (300x300 bis 1000x1000) mit vielen rechten Seiten (300 bis 1000) lösen. Zusätzlich zu diesem ersten Problem möchte ich auch verschiedene Systeme lösen, aber mit denselben Nicht-Null-Elementen (nur unterschiedlichen Werten), das sind viele spärliche Systeme mit konstantem Sparsity-Muster. Meine Matrizen sind unbestimmt. Die Leistung …
2
Wie man mit Funktionszeigern in Fortran in wissenschaftlichen Programmen arbeitet
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 aufgerufen. typedef PetscErrorCode (*TSIFunction)(TS,PetscReal,Vec,Vec,Vec,void*); PetscErrorCode …
11 petsc  fortran 
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.