Dans la simulation de semi-conducteur, il est courant que les équations soient mises à l'échelle afin qu'elles aient des valeurs normalisées. Par exemple, dans des cas extrêmes, la densité électronique dans les semi-conducteurs peut varier sur 18 ordres de grandeur, et le champ électrique peut changer de manière galbée, sur …
Je développe un code plus grand pour effectuer des calculs de valeurs propres d'énormes matrices clairsemées, dans le contexte de la physique computationnelle. Je teste mes routines contre le simple oscillateur harmonique dans une dimension, car les valeurs propres sont bien connues analytiquement. Ce faisant et en comparant mes propres …
Parfois, tout en optimisant le code, il est nécessaire de chronométrer certaines parties du code, j'utilise les éléments suivants depuis des années, mais je me demandais s'il y avait une manière plus simple / meilleure de le faire? call system_clock(count_rate=clock_rate) !Find the time rate call system_clock(count=clock_start) !Start Timer call do_something_subroutine …
J'ai entendu de façon anecdotique que lorsque l'on essaie de faire numériquement une intégrale de la forme ∫∞0F( x ) J0( x )d x∫0∞F(X)J0(X)réX\int_0^\infty f(x) J_0(x)\,\mathrm{d}x avec lisse et bien comporté (par exemple pas lui-même très oscillatoire, non singulier, etc.), alors il aidera la précision à le réécrire commeF( x …
J'ai un ensemble de données x1,x2,…,xkx1,x2,…,xkx_{1}, x_{2}, \ldots, x_{k} et je veux trouver le paramètre mmm tel qu'il minimise la somme ∑i=1k∣∣m−xi∣∣.∑i=1k|m−xi|.\sum_{i=1}^{k}\big|m-x_i\big|. C'est minm∑i=1k∣∣m−xi∣∣.minm∑i=1k|m−xi|.\min_{m}\sum_{i=1}^{k}\big|m-x_i\big|.
J'utilise HDF5 depuis des années, mais à mesure que la taille de l'ensemble de données augmente, je commence à rencontrer les mêmes problèmes répertoriés ici http://cyrille.rossant.net/moving-away-hdf5/ Pouvez-vous m'indiquer un format comme HDF5 avec - Un bon support pour les écritures parallèles - Un support pour un accès en bloc de …
Je voudrais consolider la résolution de systèmes d'ODE sur des GPU, dans un cadre «trivialement parallélisable». Par exemple, faire une analyse de sensibilité avec 512 jeux de paramètres différents. Idéalement, je veux faire une résolution ODE avec un solveur d'horodateur adaptatif intelligent comme CVODE, plutôt qu'un horodatage fixe comme Forward …
Je parallélise le code pour résoudre numériquement un modèle d'équilibre de population à 5 dimensions. Actuellement, j'ai un très bon code parallélisé MPICH2 dans FORTRAN mais à mesure que nous augmentons les valeurs des paramètres, les tableaux deviennent trop volumineux pour fonctionner en mode mémoire distribuée. J'ai accès à un …
Je fais quelques simulations de l'équation de Langevin, pour diverses forces externes. Être dit que C est rand()de stdlib.hpeut introduire un biais dans mes résultats, je suis sur un Mersenne Twister. Néanmoins, je voudrais savoir (et voir) exactement quel genre d'erreurs un générateur congruentiel linéaire peut introduire dans ma simulation. …
Je rédige actuellement ma thèse de doctorat. J'ai passé une fraction importante de mon doctorat à nettoyer et étendre le code scientifique existant, à appliquer les meilleures pratiques en génie logiciel qui n'étaient pas utilisées auparavant, et j'aimerais écrire à ce sujet dans ma thèse. Plutôt que de simplement dire …
Je suis en train de résoudre un système de deux PDE couplés dans deux dimensions spatiales et en temps de calcul. Étant donné que les évaluations de fonctions sont coûteuses, je voudrais utiliser une méthode en plusieurs étapes (initialisée à l'aide de Runge-Kutta 4-5). La méthode Adams-Bashforth utilisant cinq évaluations …
Je compare deux méthodes itératives pour inverser des matrices carrées aléatoires. Étant donné que les matrices sont aléatoires, chaque cas de test prend à la fois différentes quantités d'itérations et différents temps écoulés. Ma question est, en plus du temps CPU moyen, est la valeur moyenne des itérations prises par …
Quelqu'un ici a-t-il utilisé le calcul scientifique en double précision avec des GPU de nouvelle génération (par exemple K20) via Python? Je sais que cette technologie évolue rapidement, mais quelle est la meilleure façon de le faire actuellement? Le GPU est hors de portée des bibliothèques scientifiques Python numpy et …
Je viens de commencer à jouer avec FEniCS. Je résout Poisson avec des éléments de troisième ordre et je souhaite visualiser les résultats. Cependant, lorsque j'utilise plot (u), la visualisation n'est qu'une interpolation linéaire des résultats. J'obtiens la même chose lorsque je produis en VTK. Dans un autre code avec …
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.