Angenommen, die folgende Matrix ist mit ihrer Transponierten . Das Produkt ergibt ,[ 0,500 - 0,333 - 0,167 - 0,500 0,667 - 0,167 - 0,500 - 0,333 0,833 ] A T A T A = G [ 0,750 - 0,334 - 0,417 - 0,334 0,667 - 0,333 - 0,417 - 0,333 0,750 ]EINAA⎡⎣⎢0,500- 0,500- 0,500- 0,3330,667- 0,333- 0,167-...
(σ2(x)u′′(x))′′=f(x),0⩽x⩽1(σ2(x)u″(x))″=f(x),0⩽x⩽1 \left( \sigma^{2}(x) u ''(x) \right)'' = f(x), \;\;\; 0 \leqslant x \leqslant 1 u(0)=u(1)=0u(0)=u(1)=0u(0) = u(1) = 0u′′(0)=u′′(1)=0u″(0)=u″(1)=0u''(0) = u''(1) = 0σ(x)⩾σ0>0σ(x)⩾σ0>0\sigma(x) \geqslant \sigma_{0} > 0Au=fAu=fAu = fAAA Nach dem...
Bitte listen Sie das Python-Paket (Jewsc4py usw.) und die von ihm unterstützten spärlichen Direktlöser auf. Bitte eine (Community-Wiki) Antwort pro
(Ich hoffe, diese Frage passt zu dieser Seite; wenn nicht, akzeptiere meine Entschuldigung). Ich habe eine bestimmte Simulation ausgeführt und eine Zeitreihe y (t), t = 0, 1, ... 20 erhalten. Nachdem ich einige Funktionen ausprobiert hatte, stellte ich Folgendes fest: y(t) =~ 1 / (A t +...
Angenommen, A ist eine allgemeine Matrix mit geringer Dichte, und ich möchte die Eigenwerte berechnen. Ich weiß nicht, wie ich die Multiplizität für die Eigenwerte ermitteln soll. Soweit ich weiß, können wir für einen speziellen Fall, bei dem die Polynomwurzeln mit der Begleitmatrixmethode...
Angesichts des Systems in dem A ∈ R n × n ist , habe ich gelesen, dass bei Verwendung der Jacobi-Iteration als Löser die Methode nicht konvergiert, wenn b eine Nicht-Null-Komponente im Nullraum von A hat . Wie könnte man also formal sagen, dass die Jacobi-Methode nicht konvergent ist ,...
Wie viel langsamer ist Haustierec4py vs c / c ++ / fortran? Mir ist klar, dass dies erheblich vom ausgeführten Code abhängt, aber was ist mit etwas Einfachem wie einem
Ich implementiere das Papier " Optimaler Massentransport für Registrierung und Warping ". Mein Ziel ist es, es online zu stellen, da ich online keinen eulerschen Massentransportcode finden kann. Dies wäre zumindest für die Forschungsgemeinschaft in der Bildverarbeitung interessant. Das Papier kann...
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 habe Sätze von Mehrpunktfelddaten, jeder Punktdatensatz bezieht sich auf eine einzelne Zelle eines unstrukturierten Netzes. Ziel ist es, die Daten direkt oder indirekt auf genaueste Weise direkt in das Zellzentrum zu interpolieren. Wenn ich die inverse entfernungsgewichtete Interpolation...
Ich versuche eine Gleichung vom Typ zu lösen: ( - ∂2∂x2- f( x ) ) ψ(x)=λψ(x)(−∂2∂x2−f(x))ψ(x)=λψ(x) \left( -\tfrac{\partial^2}{\partial x^2} - f\left(x\right) \right) \psi(x) = \lambda \psi(x) Wobei einen einfachen Pol bei 0 hat , für die kleinsten N Eigenwerte und Eigenvektoren. Die...
Gibt es Algorithmen zur Community-Erkennung für zweigeteilte Graphen (2-Mode-Netzwerke), die in igraph, networkX, R oder Python usw. implementiert sind? Gibt es insbesondere eine solche Implementierung, bei der die Erkennung von Communities nur in einem der beiden Modi eingeschränkt werden...
Die folgende Matrixgleichung in Σ - für gegebene B- und C- Matrizen - erscheint in meiner Arbeit als Charakterisierung einer Kovarianzmatrix. Ich habe gelernt, dass diese Gleichung, insbesondere in der Theorie der kontinuierlichen Zeitsteuerung, als Lyapunov-Gleichung bekannt ist und dass es...
Ich habe Erfahrung in der Codierung von OpenMP für Shared Memory-Maschinen (sowohl in C als auch in FORTRAN), um einfache Aufgaben wie Matrixaddition, Multiplikation usw. auszuführen (nur um zu sehen, wie es mit LAPACK konkurriert). Ich kenne OpenMP genug, um einfache Aufgaben auszuführen, ohne die...
Die Zeitschrift Association for Computing Machinery Transactions on Mathematical Software (ACM TOMS) veröffentlicht zahlreiche Artikel zu numerischen Algorithmen, die Software-Implementierungen enthalten. Gemäß ihrer redaktionellen Richtlinie enthält die Einreichung eines Algorithmuspapiers den...
Gibt es Hoffnung, das folgende lineare System mit einer iterativen Methode effizient zu lösen? A∈Rn×n,x∈Rn,b∈Rn, with n>106A∈Rn×n,x∈Rn,b∈Rn, with n>106A \in \mathbb{R}^{n \times n}, x \in \mathbb{R}^n, b \in \mathbb{R}^n \text{, with } n > 10^6 Ax=bAx=bAx=b mit A=(Δ−K)A=(Δ−K) A=(\Delta - K) ,...
Angenommen, ist eine reelle symmetrische Matrix und ihre Eigenwertzerlegung V Λ V T ist gegeben. Es ist leicht zu erkennen, was mit den Eigenwerten der Summe A + c I passiert, wobei c eine Skalarkonstante ist (siehe diese Frage ). Können wir im allgemeinen Fall A + D, in dem D eine beliebige...
Ich produziere gerade 1d- und 2d-Bilder mit Matlab, um die Genauigkeit mit einem bestimmten Modell zu vergleichen. Ich muss meine Methoden mit dem Standard-Gaußschen .wfn-Modell vergleichen und werde dazu die Dichte eines Moleküls sowie des Laplace-Modells analysieren. Ich bin daran interessiert,...
Eine recht einfache Frage: Um ein mehrdimensionales Integral zu erstellen, gibt es einen Vorteil, den eine reguläre MC-Integration unter Verwendung von Pseudozufallszahlen gegenüber einer Quasi-Monte-Carlo-Integration unter Verwendung einer Quasirandom-Sequenz hat, wenn man entschieden hat, dass...
nVidia hat zum Beispiel CUBLAS, das eine 7-14-fache Beschleunigung verspricht. Naiv ist dies bei weitem nicht der theoretische Durchsatz einer der GPU-Karten von nVidia. Was sind die Herausforderungen bei der Beschleunigung der linearen Algebra auf GPUs, und gibt es bereits schnellere lineare...