Science computationnelle matrix

5

Calcul de la structure clairsemée pour les matrices d'éléments finis

Question: Quelles méthodes sont disponibles pour calculer avec précision et efficacité la structure de rareté d'une matrice d'éléments finis? Info: Je travaille sur un solveur d'équation de pression de Poisson, en utilisant la méthode de Galerkin avec une base de Lagrange quadratique, écrite en C, et en utilisant PETSc pour …

13 matrix finite-element petsc performance

3

Précision en virgule flottante simple ou double

Les nombres à virgule flottante simple précision occupent la moitié de la mémoire et sur les machines modernes (même sur les GPU, il semble) les opérations peuvent être effectuées avec eux à presque deux fois la vitesse par rapport à la double précision. De nombreux codes FDTD que j'ai trouvés …

13 matrix linear-solver precision

5

Résolution répétée de

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 …

12 linear-algebra matlab linear-solver matrix linear-system

1

Algorithmes pour les grandes matrices entières clairsemées

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 …

12 linear-algebra algorithms matrix

3

Quel est l'état actuel de la technique concernant les algorithmes de décomposition de valeurs singulières?

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 …

12 linear-algebra matrix svd

1

problème SVD pondéré?

É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 …

12 linear-algebra matrix reference-request

3

Tester si deux matrices 12x12 ont le même déterminant

On me donne une matrice symétrique, inversible, définie positive et dense. Je dois tester si où J est la matrice des uns.12×1212×1212 \times 12QQQJdet(Q)=det(12I−Q−J)(1)det(Q)=det(12I−Q−J)(1)\det(Q) = \det(12I-Q-J) \; \; (1)JJJ Je fais actuellement cela avec la bibliothèque de tatou mais cela s'avère trop lent. Le fait est que je dois le …

11 linear-algebra matrix c++ graph-theory c

2

Comment la SVD d'une matrice est-elle calculée en pratique

Comment MATLAB, par exemple, calcule-t-il la SVD d'une matrice donnée? Je suppose que la réponse implique probablement de calculer les vecteurs propres et les valeurs propres de A*A'. Si tel est le cas, j'aimerais également savoir comment calcule-t-on ces informations?

11 linear-algebra matrix

3

Bibliothèques pour résoudre l'équation de Lyapunov

L'équation de matrice suivante dans Σ - pour les matrices B et C données - apparaît dans mon travail comme une caractérisation d'une matrice de covariance. J'ai appris que cette équation est connue, en particulier dans la théorie du contrôle du temps continu, comme l'équation de Lyapunov , et qu'il …

11 matrix matrix-equations

1

Projection de l'espace nul de partir de dans

Étant donné le système où , j'ai lu que, dans le cas où l'itération Jacobi est utilisée comme solveur, la méthode ne convergera pas si a un non-zéro composant dans l'espace nul de . Alors, comment pourrait-on déclarer formellement que, à condition que ait une composante non nulle couvrant l'espace …

11 linear-algebra optimization matrix iterative-method linear-solver

4

Trouver la racine carrée d'une matrice laplacienne

Supposons que la matrice suivante est donnée [ 0,500 - 0,333 - 0,167 - 0,500 0,667 - 0,167 - 0,500 - 0,333 0,833 ] avec sa transposée A T . Le produit A T A = G donne [ 0,750 - 0,334 - 0,417 - 0,334 0,667 - 0,333 - …

11 linear-algebra matrix eigensystem

2

Matrice exponentielle d'une matrice hamiltonienne

Soit des matrices réelles, carrées et denses. G et Q sont symétriques. LaisserA , G , QA,G,QA, G, QgGGQQQ H= [ A- Q- G- unT]H=[A−G−Q−AT]H = \begin{bmatrix} A & -G \\ -Q &-A^T \end{bmatrix} être une matrice hamiltonienne. Je veux calculer la matrice exponentielle de . J'ai besoin de l'exponentielle …

10 linear-algebra matrix dense-matrix

1

Pourquoi SVD parle de moins que QR et LU pour une matrice clairsemée?

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 clairsemée?

10 linear-algebra sparse-matrix matrix svd eigen

4

Existe-t-il un package matriciel arithmétique quadruple double?

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 linéaires numériques sous arithmétique quadruple double. Cependant, il est …

10 sparse-matrix matrix condition-number

1

Résolution d'un système Ax = b simple en parallèle avec PETSc

Je suis nouveau dans le package PETSc. J'ai une matrice A ~ 4000x4000 au format matriciel et je veux que PETSc résolve cela en utilisant plusieurs processeurs. Je sais comment résoudre le système sur un seul processeur, mais je ne sais pas comment répartir la matrice et les vecteurs entre …

10 petsc matrix

Questions marquées «matrix»