Je comprends donc (ou du moins je crois que je comprends) comment fonctionne un cycle en V. J'ai écrit dans Matlab la version 1-D, récursive d'un V-cycle. Cependant, lorsque j'ai exécuté mon code pour FMG, ma solution ne convergeait pas. Je crois que mon problème réside dans ma compréhension de la...
Je dois résoudre le même système linéaire clairsemé (300x300 à 1000x1000) avec de nombreux côtés droits (300 à 1000). En plus de ce premier problème, je voudrais également résoudre différents systèmes, mais avec les mêmes éléments non nuls (juste des valeurs différentes), c'est-à-dire de nombreux...
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...
J'ai une collection de modèles informatiques qui pourraient être décrits comme des automates cellulaires asynchrones. Ces modèles ressemblent au modèle Ising, mais sont légèrement plus compliqués. Il semble que de tels modèles gagneraient à être exécutés sur un GPU plutôt que sur un CPU....
Je porte un code existant de MATLAB vers C ++ et j'ai un système linéaire pour résoudre (plutôt que la forme la plus typique A x = b )x A = bxA=bxA=bA x = bAx=bAx=b La matrice est dense et de forme générale, mais ne dépasse pas 1000x1000. Donc, dans MATLAB, la solution est trouvée par la fonction...
Je simule un écoulement incompressible sur un cylindre au nombre de Reynold de 500. Je résous l'équation de Navier Stokes en utilisant la méthode de correction de pression. Ma solution devient instable après un certain temps (environ 5 s). J'ai essayé d'affiner mon maillage, stepize (0,05) (en...
Il existe plusieurs familles de méthodes pour la quadrature numérique. Si j'ai une classe spécifique d'intégrandes, comment choisir la méthode idéale? Quelles sont les questions pertinentes à poser à la fois sur l'intégrande (par exemple, est-elle fluide? A-t-elle des singularités?) Et sur le...
Supposons que j'ai eu un problème de valeur limite: dud2udx2+dvdx=f in Ωd2udx2+dvdx=f in Ω\frac{d^2u}{dx^2} + \frac{dv}{dx}=f \text{ in } \Omega u=h en ∂Ωdudx+d2vdx2=g in Ωdudx+d2vdx2=g in Ω\frac{du}{dx} +\frac{d^2v}{dx^2} =g \text{ in } \Omega u=h in ∂Ωu=h in ∂Ωu=h \text{ in } \partial\Omega Mon...
Lors de la discrétisation FEM et de la résolution d'un problème de réaction-diffusion, par exemple, avec 0 < ε ≪ 1 (perturbation singulière), la solution du problème discret présentera typiquement des couches oscillatoires proches de la frontière. Avec- ε Δ u + u = 1 sur Ωu = 0 sur ...
Mon système est un problème FE symétrique avec des multiplicateurs de lagrange (par exemple le flux de Stokes incompressible): ( ABBTC)(ABTBC)\begin{pmatrix}A & B^T \\ B & C\end{pmatrix} où est le cas typique (je me suis même assuré que les équations sont numérotées pour que les multiplicateurs de...
Salut et merci d'avoir donné le temps de regarder ma question. Ceci est une version mise à jour de ma question que j'ai publiée plus tôt dans physics.stackexchange.com J'étudie actuellement un condensateur à excitons 2D Bose-Einstein Condensate et je suis curieux de connaître l'état fondamental de...
J'essaie de comprendre s'il existe un moyen plus rapide de calculer toutes les valeurs propres et vecteurs propres d'une matrice de contiguïté très grande et clairsemée que d'utiliser scipy.sparse.linalg.eigsh Pour autant que je sache, cette méthode utilise uniquement la rareté et attributs de...
Dans les équations de Navier-Stokes incompressibles, ρ(ut+(u⋅∇)u)∇⋅u=−∇p+μΔu+f=0ρ(ut+(u⋅∇)u)=−∇p+μΔu+f∇⋅u=0\begin{align*} \rho\left(\mathbf{u}_t + (\mathbf{u} \cdot \nabla)\mathbf{u}\right) &= - \nabla p + \mu\Delta\mathbf{u} + \mathbf{f}\\ \nabla\cdot\mathbf{u} &= 0 \end{align*} le terme de...
Je souhaite maximiser globalement une fonction de nombreux ( ) paramètres réels (résultat d'une simulation complexe). Cependant, la fonction en question est relativement coûteuse à évaluer, nécessitant environ 2 jours pour chaque ensemble de paramètres. Je compare différentes options et je me...
Nous travaillons sur un modèle bayésien pour un processus spatio-temporel, et utilisons un échantillonneur No-U-Turn (NUTS) qui nécessite un modèle pour la probabilité logarithmique et son gradient par rapport aux paramètres du modèle. Plus succinctement, nous avons une fonction de probabilité de...
Soit et être une fonction dans ces variables.f( → x ):[0,1]n→CX⃗ = ( x1, x2, … , Xn) ∈ [ 0 , 1 ]nX→=(X1,X2,…,Xn)∈[0,1]n\vec{x} = (x_1, x_2, \dots, x_n) \in [0,1]^nF( x⃗ ) : [ 0 , 1 ]n→ CF(X→):[0,1]n→Cf(\vec{x}): [0,1]^n \to \mathbb{C} Existe-t-il un schéma récursif pour cette intégrale itérée? ∫[ 0...
Supposons que j'ai des versions lentes et rapides d'un code et que je souhaite signaler un numéro d'accélération comparant les deux. Je lance la version lente fois et la version rapide fois, produisant des temps et . La façon la plus simple de produire une accélération est de faire la moyenne des...
Cela pourrait être une question mieux adaptée au côté des recommandations de logiciels de SE, mais je pense que les personnes qui fréquentent cette partie de SE sont plus susceptibles de pouvoir répondre à cette question. Je recherche une alternative gratuite (pas seulement en liberté) à Comsol...
Si je devais résoudre un problème relativement petit, c'est-à-dire un problème qui peut être traité par une méthode directe comme LU, le numéro de condition de l'opérateur linéaire affecte-t-il la précision de la solution? L'un des problèmes de recherche sur lesquels je travaille se concentre sur...
Je viens de lire le chapitre 3 dans "A Multigrid Tutorial" de Briggs / Henson / McCormick, lien . Le texte concerne les cycles multigrilles tels que le cycle V, le cycle mu, le FMG. Ce qui a attiré mon attention: dans la plupart des procédures itératives, on vérifie si elle a convergé vers la...