Je voudrais modéliser l'écoulement laminaire de l'eau des racines à la tige d'une plante. À la toute fin des racines, les tubes varient d'une échelle millimétrique à centimétrique en diamètre et en longueur. À mesure que nous nous rapprochons de la tige, les racines grossissent en longueur et en...
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...
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...
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 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...
Existe-t-il un algorithme SVD tronqué qui calcule les valeurs singulières une par une? Mon problème: je voudrais calculer les premières valeurs singulières (et vecteurs singuliers) d'une grande matrice dense M , mais je ne sais pas ce que serait une valeur appropriée de k . M est grand, donc pour...
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...
J'ai une question similaire à celle posée auparavant sauf en 3D, et je n'ai besoin que du volume, pas de la forme réelle de la coque. Plus précisément, on me donne un petit ensemble de points (disons, 10-15) en 3D, qui sont tous connus pour se situer sur la coque convexe de l'ensemble de points...
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...
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...
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...
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...
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 ' ) ψ k ( x ) ψ l ( x ' )⟨ij|kl⟩⟨ij|kl⟩\langle ij|kl\rangle Et il...
Étant donné une matrice symétrique définie positive, quel est l'algorithme le plus rapide pour calculer la matrice inverse et son déterminant? Pour les problèmes qui m'intéressent, la dimension de la matrice est de 30 ou moins. Une précision et une vitesse élevées sont vraiment nécessaires. (des...
J'ai une matrice carrée symétrique réelle dense. La dimension est d'environ 1000x1000. Je dois calculer le premier composant principal et me demander quel pourrait être le meilleur algorithme pour ce faire. Il semble que MATLAB utilise les algorithmes Arnoldi / Lanczos (pour eigs). Mais en lisant à...
Fermé. Cette question est hors sujet . Il n'accepte pas actuellement les réponses. Voulez-vous améliorer cette question? Mettez à jour la question afin qu'elle soit sur le sujet pour Computational Science Stack Exchange. Fermé il y a 5 ans . Existe-t-il un degré de complexité supérieur à et...
Quelles situations d'analyse numérique deviennent plus / moins stables, ont une convergence plus rapide / plus lente, ou sont autrement assez différentes lorsqu'il s'agit de fonctions de variable complexe au lieu de fonctions d'une variable réelle?
Soit et deg f > deg g . Je recherche des algorithmes asymptotiquement rapides et numériquement stables pour calculer f mod g . Dans les applications envisagées, les deux f , g sont des polynômes denses à coefficients à virgule flottante double précision. Mais, pour l'instant, je m'intéresse plus...
Je voudrais savoir s'il existe un algorithme qui, étant donné un ensemble de points o et un angle, calcule la coque convexe si l'angle est et étant donné un α > 0, calcule une enveloppe qui suit de plus près le "périmètre".α = 0α=0\alpha = 0α > 0α>0\alpha > 0 Et s'il existe une définition...
Étant donné une matrice générique clairsemée avec m << n (correction: ) éléments non nuls (typiquement ). est générique dans le sens où il n'a pas de propriétés spécifiques (par exemple, un caractère définitif positif), et aucune structure (par exemple, des bandes) n'est supposée.A ∈ Rn ×...