Vorrei eseguire un PCA su un set di dati composto da circa 40.000 campioni, ognuno dei quali visualizza circa 10.000 funzioni. L'uso della funzione princomp di Matlab richiede costantemente più di mezz'ora, a quel punto uccido il processo. Vorrei trovare un'implementazione / algoritmo che funziona in meno di 10 minuti. …
Sto producendo immagini 1d e 2d usando Matlab in questo momento per il confronto della precisione con un determinato modello. Devo confrontare i miei metodi con il modello standard .wfn gaussiano e lo farò analizzando la densità di una molecola e il laplaciano. Sono interessato a produrre immagini 2D di …
Una domanda abbastanza semplice: fare un integrale multidimensionale, dato che si è deciso che una sorta di metodo Monte Carlo è appropriato, c'è qualche vantaggio che una normale integrazione MC usando numeri pseudocasuali ha su un'integrazione quasi-Monte Carlo usando una sequenza quasirandom ? In tal caso, come riconoscerei le situazioni …
Supponiamo è una matrice reale simmetrica e la sua decomposizione autovalore è dato. È facile vedere cosa succede con gli autovalori della somma dove è una costante scalare (vedi questa domanda ). Possiamo trarre qualche conclusione nel caso generale cui è una matrice diagonale arbitraria? Grazie.V Λ V T A …
C'è qualche speranza nel risolvere efficacemente il seguente sistema lineare con un metodo iterativo? 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 con , dove Δ è una matrice molto sparsa con poche diagonali, derivante dalla …
La rivista Association for Computing Machinery Transactions on Mathematical Software (ACM TOMS) pubblica numerosi articoli su algoritmi numerici che includono implementazioni di software. Secondo la loro politica editoriale , l'invio di un documento dell'algoritmo include il codice sorgente per un'implementazione dell'algoritmo descritto in detto documento. Questo codice sorgente è soggetto …
Ho esperienza nella codifica di OpenMP per macchine a memoria condivisa (sia in C che FORTRAN) per svolgere compiti semplici come l'aggiunta di matrici, la moltiplicazione ecc. (Solo per vedere come compete con LAPACK). Conosco OpenMP abbastanza per svolgere compiti semplici senza la necessità di consultare la documentazione. Di recente, …
Esistono algoritmi per il rilevamento di comunità per grafici bipartiti (reti a 2 modalità) implementati in igraph, networkX, R o Python ecc.? In particolare, esiste una tale implementazione in cui si sarebbe in grado di limitare il rilevamento delle comunità solo in una delle due modalità?
La seguente equazione di matrice in Σ - per determinate matrici B e C - appare nel mio lavoro come una caratterizzazione di una matrice di covarianza. Ho imparato che questa equazione è nota, in particolare nella teoria del controllo del tempo continuo, come equazione di Lyapunov , e che …
Sto tentando di risolvere un'equazione del tipo: (−∂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) Dove ha un polo semplice a , per gli autovalori e gli autovettori più piccoli . Le condizioni al contorno sono: e , e sto solo guardando la funzione sopra .0 N …
Supponiamo che A sia una matrice sparsa generale e che voglio calcolare gli autovalori. Non so come rilevare la molteplicità degli autovalori. Per quanto ne so, per un caso speciale, trovando le radici polinomiali con il metodo della matrice compagna, possiamo applicare RRQR per rilevare la molteplicità delle radici.
Dato il sistema dove A ∈ R n × n , leggo che, nel caso in cui l'iterazione Jacobi sia utilizzata come solutore, il metodo non converge se b ha un componente diverso da zero nello spazio nullo di A . Quindi, come si può affermare formalmente che, a condizione …
Quanto è più lento petsc4py vs c / c ++ / fortran? Mi rendo conto che dipenderà in modo significativo dal codice in esecuzione, ma che dire di qualcosa di semplice come un prodotto a matrice vettoriale?
Ho set di dati di campo multipunto, ogni set di dati di punti si riferisce a una singola cella di una mesh non strutturata. L'obiettivo è interpolare i dati nel centro della cellula, direttamente o indirettamente, nel modo più accurato. Se utilizzo l'interpolazione ponderata a distanza inversa, nel caso in …
(Spero che questa domanda si adatti a questo sito; in caso contrario, accetta le mie scuse). Ho eseguito una certa simulazione e ho ottenuto una serie temporale y (t), t = 0, 1, ... 20. Dopo aver provato alcune funzioni, ho scoperto che: y(t) =~ 1 / (A t + …
We use cookies and other tracking technologies to improve your browsing experience on our website,
to show you personalized content and targeted ads, to analyze our website traffic,
and to understand where our visitors are coming from.
By continuing, you consent to our use of cookies and other tracking technologies and
affirm you're at least 16 years old or have consent from a parent or guardian.