Questions marquées «sparse-matrix»

J'ai plusieurs problèmes d'optimisation globale non convexe difficiles à résoudre. Actuellement, j'utilise la boîte à outils Optimization de MATLAB (en particulier, fmincon()avec algorithm = 'sqp'), ce qui est assez efficace . Cependant, la majeure partie de mon code est en Python et j'aimerais...

Les systèmes linéaires clairsemés apparaissent de plus en plus fréquemment dans les applications. On a beaucoup de routines à choisir pour résoudre ces systèmes. Au niveau le plus élevé, il existe un fossé entre les méthodes directes (par exemple, l’élimination gaussienne ou la décomposition de...

EDIT: Je teste si des valeurs propres ont une magnitude de un ou plus. J'ai besoin de trouver la plus grande valeur propre absolue d'une grande matrice clairsemée et non symétrique. J'ai utilisé la eigen()fonction de R , qui utilise l'algo QR d'EISPACK ou de LAPACK pour trouver toutes les valeurs...

J'ai des matrices et . est clairsemé et est avec très grand (peut être de l'ordre de plusieurs millions.) est une matrice haut avec plutôt petit ( ) et chaque colonne peut seulement un seul entrée avec le reste étant « s, de telle sorte que . est énorme, il est donc très difficile à inverser, et je...

Je voudrais savoir s'il existe un moyen rapide de calculer la distance euclidienne de deux vecteurs en octave. Il semble qu'il n'y ait pas de fonction spéciale pour cela, alors devrais-je simplement utiliser la formule avec

J'écris une petite bibliothèque pour des calculs matriciels clairsemés comme un moyen de m'apprendre à faire le meilleur usage de la programmation orientée objet. J'ai travaillé très dur pour avoir un joli modèle objet, où les parties (matrices clairsemées et les graphiques qui décrivent leur...

Je me demandais s'il existe une méthode rapide et efficace pour trouver à l'avance le nombre de non-zéros pour une opération de multiplication matricielle clairsemée en supposant que les deux matrices sont au format CSC ou CSR. Je sais qu'il y en a un dans le paquet smmp mais j'ai besoin de quelque...

Je développe un code plus grand pour effectuer des calculs de valeurs propres d'énormes matrices clairsemées, dans le contexte de la physique computationnelle. Je teste mes routines contre le simple oscillateur harmonique dans une dimension, car les valeurs propres sont bien connues analytiquement....

J'essaie de résoudre une équation de Poisson 2D par des différences finies. Dans le processus, j'obtiens une matrice clairsemée avec seulement variables dans chaque équation. Par exemple, si les variables étaient U , alors la discrétisation donnerait:555UUU Ui−1,j+Ui+1,j−4Ui,j+Uje,j−1+Uje,j+1=fi ,...

J'ai deux graphiques avec près de n ~ 100 000 nœuds chacun. Dans les deux graphiques, chaque nœud est connecté à exactement 3 autres nœuds, de sorte que la matrice d'adjacence est symétrique et très clairsemée. La partie difficile est que j'ai besoin de toutes les valeurs propres de la matrice...

Je me demande si l'algorithme de Thomas est le moyen le plus rapide (de manière probable?) De résoudre un système tridiagonal clairsemé à dominante diagonale symétrique en termes de complexité algorithmique (ne pas chercher de packages d'implémentation comme LAPACK, etc.). Je sais que l'algorithme...

Je dois résoudre le même système linéaire clairsemé (300x300 à 1000x1000) avec de nombreux côtés droits (300 à 1000). En plus de ce premier problème, je voudrais également résoudre différents systèmes, mais avec les mêmes éléments non nuls (juste des valeurs différentes), c'est-à-dire de nombreux...

Veuillez lister le paquet Python (petsc4py, etc ...) et les solveurs directs clairsemés qu'il prend en charge. Une réponse (wiki communautaire) par paquet, s'il vous

La multiplication matricielle (Mat * Mat et Mat * Vec) est-elle proportionnelle au nombre de non-zéros ou à la taille de la matrice? Ou une combinaison des deux. Et la forme. Par exemple, j'ai une matrice 100 x 100 contenant 100 valeurs ou une matrice 1000 x 1000 contenant 100 valeurs. Lors de la...

La compilation de PETSc avec une bibliothèque BLAS / LAPACK externe affecte-t-elle de manière significative les performances sur des matrices clairsemées, ou utilise-t-elle uniquement ces bibliothèques pour les mathématiques matricielles

Je travaille actuellement sur la résolution de très grands systèmes symétriques (mais pas définis positifs), générés par certains algorithmes. Ces matrices ont une belle densité de blocs qui peut être utilisée pour la résolution parallèle. Mais je ne peux pas décider si je dois utiliser une...

Je travaille sur un grand système d'équations linéaire clairsemé et mal conditionné. Je veux utiliser l'arithmétique double-double ou l'arithmétique quad-double pour les résoudre. Je sais qu'il existe un package nommé MPACK développé par Nakata, Maho, qui peut effectuer des calculs algébriques...

Par exemple, les bibliothèques de matrice clairsemée C ++ que j'ai utilisées - Eigen et SuiteSparse, elles ne semblent pas avoir de fonctionnalité SVD pour la matrice clairsemée. Donc, juste curieux, SVD est-il plus difficile que QR / LU pour une matrice

Supposons que j'ai le grand système linéaire clairsemé d'origine: . Maintenant, je n'ai pas A - 1 car A est trop grand pour factoriser ou toute sorte de décomposition de A , mais supposons que j'ai la solution x 0 trouvée avec une résolution itérative.A x0= b0Ax0=b0A\textbf{x}_0=\textbf{b}_0UNE-...

J'ai l'intention de résoudre Ax = b où A est une matrice carrée ou rectangulaire complexe, clairsemée, asymétrique et très mal conditionnée (nombre de conditions ~ 1E + 20). J'ai pu résoudre le système avec ZGELSS dans LAPACK avec précision. Mais à mesure que les degrés de liberté dans mon système...