Étant donné un état souhaité et un paramètre de régularisation , considérons le problème de trouver un état et un contrôle pour minimiser une soumis à la contrainte \ begin {equation} Ay = u. \ end {equation} où pour simplifier on peut penser à y, y_0, u \ in \ mathbb R ^ n et A \ in \ mathbb R ^...
n × n G n n GUNEUNEA est une matrice clairsemée symétrique positive définie (SPD). est une matrice diagonale clairsemée. est grand ( > 10000) et le nombre de nonzeros dans le est généralement de 100 à 1000.n × nn×nn \times nggGnnnnnnggG UNEUNEA a été factorisée sous la forme de Cholesky comme .L...
La recherche reproductible en calcul vise à mettre le code nécessaire pour générer les résultats dans un document de calcul à la disposition d'autres chercheurs afin qu'ils puissent exécuter ce code pour reproduire les résultats dans ce document. Je voudrais que toutes mes recherches soient...
Je cherche à utiliser OpenFOAM pour résoudre les flux internes de base dans CFD. Quelle est la meilleure façon de commencer, et est-ce que quelqu'un pourrait m'indiquer une bonne référence en ligne pour répondre à mes questions une fois que j'aurai plongé? J'ai entendu dire que c'est un créneau...
Je voudrais savoir s'il existe un moyen rapide de calculer la distance euclidienne de deux vecteurs en octave. Il semble qu'il n'y ait pas de fonction spéciale pour cela, alors devrais-je simplement utiliser la formule avec
Supposons que j'ai une équation matricielle définie récursivement comme A[n] = inverse([1 - b[n]A[n+1]]) * a[n] L'équation pour A [1] ressemble alors à une fraction continue, pour laquelle il existe des méthodes très efficaces qui évitent un recalcul fastidieux (voir «Recettes numériques» pour...
En lisant de la littérature sur les solveurs PDE, je suis tombé sur le terme pseudo pas de temps aujourd'hui. Il semble que ce soit un terme courant, mais je n'ai pas réussi à trouver une bonne définition ou un article d'introduction. Par conséquent: Qu'est-ce qu'un pseudo pas de temps et comment...
Les quantités physiques comme la pression, la densité, l'énergie, la température et la concentration doivent toujours être positives, mais les méthodes numériques calculent parfois des valeurs négatives au cours du processus de solution. Ce n'est pas correct car les équations calculent des valeurs...
Vous avez enfin réussi à découvrir comment les atomes sont disposés spatialement sur votre entité moléculaire nouvellement découverte. Par exemple, grâce à des moyens spectroscopiques, vous êtes maintenant en possession d'un tas de coordonnées atomiques, de types d'atomes, de longueurs de liaisons,...
Je viens d'installer une carte graphique Nvidia GT660 sur mon bureau et, après quelques difficultés, j'ai réussi à l'interfacer avec R. J'ai joué avec plusieurs packages R qui utilisent des GPU, en particulier gputools, et je comparais le temps pris par mon GPU et mon processeur pour effectuer...
Dans la plupart (toutes?) Des implémentations de la méthode multipolaire rapide (FMM), les octrees sont utilisés pour décomposer le domaine concerné. Théoriquement, les octrees fournissent une borne volumétrique simple, qui est utile pour prouver le temps d'exécution O (n) d'un FMM. Au-delà de...
Quels sont les meilleurs algorithmes (évolutivité et efficacité) pour générer des maillages quadruples non structurés en 2D? Où puis-je trouver un bon générateur de maillage quadruple non structuré? (open-source
Je voudrais utiliser f2pyavec le Fortran moderne. En particulier, j'essaie de faire fonctionner l'exemple de base suivant. Ceci est le plus petit exemple utile que j'ai pu générer. ! alloc_test.f90 subroutine f(x, z) implicit none ! Argument Declarations ! real*8, intent(in) :: x(:) real*8,...
Je souhaite développer un logiciel de calcul scientifique parallèle à partir de zéro. Je veux quelques réflexions sur quelle langue commencer. Le programme consiste à lire / écrire des données dans des fichiers txt et à effectuer des calculs lourds en parallèle, avec de nombreuses factorisations LU...
J'essaie d'implémenter la fonction suivante en virgule flottante double précision avec une faible erreur relative : logsum(x,y)=log(exp(x)+exp(y))logsum(x,y)=log(exp(x)+exp(y))\mathrm{logsum}(x,y) = \log(\exp(x) + \exp(y)) Ceci est largement utilisé dans les applications statistiques pour...
Voulez-vous améliorer ce post? Fournissez des réponses détaillées à cette question, y compris des citations et une explication de la raison pour laquelle votre réponse est correcte. Les réponses sans suffisamment de détails peuvent être modifiées ou supprimées. Le raffinement adaptatif du maillage...
J'essaie de résoudre certains problèmes d'optimisation non linéaire sans contrainte sur le GPU (CUDA). La fonction objectif est une fonction non linéaire lisse, et son gradient est relativement bon marché à calculer analytiquement, donc je n'ai pas besoin de me soucier de l'approximation numérique....
Je voudrais savoir comment les méthodes d'ondelettes peuvent être appliquées à la PDE, mais malheureusement je ne connais pas de bonne ressource pour en savoir plus sur ce sujet. Il semble que de nombreuses introductions aux ondelettes se concentrent sur la théorie de l'interpolation, par exemple,...
Considérons un IVP pour le système ODE , y ( x 0 ) = y 0 . Le plus souvent, ce problème est considéré comme rigide lorsque la matrice de Jacobi ∂ fy′= f( x , y)y′=F(X,y)y'=f(x,y)y( x0) = y0y(X0)=y0y(x_0)=y_0a à lafois desvaleurs propres avec une très grande partie réelle négative et des valeurs...
Dans la littérature sur les méthodes numériques, de nombreux articles de recherche consistent en une description d'une nouvelle variation algorithmique, suivie de quelques problèmes de test comparant la nouvelle méthode avec une ou deux méthodes existantes. Il est donc difficile de déterminer...