Science computationnelle

3

Pourquoi la dimension temporelle est-elle spéciale?

De manière générale, j'ai entendu des analystes numériques exprimer l'opinion que "Bien sûr, mathématiquement parlant, le temps n'est qu'une autre dimension, mais quand même, le temps est spécial" Comment justifier cela? Dans quel sens le temps est-il spécial pour la science informatique? De plus, pourquoi préférons-nous si souvent utiliser des …

24 pde discretization

5

Quelles sont les principales différences entre PETSc et Trilinos?

Pour autant que je sache, les deux grands cadres génériques des logiciels de science informatique du Département américain de l'énergie sont PETSc et Trilinos . Ils semblent similaires à première vue, au-delà des différences de langage (C versus C ++). Quelles sont les principales différences entre les deux cadres et …

24 petsc software trilinos

8

Quel logiciel est bon à utiliser pour le débogage parallèle?

Je n'exécute aucun code parallèle pour le moment, mais je prévois d'exécuter du code parallèle à l'avenir en utilisant un hybride d'OpenMP et de MPI. Les débogueurs ont été des outils précieux pour moi lors de l'exécution de projets en série. Quelqu'un peut-il recommander un débogueur parallèle (ou plusieurs débogueurs) …

24 software parallel-computing

4

Comment ajouter de gros termes exponentiels de manière fiable sans erreurs de débordement?

Un problème très courant dans Markov Chain Monte Carlo implique le calcul de probabilités qui sont la somme de grands termes exponentiels, eune1+ eune2+ . . .eune1+eune2+... e^{a_1} + e^{a_2} + ... où les composants d' peut varier de très petit à très grand. Mon approche a été de factoriser …

24 monte-carlo floating-point statistics

5

Pourquoi les points équidistants se comportent-ils mal?

Description de l'expérience: Dans l'interpolation de Lagrange, l'équation exacte est échantillonnée en NNN points (ordre polynomial N−1N−1N - 1 ) et interpolée en 101 points. Ici, NNN varie de 2 à 64. Chaque fois que L1L1L_1 , L2L2L_2 et L∞L∞L_\infty des tracés d'erreur sont préparés. On voit que, lorsque la …

24 finite-difference interpolation error-estimation polynomials

3

Quel est le principe derrière la convergence des méthodes du sous-espace de Krylov pour résoudre des systèmes d'équations linéaires?

Si je comprends bien, il existe deux grandes catégories de méthodes itératives pour résoudre des systèmes linéaires d'équations: Méthodes stationnaires (Jacobi, Gauss-Seidel, SOR, Multigrid) Méthodes Krylov Subspace (Gradient Conjugué, GMRES, etc.) Je comprends que la plupart des méthodes stationnaires fonctionnent en relaxant de manière itérative (lissage) les modes de Fourier …

24 linear-algebra convergence krylov-method conjugate-gradient

4

Quand une méthode d'ordre élevé est-elle utile pour les simulations de dynamique des fluides computationnelles?

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 donné?

23 fluid-dynamics accuracy

5

Quel matériel dois-je inclure dans un article de revue (ou publier en ligne) afin de rendre mes recherches informatiques reproductibles?

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 devrais inclure avec un …

23 software publications documentation reproducibility

5

Quand est-il moins cher de créer un cluster dans le cloud que d'en créer un dans mon laboratoire pour des simulations MD?

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 …

23 cloud-computing hpc

3

Quel est l'état de l'art du calcul intégral hautement oscillatoire?

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 arbitraire?

23 quadrature precision numerical-analysis

1

Existe-t-il un algorithme numérique pour trouver une pente asymptotique?

( 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 …

23 algorithms

12

Est-il possible d'utiliser Octave pour apprendre la programmation MATLAB?

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 …

23 software matlab languages octave

3

Utilisation de langages de programmation non conventionnels pour le calcul scientifique [fermé]

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 discussion approfondie. …

22 languages

10

Quel algorithme est le plus précis pour calculer la somme d'un tableau trié de nombres?

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} ; …

22 algorithms floating-point

1

Pourquoi la méthode de Newton ne converge-t-elle pas?

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 équations?

22 petsc pde implicit-methods