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 voudrais prédire les temps d'exécution des opérations d'algèbre linéaire dense sur une architecture spécifique en utilisant une bibliothèque spécifique. Je voudrais apprendre un modèle qui se rapproche de la fonction Fo p: :Fop::F_{op} \;::\; tailles d'entrée runtime→→ \rightarrow pour des opérations telles que multiplication matricielle, ajout par élément, résolution …
Je fais des recherches sur la structure des compléments Schur et trouve un phénomène intéressant: Supposons que A est du laplacien 5 pt. Si j'utilise l'ordre de dissection imbriqué et la méthode multifrontale pour calculer la factorisation LU, puis vérifier le dernier bloc de complément de Schur, il a un …
VVV∥⋅∥‖⋅‖\|\cdot\|F:V→RF:V→RF : V \rightarrow \mathbb R Je voudrais estimer la norme de (d'en haut et d'en bas). Comme est une boîte noire, la seule façon de le faire est de le tester avec des vecteurs unitaires de et, en fonction du résultat, de trouver qui maximise.FFFFFFVVVv∈S1Vv∈S1Vv \in S^1 V|F(v)||F(v)||F(v)| Connaissez-vous …
Je sais que pour résoudre le problème des valeurs propres symétriques , nous pouvons utiliser la loi d'inertie de Sylvester, c'est-à-dire que le nombre de valeurs propres de A inférieur à a est égal au nombre d'entrées négatives de D où la matrice diagonale D provient de la factorisation de …
Étant donné une matrice générique clairsemée avec m << n (correction: ) éléments non nuls (typiquement ). est générique dans le sens où il n'a pas de propriétés spécifiques (par exemple, un caractère définitif positif), et aucune structure (par exemple, des bandes) n'est supposée.A ∈ Rn × nUNE∈Rn×nA \in \mathbb{R}^{n\times …
Il est montré (Yousef Saad, Méthodes itératives pour les systèmes linéaires clairsemés , p. 260) quec o n d( A′A ) ≈ c o n d( A )2cond(A′A)≈cond(A)2cond(A'A) \approx cond(A)^2 Est-ce également vrai pour les ?UNEA′AA′AA' Dans le cas où est avec , j'observe queN × M N ≪ M …
Supposons que le système linéaire suivant soit donné où est le Laplacien pondéré connu pour être défini positif avec un espace nul unidimensionnel par , et la variance de translation de , c'est-à-dire que ne change pas la valeur de la fonction (dont la dérivée est ). Les seules entrées …
J'ai quelques questions concernant les points suivants: J'essaie de résoudre l'équation de Schrodinger en 1D en utilisant la discrétisation de manivelle nicolson suivie de l'inversion de la matrice tridiagonale résultante. Mon problème a maintenant évolué vers un problème avec des conditions aux limites périodiques et j'ai donc modifié mon code …
Des systèmes indéfinis de matrices apparaissent par exemple dans la discrétisation des problèmes de point de selle par des éléments finis mixtes. La matrice du système peut alors être mise sous la forme ( ABBtC)(UNEBtBC)\begin{pmatrix} A & B^t \\ B & C\end{pmatrix} où est négatif (semi) défini, C est positif …
Lors de la résolution de systèmes linéaires clairsemés à l'aide de méthodes de factorisation directe, la stratégie de classement utilisée a un impact significatif sur le facteur de remplissage des éléments non nuls dans les facteurs. Une telle stratégie d'ordonnancement est la dissection imbriquée. Je me demande s'il est possible …
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.