J'essaie de résoudre un système linéaire d'équations qui est clairsemé, mais qui manque de tout type de structure en bandes. J'ai entendu dire qu'il existe un moyen d'étendre les principes d'un solveur multigrille pour les schémas de différences finies implicites à un problème linéaire général (si je ne me trompe …
Après quelques étapes de mise à jour bayésienne, il me reste une distribution postérieure de la forme d'un mélange de distributions normales, C'est-à-dire que le paramètre θ est tiré d'une distribution dont le PDF est donné comme un mélange pondéré de PDF normaux, et n'est pas une somme de RV …
Si je comprends bien, la relaxation successive fonctionne en choisissant un paramètre 0 ≤ ω ≤ 20≤ω≤20\leq\omega\leq2 et en utilisant une combinaison linéaire d'une itération (quasi) de Gauss-Seidel et de la valeur au pas de temps précédent ... c'est-à-dire uk + 1= ( ω ) ugsk + 1+ ( 1 …
J'ai une ODE: u ( 0 ) = - 1u′=−1000u+sin(t)u′=−1000u+sin(t)u'=-1000u+sin(t) u(0)=−11000001u(0)=−11000001u(0)=-\frac{1}{1000001} Je sais que cette ODE particulière est rigide, analytiquement. Je sais également que si nous utilisons une méthode de pas de temps explicite (vers l'avant) (Euler, Runge-Kutta, Adams, etc.), la méthode devrait renvoyer de très grandes erreurs si le …
Quelles situations d'analyse numérique deviennent plus / moins stables, ont une convergence plus rapide / plus lente, ou sont autrement assez différentes lorsqu'il s'agit de fonctions de variable complexe au lieu de fonctions d'une variable réelle?
On m'a toujours dit que le diagramme de Voronoi est le double du problème de triangulation de Delaunay. Dans quel sens peuvent-ils être des duels les uns des autres? Je pensais que les problèmes doubles (c'est-à-dire en programmation linéaire) sont censés produire la même réponse. De toute évidence, les deux …
Il existe au moins une encyclopédie assez complète de règles de quadrature qui ne semble pas avoir été mise à jour depuis longtemps et dont l'accès est restreint. Cette source fait référence à plusieurs sources classiques et modernes, et est généralement bien constituée. Cependant, il aborde la construction de règles …
Est-il possible d'obtenir une précision de second ordre pour les éléments finis hexaédriques avec moins de 8 points de Gauss sans introduire de modes non physiques? Un seul point central de Gauss introduit un mode de cisaillement non physique, et l'agencement symétrique standard de 8 points de Gauss est coûteux …
Divers ensembles importants de polynômes (Legendre, Chebyshev, etc.) sont orthogonaux sur un certain intervalle réel avec une certaine pondération. Existe-t-il des familles connues de polynômes orthogonaux par rapport à d'autres courbes dans le plan complexe? Par exemple, je voudrais une base pour les polynômes de degré n qui sont orthogonaux …
En ce moment, j'ai un code qui utilise l'algorithme Crank-Nicholson, mais je pense que je voudrais passer à un algorithme d'ordre supérieur pour l'horodatage. Je sais que l'algorithme de Crank-Nicholson est stable dans le domaine dans lequel je veux travailler, mais je crains que certains autres algorithmes ne le soient …
Lequel est calculé plus rapidement, ou ou ? , et sont des réels positifs avec .log a c b √ababa^blogaclogac\log_a c abcb>1c√bcb\sqrt[b]{c}aaabbbcccb>1b>1b>1 Quels types d'algorithmes utiliserez-vous dans la comparaison? Quelles sont leurs complexités? Par exemple, lorsque ou c ≈ a bc≡abc≡abc \equiv a^bc≈abc≈abc \approx a^b Cette question a été inspirée …
J'ai utilisé le SVD d'Intel MKL ( dgesvdvia SciPy) et j'ai remarqué que les résultats sont considérablement différents lorsque je change de précision entre float32et float64lorsque ma matrice est mal conditionnée / n'est pas pleine. Existe-t-il un guide sur le montant minimum de régularisation à ajouter pour rendre les résultats …
Je suis nouveau en science informatique et j'ai déjà appris des méthodes de base pour l'intégration, l'interpolation, des méthodes comme RK4, Numerov etc. sur c ++ mais récemment mon professeur m'a demandé d'apprendre à utiliser LAPACK pour résoudre des problèmes liés aux matrices. Comme par exemple trouver des valeurs propres …
J'ai besoin de calculer beaucoup d' inverses matricielles (pour la décomposition polaire d'itération de Newton), avec un très petit nombre de cas dégénérés ( ).3×33×33\times3<0.1%<0.1%<0.1\% L'inverse explicite (via les matrices mineures divisées par le déterminant) semble fonctionner, et est d'environ ~ 32 ~ 40 flops fusionnés (selon la façon dont …
Je suis à la recherche de références crédibles indiquant combien de ressources les superordinateurs dépensent pour la coordination par rapport au travail réel lié aux tâches. Les ressources pourraient être une puissance de traitement disponible, mais même les watts semblent être une unité valide. Je crois qu'un de mes professeurs …
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.