Dans mon projet, je dois résoudre quelques matrices tridiagonales à chaque pas de temps, il est donc crucial d'avoir un bon solveur pour celles-ci. J'ai fait ma propre implémentation, juste la manière classique de le faire décrite sur Wikipedia. J'ai ensuite essayé d'utiliser Lapack à la place, et à ma …
J'ai fait une observation étrange en calculant la solution d'une simple équation de diffusion de réaction 1D: ∂∂ta=∂2∂x2a - a b∂∂ta=∂2∂x2a−ab\frac{\partial}{\partial t}a=\frac{\partial^2}{\partial x^2}a-ab ∂∂tb = - a b∂∂tb=−ab\frac{\partial}{\partial t}b=-ab ∂∂tc = a∂∂tc=a\frac{\partial}{\partial t}c = a La valeur initiale de est une constante ( b ( 0 , x ) = …
J'ai un maillage, des faces , des arêtes E et des sommets V , et j'ai une liste de contours de niveaux prédéfinis.FFFEEEVVV Quel algorithme puis-je utiliser pour construire les contours de la manière la plus efficace? Un tracé du contour est illustré ci-dessus. Les lignes de même couleur ont …
J'ai un système d'équations différentielles ordinaires - 7 équations et ~ 30 paramètres régissant leur comportement dans le cadre d'un modèle mathématique de transmission de maladies. Je voudrais trouver les états stables pour les équations Modification dx/dt = rest of the equationà 0 = equationchacun des équations en fait un …
Je voudrais comprendre plus de détails sur la mise en œuvre des méthodes multigrilles algébriques (AMG). J'ai lu "A Multigrid Tutorial", qui est assez bon et explique tous les détails de l'interpolation, de l'opérateur de grille grossière et de la sélection de grille grossière pour AMG. Cependant, je pense qu'il …
Il est montré (Yousef Saad, Méthodes itératives pour les systèmes linéaires clairsemés , p. 260) quec o n d( A′A ) ≈ c o n d( A )2cond(A′A)≈cond(A)2cond(A'A) \approx cond(A)^2 Est-ce également vrai pour les ?UNEA′AA′AA' Dans le cas où est avec , j'observe queN × M N ≪ M …
J'ai entendu dire qu'une transformée de Fourier rapide peut être utilisée pour résoudre le problème du poisson lorsque les conditions aux limites sont toutes d'un même type ... Série sinus pour dirichlet, cosinus pour neumann et les deux pour périodique. En considérant un domaine rectangulaire 2D, supposons que deux côtés …
Supposons que le système linéaire suivant soit donné où est le Laplacien pondéré connu pour être défini positif avec un espace nul unidimensionnel par , et la variance de translation de , c'est-à-dire que ne change pas la valeur de la fonction (dont la dérivée est ). Les seules entrées …
J'ai quelques questions concernant les points suivants: J'essaie de résoudre l'équation de Schrodinger en 1D en utilisant la discrétisation de manivelle nicolson suivie de l'inversion de la matrice tridiagonale résultante. Mon problème a maintenant évolué vers un problème avec des conditions aux limites périodiques et j'ai donc modifié mon code …
Par exemple, je voudrais calculer numériquement la norme de u = 1L2L2L^2 danscertain domaine qui comprend zéro, j'essayé quadratureGauss et d'échec, il estpeu loin de la vraieL2-norme sur la boule unité en coordonnées sphériques d'intégration, est existe-t-il un bon moyen de le faire? Ce problème est souvent observé dans les …
Des systèmes indéfinis de matrices apparaissent par exemple dans la discrétisation des problèmes de point de selle par des éléments finis mixtes. La matrice du système peut alors être mise sous la forme ( ABBtC)(UNEBtBC)\begin{pmatrix} A & B^t \\ B & C\end{pmatrix} où est négatif (semi) défini, C est positif …
Lors de la résolution de systèmes linéaires clairsemés à l'aide de méthodes de factorisation directe, la stratégie de classement utilisée a un impact significatif sur le facteur de remplissage des éléments non nuls dans les facteurs. Une telle stratégie d'ordonnancement est la dissection imbriquée. Je me demande s'il est possible …
J'essaye d'intégrer ∫10t2 n + 2exp( α r0t) dt∫01t2n+2exp(αr0t)rét\int^1_0 t^{2n+2}\exp\left({\frac{\alpha r_0}{t}}\right)dt qui est une simple transformation de ∫∞1X2 nexp( - α r0x ) dX∫1∞X2nexp(-αr0X)réX\int^{\infty}_1 x^{2n}\exp(-\alpha r_0 x)dx en utilisant car il est difficile d'approximer numériquement les intégrales incorrectes. Cela pose cependant le problème de l'évaluation de la nouvelle intégrale proche …
Existe-t-il de bons papiers et / ou codes qui couplent les solveurs à éléments finis Galerkin discontinus aux solveurs Riemann? J'ai besoin d'explorer le couplage des problèmes elliptiques et hyperboliques, mais la plupart des méthodes de fractionnement sont au mieux ad hoc. Comme j'ai une grande quantité de code FEniCS, …
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.