Angenommen, ich habe zwei Matrizen Nx2, Mx2, die jeweils N, M 2d Vektoren darstellen. Gibt es eine einfache und gute Möglichkeit, die Abstände zwischen den einzelnen Vektorpaaren (n, m) zu berechnen? Der einfache, aber ineffiziente Weg ist natürlich: d = zeros(N, M); for i = 1:N, for j = 1:M,...
Ich versuche herauszufinden, ob es einen schnelleren Weg gibt, alle Eigenwerte und Eigenvektoren einer sehr großen und spärlichen Adjazenzmatrix zu berechnen, als mit scipy.sparse.linalg.eigsh Soweit ich weiß, verwendet diese Methode nur die Spärlichkeit und Symmetrieattribute der Matrix. Eine...
Sowohl bei der Domänenzerlegung (DD) als auch bei der Multigrid-Methode (MG) kann man die Anwendung der Blockaktualisierungen oder der Grobkorrekturen entweder additiv oder multiplikativ zusammenstellen . Für Punktlöser ist dies der Unterschied zwischen der Jacobi- und der Gauß-Seidel-Iteration....
Betrachten Sie zwei Computer mit unterschiedlichen Hardware- und Softwarekonfigurationen. Wenn auf jeder Plattform genau derselbe serielle Navier-Stokes-Code ausgeführt wird, dauert es x- und y-Zeit, um eine Iteration für Computer 1 bzw. 2 auszuführen. In diesem Fall ist die Iterationszeitdifferenz...
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...
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 ' ) ψ k ( x ) ψ L ( x ' )⟨ij|kl⟩⟨ij|kl⟩\langle ij|kl\rangle...
Ich suche nach der wichtigsten Literatur und Folienreferenzen zur Modellierung der aktuellen und zukünftigen Energiekosten von Gleitkommaoperationen und Datenübertragungen über CPU, Speicher, Netzwerk und Speicher. Ich habe diese Frage als Community-Wiki markiert und möchte, dass Sie jede Antwort...
Skaliert die Matrixmultiplikation (sowohl Mat * Mat als auch Mat * Vec) mit der Anzahl der Nicht-Nullen oder mit der Größe der Matrix? Oder eine Kombination aus beiden. Was ist mit der Form? Zum Beispiel habe ich eine 100 x 100-Matrix mit 100 Werten oder eine 1000 x 1000-Matrix mit 100 Werten. Wird...
Ich muss eine Koordinatentransformation zwischen zwei Referenzsystemen (Achsen) durchführen. Dafür müssen drei Matrizen ( ) multipliziert werden, da einige Zwischenachsen verwendet werden. Ich habe über zwei Lösungsansätze nachgedacht:3 × 33×33\times3 Methode 1 : Die Multiplikation direkt...
C ++ 11 führt eine Verschiebungssemantik ein, die beispielsweise die Codeleistung in Situationen verbessern kann, in denen C ++ 03 eine Kopierkonstruktion oder eine Kopierzuweisung durchführen müsste. Dieser Artikel berichtet, dass der folgende Code beim Kompilieren mit C + 11 eine 5-fache...
Ich versuche den Unterschied zwischen diesen beiden Grafikkarten für das akademische Rechnen zu verstehen, speziell für die DGEMM-Komponente. Wenn wir uns die Rohstatistik ansehen, haben beide den gleichen GK110-Chip, vergleichbare Statistiken in praktisch jeder Kategorie und, glaube ich, die...
Ich habe ein Problem, bei dem ich alle positiven (wie im Eigenwert positiv ist) Eigenpaare einer kleinen (normalerweise kleiner als 60x60) unsymmetrischen Matrix finden muss. Ich kann aufhören zu berechnen, wenn der Eigenwert kleiner als ein bestimmter Schwellenwert ist. Ich weiß, dass die...
Der Titel sagt das meiste davon. Ich suche eine leichte und benutzerfreundliche Bibliothek, die ich für Android (NDK) -Projekte verwenden kann. Für dichtes Material verwende ich gerne Eigen, aber ich habe nicht viele umfassende (und dokumentierte!) Bibliotheken für spärliches Material gefunden, die...
double trap(double func(double), double b, double a, double N) { double j; double s; double h = (b-a)/(N-1.0); //Width of trapezia double func1 = func(a); double func2; for (s=0,j=a;j<b;j+=h){ func2 = func(j+h); s = s + 0.5*(func1+func2)*h; func1 = func2; } return s; } Das Obige...
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...
Ich habe viele Veröffentlichungen in Fachzeitschriften für Computerphysik gesehen, die unterschiedliche Metriken für die Leistung ihres Codes verwenden. Insbesondere für GPGPU-Code scheint es eine Vielzahl von Timing-Ergebnissen zu geben, die von Menschen veröffentlicht werden. Insbesondere habe...
Es ist allgemein üblich zu sagen, dass Linux aus guten Gründen schneller ist. Aber wie wichtig ist die Auswahl des Betriebssystems für die Leistung des wissenschaftlichen Computercodes, wie im Titel angegeben? Bei Dingen, an denen das Betriebssystem direkter beteiligt ist (z. B. Spielen oder...
Ich habe CG in FORTRAN implementiert, indem ich es mit Intel MKL verknüpft habe. Wenn es Aussagen gibt wie: ( Siehe Wikipedia ) p=r; x=x+alpha*p r=r-alpha*Ap; oder ähnliche in QMR (in viel größerer Menge) v_tld = r; y = v_tld; rho = norm( y ); w_tld = r; z = w_tld; xi = norm( z ); (and more) Ist es...
Ich stelle die Leistung der linearen Löser von PETSc vor. So wie ich es verstehe, speedup=Sequential TimeParallel Time.speedup=Sequential TimeParallel Time.\text{speedup}=\frac{\text{Sequential Time}}{\text{Parallel Time}}. Ich weiß, dass das Ausführen des parallelen Codes auf einem Prozessor als...
Ich möchte künstliche neuronale Netze (ANN) verwenden, um die Reaktionsgeschwindigkeiten in meiner Flüssigkeit vorherzusagen, anstatt das gesamte System steifer ODEs zu lösen. Einige Leute aus meinem Labor haben bereits daran gearbeitet, damit ich nicht bei Null anfange, sondern Probleme mit meinen...