Science computationnelle

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

3
Alternatives à Comsol Multiphysics
Cela pourrait être une question mieux adaptée au côté des recommandations de logiciels de SE, mais je pense que les personnes qui fréquentent cette partie de SE sont plus susceptibles de pouvoir répondre à cette question. Je recherche une alternative gratuite (pas seulement en liberté) à Comsol Multiphysics . Voici …
2
En ce qui concerne la différenciation automatique, la transformation du code source (STC) est-elle plus efficace que la surcharge d'opérateur (OO)?
Nous travaillons sur un modèle bayésien pour un processus spatio-temporel, et utilisons un échantillonneur No-U-Turn (NUTS) qui nécessite un modèle pour la probabilité logarithmique et son gradient par rapport aux paramètres du modèle. Plus succinctement, nous avons une fonction de probabilité de log assez compliquée , impliquant des distributions statistiques, …
1
La pression comme multiplicateur de Lagrange
Dans les équations de Navier-Stokes incompressibles, ρ(ut+(u⋅∇)u)∇⋅u=−∇p+μΔu+f=0ρ(ut+(u⋅∇)u)=−∇p+μΔu+f∇⋅u=0\begin{align*} \rho\left(\mathbf{u}_t + (\mathbf{u} \cdot \nabla)\mathbf{u}\right) &= - \nabla p + \mu\Delta\mathbf{u} + \mathbf{f}\\ \nabla\cdot\mathbf{u} &= 0 \end{align*} le terme de pression est souvent appelé multiplicateur de Lagrange imposant la condition d'incompressibilité. Dans quel sens est-ce vrai? Existe-t-il une formulation des équations de Navier-Stokes …
2
Quel est le moyen le plus rapide pour calculer toutes les valeurs propres d'une matrice d'adjacence très grande et clairsemée en python?
J'essaie de comprendre s'il existe un moyen plus rapide de calculer toutes les valeurs propres et vecteurs propres d'une matrice de contiguïté très grande et clairsemée que d'utiliser scipy.sparse.linalg.eigsh Pour autant que je sache, cette méthode utilise uniquement la rareté et attributs de symétrie de la matrice. Une matrice d'adjacence …
3
Choix de la taille des pas à l'aide des ODE dans matlab
Salut et merci d'avoir donné le temps de regarder ma question. Ceci est une version mise à jour de ma question que j'ai publiée plus tôt dans physics.stackexchange.com J'étudie actuellement un condensateur à excitons 2D Bose-Einstein Condensate et je suis curieux de connaître l'état fondamental de ce système. La méthode …
2
Quels préconditionneurs (et solveurs) en PETSc pour les systèmes symétriques indéfinis dois-je utiliser?
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 Lagrange apparaissent en …
1
Utiliser l'itération à virgule fixe pour découpler un système de pde
Supposons que j'ai eu un problème de valeur limite: dud2udx2+dvdx=f in Ωd2udx2+dvdx=f in Ω\frac{d^2u}{dx^2} + \frac{dv}{dx}=f \text{ in } \Omega u=h en ∂Ωdudx+d2vdx2=g in Ωdudx+d2vdx2=g in Ω\frac{du}{dx} +\frac{d^2v}{dx^2} =g \text{ in } \Omega u=h in ∂Ωu=h in ∂Ωu=h \text{ in } \partial\Omega Mon objectif est de décomposer la solution de …
2
résoudre
Je porte un code existant de MATLAB vers C ++ et j'ai un système linéaire pour résoudre (plutôt que la forme la plus typique A x = b )x A = bxA=bxA=bA x = bAx=bAx=b La matrice est dense et de forme générale, mais ne dépasse pas 1000x1000. Donc, dans …
1
Sélection de la méthode pour la quadrature numérique
Il existe plusieurs familles de méthodes pour la quadrature numérique. Si j'ai une classe spécifique d'intégrandes, comment choisir la méthode idéale? Quelles sont les questions pertinentes à poser à la fois sur l'intégrande (par exemple, est-elle fluide? A-t-elle des singularités?) Et sur le problème de calcul (par exemple, tolérance aux …
12 quadrature 
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.