Science computationnelle

Questions-réponses pour les scientifiques utilisant des ordinateurs pour résoudre des problèmes scientifiques

3
Calcul d'une série légèrement oscillatoire à haute précision?
Supposons que j'ai la fonction intéressante suivante: f(x)=∑k≥1coskxk2(2−coskx).F(X)=∑k≥1cos⁡kXk2(2-cos⁡kX). f(x) = \sum_{k\geq1} \frac{\cos k x}{k^2(2-\cos kx)}. Il a des propriétés désagréables, comme sa dérivée qui n'est pas continue à des multiples rationnels deππ\pi . Je soupçonne qu'il n'existe pas de formulaire fermé. Je peux le calculer en calculant des sommes partielles …




5
Calculer
La fonction a une singularité proche de . Cette singularité peut être levée, cependant: pour , on devrait avoir , car Et donc Cependant, la forme n'est pas seulement non définie en , il est également numériquement instable au voisinage de ce point; afin d'évaluer pour de très petits numériquement, …
13 c++  c 


4
Tester des suites pour des applications numériques en C ++?
Récemment, j'ai poussé mon groupe à inclure plus de tests lors de l'écriture de leur code. Il y avait plusieurs bugs majeurs qui ont pris beaucoup plus de temps à attraper que probablement parler était nécessaire, car nous n'avions pas un bon régime de test en place. Cependant, je pense …
13 testing 

3
Quels sont les principes de base derrière la génération d'un maillage mobile?
Je m'intéresse à l'implémentation d'un maillage mobile pour un problème d'advection-diffusion. Les méthodes de maillage mobile adaptatif donnent un bon exemple de la façon de procéder pour l'équation de Burger en 1D en utilisant la différence finie. Quelqu'un pourrait-il offrir un exemple concret de résolution de l'équation d'advection-diffusion 1D en …

2
Confusion au sujet de la règle d'Armijo
J'ai cette confusion sur la règle Armijo utilisée dans la recherche en ligne. Je relisais la recherche de ligne de suivi, mais je n'ai pas compris de quoi s'agissait cette règle Armijo. Quelqu'un peut-il expliquer ce qu'est la règle d'Armijo? Le wikipedia ne semble pas bien expliquer. Merci

3
Confusion au sujet du problème de détection compressé
J'ai lu quelques références dont celle-ci . Je suis un peu confus quel problème d'optimisation compressé détecte les builds et essaie de résoudre. Est-ce minimizesubject to∥x∥1Ax=bminimize‖x‖1subject toAx=b\begin{array}{ll} \text{minimize} & \|x\|_1\\ \text{subject to} & Ax=b\end{array} ou et minimizesubject to∥x∥0Ax=bminimize‖x‖0subject toAx=b\begin{array}{ll} \text{minimize} & \|x\|_0\\ \text{subject to} & Ax=b\end{array} ou / et autre …



3
Comprendre le «taux de convergence» des méthodes itératives
Selon Wikipedia, le taux de convergence est exprimé comme un rapport spécifique de normes vectorielles. J'essaie de comprendre la différence entre les taux "linéaires" et "quadratiques", à différents moments (en gros, "au début" de l'itération et "à la fin"). Pourrait-on dire que: ek + 1ek+1e_{k+1}Xk + 1Xk+1x_{k+1}∥ ek∥‖ek‖\|e_k\| avec convergence …


2
Alternatives à l'analyse de stabilité de von neumann pour les méthodes de différences finies
Je travaille sur la résolution des équations de poroélasticité unidimensionnelles couplées (modèle de Biot), étant donné que: ∂- ( λ + 2 μ ) ∂2u∂X2+ ∂p∂X= 0-(λ+2μ)∂2u∂X2+∂p∂X=0-(\lambda+ 2\mu) \frac{\partial^2 u}{\partial x^2} + \frac{\partial p}{\partial x} = 0 ∂∂t[ γp + ∂u∂X] - κη[ ∂2p∂X2] =q( x , t )∂∂t[γp+∂u∂X]-κη[∂2p∂X2]=q(X,t)\frac{\partial}{\partial t} …

En utilisant notre site, vous reconnaissez avoir lu et compris notre politique liée aux cookies et notre politique de confidentialité.
Licensed under cc by-sa 3.0 with attribution required.