Dans l'article http://www.sciencedirect.com/science/article/pii/S0045782509003521 , une équation locale d'élément HDG est décrite à la page 584 équation (4), l'une des équations prenant la forme suivante - ( uh, ∇ q)K= - ⟨ u^h⋅ n , q- q¯⟩∂K−(uh,∇q)K=−⟨u^h⋅n,q−q¯⟩∂K-(u_h,\nabla q)_K =...
Dans la recherche sur les problèmes inverses, il est courant de construire un ensemble de données synthétiques à partir d'un ensemble connu de paramètres, puis de tester si la technique d'inversion peut reconstruire ces paramètres. Pour ce faire, il est important d'ajouter des niveaux appropriés de...
Soit des matrices réelles, carrées et denses. G et Q sont symétriques. LaisserA , G , QA,G,QA, G, QgGGQQQ H= [ A- Q- G- unT]H=[A−G−Q−AT]H = \begin{bmatrix} A & -G \\ -Q &-A^T \end{bmatrix} être une matrice hamiltonienne. Je veux calculer la matrice exponentielle de . J'ai besoin de l'exponentielle...
Chaque matrice réelle peut être réduite à la forme réelle de Schur T = U T A U en utilisant une transformation orthogonale similiary U . Ici, la matrice T est de forme quasi triangulaire avec 1 par 1 ou 2 par 2 blocs sur la diagonale principale. Chaque 1 de 1 correspond à un bloc de valeur propre...
Les algorithmes existants pour résoudre les ODE gèrent les fonctions , oùy∈Rn. Mais dans de nombreux systèmes physiques, l'équation différentielle est autonome, doncdyréyrét= f( y, t )dydt=f(y,t)\frac{dy}{dt} = f(y, t)y∈ Rny∈Rny \in \mathbb R^n,y∈Rn, sanst. Avec cette hypothèse simplificatrice,...
J'ai constaté que beaucoup de mes programmes de science informatique ont des exigences de test qui ne sont pas couvertes par les cadres de test standard: Test de temps de calcul Pour vous assurer que les algorithmes ne ralentissent pas. Je pourrais faire quelque chose comme
Je résous des équations de Navier-Stokes compressibles multi-espèces sur une grille structurée 3D. J'ai obtenu une solution sur une grille donnée (disons une relativement grossière). Je veux maintenant affiner ma grille et interpoler ma solution précédente sur ma nouvelle grille avant de redémarrer...
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...
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....
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...
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...
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...
J'ai entendu dire que certaines revues sont mieux notées que d'autres. Est-ce vrai? Et si oui, quels sont les critères pour juger de la valeur d'une revue évaluée par des pairs par rapport à une autre? Comment connaître sa note? Ma publication aura-t-elle moins de "valeur" si elle est acceptée dans...
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,...
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...
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,...
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 la...
Existe-t-il un format standard pour les maillages par éléments finis qui est largement utilisé dans l'industrie?
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 positives de...
Soit et deg f > deg g . Je recherche des algorithmes asymptotiquement rapides et numériquement stables pour calculer f mod g . Dans les applications envisagées, les deux f , g sont des polynômes denses à coefficients à virgule flottante double précision. Mais, pour l'instant, je m'intéresse plus...