Différences
Ci-dessous, les différences entre deux révisions de la page.
| Les deux révisions précédentes Révision précédente Prochaine révision | Révision précédente | ||
|
developpement:productions:logiciels:aquasol:optimisation [2013/06/17 09:55] cicaluga [Optimisation des méthodes d’approximation numérique déjà implémentées. Mise en œuvre d’autres méthodes et de nouveaux algorithmes de discrétisation] |
developpement:productions:logiciels:aquasol:optimisation [2015/01/07 10:04] (Version actuelle) |
||
|---|---|---|---|
| Ligne 3: | Ligne 3: | ||
| ====== Optimisation des méthodes d’approximation numérique déjà implémentées. Mise en œuvre d’autres méthodes et de nouveaux algorithmes de discrétisation ====== | ====== Optimisation des méthodes d’approximation numérique déjà implémentées. Mise en œuvre d’autres méthodes et de nouveaux algorithmes de discrétisation ====== | ||
| - | * **Contact** : Ralf Everaers, Cerasela Calugaru, Sam Meyer\\ | + | {{:developpement:productions:figureaquasol.png?200 |}} |
| - | * **Objectif** : Optimisation (nouveaux algorithmes, parallélisation), implémentation des nouvelles fonctionnalités, validation et tests | + | * **Contact** : Cerasela Calugaru, Ralf Everaers, Sam Meyer\\ |
| + | * **Objectif** : Optimisation (nouveaux algorithmes, parallélisation), implémentation des nouvelles fonctionnalités, analyse du comportement numérique et informatique du code AquaSol | ||
| ➢ **Approximation des dérivées partielles premières** | ➢ **Approximation des dérivées partielles premières** | ||
| Ligne 13: | Ligne 13: | ||
| ➢ **Méthode « focusing »** | ➢ **Méthode « focusing »** | ||
| - | + | {{:developpement:productions:logiciels:aquasol:figfocusing2.png?200 |}} | |
| Cette méthode permet de « zoomer » successivement la résolution du problème de type Poisson-Boltzmann autour de la molécule. Dans un premier temps, elle a été implémentée de la manière la plus non intrusive possible, à savoir en minimisant les modifications apportées au code. Les domaines de simulations successifs sont pilotés via un script shell qui les traitent l’un après l’autre. Le passage d’une simulation plus « grossière » à une simulation plus « fine » se fait par une transmission des valeurs du potentiel obtenues à l’intérieur du domaine grossier (correspondant à la frontière du domaine fin) comme conditions aux limites pour le potentiel pour le domaine fin. Cela a nécessité l’implémentation des procédures pour déterminer la correspondance des points de cette frontière entre les deux maillages. Ainsi, deux techniques ont été implémentées : une première basée sur la recherche globale du plus proche point qui peut traiter tout type de maillage mais qui s’avère pénalisante en terme de temps de calcul et une seconde technique basée sur une méthode de splitting et que nous avons développée pour traiter cette recherche en deux étapes - une recherche 1D suivie d’une autre 2D et qui s’avère très efficace dans le cas particulier des maillages cartésiens (seuls qui sont actuellement pris en compte par AquaSol). Les résultats obtenus jusqu’à maintenant avec cette méthode sont très encourageants et incitent la poursuite des travaux dans cette direction avec une généralisation des ces techniques et un rajout des autres fonctionnalités. | Cette méthode permet de « zoomer » successivement la résolution du problème de type Poisson-Boltzmann autour de la molécule. Dans un premier temps, elle a été implémentée de la manière la plus non intrusive possible, à savoir en minimisant les modifications apportées au code. Les domaines de simulations successifs sont pilotés via un script shell qui les traitent l’un après l’autre. Le passage d’une simulation plus « grossière » à une simulation plus « fine » se fait par une transmission des valeurs du potentiel obtenues à l’intérieur du domaine grossier (correspondant à la frontière du domaine fin) comme conditions aux limites pour le potentiel pour le domaine fin. Cela a nécessité l’implémentation des procédures pour déterminer la correspondance des points de cette frontière entre les deux maillages. Ainsi, deux techniques ont été implémentées : une première basée sur la recherche globale du plus proche point qui peut traiter tout type de maillage mais qui s’avère pénalisante en terme de temps de calcul et une seconde technique basée sur une méthode de splitting et que nous avons développée pour traiter cette recherche en deux étapes - une recherche 1D suivie d’une autre 2D et qui s’avère très efficace dans le cas particulier des maillages cartésiens (seuls qui sont actuellement pris en compte par AquaSol). Les résultats obtenus jusqu’à maintenant avec cette méthode sont très encourageants et incitent la poursuite des travaux dans cette direction avec une généralisation des ces techniques et un rajout des autres fonctionnalités. | ||