Computational Science

Fragen und Antworten für Wissenschaftler, die Computer verwenden, um wissenschaftliche Probleme zu lösen

1
Gibt es Open Source-ILU-Implementierungen mit mehreren Ebenen auf Invers-Basis?
Ich bin sehr beeindruckt von der Serienleistung von mehrstufigen inversen ILU-Vorkonditionierern , insbesondere für heterogenes Helmholtz , aber ich bin überrascht, keine Open-Source-Implementierungen finden zu können. Insbesondere stellt ILUPACK den Wissenschaftlern Binärdateien frei zur Verfügung, aber es scheint nicht, dass sie ihren Quellcode veröffentlichen. Ist es wirklich so, dass niemand …
3
Warum wird meine Matrix-Vektor-Multiplikationsskalierung nicht durchgeführt?
Entschuldigung für den langen Beitrag, aber ich wollte alles, was ich für relevant hielt, gleich mit einbeziehen. Was ich möchte Ich implementiere eine parallele Version der Krylov-Subraummethoden für dichte Matrizen. Hauptsächlich GMRES, QMR und CG. Mir wurde (nach der Profilerstellung) klar, dass meine DGEMV-Routine erbärmlich war. Also beschloss ich, mich …
1
Kann eine Krylov-Subraummethode als Glättungsfaktor für Mehrgitter verwendet werden?
Nach meinem Kenntnisstand verwenden Multigrid-Löser iterative Glätter wie Jacobi, Gauss-Seidel und SOR, um den Fehler bei verschiedenen Frequenzen zu dämpfen. Könnte stattdessen eine Krylov-Subraummethode (wie Konjugatgradient, GMRES usw.) verwendet werden? Ich glaube nicht, dass sie als "Glätteisen" klassifiziert sind, aber sie können verwendet werden, um die Grobgitterlösung anzunähern. Können wir …
4
Welche Anwendungen erfordern Intervallarithmetik?
Ich habe eine sehr grundlegende Vorstellung von Intervallarithmetik (IA), aber es scheint sowohl theoretisch als auch praktisch ein sehr interessanter Zweig der Computerwissenschaft zu sein. Es ist klar, dass es sich bei den offensichtlichen Anwendungen um verifizierte Computer- und unsachgemäße Probleme handelt, aber dies ist zu abstrakt. Da hier viele …
4
Wie kann man Variablen neu anordnen, um eine gebänderte Matrix mit minimaler Bandbreite zu erzeugen?
Ich versuche, eine 2D-Poisson-Gleichung durch endliche Differenzen zu lösen. Dabei erhalte ich eine spärliche Matrix mit nur Variablen in jeder Gleichung. Wenn die Variablen beispielsweise wären, würde die Diskretisierung ergeben:555UUU Ui−1,j+Ui+1,j−4Ui,j+Ui,j−1+Ui,j+1=fi,jUi−1,j+Ui+1,j−4Ui,j+Ui,j−1+Ui,j+1=fi,jU_{i-1,j} + U_{i+1,j} -4U_{i,j} + U_{i,j-1} + U_{i,j+1} = f_{i,j} Ich weiß, dass ich dieses System durch eine iterative Methode …
1
Wie kann die schwache Formulierung einer partiellen Differentialgleichung für die Finite-Elemente-Methode abgeleitet werden?
Ich habe eine grundlegende Einführung in die Finite-Elemente-Methode gegeben, in der ein differenziertes Verständnis einer „schwachen Formulierung“ nicht betont wurde. Ich verstehe, dass wir mit der Galerkin-Methode beide Seiten der (elliptischen) PDE mit einer Testfunktion multiplizieren und dann integrieren (nach Teilen oder nach dem Divergenzsatz). Manchmal musste ich Teile zweimal …
3
Was ist ein skalierbares Vorkonditionierer für Hochfrequenz Helmholtz?
Standard Mehrgitter- und Gebietszerlegungsmethoden nicht funktionieren, aber ich habe große 3D-Probleme und direkte Löser sind keine Option. Welche Methoden sollte ich versuchen? Wie werden meine Entscheidungen von folgenden Überlegungen betroffen? Koeffizienten variieren über mehrere Größenordnungen oder Finite-Elemente im Vergleich zu finite verschiedenen Methoden werden verwendet,
15 pde 
3
Wie wird B3LYP in Gaussin 0 *, GAMESS-US, Molpro usw. implementiert?
Insbesondere möchte ich die Arbeit mit B3LYP erweitern, die mit Gaussian 03 begann, aber mit GAMESS-US fortgesetzt wurde. Die von den Standard-B3LYP-Methoden bereitgestellten Energien sind nicht dieselben. Im Handbuch zu GAMESS-US (Abschnitt Weitere Informationen) wird darüber diskutiert: Beachten Sie, dass B3LYP in GAMESS teilweise auf der VWN5-Elektronengaskorrelationsfunktion basiert. Da es …
1
Wie sollen Randbedingungen bei Verwendung der Finite-Volumen-Methode angewendet werden?
Ausgehend von meiner vorherigen Frage versuche ich, Randbedingungen auf dieses ungleichmäßige endliche Volumen-Netz anzuwenden. Ich möchte eine Robin-Typ-Randbedingung auf die lhs der Domäne anwenden ( x=xL)x=xL)x=x_L) , so dass σL=(dux+au)∣∣∣x=xLσL=(dux+au)|x=xL \sigma_L = \left( d u_x + a u \right) \bigg|_{x=x_L} wobei σLσL\sigma_L der Grenzwert ist; a,da,da, d sind an der …
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.