Je m'intéresse aux cas où le gradient conjugué fonctionne bien mieux que la méthode GMRES. En général, le CG est le choix préférable dans de nombreux cas de SPD (symétrique-positif-défini) car il nécessite moins de stockage et le taux théorique de convergence pour le CG est le double de celui du...
J'écris une petite bibliothèque pour des calculs matriciels clairsemés comme un moyen de m'apprendre à faire le meilleur usage de la programmation orientée objet. J'ai travaillé très dur pour avoir un joli modèle objet, où les parties (matrices clairsemées et les graphiques qui décrivent leur...
Une réduction parallèle suppose que l'opération correspondante est associative. Cette hypothèse est violée pour l'ajout de nombres à virgule flottante. Vous pourriez vous demander pourquoi je m'en soucie. Eh bien, cela rend les résultats moins reproductibles. Et cela empire lorsque le recuit simulé...
Le développement d'algorithmes scientifiques est un processus hautement itératif impliquant souvent la modification de nombreux paramètres que je souhaiterais faire varier soit dans le cadre de ma conception expérimentale, soit dans le cadre de l'amélioration des performances de l'algorithme....
Une grande partie de mon propre travail consiste à améliorer la mise à l'échelle des algorithmes, et l'une des façons préférées de montrer la mise à l'échelle parallèle et / ou l'efficacité parallèle est de tracer les performances d'un algorithme / code sur le nombre de cœurs, par exemple où l' axe...
L'autre jour, mon professeur de dynamique des fluides numérique était absent et il a envoyé son doctorant pour le remplacer. Dans la conférence qu'il a donnée, il a semblé indiquer plusieurs inconvénients associés à divers schémas de discrétisation pour les simulations d'écoulement de fluide:...
Ma thèse porte sur le développement de méthodes numériques pour la réduction des modèles de combustion. J'exécute mes méthodes uniquement sur la partie chimie des simulations de combustion, et j'ai de nombreuses études de cas pour les simulations 0-D (sans débit). Ce que j'aimerais, c'est exécuter...
Je prévois d'utiliser le multigrille pour calculer certaines valeurs et vecteurs propres, et j'ai remarqué que PETSc a un support de haut niveau pour le multigrille. La documentation PETSc indique que cette partie de PETSc ne doit pas être utilisée, car elle sera bientôt remplacée. Quelles autres...
J'ai besoin d'écrire quelques données à partir d'un calcul, qui seront lues plus tard par Paraview (fichier .vtu ou vtk). En ce qui concerne la taille du fichier, dois-je opter pour le format ASCII ou le format
Supposons que j'ai une fonction FFf et que je veuille trouver XXx tel que F( x ) ≈ 0F(X)≈0f(x)\approx 0 . Je pourrais utiliser la méthode Newton-Raphson. Mais cela nécessite que je connaisse la fonction dérivée F′( x )F′(X)f'(x) . Une expression analytique pour FFf peut ne pas être disponible. Par...
À l'aube d'OS X, il semblait y avoir beaucoup de brouhaha, du moins dans le monde Mac (j'étais loin d'être du calcul scientifique à l'époque) à propos du Mac OS en tant que plate-forme pour le calcul scientifique et les applications HPC. XGrid est sorti de la boîte, Virginia Tech avait son cluster...
Le langage C ++ fournit une programmation générique et une métaprogrammation via des modèles. Ces techniques ont trouvé leur chemin dans de nombreux progiciels de calcul scientifique à grande échelle (par exemple, MPQC , LAMMPS , CGAL , Trilinos ). Mais qu'est-ce qu'ils ont réellement contribué au...
Lequel est le meilleur: FORTRAN ou Python? Et je suppose que dans les deux cas, vous avez besoin de Gnuplot, ai-je raison? Je travaille actuellement sur une machine Windows. J'aimerais l'utiliser pour obtenir des solutions numériques pour des problèmes de physique, y compris des simulations...
Je recherche un logiciel de génération de maillage qui est gratuit et open source, fournit une interface de script saine pour la spécification de domaine, travaille pour des géométries complexes, peut générer des maillages 2D et 3D, Quelles options
Je suis un étudiant en génie mécanique intéressé par le domaine de l'ingénierie aérospatiale où, me dit-on, Fortran est encore couramment utilisé. Quelle version de Fortran dois-je investir mon temps pour
Je sais que la méthode de Nitsche est une méthode très intéressante car elle permet de prendre en compte les conditions aux limites de type Dirichlet ou le contact avec les conditions aux limites de frottement de façon faible sans utiliser de multiplicateurs de Lagrange. Et son avantage, qui est de...
Considérons un IVP pour le système ODE , y ( x 0 ) = y 0 . Le plus souvent, ce problème est considéré comme rigide lorsque la matrice de Jacobi ∂ fy′= f( x , y)y′=F(X,y)y'=f(x,y)y( x0) = y0y(X0)=y0y(x_0)=y_0a à lafois desvaleurs propres avec une très grande partie réelle négative et des valeurs...
Dans la littérature sur les méthodes numériques, de nombreux articles de recherche consistent en une description d'une nouvelle variation algorithmique, suivie de quelques problèmes de test comparant la nouvelle méthode avec une ou deux méthodes existantes. Il est donc difficile de déterminer...
Je me demandais s'il existe une méthode rapide et efficace pour trouver à l'avance le nombre de non-zéros pour une opération de multiplication matricielle clairsemée en supposant que les deux matrices sont au format CSC ou CSR. Je sais qu'il y en a un dans le paquet smmp mais j'ai besoin de quelque...
Il existe deux approches générales pour représenter des solutions dans la méthode de Galerkin discontinue: nodale et modale. Modal : les solutions sont représentées par des sommes de coefficients modaux multipliées par un ensemble de polynômes, par exemple où est généralement des polynômes...