Ich habe begonnen, Python als Programmiersprache für alle meine Aufgaben in CFD zu verwenden. Ich habe sehr wenig Erfahrung in der Programmierung. Ich komme aus dem Maschinenbau und studiere Luft- und Raumfahrttechnik. Manchmal wird der Rechenaspekt von CFD langwieriger als das Manipulieren der Gleichungen oder das Berechnen. Welche allgemeinen Richtlinien …
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 …
Ich löse ein System von zwei gekoppelten PDEs in zwei räumlichen Dimensionen und in der Zeit rechnerisch. Da die Funktionsauswertungen teuer sind, würde ich gerne eine mehrstufige Methode verwenden (initialisiert mit Runge-Kutta 4-5). Die Adams-Bashforth-Methode unter Verwendung von fünf vorherigen Funktionsbewertungen hat einen globalen Fehler von (dies ist der Fall, …
Ich habe einen Datensatz, der in Millionen von Datenpunkten in 3D läuft. Für die Berechnung, die ich mache, muss ich Nachbarn (Entfernungssuche) für jeden Datenpunkt in einem Radius berechnen, versuchen, eine Funktion anzupassen, den Fehler für die Anpassung berechnen, dies für den nächsten Datenpunkt wiederholen und so weiter. Mein Code …
Viele Algorithmen, die beim wissenschaftlichen Rechnen verwendet werden, haben eine andere inhärente Struktur als Algorithmen, die üblicherweise in weniger mathematisch intensiven Formen des Software-Engineerings berücksichtigt werden. Insbesondere einzelne mathematische Algorithmen sind in der Regel hochkomplex und umfassen häufig Hunderte oder Tausende von Codezeilen, enthalten jedoch keinen Zustand (dh sie wirken …
Bei Anwendung der klassischen Formel für den Winkel zwischen zwei Vektoren: α=arccosv1⋅v2∥v1∥∥v2∥α=arccosv1⋅v2‖v1‖‖v2‖\alpha = \arccos \frac{\mathbf{v_1} \cdot \mathbf{v_2}}{\|\mathbf{v_1}\| \|\mathbf{v_2}\|} man stellt fest, dass bei sehr kleinen / spitzen Winkeln ein Präzisionsverlust auftritt und das Ergebnis nicht genau ist. Wie in dieser Antwort zum Stapelüberlauf erläutert , besteht eine Lösung darin, stattdessen …
Ich bin auf der Suche nach einer Bibliothek, die Matrixoperationen mit großen, spärlichen Matrizen ausführt, ohne die numerische Stabilität zu beeinträchtigen. Matrizen werden 1000+ mal 1000+ sein und die Werte der Matrix werden zwischen 0 und 1000 liegen. Ich werde den Indexberechnungsalgorithmus ausführen , so dass ich (spärliche) Zeilenvektoren der …
Gibt es einen Algorithmus, der die Graphen auflistet, die einer Delaunay-Tessellation von Punkten in 3D entsprechen? Wenn ja, gibt es eine effiziente Parametrisierung von Geometrien, die einem "Delaunay-Graphen" entsprechen? Ich versuche, systematisch alle stabilen Geometrien von Molekülen einer bestimmten Zusammensetzung ohne vorherige Kenntnis der Bindung usw. aufzulisten. BEARBEITEN: Sei die …
Ich versuche, die folgende integrale Transformation numerisch zu invertieren: F(y)=∫∞0yexp[−12(y2+x2)]I0(xy)f(x)dxF(y)=∫0∞yexp[−12(y2+x2)]I0(xy)f(x)dxF(y) = \int_{0}^{\infty} y\exp{\left[-\frac{1}{2}(y^2 + x^2)\right]} I_0\left(xy\right)f(x)\;\mathrm{d}x Für ein gegebenes ich also approximieren, wobei:f ( x )F(y)F(y)F(y)f(x)f(x)f(x) f(x)f(x)f(x) und sind real und positivF(y)F(y)F(y) (sie sind kontinuierliche Wahrscheinlichkeitsverteilungen) x,yx,yx,y sind real und positiv (sie sind Größen) Ich habe momentan eine sehr unordentliche …
Der Satz der Cholesky-Zerlegung besagt also, dass jede reelle symmetrische positiv-definitive Matrix eine Cholesky-Zerlegung M = L L ⊤ hat, wobei L eine untere Dreiecksmatrix ist.MMMM=LL⊤M=LL⊤M= LL^\topLLL Angesichts von wissen wir bereits, dass es schnelle Algorithmen gibt, um den Cholesky-Faktor L zu berechnen .MMMLLL Angenommen, ich habe eine rechteckige Matrix …
Hallo Scicomp-Community, Ich habe im Bereich der Graph-Algorithmen mit Frameworks wie NetworkX (Python), JUNG und YFiles (Java) gearbeitet. Ich betrete jetzt den Bereich des parallelen und hochleistungsfähigen Rechnens. Für ein neues Projekt suche ich eine C ++ - Grafikbibliothek mit folgenden Funktionen: verfügt über eine intuitive Benutzeroberfläche, die die Entwicklung …
Dies ist ein einfaches Problem der Symmetrieaufzählung. Ich gebe hier den vollständigen Hintergrund an, aber es sind keine Kenntnisse der Quantenchemie erforderlich. Die beiden Teilchen integral ist: ⟨ i j | k l ⟩ = ∫ ψ * i ( x ) ψ * j ( x ' ) ψ …
Ich arbeite an einem Problem, das mithilfe einer einzelnen mpi_allgather-Operation oder einer mpi_scatter- und einer mpi_gather-Operation parallelisiert werden kann. Diese Operationen werden innerhalb einer while-Schleife aufgerufen, sodass sie häufig aufgerufen werden können. In der Implementierung mit einem MPI_allgather-Schema sammle ich einen verteilten Vektor auf alle Prozesse zur Lösung doppelter Matrix. …
Ich mache eine Lanczos-Diagonalisierung einer großen Matrix mit geringer Dichte (~ 2 Millionen Elemente). Fast alle Schritte im Lanzcos-Algorithmus werden parallel auf der GPU ausgeführt, mit Ausnahme der Diagonalisierung der Lanczos-Matrix, um die Konvergenz zu überprüfen. Dafür habe ich den TQLI-Algorithmus von Numerical Recipes verwendet. Gibt es Methoden, um das …
Frage Wie würden Sie eine Punktwolke in Bezug auf ein unstrukturiertes Netz hexaedrischer Zellen sortieren? Jede Zelle hat ein Zentrum und eine eindeutige Bezeichnung, um sie darzustellen. Grundsätzlich gibt es zwei Trübungspunkte (ursprüngliche Punktwolke und eine Punktwolke der Zellzentren), aber die Informationen zur Zellgeometrie (Begrenzungsrahmen) können von Nutzen sein, da …
We use cookies and other tracking technologies to improve your browsing experience on our website,
to show you personalized content and targeted ads, to analyze our website traffic,
and to understand where our visitors are coming from.
By continuing, you consent to our use of cookies and other tracking technologies and
affirm you're at least 16 years old or have consent from a parent or guardian.