De nombreuses approches numériques de la CFD peuvent être étendues à un ordre arbitrairement élevé (par exemple, les méthodes Galerkin discontinues, les méthodes WENO, la différenciation spectrale, etc.). Comment choisir un ordre de précision approprié pour un problème...
La reproductibilité est devenue de plus en plus importante dans la recherche en science informatique. (Par exemple, consultez cet article de Roger Peng dans Science ; je connais également d'autres articles et sites Web de ce type.) Cependant, je ne sais pas exactement combien d'informations je...
Quel est l'état de l'art dans l'approximation d'intégrales hautement oscillatoires à la fois dans une dimension et dans des dimensions supérieures avec une précision
( xje, yje)(Xje,yje)(x_i,y_i)y( x )y(X)y(x)XXxF( x ) ≡ y( x ) - ( a x + b )F(X)≡y(X)-(uneX+b)f(x) \equiv y(x) - (ax + b)x → ∞X→∞x \to \inftyF′( x )F′(X)f'(x)F′ ′( x )F″(X)f''(x), etc. Mais je ne sais pas quelle est la forme fonctionnelle de , si elle en a même une qui peut être décrite en termes de...
Cette question a été migrée à partir de Cross Validated car il est possible d'y répondre sur Computational Science Stack Exchange. Migré il y a 7 ans . Je veux apprendre la programmation MATLAB afin de pouvoir effectuer des recherches / analyses par moi-même et aussi, afin de pouvoir étudier /...
Une Amazon EC2 Compute les coûts de cluster environ $ 800- $ 1000 (selon le cycle de service) par cœur de processeur physique au cours de 3 ans. Lors de notre dernier cycle d'acquisition de matériel, mon laboratoire a récupéré 48 cœurs de matériel très similaire à celui des clusters d'Amazon pour...
Comme nous le savons, les méthodes de transformations orthogonales (rotations de Givens et réflexions de Housholder) pour les systèmes d'équations linéaires sont plus coûteuses que l'élimination gaussienne, mais ont théoriquement de plus bonnes propriétés de stabilité dans le sens où elles ne...
J'utilise le package de résolution non linéaire SNES de PETSc pour résoudre un système d'équations non linéaires obtenues en discrétisant une équation différentielle partielle. Comment puis-je déterminer pourquoi le solveur ne converge pas et que puis-je faire pour résoudre avec succès mes...
Etant donné une suite finie croissante de nombres positifs . Lequel des deux algorithmes suivants est le meilleur pour calculer la somme des nombres?z1, z2, . . . . . znz1,z2,.....znz_{1} ,z_{2},.....z_{n} s=0; for \ i=1:n s=s + z_{i} ; end Ou: s=0; for \ i=1:n s=s + z_{n-i+1} ; end À mon avis, il...
Aller enseigner aux étudiants de premier cycle un cours intitulé Introduction à la programmation informatique. Je suis un peu confus. En physique computationnelle, les scientifiques utilisent C / C ++ ou Python ou Fortran, CUDA etc ..... il est temps de construire leur base. Que dois-je utiliser?...
Quelle est l'approche préférée et efficace pour interpoler des données multidimensionnelles? Choses qui m'inquiètent: performances et mémoire pour la construction, évaluation unique / batch dimensions de manipulation de 1 à 6 linéaire ou d'ordre supérieur pouvoir obtenir des gradients (sinon...
Dans l'état actuel des choses, cette question ne convient pas à notre format de questions / réponses. Nous nous attendons à ce que les réponses soient étayées par des faits, des références ou une expertise, mais cette question suscitera probablement un débat, des arguments, des sondages ou une...
Le paradigme de programmation OpenCL promet d'être un standard ouvert sans redevance pour l'informatique hétérogène. Faut-il investir notre temps dans le développement de logiciels basés sur OpenCL? Avantages /
J'ai une fonction bidimensionnelle dont je voudrais échantillonner les valeurs. La fonction est très coûteuse à calculer et elle a une forme complexe, j'ai donc besoin de trouver un moyen d'obtenir le plus d'informations sur sa forme en utilisant le moins de points d'échantillonnage.F( x ,...
J'ai des matrices et . est clairsemé et est avec très grand (peut être de l'ordre de plusieurs millions.) est une matrice haut avec plutôt petit ( ) et chaque colonne peut seulement un seul entrée avec le reste étant « s, de telle sorte que . est énorme, il est donc très difficile à inverser, et je...
Quelle serait une bonne discrétisation par différence finie pour l'équation suivante: ∂ρ∂t+∇⋅(ρu)=0∂ρ∂t+∇⋅(ρu)=0\frac{\partial \rho}{\partial t} + \nabla \cdot \left(\rho u\right)=0 ? Nous pouvons prendre le cas 1D: ∂ρ∂t+ddx(ρu)=0∂ρ∂t+ddx(ρu)=0\frac{\partial \rho}{\partial t} +...
Je recherche des livres ou des articles, ou des articles de blog, ou tout matériel publié en général, qui traitent spécifiquement des utilisations des fonctionnalités modernes C ++ (déplacer la sémantique, la STL, les itérateurs, l'évaluation paresseuse, etc.) dans le calcul scientifique....
Je recherche des algorithmes pour dessiner des graphiques 2D standard pour des fonctions qui peuvent ou non avoir des singularités. Le but est d'écrire un "Mini-CAS", donc je n'ai pas de connaissance a priori sur les types de fonctions, les utilisateurs veulent faire un graphique. Ce problème est...
Je vais commencer par mon expérience personnelle dans notre laboratoire. Dans les ifort 9 et 10 jours, nous avions l'habitude d'être assez agressifs avec les optimisations, en compilant avec -O3 et des drapeaux spécifiques au processeur (-xW -xSSE4.2 par exemple). Mais à partir d'ifort 11, nous...
Je cherche à résoudre un problème d'optimisation contraint où je connais les limites de certaines variables (en particulier une contrainte encadrée). argminuf(u,x)argminuf(u,x) \arg \min_u f(u,x) sujet à c(u,x)=0c(u,x)=0 c(u,x) = 0 a≤d(u,x)≤ba≤d(u,x)≤b a \le d(u,x) \le b où uuu est un vecteur de...