Existe-t-il un algorithme qui énumère les graphiques qui correspondent à une tessellation Delaunay de points en 3D? Dans l'affirmative, existe-t-il une paramétrisation efficace des géométries correspondant à tout "graphe de Delaunay"? Je cherche à énumérer systématiquement toutes les géométries stables de molécules d'une composition spécifiée sans aucune connaissance a priori …
Pour utiliser la transformée de Fourier rapide (FFT) sur des données échantillonnées uniformément, par exemple en relation avec des solveurs PDE, il est bien connu que la FFT est un algorithme ). Dans quelle mesure l'échelle FFT est-elle traitée en parallèle pour n → ∞ (c'est-à-dire très grande)?O (nlog( n …
De nombreux problèmes importants peuvent être exprimés sous la forme d'un programme linéaire entier mixte . Malheureusement, le calcul de la solution optimale à cette classe de problèmes est NP-Complete. Heureusement, il existe des algorithmes d'approximation qui peuvent parfois fournir des solutions de qualité avec seulement des quantités modérées de …
Je travaille avec la bibliothèque OpenFOAM C ++ Computational Continuum Mechanics (elle peut gérer les interactions fluide-solide, les flux MHD ...) qui utilise des maillages arbitraires non structurés. Cela a été motivé par l'idée d'utiliser l'avantage de la génération rapide (automatique généralement) de maillages non structurés pour simuler des problèmes …
Dans les deux méthodes de décomposition de domaine (DD) et multigrille (MG), on peut composer l'application des mises à jour de bloc ou des corrections grossières comme additive ou multiplicative . Pour les solveurs ponctuels, c'est la différence entre les itérations de Jacobi et de Gauss-Seidel. Le lisseur multiplicatif pour …
Dans une classe particulière de détecteurs, nos données se présentent sous la forme de paires de points en deux dimensions, et nous voulons enchaîner ces points en lignes. Les données sont bruyantes et regroupées dans un sens mais pas dans l'autre. Nous ne pouvons pas garantir un coup dans chaque …
La différenciation automatique nous permet d'évaluer numériquement la dérivée d'un programme sur une entrée particulière. Il existe un théorème selon lequel ce calcul peut être effectué à un coût inférieur à cinq fois le coût d'exécution du programme d'origine. Ce facteur de cinq est une borne supérieure. Dans quelles situations …
Je résous un problème physique en utilisant un schéma numérique implicite. Cela m'amène à résoudre une équation linéaire avec une matrice tridiagonale. J'ai codé cet algorithme à partir de Wikipedia. Je me demande s'il existe une bibliothèque efficace qui permet de résoudre ce type d'équation de manière optimisée. Une note …
J'ai utilisé GSL comme base de plusieurs de mes simulations, mais c'est un peu exagéré pour mes besoins et il définit son propre type complexe pour des raisons héritées. Plutôt que de coder mon propre solveur Runge-Kutta ODE, qui ne serait probablement pas très efficace, existe-t-il des solveurs ODE open …
Les flux à nombre élevé de Reynolds produisent des couches limites très minces. Si la résolution du mur est utilisée dans la simulation de grands tourbillons, le rapport d'aspect peut être de l'ordre de . De nombreuses méthodes deviennent instables dans ce régime car la constante inf-sup se dégrade en …
Je souhaite maximiser globalement une fonction de nombreux ( ) paramètres réels (résultat d'une simulation complexe). Cependant, la fonction en question est relativement coûteuse à évaluer, nécessitant environ 2 jours pour chaque ensemble de paramètres. Je compare différentes options et je me demandais si quelqu'un avait des suggestions.≈ 30≈30\approx 30 …
Je voudrais connaître une bonne méthode pour interpoler des données entre deux grilles non structurées, où une grille est une version plus grossière de l'autre. L'efficacité est très importante pour moi car je résous un problème PDE transitoire où je dois transférer des données entre les grilles à chaque étape …
L'essence de ma question est la suivante: j'ai un système de deux ODE. L'un a une contrainte de valeur initiale et l'autre une contrainte de valeur finale. Cela peut être considéré comme un système unique avec une contrainte de valeur initiale sur certaines variables et une contrainte de valeur finale …
Étant donné deux matrices et B , je voudrais trouver les vecteurs x et y , tels que, min ∑ i j ( A i j - x i y j B i j ) 2 . Sous forme de matrice, j'essaie de minimiser la norme de Frobenius de A …
Je tente actuellement de calculer à peu de frais une bonne estimation de rang pour une matrice . Par conséquent, je calcule une décompostion QR pivotante en utilisantUNEAA [Q,R,E]=qr(A) dans Matlab. J'évalue le rang de utilisantUNEAA tol = size(A,n)*eps*norm(A,'fro'); r = sum(abs(diag(R))>tol) Cela fonctionne bien et un tracé sur toutes …
We use cookies and other tracking technologies to improve your browsing experience on our website,
to show you personalized content and targeted ads, to analyze our website traffic,
and to understand where our visitors are coming from.
By continuing, you consent to our use of cookies and other tracking technologies and
affirm you're at least 16 years old or have consent from a parent or guardian.