linear-algebra

5calore

2risposta

Avvolgere una funzione LAPACKE utilizzando Cython

Sto provando a racchiudere la funzione LAPACK dgtsv (un risolutore per i sistemi di equazioni tridiagonali) utilizzando Cython. mi sono imbattuto this previous answer, ma dal momento che dgtsv non è u

6calore

2risposta

colonna-saggio Eigen C++

C'è un modo semplice per valutare colonna prodotto scalare saggio di 2 matrici (lascia chiamarli A e B, di tipo Eigen::MatrixXd) che hanno dimensioni mxn, senza valutare A*B o senza dover ricorrere ai

5calore

1risposta

Risolvi sistema triangolare superiore spoglio

Sto cercando di capire come risolvere un sistema triangolare sparse in modo efficiente, Au * x = b in scarse sparse. Ad esempio, si può costruire una rada matrice triangolare superiore, Au, e un lato

5calore

1risposta

Come posso proiettare i vettori su un piano definito dal suo vettore ortogonale in Python?

Ho un piano, plane A, definito dal suo vettore ortogonale, ad esempio (a, b, c). (cioè il vettore (a, b, c) ortogonale al plane A) desidero proiettare un vettore (d, e, f) su plane A. Come posso farlo

8calore

2risposta

Perché l'einsum di Numpy è più lento delle funzioni incorporate di numpy?

Di solito ho ottenuto una buona performance dalla funzione di einsum di Numpy (e mi piace che sia sintassi). @ La risposta di Ophion a this question mostra che - per i casi testati - einsum supera cos

6calore

3risposta

Algoritmo più veloce per calcolare il determinante di una matrice?

Per un documento di ricerca, sono stato assegnato alla ricerca l'algoritmo più veloce per calcolare il determinante di una matrice. so già su decomposizione LU e algoritmo Bareiss che sia eseguito in

18calore

1risposta

Come implementare la mldivide di Matlab (ovvero l'operatore backslash "\")

Attualmente sto cercando di sviluppare una piccola libreria matematica orientata alla matrice (sto usando Eigen 3 per le strutture e le operazioni di dati matriciali) e volevo implementare alcuni prat

22calore

4risposta

Matrice triangolare superiore dell'indice lineare

Se si dispone della porzione triangolare superiore di una matrice, spostata sopra la diagonale, memorizzata come matrice lineare, come si possono estrarre gli indici di un elemento di matrice dall'ind

5calore

1risposta

Traccia di una matrice di crossproduct - un calcolo più veloce?

Sto cercando un calcolo veloce nella R della traccia (trace (A)) di una matrice A = B' C. Il modo più veloce che posso pensare è il seguente: set.seed(123) n <- 10^6 B <- matrix(rnorm(n), ncol=sqrt(

6calore

2risposta

Numpy: Generalized autovalori Problema

Sto cercando di risolvere un problema del tipo: dove Aw = xBwx è uno scalare (autovalore), w è un autovettore e A e B sono simmetriche, matrici NumPy quadrati di uguale dimensione. Dovrei essere in gr