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...
Voici une utilisation typique des pointeurs de fonction en C. Je voudrais faire quelque chose de similaire dans Fortran. J'ai quelques idées, mais j'aimerais savoir s'il existe un moyen canonique de le faire. Les pointeurs de fonction et les contextes transmis par l'utilisateur sont stockés, puis...
Le calcul parallèle est fréquemment modélisé à l'aide d'un taux de calcul local déterministe, d'une surcharge de latence et d'une bande passante réseau. En réalité, ceux-ci sont spatialement variables et non déterministes. Des études telles que Skinner et Kramer (2005) observent des distributions...
Combien est plus lent petsc4py vs c / c ++ / fortran? Je me rends compte que cela dépendra considérablement du code exécuté, mais qu'en est-il de quelque chose de simple comme un produit
On me donne une matrice symétrique, inversible, définie positive et dense. Je dois tester si où J est la matrice des uns.12×1212×1212 \times 12QQQJdet(Q)=det(12I−Q−J)(1)det(Q)=det(12I−Q−J)(1)\det(Q) = \det(12I-Q-J) \; \; (1)JJJ Je fais actuellement cela avec la bibliothèque de tatou mais cela...
La compilation de PETSc avec une bibliothèque BLAS / LAPACK externe affecte-t-elle de manière significative les performances sur des matrices clairsemées, ou utilise-t-elle uniquement ces bibliothèques pour les mathématiques matricielles
Je voudrais utiliser la méthode Runge-Kutta de 8e ordre (89) dans une application de mécanique / astrodynamique céleste, écrite en C ++, à l'aide d'une machine Windows. Par conséquent, je me demande si quelqu'un connaît une bonne bibliothèque / implémentation documentée et gratuite à utiliser?...
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...
Il existe de nombreuses méthodes numériques bien connues pour résoudre des équations de type par exemple la méthode de la bissection, la méthode de Newton, etc.f(x)=0,x∈Rn,f(x)=0,x∈Rn, f(x) = 0, \quad x \in \mathbb{R}^n, Dans mon application, est calculé avec une méthode stochastique (le résultat...
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 suis nouveau dans le package PETSc. J'ai une matrice A ~ 4000x4000 au format matriciel et je veux que PETSc résolve cela en utilisant plusieurs processeurs. Je sais comment résoudre le système sur un seul processeur, mais je ne sais pas comment répartir la matrice et les vecteurs entre...
J'ai un modèle non négligeable (~ 5000 lignes) écrit en C. C'est un programme en série, sans génération de nombre aléatoire nulle part. Il utilise la bibliothèque FFTW pour les fonctions utilisant FFT - Je ne connais pas les détails de l'implémentation FFTW, mais je suppose que les fonctions y sont...
Au cours des dernières années, plusieurs bibliothèques / projets logiciels sont apparus qui offrent une forme ou une autre de parallélisme à mémoire partagée basé sur des données à usage général. L'idée principale est qu'au lieu d'écrire un code explicitement threadé, les programmeurs implémentent...
Je suis très nouveau dans la programmation GPGPU, alors pardonnez-moi si la question n'est pas particulièrement appropriée. D'après ce que je comprends, la programmation GPU est un travail d'ingénierie très complexe par rapport à la programmation CPU habituelle. Il faut faire très attention aux...
Fondamentalement, FEM semble être un problème à peu près «résolu». Il existe de nombreux frameworks puissants, comme Trilinos, PETSc, FEniCS, Libmesh ou MOOSE. Une chose qu'ils ont en commun: ils sont extrêmement "lourds". Tout d'abord, l'installation est normalement super douloureuse....
Je plonge dans le monde fascinant de l'analyse par éléments finis et je voudrais résoudre un gros problème thermo-mécanique (uniquement thermique mécanique à , pas de feedback).→→\rightarrow Pour le problème mécanique, j'ai déjà compris de la réponse de Geoff , que je devrai utiliser un solveur...
J'ai besoin de résoudre Ax = b, mais je me rends compte que même s'il est rare, le stockage des coefficients matriciels de mon problème prendra trop de mémoire. Alors maintenant, j'envisage d'utiliser une méthode sans matrice, car les mêmes coefficients apparaissent beaucoup de temps dans la...
La plupart de mes programmes sont des codes de recherche uniques en C pour mon propre usage. Je n'ai jamais distribué de code à des collaborateurs proches. J'ai développé un algorithme que je publie dans une revue scientifique. Je veux fournir le code source et peut-être le code exécutable dans le...
Comme l'illustre la réponse de Jed Brown aux coûts des recherches par rapport aux calculs , l'utilisation d'opérations à virgule flottante vectorisées vs non vectorisées résulte en un code beaucoup plus rapide. De nombreux compilateurs modernes affirment qu'ils peuvent effectuer une vectorisation...