Ich möchte eine PCA für einen Datensatz durchführen, der aus ungefähr 40 000 Proben besteht, wobei jede Probe ungefähr 10 000 Merkmale aufweist. Die Verwendung der Matlab-Princomp-Funktion dauert durchweg über eine halbe Stunde. An diesem Punkt beende ich den Prozess. Ich möchte eine...
CVXOPT: http://abel.ee.ucla.edu/cvxopt/index.html OpenOpt: http://openopt.org/Welcome Wie ist die Beziehung zwischen ihnen? Was sind die Vor- und Nachteile von ihnen? Übrigens, gibt es eine andere hochwertige konvexe Allzweck-Optimierungsbibliothek für Python / C ++, die es wert ist, erwähnt zu...
Gibt es eine veröffentlichte Liste bewährter Verfahren zur Gewährleistung der Langlebigkeit von Code mit Blick auf reproduzierbare wissenschaftliche Ergebnisse? (z. B. Open Source, Dokumentationspraktiken, Auswahl von Abhängigkeiten, Auswahl einer Sprache, virtueller Maschinen usw.). Kennen Sie...
Die meisten Mathematikbibliotheken verfügen über eine Reihe von Versionen von Logarithmusfunktionen. Die meiste Zeit gehen wir davon aus, dass sie perfekt sind, aber tatsächlich bieten ziemlich viele von ihnen nur eine bestimmte Anzahl von Stellen mit Präzision. Für einige Funktionen gibt es...
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...
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...
Ich löse für eine riesige, spärlich positive, definitive Matrix Verwendung der Methode des konjugierten Gradienten (CG). Es ist möglich, die Determinante von Verwendung der während der Lösung erzeugten Informationen zu berechnen .A A.A x =
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...
Möchten Sie diesen Beitrag verbessern? Geben Sie detaillierte Antworten auf diese Frage, einschließlich Zitaten und einer Erklärung, warum Ihre Antwort richtig ist. Antworten ohne ausreichende Details können bearbeitet oder gelöscht werden. Ich bin völlig neu in der...
Angenommen, ich habe eine Funktion , die ich über eine Tetraeder- integrieren möchte . Wenn willkürlich wäre, wäre die Gauß-Quadratur eine gute Lösung, aber ich weiß zufällig, dass harmonisch ist. Wie viel kann die Gauß-Quadratur mit diesen Informationen beschleunigt werden? T ⊂ R 3 f...
Wie skaliert theoretisch die Zeit für eine Berechnung der Dichtefunktionaltheorie (DFT) mit der Anzahl der Elektronen? Ich interessiere mich für "typische" DFT-Implementierungen wie VASP, ABINIT usw., nicht für O (N)
Ich habe ein inhomogenes lineares System Ax=bAx=b Ax=b wobei AAA eine reelle n×nn×nn\times n Matrix mit n≤4n≤4n\leq 4 . AAA wird garantiert, dass der Nullraum von A eine Dimension von Null hat, sodass die Gleichung eine eindeutige Umkehrung von x=A−1bx=A−1bx=A^{-1} b . Da das Ergebnis auf der...
Bei einer unbekannten Funktion können wir ihren Wert an jedem Punkt in ihrer Domäne bewerten, aber wir haben keinen Ausdruck. Mit anderen Worten, f ist für uns wie eine Black Box.f:Rd→Rf:Rd→Rf:\mathbb R^d \to \mathbb Rfff Wie heißt das Problem, den Minimierer von ? Welche Methoden gibt es da...
Wie berechnet MATLAB beispielsweise die SVD einer bestimmten Matrix? Ich gehe davon aus, dass die Antwort wahrscheinlich die Berechnung der Eigenvektoren und Eigenwerte von beinhaltet A*A'. Wenn das der Fall ist, würde ich auch gerne wissen, wie es diese
Ich habe ein physikalisches Problem, das durch die Poisson-Gleichung in zwei Dimensionen bestimmt wird Ich habe Messungen der beiden Gradientenkomponenten ∂ u / ∂ x und ∂ u / ∂ y entlang eines Teils der Grenze Γ m , möchte also ∂ u auferlegen −∇2u=f(x,y),inΩ−∇2u=f(x,y),inΩ -\nabla^2 u = f(x,y), \;...
Bei der Arbeit arbeite ich im Wesentlichen als unabhängiger Berater. Für das Management und die Kunden muss ich die Zeit abschätzen, die für die Entwicklung von Software im Rahmen meiner rechnerwissenschaftlichen Forschung benötigt wird. Meine Zeitschätzungen sind jedoch normalerweise nicht...
Ich habe einen Datensatz mit dreidimensionalen Punkten, für die ich mit Python ein Netz erstellen möchte. Für die gesamte Software, die ich gesehen habe, müssen Sie die Kanten bereitstellen. Gibt es ein Programm in Python, das eine Reihe von Punkten in 3D als Eingabe verwendet und ein Netz ausgibt?...
Beginnend mit der Advektionsgleichung in Erhaltungsform. ut= ( a ( x ) u )xut=(a(x)u)x u_t = (a(x)u)_x wobei eine Geschwindigkeit ist, die vom Raum abhängt, und u eine Konzentration einer Spezies ist, die konserviert ist.a ( x )a(x)a(x)uuu Die Diskretisierung des Flusses (wobei der Fluss an den...
Ich entwickle derzeit eine Domänenzerlegungsmethode zur Lösung des Streuproblems. Grundsätzlich löse ich iterativ ein System von Helmholtz-BVPs. Ich diskretisiere die Gleichungen mit der Finite-Elemente-Methode über Dreiecks- oder Tetraedernetzen. Ich entwickle den Code für meine Doktorarbeit. Ich...
Ich bekomme eine Matrix , die symmetrisch, invertierbar, positiv definitiv und dicht ist. Ich muss testen, ob wobei J die All-One-Matrix ist.Q det ( Q ) = det ( 12 I - Q - J )12×1212×1212 \times 12QQQdet(Q)=det(12I−Q−J)(1)det(Q)=det(12I−Q−J)(1)\det(Q) = \det(12I-Q-J) \; \; (1)JJJ Ich mache das...