Ich möchte für meine wissenschaftliche Arbeit unabhängig von kommerzieller Software sein. Ich finde eine Abhängigkeit von kommerziellen Paketen wie Matlab und seinen Toolboxen unbefriedigend, weil ich nicht weiß, ob ich in Zukunft Zugang zu Matlab haben werde und weil mir die Sprache nicht gefällt....
Betrachte mit nahezu singulär, was bedeutet, dass es einen sehr kleinen Eigenwert von gibt. Das übliche einer iterativen Methode basiert auf dem Residuum und die Iterationen anhalten können, wenn mit der Iterationsnummer. Aber in dem Fall, den wir betrachten, könnte es einen großen Fehler , der in...
Angenommen , ich habe eine 2D - Gitter , bestehend aus nicht - überlappende Dreiecken , und eine Menge von Punkten { P i } M i = 1 ⊂ ∪ N k = 1 T K . Wie lässt sich am besten bestimmen, in welchem Dreieck jeder der Punkte liegt?{ Tk}Nk = 1{Tk}k=1N\{T_k\}_{k=1}^N{ pich}Mi = 1⊂ ∪Nk =...
Ich bin ein großer Verfechter der testgetriebenen Entwicklung im wissenschaftlichen Rechnen. Die Nützlichkeit in der Praxis ist einfach umwerfend und lindert die klassischen Probleme, die Codeentwickler kennen. Das Testen von wissenschaftlichen Codes, die bei der allgemeinen Programmierung nicht...
Ich würde gerne wissen, ob es einen schnellen Weg gibt, den euklidischen Abstand zweier Vektoren in Oktave zu berechnen. Es scheint, dass es dafür keine spezielle Funktion gibt. Soll ich also einfach die Formel mit verwenden
Ich führe eine Textklassifizierungsaufgabe mit R durch und erhalte eine Dokumentterm-Matrix mit einer Größe von 22490 mal 120.000 (nur 4 Millionen Einträge ungleich Null, weniger als 1% Einträge). Jetzt möchte ich die Dimensionalität mit PCA (Principal Component Analysis) reduzieren. Leider kann R...
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...
Was ich an MPI am wenigsten mag, ist der Umgang mit Datentypen (dh Datenmaps / Masken), weil sie nicht so gut zu objektorientiertem C ++ passen. boost::mpiunterstützt jedoch nur MPI 1.1 von ihrer Website: boost :: mpi ist eine C ++ - freundliche Schnittstelle zur Standardschnittstelle für die...
Ich habe in Wikipedia über Computerwissenschaften gelesen, aber mein Verständnis ist nicht sehr klar. Umfasst die Computerwissenschaft das Programmieren? Wie unterschiedlich ist die Computerwissenschaft von der Computerwissenschaft _ ____ , bei der es sich um eine beliebige Disziplin handeln kann...
Bei der Berechnung der QR-Faktorisierung in der Praxis werden Householder-Reflexionen verwendet, um den unteren Teil einer Matrix auf Null zu setzen. Ich weiß, dass für die Berechnung von Eigenwerten symmetrischer Matrizen das Beste, was Sie mit Householder-Reflexionen tun können, darin besteht,...
Gibt es beim Programmieren dichter Matrixberechnungen einen Grund, ein Zeilen-Hauptlayout des über dem Spalten-Hauptlayout liegenden zu wählen? Ich weiß, dass wir abhängig vom Layout der gewählten Matrix den entsprechenden Code schreiben müssen, um die Cache-Speicher aus Geschwindigkeitsgründen...
Es ist wahrscheinlich eine Frage auf Schülerebene, aber ich kann sie nicht genau für mich selbst tun. Warum ist es genauer, in den numerischen Methoden ungleichmäßige Gitter zu verwenden? Ich denke im Kontext einer Finite-Differenzen-Methode für die PDE der Form . Und nehme an, ich interessiere...
Da sich diese Frage auf die Berechnung bezieht, habe ich mich entschlossen, hier zu posten. Hoffentlich wird es als angemessen angesehen. Ich habe gerade angefangen, atmosphärische und ozeanische Modelle zu verwenden, und mir ist klar, dass ich mehr Kerne, Speicher und Speicherplatz benötige als...
Annehmen min Av e c (U)vorbehaltlich Uich , j≤ max { Uich , k, Uk , j} ,i , j , k = 1 , ... , nMindestEINvec(U)unterliegen Uich,j≤max{Uich,k,Uk,j},ich,j,k=1,…,n\begin{align*} \min A &\mathrm{vec}(U) \\ &\text{subject to } U_{i,j} \leq \max\{U_{i,k}, U_{k,j}\}, \quad i,j,k = 1, \ldots,...
Gibt es Gründe, warum man implizite Runge Kutta (IMRK) höherer Ordnung anstelle von BDF-Zeitschritten wählen sollte? BDF scheint mir viel einfacher zu sein, da Stufe IMRK q lineare Lösungen pro Zeitschritt benötigt. Stabilität für BDF und IMRK scheint ein strittiger Punkt zu sein. Ich kann keine...
Ich benutze Valgrind + Callgrind, um einen von mir geschriebenen Solver zu profilieren. Wie im Valgrind-Benutzerhandbuch angegeben, habe ich meinen Code mit den Debugging-Optionen für den Compiler kompiliert: "Ohne Debugging-Informationen können die Valgrind-Tools nur raten, zu welcher Funktion...
Ich frage mich, was mit Polynom-Vorkonditionierern passiert ist. Ich interessiere mich für sie, weil sie aus mathematischer Sicht vergleichsweise elegant erscheinen, aber soweit ich in Umfragen zu Krylov-Methoden gelesen habe, schneiden sie im Allgemeinen als Vorkonditionierer sehr schlecht ab....
Ich habe einen Datensatz x1,x2,…,xkx1,x2,…,xkx_{1}, x_{2}, \ldots, x_{k} und möchte den Parameter mmm so finden, dass er die Summe minimiert m - x i | . ∑i=1k∣∣m−xi∣∣.∑i=1k|m−xi|.\sum_{i=1}^{k}\big|m-x_i\big|. das
Nach meinem Kenntnisstand verwenden Multigrid-Löser iterative Glätter wie Jacobi, Gauss-Seidel und SOR, um den Fehler bei verschiedenen Frequenzen zu dämpfen. Könnte stattdessen eine Krylov-Subraummethode (wie Konjugatgradient, GMRES usw.) verwendet werden? Ich glaube nicht, dass sie als...
Ich versuche, eine 2D-Poisson-Gleichung durch endliche Differenzen zu lösen. Dabei erhalte ich eine spärliche Matrix mit nur Variablen in jeder Gleichung. Wenn die Variablen beispielsweise wären, würde die Diskretisierung ergeben:555UUU