Lorsque nous avons un modèle numérique qui représente un véritable système physique et qui présente le chaos (par exemple, les modèles de dynamique des fluides, les modèles climatiques), comment pouvons-nous savoir que le modèle fonctionne comme il se doit? Nous ne pouvons pas comparer directement...
Je recherche la documentation et les références de diapositives les plus importantes pour modéliser les coûts énergétiques actuels et futurs des opérations en virgule flottante et des transferts de données à travers le CPU, la mémoire, le réseau et le stockage. J'ai marqué cette question comme un...
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...
maxf(x) subject to Ax=bmaxf(x) subject to Ax=b\max f(\mathbf{x}) \text{ subject to } \mathbf{Ax} = \mathbf{b} où , \ mathbf {x} = [x_1, x_2, ..., x_N] ^ T \ in \ mathbb {R} ^ {N \ times 1} , et \ mathbf {A} \ in \ mathbb {R} ^ {M \ fois
En lisant aujourd'hui sur la théorie des formes différentielles, j'ai été impressionné à quel point cela me rappelait la méthode du volume fini (FVM) de second ordre. J'ai du mal à comprendre est de penser que cette façon de faire est tout simplement trivial ou y a-t-il une connexion plus profonde....
Je vois très peu de bibliothèques / packages informatiques à virgule non flottante. Étant donné les diverses inexactitudes de la représentation en virgule flottante, la question se pose de savoir pourquoi il n'y a pas au moins certains domaines où cette précision accrue pourrait valoir la...
Je travaille sur une simulation fluide 2D utilisant des particules / vortons vortex comme décrit dans Simulation fluide pour les jeux vidéo . Ce que je pense, c'est la même chose que la "méthode du vortex discret". Fondamentalement, vous représentez le fluide avec une collection de particules avec...
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...
dans matlab, linsolve et mldivide sont utilisés pour résoudre un système d'équations linéaires, dans tous les cas déterminés, surdéterminés et sous-déterminés. En lisant leurs documents, je me demandais quelles sont les différences entre eux? Utilisent-ils presque les mêmes algorithmes de...
Fermé. Cette question est hors sujet . Il n'accepte pas actuellement les réponses. Voulez-vous améliorer cette question? Mettez à jour la question afin qu'elle soit sur le sujet pour Computational Science Stack Exchange. Fermé il y a 2 ans . Comme le dit le titre, je veux pouvoir arrêter le code à...
C ++ 11 introduit la sémantique de déplacement qui peut, par exemple, améliorer les performances du code dans les situations où C ++ 03 devrait effectuer une construction de copie ou une affectation de copie. Cet article signale que le code suivant connaît une accélération de 5 fois lorsqu'il est...
Nous connaissons tous les nombreuses méthodes de calcul pour résoudre le système linéaire standard Ax=b.Ax=b. Ax=b. Cependant, je suis curieux de savoir s'il existe des méthodes de calcul "standard" pour résoudre un système linéaire plus général (de dimension finie) de la forme LA=B,LA=B, LA=B, où,...
Comme tout le monde sait que le flux incompressible n'existe pas dans la réalité, c'est une hypothèse introduite pour simplifier les équations gouvernantes. Nous ne pouvons pas appliquer cette hypothèse directement. Généralement, le nombre de Mach (M <0,3 pour un écoulement incompressible), la...
Dans la méthode des solutions fabriquées (MMS), on postule une solution exacte, la substitue dans les équations et calcule le terme source correspondant. La solution est ensuite utilisée pour la vérification du code. Pour les équations de Navier-Stokes incompressibles, le MMS conduit facilement à...
J'utilise le schéma de différence finie Crank-Nicolson pour résoudre une équation de chaleur 1D. Je me demande si le principe maximum / minimum de l'équation de la chaleur (c'est-à-dire que le maximum / minimum se produit à la condition initiale ou aux limites) est également valable pour la...
Je cherche à porter du code qui résout un ensemble d'équations différentielles partielles (PDE) par la méthode des volumes finis sous forme IMPLICITE (pour la discrétisation temporelle). Il en résulte qu'il existe un système tridiagonal d'équations dans les directions x, y, z qui est géré par le...
Je suis nouveau dans le calcul scientifique. Je recherche une implémentation Fortran (de préférence f90) d'un Octree. Mon problème nécessite un octree qui divise mon domaine jusqu'à ce qu'il n'y ait plus que quelques particules N (ou sources où je connais la valeur de densité qui peut être branchée...
J'essaie de résoudre l'équation de Poisson 2D en utilisant la méthode de Galerkin discontinu (DG) et la discrétisation suivante (j'ai un fichier png mais je ne suis pas autorisé à le télécharger, désolé): Équation: ∇ ⋅ ( κ ∇ T) + f= 0∇⋅(κ∇T)+F=0\nabla \cdot( \kappa \nabla T) + f = 0 Nouvelles...
J'utilise la fonction de noyau RBF pour implémenter un algorithme d'apprentissage machine basé sur le noyau (KLPP), la matrice de noyau résultante KKK K(i,j)=exp(−(xi−xj)2σ2m)K(i,j)=exp(−(xi−xj)2σm2)K(i,j)= \exp\left({\frac{-(x_{i}-x_{j})^2}{ \sigma_{m}^2}}\right) est extrêmement mal...
Pour la réduction du modèle, je veux calculer les vecteurs singuliers de gauche associés aux - disons 20 - plus grandes valeurs singulières d'une matrice , où N ≈ 10 6 et k ≈ 10 3 . Malheureusement, ma matrice A sera dense sans aucune structure.A∈RN,kA∈RN,kA \in \mathbb R^{N,k}N≈106N≈106N\approx...