Je veux résoudre un problème non linéaire avec des contraintes d'égalité non linéaires et j'utilise un lagrangien augmenté avec un terme de régularisation de pénalité qui, comme on le sait, gâche le numéro de condition de mes systèmes linéarisés (à chaque itération de Newton, je veux dire) . Plus le …
Nous comparons les performances de diverses méthodes numériques qui peuvent être utilisées pour résoudre l'équation de Schrodinger pour l'atome d'hydrogène interagissant avec une forte impulsion laser (trop forte pour utiliser des méthodes de perturbation). Lorsque vous utilisez des schémas de discrétisation pour la partie radiale, il semble que la plupart …
J'ai une intégrale incorrecte (bidimensionnelle) je= ∫UNEW( x , y)F( x , y)d x d yI=∫AW(x,y)F(x,y)dxdyI=\int_A \frac{W(x,y)}{F(x,y)}\,\mbox{d}x\mbox{d}y où le domaine d'intégration est plus petit que , mais encore restreint par . Puisque et sont lisses etx = [ - 1 , 1 ] y = [ - 1 , 1 …
J'ai besoin de résoudre s.t.minx∥Ax−b∥22,∑ixi=1,xi≥0,∀i.minx‖Ax−b‖22,s.t.∑ixi=1,xi≥0,∀i.\begin{alignat}{1} & \min_{x}\|Ax - b\|^2_{2}, \\ \mathrm{s.t.} & \quad\sum_{i}x_{i} = 1, \\ & \quad x_{i} \geq 0, \quad \forall{i}. \end{alignat} Je pense que c'est un problème quadratique qui devrait être résolu avec CVXOPT , mais je ne sais pas comment.
Je souhaite configurer des calculs pour vérifier si un critère de distance est satisfait: c'est-à-dire que la distance entre un vecteur et un vecteur d'anthère doit être inférieure à une valeur . Mes données sont partitionnées selon une grille orthogonale de coordonnées. Étant donné que ma coupure est plus petite …
Je me demandais si quelqu'un avait des suggestions de textes ou d'articles d'enquête sur les méthodes de décomposition (par exemple les décompositions primaires, doubles, Dantzig – Wolfe) pour résoudre de gros problèmes de programmation mathématique. J'ai aimé les "Notes sur les méthodes de décomposition" de Stephen Boyd , et ce …
L'article «Expression Templates Revisited: A Performance Analysis of Current Methodologies» du SIAM Journal of Scientific Computing fait référence à la bibliothèque d'algèbre linéaire «Blaze». Je n'en ai jamais entendu parler auparavant et je n'arrive pas à trouver de références en ligne. (Les recherches google évidentes rendent le papier ci-dessus.) Alors, …
J'ai demandé des éclaircissements sur une question récente à propos de minpack et j'ai obtenu le commentaire suivant: Tout système d'équations équivaut à un problème d'optimisation, c'est pourquoi les méthodes basées sur Newton en optimisation ressemblent beaucoup aux méthodes basées sur Newton pour résoudre des systèmes d'équations non linéaires. Ce …
Je suis nouveau dans la programmation et j'essaie de résoudre mon premier gros problème et d'écrire mon premier gros programme. J'ai cherché des exemples de code open source pour apprendre, mais jusqu'à présent, je n'ai trouvé que du code dans des langues que je ne comprends pas complètement ou qui …
Je travaille sur une bibliothèque de matrices uniquement en-tête pour fournir un certain degré raisonnable de capacité d'algèbre linéaire dans un package aussi simple que possible, et j'essaie de faire le point sur l'état actuel de la technique: le calcul de la SVD d'un matrice complexe. Je fais une décomposition …
Les protéines sont généralement représentées sous forme de dessin animé, avec des feuilles β sous forme de flèches et des hélices α sous forme de bobines: Je me demande, y a-t-il quelque part une référence qui décrit la construction de cette représentation? Autrement dit, quels objets mathématiques sont utilisés pour …
Supposons que j'ai deux matrices Nx2, Mx2 représentant respectivement N, M 2d vecteurs. Existe-t-il un moyen simple et efficace de calculer les distances entre chaque paire de vecteurs (n, m)? La manière simple mais inefficace est bien sûr: d = zeros(N, M); for i = 1:N, for j = 1:M, …
Habituellement, on consulte un document ou un livre pour trouver des points d'intégration et des poids pour le triangle unitaire et les tétraèdres. Je recherche une méthode pour calculer automatiquement ces points et ces poids. L' exemple de code Mathematica suivant calcule les poids et points d'intégration pour les éléments …
J'ai un programme Mathematica qui effectue des intégrales en 3 ou 4 dimensions en utilisant la QuasiMonteCarlométhode. Le problème est qu'il prend un temps énormément long pour s'exécuter, au point que certains de ces calculs ne peuvent pas se terminer dans le temps de travail maximum disponible sur notre cluster …
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.