Computational Science

Die Reproduzierbarkeit wird in der rechnerwissenschaftlichen Forschung immer wichtiger. (Siehe z. B. diesen Artikel von Roger Peng in Science ; mir sind auch andere Artikel und Websites dieser Art bekannt.) Es ist mir jedoch nicht klar, wie viele Informationen ich in einem Zeitschriftenartikel...

Ein Amazon EC2-Computing-Cluster kostet im Verlauf von 3 Jahren etwa 800 bis 1000 US- Dollar (abhängig vom Arbeitszyklus) pro physischem CPU-Kern. In unserer letzten Runde der Hardware-Akquisition hat mein Labor Hardware im Wert von 48 Kernen für ca. 300 US-Dollar pro Kern gefunden, die denen von...

Viele numerische CFD-Ansätze können auf beliebig hohe Ordnungen erweitert werden (z. B. diskontinuierliche Galerkin-Methoden, WENO-Methoden, Spektraldifferenzierung usw.). Wie soll ich für ein bestimmtes Problem eine geeignete Genauigkeitsreihenfolge

Wie ist der Stand der Technik bei der Näherung hochschwingender Integrale sowohl in einer Dimension als auch in höheren Dimensionen mit willkürlicher

Ich habe eine Reihe von Datenpunkten denen ich erwarte, dass sie (ungefähr) einer Funktion folgen , die zu einer Linie im großen asymptotisch ist . Im Wesentlichen nähert sich Null, wenn , und dasselbe kann wahrscheinlich von allen Ableitungen , gesagt werden. usw. Aber ich weiß nicht, wie die...

Diese Frage wurde von Cross Validated migriert, da sie über Computational Science Stack Exchange beantwortet werden kann. Vor 7 Jahren migriert . Ich möchte die MATLAB-Programmierung erlernen, damit ich einige Recherchen / Analysen selbst durchführen und auch einige MATLAB-Skripte studieren /...

Ich bin auf der Suche nach Büchern, Artikeln, Blog-Posts oder allgemein veröffentlichtem Material, das sich speziell mit der Verwendung moderner C ++ - Funktionen (Bewegungssemantik, STL, Iteratoren, verzögerte Auswertung usw.) im wissenschaftlichen Rechnen befasst. Kannst du welche...

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...

Ich habe Matrizen EINEINAGGGEINEINAn × nn×nn\times nnnnGGGn × mn×mn\times mmmm1 < m < 10001<m<10001 \lt m \lt 1000111000GTG = ichGTG=ichG^TG = IEINEINAA x = bEINx=bAx = bB i C G S t a b (l)BichCGSteinb(l)\mathrm{BiCGStab}(l)EIN- 1EIN-1A^{-1} Ich möchte ein System der Form lösen: (...

Aus heutiger Sicht passt diese Frage nicht zu unserem Q & A-Format. Wir erwarten, dass die Antworten durch Fakten, Referenzen oder Fachwissen gestützt werden, aber diese Frage wird wahrscheinlich Debatten, Argumente, Abstimmungen oder erweiterte Diskussionen hervorrufen. Wenn...

Was wäre eine gute Finite-Differenzen-Diskretisierung für die folgende Gleichung: ∂ρ∂t+∇⋅(ρu)=0∂ρ∂t+∇⋅(ρu)=0\frac{\partial \rho}{\partial t} + \nabla \cdot \left(\rho u\right)=0? Wir können den 1D Fall nehmen: ∂ρ∂t+ddx(ρu)=0∂ρ∂t+ddx(ρu)=0\frac{\partial \rho}{\partial t} + \frac{d}{dx}\left(\rho...

Gegeben ist eine zunehmende endliche Folge von positiven Zahlen . Welcher der beiden folgenden Algorithmen ist besser für die Berechnung der Summe der Zahlen?z1,z2,.....znz1,z2,.....znz_{1} ,z_{2},.....z_{n} s=0; for \ i=1:n s=s + z_{i} ; end Oder: s=0; for \ i=1:n s=s + z_{n-i+1} ;...

Was ist der bevorzugte und effiziente Ansatz für die Interpolation mehrdimensionaler Daten? Dinge, über die ich mir Sorgen mache: Leistung und Speicher für Konstruktion, Einzel- / Chargenbewertung Handhabungsabmessungen von 1 bis 6 linear oder höherer Ordnung Fähigkeit, Gradienten zu erhalten (wenn...

Unterrichten von Studenten im Grundstudium in einem Kurs mit dem Titel Einführung in die Computerprogrammierung. Ich bin ein bisschen verwirrt. In der Computerphysik verwenden Wissenschaftler C / C ++ oder Python oder Fortran, CUDA usw. ... dies ist die Zeit, um ihre Basis aufzubauen. Was soll ich...

Das OpenCL-Programmierparadigma verspricht, ein lizenzfreier Open-Standard für heterogenes Computing zu sein. Sollen wir unsere Zeit in die Entwicklung von OpenCL-basierter Software investieren? Für und

Wie wir wissen, sind orthogonale Transformationsmethoden (Givens-Rotationen und Housholder-Reflexionen) für lineare Gleichungssysteme teurer als die Gauß-Elimination, haben jedoch theoretisch bessere Stabilitätseigenschaften in dem Sinne, dass sie die Bedingungszahl des Systems nicht ändern. Obwohl...

Ich habe eine zweidimensionale Funktion deren Werte ich abtasten möchte. Die Berechnung der Funktion ist sehr aufwendig und weist eine komplexe Form auf. Ich muss daher einen Weg finden, um mit der geringsten Anzahl von Abtastpunkten die meisten Informationen über ihre Form zu erhalten.f( x ,...

Ich suche nach Algorithmen, um Standard-2D-Graphen für Funktionen zu zeichnen, die Singularitäten haben oder nicht. Der Zweck ist es, ein "Mini-CAS" zu schreiben, so dass ich keine A-priori-Kenntnisse über die Arten von Funktionen habe, die die Benutzer grafisch darstellen möchten. Dieses Problem...

Ich werde mit meiner persönlichen Erfahrung in unserem Labor beginnen. In den letzten 9 und 10 Tagen waren wir ziemlich aggressiv bei den Optimierungen. Wir kompilierten mit -O3 und prozessorspezifischen Flags (-xW -xSSE4.2 zum Beispiel). Ab ifort 11 bemerkten wir jedoch: 1. einige Inkonsistenzen...

Wie skalieren Python / Numpy-Arrays mit zunehmenden Array-Dimensionen? Dies basiert auf einem Verhalten, das ich beim Benchmarking von Python-Code für diese Frage festgestellt habe: Wie kann man diesen komplizierten Ausdruck mit numpy-Slices ausdrücken? Das Problem bestand hauptsächlich in der...