Questions sur les aspects algorithmiques / informatiques de l'algèbre linéaire, y compris la solution des systèmes linéaires, les problèmes des moindres carrés, les problèmes propres et d'autres questions de ce type.
Je travaille actuellement avec "A Multigrid Tutorial" de Briggs et al, Chapter 8. La construction de l'opérateur d'interpolation est donnée comme suit: La construction de l'opérateur de restriction et de l'opérateur de grille fine est donnée comme suit: Supposons que nous ayons trois points de grille x0, x1, x2 avec …
J'essaie de comprendre s'il existe un moyen plus rapide de calculer toutes les valeurs propres et vecteurs propres d'une matrice de contiguïté très grande et clairsemée que d'utiliser scipy.sparse.linalg.eigsh Pour autant que je sache, cette méthode utilise uniquement la rareté et attributs de symétrie de la matrice. Une matrice d'adjacence …
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 systèmes clairsemés avec un motif …
J'utilise MATLAB pour résoudre un problème qui implique de résoudre à chaque pas de temps, où change avec le temps. En ce moment, j'accomplis ceci en utilisant MATLAB :bA x = bAx=b\mathbf{A} \mathbf{x}=\mathbf{b}bb\mathbf{b}mldivide x = A\b J'ai la flexibilité de faire autant de précalculs que nécessaire, donc je me demande …
Je recherche une bibliothèque qui effectue des opérations matricielles sur de grandes matrices clairsemées sans sacrifier la stabilité numérique. Les matrices seront 1000+ par 1000+ et les valeurs de la matrice seront comprises entre 0 et 1000. Je vais exécuter l' algorithme de calcul d'index donc je vais générer des …
Je me demande: quel est le meilleur algorithme pour résoudre réurét= A uréurét=UNEu\begin{equation} \frac{du}{dt} = Au \end{equation} OùUNEUNEAest unematricen × nn×nn\times nréelle. A n'est pas explicitement dépendant du temps, généralement clairsemé mais pas nécessairement en bandes. Ses valeurs propres ont des parties réelles non positives. A est également diagonalisable mais …
L'article «Expression Templates Revisited: A Performance Analysis of Current Methodologies» du SIAM Journal of Scientific Computing fait référence à la bibliothèque d'algèbre linéaire «Blaze». Je n'en ai jamais entendu parler auparavant et je n'arrive pas à trouver de références en ligne. (Les recherches google évidentes rendent le papier ci-dessus.) Alors, …
Je travaille sur une bibliothèque de matrices uniquement en-tête pour fournir un certain degré raisonnable de capacité d'algèbre linéaire dans un package aussi simple que possible, et j'essaie de faire le point sur l'état actuel de la technique: le calcul de la SVD d'un matrice complexe. Je fais une décomposition …
Dans les deux méthodes de décomposition de domaine (DD) et multigrille (MG), on peut composer l'application des mises à jour de bloc ou des corrections grossières comme additive ou multiplicative . Pour les solveurs ponctuels, c'est la différence entre les itérations de Jacobi et de Gauss-Seidel. Le lisseur multiplicatif pour …
Je résous un problème physique en utilisant un schéma numérique implicite. Cela m'amène à résoudre une équation linéaire avec une matrice tridiagonale. J'ai codé cet algorithme à partir de Wikipedia. Je me demande s'il existe une bibliothèque efficace qui permet de résoudre ce type d'équation de manière optimisée. Une note …
Étant donné deux matrices et B , je voudrais trouver les vecteurs x et y , tels que, min ∑ i j ( A i j - x i y j B i j ) 2 . Sous forme de matrice, j'essaie de minimiser la norme de Frobenius de A …
Supposons que A∈Rn×nA∈Rn×nA\in\mathbb{R}^{n\times n} est une matrice définie symétrique positive. AAA est suffisamment grand pour qu'il soit coûteux de résoudre directement .Ax=bAx=bAx=b Existe-t-il un algorithme itératif pour trouver la plus petite valeur propre de qui n'implique pas l'inversion de à chaque itération?AAAAAA Autrement dit, je devrais utiliser un algorithme itératif …
Dans le livre de Nocedal & Wright sur l'optimisation numérique, il y a une déclaration dans la section 2.2 (page 27), "D'une manière générale, il est plus facile de conserver l'invariance d'échelle pour les algorithmes de recherche de ligne que pour les algorithmes de région de confiance". Dans cette même …
Le produit intérieur complexe a deux définitions différentes décidées par des conventions: ou . Dans BLAS, j'ai trouvé les routines cdotu, zdotu et cdotc, zdotc. Les deux anciennes routines calculent en fait (un faux produit intérieur!) Et les deux dernières routines conjuguent le premier vecteur du produit intérieur. Aussi, par …
J'ai une liste LL{\cal L} de matrices symétriques dont j'ai besoin pour vérifier la semi-définition positive (c'est-à-dire que leurs valeurs propres ne sont pas négatives.) Le commentaire ci-dessus implique que l'on pourrait le faire en calculant les valeurs propres respectives et en vérifiant si elles ne sont pas négatives (peut-être …
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.