J'ai commencé à utiliser python comme langage de programmation pour faire toutes mes tâches en CFD. J'ai très peu d'expérience en programmation. Je suis issu du génie mécanique et je poursuis des études supérieures en génie aérospatial. Parfois, l'aspect informatique du CFD devient plus fastidieux que la manipulation des équations …
Je suis très impressionné par les performances en série des préconditionneurs ILU multiniveaux , en particulier pour les Helmholtz hétérogènes , mais je suis surpris de ne pouvoir trouver aucune implémentation open source. En particulier, ILUPACK met les binaires à la disposition des universitaires, mais il ne semble pas qu'ils …
Je suis en train de résoudre un système de deux PDE couplés dans deux dimensions spatiales et en temps de calcul. Étant donné que les évaluations de fonctions sont coûteuses, je voudrais utiliser une méthode en plusieurs étapes (initialisée à l'aide de Runge-Kutta 4-5). La méthode Adams-Bashforth utilisant cinq évaluations …
J'ai un ensemble de données fonctionnant sur des millions de points de données en 3D. Pour le calcul que je fais, je dois calculer le voisin (recherche de plage) de chaque point de données dans un rayon, essayer d'ajuster une fonction, calculer l'erreur pour l'ajustement, répéter cela pour le prochain …
Lors de l'application de la formule classique de l'angle entre deux vecteurs: α=arccosv1⋅v2∥v1∥∥v2∥α=arccosv1⋅v2‖v1‖‖v2‖\alpha = \arccos \frac{\mathbf{v_1} \cdot \mathbf{v_2}}{\|\mathbf{v_1}\| \|\mathbf{v_2}\|} on constate que, pour des angles très petits / aigus, il y a une perte de précision et le résultat n'est pas précis. Comme expliqué dans cette réponse Stack Overflow , …
De nombreux algorithmes utilisés dans le calcul scientifique ont une structure inhérente différente de celle des algorithmes couramment considérés dans les formes de génie logiciel moins gourmandes en mathématiques. En particulier, les algorithmes mathématiques individuels ont tendance à être très complexes, impliquant souvent des centaines ou des milliers de lignes …
Je recherche une bibliothèque qui effectue des opérations matricielles sur de grandes matrices clairsemées sans sacrifier la stabilité numérique. Les matrices seront 1000+ par 1000+ et les valeurs de la matrice seront comprises entre 0 et 1000. Je vais exécuter l' algorithme de calcul d'index donc je vais générer des …
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 …
J'essaie d'inverser numériquement la transformation intégrale suivante: F(y)=∫∞0yexp[−12(y2+x2)]I0(xy)f(x)dxF(y)=∫0∞yexp[−12(y2+x2)]I0(xy)f(x)dxF(y) = \int_{0}^{\infty} y\exp{\left[-\frac{1}{2}(y^2 + x^2)\right]} I_0\left(xy\right)f(x)\;\mathrm{d}x Donc, pour un j'ai besoin d'approximer où:F(y)F(y)F(y)f(x)f(x)f(x) f(x)f(x)f(x) et sont réels et positifsF(y)F(y)F(y) (ce sont des distributions de probabilité continues) x,yx,yx,y sont réels et positifs (ce sont des grandeurs) J'ai une méthode très compliquée et brutale …
Ainsi , le théorème de décomposition de Cholesky que que toute réel symétrique définie positive matrice a une décomposition de Cholesky où est une matrice triangulaire inférieure.MMMM=LL⊤M=LL⊤M= LL^\topLLL Compte tenu , nous savons déjà , il y a des algorithmes rapides pour calculer le facteur de Cholesky .MMMLLL Supposons maintenant …
Bonjour communauté scicomp, J'ai travaillé dans le domaine des algorithmes de graphes en utilisant des frameworks tels que NetworkX (Python), JUNG et YFiles (Java). J'entre maintenant dans le domaine de l'informatique parallèle et haute performance. Pour un nouveau projet, je recherche une bibliothèque de graphes C ++ avec les fonctionnalités …
Il s'agit d'un simple problème d'énumération de symétrie. Je donne ici le contexte complet, mais aucune connaissance de la chimie quantique n'est nécessaire. L'intégrale de deux particules est: ⟨ i j | k l ⟩ = ∫ ψ * i ( x ) ψ * j ( x ' ) …
Je travaille sur un problème qui peut être mis en parallèle en utilisant une seule opération mpi_allgather ou une opération mpi_scatter et une opération mpi_gather. Ces opérations sont appelées dans une boucle while, elles peuvent donc être appelées plusieurs fois. Dans l'implémentation avec un schéma MPI_allgather, je rassemble un vecteur …
Je fais une diagonalisation Lanczos d'une grande matrice clairsemée (~ 2 millions d'éléments). Presque toutes les étapes de l'algorithme de Lanzcos sont effectuées en parallèle sur le GPU, à l'exception de la diagonalisation de la matrice de Lanczos pour vérifier la convergence. Pour cela, j'ai utilisé l'algorithme TQLI de Recettes …
Question Comment classeriez-vous un nuage de points par rapport à un maillage non structuré de cellules hexaédriques? Chaque cellule a un centre et une étiquette unique pour la représenter. Il y a essentiellement deux points de nuage (nuage de points d'origine et un nuage de points des centres cellulaires), mais …
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.