Quels seraient, selon vous, les avantages / inconvénients de deux approches de codage d'une bibliothèque générale (volumes finis, fem, dg) pour la mécanique du continu informatique? C'est comme ça que je vois les choses en ce moment, alors s'il vous plaît fournissez vos propres expériences et ne me faites pas …
Lors de l'évaluation du nombre de FLOP dans une fonction simple, on peut souvent simplement descendre l'expression en comptant les opérateurs arithmétiques de base. Cependant, dans le cas d'instructions mathématiques impliquant une division paire, on ne peut pas le faire et s'attendre à pouvoir comparer avec les nombres FLOP à …
Les termes sources, tels que ceux dus à la bathymétrie dans les équations des eaux peu profondes, doivent être intégrés de manière spéciale afin de préserver les états physiques stables. Existe-t-il un moyen général de construire des méthodes bien équilibrées ou nécessite-t-il des techniques spéciales pour chaque équation?
Je sais que les gens remplacent souvent des matrices de masse cohérentes par des matrices diagonales groupées. Dans le passé, j'ai également implémenté un code dans lequel le vecteur de charge est assemblé de manière groupée plutôt que de manière cohérente FEM. Mais je n'ai jamais cherché à savoir pourquoi …
Je recherche des ressources sur les modèles de conception et les principes des logiciels numériques, potentiellement mais pas nécessairement en mettant l'accent sur les approches orientées objet des codes numériques. Je connais le texte, Writing Scientific Software: A Guide for Good Style , mais cela traite encore beaucoup de choses …
Dans une méthode comme le gmres ou le bicgstab, il pourrait être intéressant d'utiliser une autre méthode krylov comme préconditionneur. Après tout, ils sont faciles à mettre en œuvre sans matrice et dans un environnement parallèle. Par exemple, on peut utiliser quelques (disons ~ 5) itérations de bigcstab non conditionné …
Je ne connais pas le format de fichier utilisé dans HDF5, mais je me demande si les fichiers HDF5 conviennent au contrôle de révision avec git (ou par exemple Mercurial ou Subversion)? Je suppose que ce que je veux dire est: les fichiers HDF5 sont-ils adaptés à la diff-line en …
Une équation de Poisson avec toutes les conditions aux limites de Neumann a un seul espace nul dimensionnel constant. Lors de la résolution via une méthode de Krylov, l'espace nul peut être supprimé soit en soustrayant la moyenne de la solution à chaque itération, soit en épinglant la valeur d'un …
Dans un système où l'énergie devrait théoriquement être conservée, la simulation la plus précise permettrait de conserver l'énergie (tout en donnant des positions, des vitesses, etc. précises). Le RK4 est plus précis que le saute-mouton, mais le saute-mouton économise l'énergie et pas le RK4. Pourquoi est-ce?
J'ai une question concernant l'ajustement quadrique à un ensemble de points et de normales correspondantes (ou de manière équivalente, des tangentes). L'ajustement de surfaces quadriques à des données ponctuelles est bien exploré. Certains travaux sont les suivants: Ajustement direct contraint par type de surfaces quadriques , James Andrews, Carlo H. …
La forme forte d'une PDE nécessite que la solution inconnue appartienne à . Mais la forme faible nécessite seulement que la solution inconnue appartienne à .H2H2H^2H1H1H^1 Comment conciliez-vous cela?
J'essaie de comprendre un peu mieux l'équation d'advection avec un coefficient de vitesse variable. En particulier, je ne comprends pas comment l'équation peut être conservatrice. L' équation d'advection , ∂u∂t+ ∂∂X( v u ) = 0∂u∂t+∂∂x(vu)=0 \frac{\partial u}{\partial t} + \frac{\partial}{\partial x}(\boldsymbol{v}u) = 0 Interprétons comme étant la concentration de …
J'ai lu que si, par exemple, j'ai une double forboucle qui court sur les index d'une matrice, alors mettre l'index courant de la colonne dans la boucle externe est plus efficace. Par exemple: a=zeros(1000); for j=1:1000 for i=1:1000 a(i,j)=1; end end Quelle est la manière la plus efficace de la …
Je travaille sur un problème inverse pour mon doctorat. la recherche, qui pour simplifier, nous dirons est de déterminer dansββ\beta L(β)u≡−∇⋅(k0eβ∇u)=fL(β)u≡−∇⋅(k0eβ∇u)=fL(\beta)u \equiv -\nabla\cdot(k_0e^\beta\nabla u) = f à partir de quelques observations ; k 0 est une constante et f est connu. Ceci est généralement formulé comme un problème d'optimisation pour …
Seules les deux premières sections de cette longue question sont essentielles. Les autres sont juste à titre d'illustration. Contexte Des quadratures avancées telles que Newton – Cotes composites de degré supérieur, Gauß – Legendre et Romberg semblent être principalement destinées aux cas où l'on peut finement échantillonner la fonction sans …
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.