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 - 0,417 - 0,333 0,750 ]
Je résous l'équation différentielle avec les conditions initiales u (0) = u (1) = 0 , u '' (0) = u '' (1) = 0 . Ici \ sigma (x) \ geqslant \ sigma_ {0}> 0 est le paramètre. Sous forme d'opérateur, nous pouvons réécrire l'équation différentielle comme Au = f , où l'opérateur A est défini...
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
(J'espère que cette question convient à ce site; sinon, acceptez mes excuses). J'ai exécuté une certaine simulation et obtenu une série chronologique y (t), t = 0, 1, ... 20. Après avoir essayé certaines fonctions, j'ai constaté que: y(t) =~ 1 / (A t + B) Où A et B sont des coefficients I calculés...
Supposons que A soit une matrice générale clairsemée et que je souhaite calculer les valeurs propres. Je ne sais pas comment détecter la multiplicité des valeurs propres. Pour autant que je sache, pour un cas particulier, trouver les racines polynomiales par la méthode de la matrice compagnon, nous...
É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...
Combien est plus lent petsc4py vs c / c ++ / fortran? Je me rends compte que cela dépendra considérablement du code exécuté, mais qu'en est-il de quelque chose de simple comme un produit
J'implémente le document " Transport de masse optimal pour l'enregistrement et la déformation ", mon objectif étant de le mettre en ligne car je ne trouve aucun code de transport de masse eulérien en ligne et cela serait intéressant au moins pour la communauté des chercheurs en traitement d'images....
Il s'agit d'un simple problème d'énumération de symétrie. Je donne ici le contexte complet, mais aucune connaissance de la chimie quantique n'est nécessaire. L'intégrale de deux particules est: ⟨ i j | k l ⟩ = ∫ ψ * i ( x ) ψ * j ( x ' ) ψ k ( x ) ψ l ( x ' )⟨ij|kl⟩⟨ij|kl⟩\langle ij|kl\rangle Et il...
J'ai des ensembles de données de champ multipoint, chaque ensemble de données de point concerne une seule cellule d'un maillage non structuré. Le but est d'interpoler les données au centre de la cellule, directement ou indirectement, de la manière la plus précise. Si j'utilise l'interpolation...
J'essaie de résoudre une équation du type: (−∂2∂x2−f(x))ψ(x)=λψ(x)(−∂2∂x2−f(x))ψ(x)=λψ(x) \left( -\tfrac{\partial^2}{\partial x^2} - f\left(x\right) \right) \psi(x) = \lambda \psi(x) Où a un pôle simple à 0 , pour les N valeurs propres et vecteurs propres les plus petits . Les conditions aux...
Existe-t-il des algorithmes de détection de communauté pour les graphes bipartites (réseaux à 2 modes) implémentés dans igraph, networkX, R ou Python etc.? En particulier, existe-t-il une telle mise en œuvre dans laquelle on pourrait restreindre la détection des communautés uniquement sur l'un des...
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...
J'ai de l'expérience dans le codage d'OpenMP pour les machines à mémoire partagée (en C et FORTRAN) pour effectuer des tâches simples comme l'ajout de matrice, la multiplication, etc. (Juste pour voir comment il est en concurrence avec LAPACK). Je connais suffisamment OpenMP pour effectuer des...
La revue Association for Computing Machinery Transactions on Mathematical Software (ACM TOMS) publie de nombreux articles sur les algorithmes numériques qui incluent les implémentations logicielles. Selon leur politique éditoriale , la soumission d'un article d'algorithme inclut le code source...
Y a-t-il un espoir de résoudre efficacement le système linéaire suivant avec une méthode itérative? A∈Rn×n,x∈Rn,b∈Rn, with n>106A∈Rn×n,x∈Rn,b∈Rn, with n>106A \in \mathbb{R}^{n \times n}, x \in \mathbb{R}^n, b \in \mathbb{R}^n \text{, with } n > 10^6 Ax=bAx=bAx=b avec , où Δ est une matrice...
Supposons que soit une matrice symétrique réelle et que sa décomposition en valeurs propres V Λ V T soit donnée. Il est facile de voir ce qui se passe avec les valeurs propres de la somme A + c I où c est une constante scalaire (voir cette question ). Peut-on tirer une conclusion dans le cas...
Je produis actuellement des images 1d et 2d en utilisant Matlab pour comparer la précision avec un modèle donné. Je dois comparer mes méthodes avec le modèle .wfn gaussien standard et je vais le faire en analysant la densité d'une molécule ainsi que celle du laplacien. Je suis intéressé à produire...
Une question assez simple: pour faire une intégrale multidimensionnelle, étant donné que l'on a décidé qu'une sorte de méthode de Monte Carlo est appropriée, y a-t-il un avantage qu'une intégration MC régulière utilisant des nombres pseudo-aléatoires a sur une intégration quasi-Monte Carlo...
nVidia, par exemple, a CUBLAS, qui promet une accélération de 7 à 14 fois. Naïvement, cela est loin du débit théorique de l'une des cartes GPU de nVidia. Quels sont les défis liés à l'accélération de l'algèbre linéaire sur les GPU, et existe-t-il déjà des routages d'algèbre linéaire plus...