Science computationnelle

1

Quels sont les schémas numériques possibles pour une équation de diffusion avec un terme de réaction non linéaire?

Pour un domaine convexe simple en 2D, nous avons certains u ( x ) satisfaisant l'équation suivante: - d i v ( A ∇ u ) + c u n = f avec certaines conditions aux limites de Dirichlet et / ou Neumann. À ma connaissance, l'application de la méthode …

11 pde finite-element

6

Quels outils open source sont disponibles pour visualiser les vibrations moléculaires?

Je voudrais visualiser une vibration moléculaire qui n'est pas un mode normal. Je voudrais présenter une représentation vectorielle statique du mouvement et je voudrais une certaine flexibilité dans le style vectoriel (taille, couleur, etc.). Je suis également intéressé à produire une vidéo de la vibration. Quelles sont les bonnes ressources …

11 computational-chemistry

6

Différences finies sur les domaines aux frontières irrégulières

Quelqu'un peut-il m'aider à trouver les livres sur les solutions numériques (différences finies et méthodes de Crank-Nicolson) des équations de Poisson et de diffusion, y compris des exemples sur la géométrie irrégulière, comme un domaine composé de l'aire entre un rectangle et un cercle (en particulier des livres ou des …

11 pde finite-difference

2

Comment travailler avec des pointeurs de fonction dans Fortran dans les programmes scientifiques

Voici une utilisation typique des pointeurs de fonction en C. Je voudrais faire quelque chose de similaire dans Fortran. J'ai quelques idées, mais j'aimerais savoir s'il existe un moyen canonique de le faire. Les pointeurs de fonction et les contextes transmis par l'utilisateur sont stockés, puis appelés ultérieurement. typedef PetscErrorCode …

11 petsc fortran

5

Est-il préférable de se concentrer sur l'étude des mathématiques ou du calcul?

Parallèlement à mes recherches sur les méthodes de sous-espace de Krylov, j'ai la possibilité d'explorer les mathématiques derrière HPC un pas en avant ou la théorie du calcul (matériel, OS, compilateurs etc.). Actuellement, je connais suffisamment les deux pour me débrouiller. Par exemple, je sais comment dériver les équations de …

11 parallel-computing education

1

Comment formuler une matrice de masse localisée dans FEM

Lors de la résolution d'EDP dépendant du temps en utilisant la méthode des éléments finis, par exemple, l'équation de la chaleur, si nous utilisons un pas de temps explicite, nous devons résoudre un système linéaire en raison de la matrice de masse. Par exemple, si nous nous en tenons à …

11 finite-element time-integration navier-stokes

2

Limitations de Pandas et ses alternatives en Python

J'ai lu quelque part que Pandas a été développé pour la première fois pour le monde financier, du moins pas spécialement pour les sciences naturelles (physique, biologie, etc.). Je viens de commencer à utiliser Pandas et j'ai déjà rencontré deux problèmes que je ne pouvais pas résoudre sans utiliser d'autres …

11 python data-analysis data-management

2

calculer la SVD tronquée, une valeur / vecteur singulier à la fois

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 …

11 linear-algebra algorithms svd

3

Volume de coque convexe 3D de petits ensembles de points tout sur la coque

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 …

11 algorithms reference-request computational-geometry

2

Descente en gradient et descente en gradient conjugué

Pour un projet, je dois implémenter ces deux méthodes et comparer leurs performances sur différentes fonctions. Il semble que la méthode du gradient conjugué soit destinée à résoudre des systèmes d'équations linéaires du for A x = bUNEX=b A\mathbf{x} = \mathbf{b} Où est une matrice n par n symétrique, définie …

11 optimization conjugate-gradient

1

Mise en œuvre optimale de la déformation de transport dans Matlab

J'implémente le document " Transport de masse optimal pour l'enregistrement et la déformation ", mon objectif étant de le mettre en ligne car je ne trouve aucun code de transport de masse eulérien en ligne et cela serait intéressant au moins pour la communauté des chercheurs en traitement d'images. Le …

11 pde finite-difference matlab fluid-dynamics image-processing

1

Problèmes de référence pour les algorithmes de réorganisation des valeurs propres recherchés

Chaque matrice réelle peut être réduite à la forme réelle de Schur T = U T A U en utilisant une transformation orthogonale similiary U . Ici, la matrice T est de forme quasi triangulaire avec 1 par 1 ou 2 par 2 blocs sur la diagonale principale. Chaque 1 …

10 eigenvalues eigensystem lapack numerics scalapack

1

Existe-t-il des raccourcis pour les systèmes d'approximation numérique d'équations différentielles ordinaires lorsqu'ils sont autonomes?

Les algorithmes existants pour résoudre les ODE gèrent les fonctions , oùy∈Rn. Mais dans de nombreux systèmes physiques, l'équation différentielle est autonome, doncdyréyrét= f( y, t )dydt=f(y,t)\frac{dy}{dt} = f(y, t)y∈ Rny∈Rny \in \mathbb R^n,y∈Rn, sanst. Avec cette hypothèse simplificatrice, quelles améliorations peuvent être observées dans les méthodes numériques existantes? Par …

10 ode numerics

1

Existe-t-il des cadres de test pour le développement de logiciels numériques

J'ai constaté que beaucoup de mes programmes de science informatique ont des exigences de test qui ne sont pas couvertes par les cadres de test standard: Test de temps de calcul Pour vous assurer que les algorithmes ne ralentissent pas. Je pourrais faire quelque chose comme assureSmallerEqual(RuntimeWrapper(algorithm),53)mais j'aimerais que le …

10 testing

1

Première apparition de l'expression «crime inverse»

Dans la recherche sur les problèmes inverses, il est courant de construire un ensemble de données synthétiques à partir d'un ensemble connu de paramètres, puis de tester si la technique d'inversion peut reconstruire ces paramètres. Pour ce faire, il est important d'ajouter des niveaux appropriés de bruit aléatoire aux données …

10 reference-request inverse-problem