Considérons avec presque singulier, ce qui signifie qu'il existe une valeur propre de qui est très petite. Le critère d'arrêt habituel d'une méthode itérative est basé sur le résiduel et concerne les itérations peuvent s'arrêter lorsque avec n le numéro d'itération. Mais dans le cas que nous...
Je suis un ardent défenseur du développement piloté par les tests en informatique scientifique. Son utilité dans la pratique est tout simplement stupéfiante et atténue vraiment les problèmes classiques que les développeurs de code connaissent. Cependant, il y a des difficultés inhérentes à tester...
À des fins d'enseignement, j'aurais besoin d'une fonction continue d'une seule variable qui est "difficile" à approximer avec des polynômes, c'est-à-dire qu'il faudrait des puissances très élevées dans une série de puissances pour "bien" s'adapter à cette fonction. J'ai l'intention de montrer à mes...
La chose que je n'aime pas le plus à propos de MPI concerne les types de données (c'est-à-dire les cartes / masques de données) car ils ne correspondent pas très bien au C ++ orienté objet. boost::mpine prend en charge que MPI 1.1, cependant, à partir de leur site Web: boost :: mpi est une...
Je fais une tâche de classification de texte avec R, et j'obtiens une matrice de termes de document avec la taille 22490 par 120 000 (seulement 4 millions d'entrées non nulles, moins de 1% d'entrées). Maintenant, je veux réduire la dimensionnalité en utilisant PCA (analyse en composantes...
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...
Je suis curieux de savoir quels sont les bons algorithmes numériques pour évaluer la fonction hypergéométrique généralisée (ou série), définie comme pFq( un1, … , Unp; b1, … , Bq; z) = ∑k = 0∞( un1)k⋯ ( ap)k( b1)k⋯ ( bq)kzkk
Lors du calcul de la factorisation QR dans la pratique, on utilise les réflexions de Householder pour mettre à zéro la partie inférieure d'une matrice. Je sais que pour calculer les valeurs propres des matrices symétriques, le mieux que vous puissiez faire avec les réflexions de Householder est de...
Étant donné une matrice dense quelle est la meilleure façon de trouver sa base d'espace nul dans une certaine tolérance ?A∈Rm×n,m>>n;max(m)≈100000A∈Rm×n,m>>n;max(m)≈100000A \in R^{m \times n}, m >> n; max(m) \approx 100000 ϵϵ\epsilon Sur cette base, puis-je dire que certains cols...
Supposer minAvec(U)subject to Ui,j≤max{Ui,k,Uk,j},i,j,k=1,…,nminAvec(U)subject to Ui,j≤max{Ui,k,Uk,j},i,j,k=1,…,n\begin{align*} \min A &\mathrm{vec}(U) \\ &\text{subject to } U_{i,j} \leq \max\{U_{i,k}, U_{k,j}\}, \quad i,j,k = 1, \ldots, n \end{align*} où UUU est une matrice n \ fois n symétrique...
Puisque cette question est liée au calcul, j'ai décidé de poster ici. Espérons que cela sera considéré comme approprié. Je viens de commencer à exécuter des modèles atmosphériques et océaniques, et je me rends compte que j'ai besoin de plus de cœurs, de mémoire et d'espace disque que mon bureau...
J'utilise Valgrind + Callgrind pour profiler un solveur que j'ai écrit. Comme l'indique le manuel de l'utilisateur de Valgrind, j'ai compilé mon code avec les options de débogage du compilateur: "Sans informations de débogage, les meilleurs outils Valgrind seront capables de deviner à quelle...
J'aimerais être indépendant des logiciels commerciaux pour mon travail scientifique. Je trouve une dépendance à des packages commerciaux tels que Matlab et ses boîtes à outils insatisfaisante, car je ne sais pas si j'aurai accès à Matlab à l'avenir, et parce que je n'aime pas la langue. Par...
J'ai appris la méthode des éléments finis (également un peu sur d'autres méthodes numériques) mais je ne sais pas quelle est exactement la définition de ces deux erreurs et les différences entre
J'ai lu sur la science informatique sur Wikipédia, mais ma compréhension n'est pas très claire. La science informatique implique-t-elle une programmation? Dans quelle mesure la science informatique est-elle différente de l'informatique _ ____ , où le blanc peut être n'importe quelle discipline...
Je suis conscient du fait qu'inverser une matrice pour résoudre un système linéaire n'est pas une bonne idée, car il n'est pas aussi précis et aussi efficace que de résoudre directement le système ou d'utiliser la décomposition LU, Cholesky ou QR. Cependant, je n'ai pas pu vérifier cela avec un...
Supposons que j'ai un maillage 2D constitué de triangles non chevauchants , et un ensemble de points { p i } M i = 1 ⊂ ∪ N k = 1 T K . Quelle est la meilleure façon de déterminer dans quel triangle chacun des points se
Y a-t-il des raisons pour lesquelles on devrait choisir Runge Kutta implicite d'ordre élevé (IMRK) plutôt que le pas de temps BDF? BDF me semble beaucoup plus facile car étape IMRK a besoin de q résolutions linéaires par pas de temps. La stabilité pour BDF et IMRK semble être un point discutable....
Quelles approches sont utilisées dans la pratique pour estimer le nombre de conditions de grandes matrices
Il existe un système de contraintes linéaires Ax≤bAx≤b{\bf Ax} \leq {\bf b} . Je souhaite trouver un vecteur strictement positif x>0x>0{\bf x} > 0 qui satisfait ces contraintes. Cela signifie que est requis pour chaque composant de . Comment puis-je utiliser un solveur LP pour trouver un...