J'ai un système de équations non linéaires que je veux résoudre numériquement:nnn f = ( f 1 , … , f n )F( x ) = af(x)=a\mathbf{f}(\mathbf{x})=\mathbf{a} F= ( f1, … , Fn)x =( x1, … , Xn)f=(f1,…,fn)x=(x1,…,xn)\mathbf{f}=(f_1,\dots,f_n)\quad\mathbf{x}=(x_1,\dots,x_n) Ce système présente un certain...
J'ai utilisé le SVD d'Intel MKL ( dgesvdvia SciPy) et j'ai remarqué que les résultats sont considérablement différents lorsque je change de précision entre float32et float64lorsque ma matrice est mal conditionnée / n'est pas pleine. Existe-t-il un guide sur le montant minimum de régularisation à...
Il y a quelque temps, j'avais écrit un code qui tentait de calculer sans utiliser les fonctions de bibliothèque. Hier, je passais en revue l'ancien code et j'ai essayé de le rendre aussi rapide que possible (et correct). Voici ma tentative jusqu'à présent:l o g( x )log(X)log(x) const double ee =...
Je ne peux pas imaginer que je suis le premier à penser au problème suivant, donc je serai satisfait d'une référence (mais une réponse complète et détaillée est toujours appréciée): Disons que vous avez une symétrique définie positive . n est considéré comme très grand, donc garder Σ en mémoire est...
Je lis un article [1] où ils résolvent l'équation non linéaire suivante utilisant des méthodes de différences finies. Ils analysent également la stabilité des schémas à l'aide de l'analyse de stabilité de Von Neumann. Cependant, comme le réalisent les auteurs, cela ne s'applique qu'aux PDE...
Je m'intéresse au calcul de la solution d'un système de lage d'OD en utilisant une méthode krylov comme dans [1]. Une telle méthode implique des fonctions liées à l'exponentielle (les fonctions dites ). Elle consiste essentiellement à calculer l'action de la fonction matricielle en construisant un...
J'effectue une recherche de ligne dans le cadre d'un algorithme BFGS quasi-Newton. Dans une étape de la recherche de ligne, j'utilise une interpolation cubique pour me rapprocher du minimiseur local. Soit la fonction d'intérêt. Je veux trouver un tel que .x ∗ f ′ ( x ∗ ) ≈ 0f:R→R,f∈C1f:R→R,f∈C1f :...
Quelles sont les principales différences entre FEM et XFEM? Quand devrions-nous (ne pas) utiliser XFEM à la place de FEM et vice versa? En d'autres termes, lorsque je rencontre un nouveau problème, comment puis-je savoir lequel