Home | Project | Software | Gallery | News | Development | Use | Documentation | Search/Ask
 
News archive
Developers
 
 
© 1999-2016 ONERA Terms of use
  

Development news archive



30/04/02 : Bilan Apps de la version V2.2.13  

From Anne Jollès.
La base des tests de developpement Apps a ete mise a jour avec la version v2.2.13. (cvs rtag: I2213e) Passage de la base pour les modes de production: sgiet sgi_dbg. Le fichier apps2213_sgi.gz est le resultat du passage de Apps pour le mode sgi . Le fichier apps2213_sgi_dbg.gz est le resultat de son passage en mode sgi_dbg. Tests corriges depuis le bilan Apps de la version v2.2.08 : --------------------------------------------------------- TmoSquaredNozzle3D.C (sgi et sgi_dbg), GlobBorderIsConnexTest.C (sgi et sgi_dbg), tous les tests de Lhs pour le mode sgi, les tests de Lhs de type LhsLuRelaxMat pour le mode sgi_dbg. Tests toujours non corriges depuis la v2.2.08 (sgi et sgi_dbg) : -------------------------------------------------------------- SouDts.C, OperGradMultiDom.C. Tests toujours non corriges depuis la v2.2.08 (sgi_dbg) : ------------------------------------------------------- Les tests de Lhs couples 7 equations. ------------------------------------------------------------------ HISTORIQUE DES DERNIERS CAS TESTS INTRODUITS DANS APPS ------------------------------------------------------ Cas tests introduits depuis la v2.2.08 : -------------------------------------- BndStageRedCoinc.C : S. Plot (Interface instationnaire avec reduction du nombre de canaux) BndStageRedNoCoinc.C : S. Plot (Interface instationnaire avec reduction du nombre de canaux) PcmGraphTest.C : X. Juvigny (Graphe de communication) Cas tests introduits dans la v2.2.08 : ---------------------------- ------- BndVanne.C : A. Roure, L. Castillon (C.L. de vanne) JoinNoMatch.C : M. Montagnac (Patch grid) JoinNoMatch2.C : M. Montagnac (Patch grid) TurInterpolWallLine.C : R. Houdeville (Criteres de transition non locaux) TurWallLine2.C : R. Houdeville (Criteres de transition non locaux) TurTuEffect.C : R. Houdeville (Criteres de transition locaux) Cas tests introduits dans la v2.2.05 : ------------------------------------ BlkInterpCyl.C : C. Benoit (Interpolation grilles cylindriques) GlobBorderConnexCompFrames.C : X. Juvigny (Super-frontieres connexes) GlobBorderConnexIteratorTest.C : X. Juvigny (Super-frontieres connexes) GlobBorderIsConnexTest.C : X. Juvigny (Super-frontieres connexes)


29/04/02 : Non-regression pour la version V2.2.13  

From Anne Jollès.
La base de non-regression se trouve maintenant dans : /data/Public/NoReg. Vous pouvez acceder au Makefile a utiliser. Les references v2.2.13 de la non-regression se trouvent dans: /data/Public/NoReg/Ref2213_r8 pour le mode de production sgi_r8, et /data/Public/NoReg/Ref2213 pour le mode de production sgi, sauf pour les cas tests des repertoires suivants: - Join - OAT15A - 2Cell - RAE2822 - caimenWingTurb - shockTube - bluntCyl qui correspondent malencontreusement a un passage en mode sgi_r8.

26/04/02 : Mise à jour de la base de validation  

De la part de Anne Jollès

La base de validation a ete mise a jour avec la v2.2.13 et un tag
lui a ete appose. Ce tag est le V2213.


26/04/02 : Production de la version V2.2.13  

From Michel Gazaix.
Une nouvelle production elsA est disponible : ELSAVERSION = v2.2.13 cvs rtag : I2213c Date : 26 Avril 2002 Cette production a ete compilee avec : ELSAPROD = sgi, sgi_dbg, sgi_r8, sgi_mpi, nec_r8, nec_r8_mpi, linux, hp_mpi ---------------------------- NOUVELLES FONCTIONNALITES ---------------------------- 1. M. Lazareff: Evolutions importantes de l'interface utilisateur [python-Utilisateur, interface graphique] ------------ rtag I2213a ------------- ---------------------------- CORRECTIONS BUGS, PORTAGE ---------------------------- 1. M. Gazaix: Verification fenetre dans DesWindow::createWindow (D02-11-10) ---------------------------- DIVERS (Taches realisees par les membres de l'unite elsA) ---------------------------- 1. : Generation documentation Doc++ 2. : Mise a jour Apps (test unitaires) 3. : Passage Base non regression (sgi, sgi_dbg, sgi_r8) 4. : Passage Base de Validation (sgi, nec_r8) et realisation Document de Validation (PS et HTML) 5. : Mise a jour du Manuel Utilisateur


19/04/02 : Production de la version V2.2.12  

From Michel Gazaix.
Une nouvelle production elsA est disponible : ELSAVERSION = v2.2.12 cvs rtag : I2212d Date : 19 Avril 2002 Cette production a ete compilee avec : ELSAPROD = sgi, sgi_dbg, sgi_r8, sgi_mpi, nec_r8, linux, hp, hp_r8_mpi ---------------------------- CORRECTIONS BUGS, PORTAGE ---------------------------- 1. J.C. Boniface: Decentre a 5+n equations (FxcCompDlF.for)(D02-15-02) 2. J. Peter: Lurelax + laminaire (D02-13-18) 3. M. Gazaix: Portage FUJITSU (SioApi.C et Make_fuji.mk) (D02.14.01) 4. M. Gazaix: Extraction temperature en Euler (D02-13-15) 5. M. Montagnac: Extraction fmt_tp en linux_r8 (D02-13-12) ------------ rtag I2212a ------------- 6. S. Heib: Correction restrict_vfield=2 (D02-15-08) 7. S. Heib: Regression CPU STEP-ASM-TL (D02-13-20) 8. M. Gazaix: printf dans DefNumError (D02-11-11) 9. S. Heib: Bug GeoConnect::getIJK 10. S. Heib: Bug createLevel.C et query (D02-11-01) ------------ rtag I2212b ------------- 11. S. Beaugrand:Bug extraction aux noeuds du frottement, sur "certains" plans (D02-16-02) 12. S. Beaugrand:Extraction mut/mu aux noeuds (D02-13-08) 13. S. Beaugrand:Priorite dans les extractions aux conditions de parois pour les moyennes entre CL (D02-13-05) ---------------------------- MODIFICATIONS INTERFACES ---------------------------- 1. 'period', 'nbcoarsegrid': Suppression valeur par defaut 2. 'comp_bary': Remplacer par jcomp_bary 4. 'minmod', valbada : au lieu de 'min', 'val' 4. 'superbee' : Nouveau limiteur 4. 'coquel_d', 'coquel_i' : Nouveaux flux convectifs (decentres) 5. 'config' : Calcul 2D/3D/Axi (possibilite un seul plan en 2D et en Axi) (D01-48-06) 6. 'cell2node_coarse': : Remplace 'is_geometric' ---------------------------- OPTIMISATION CPU ---------------------------- 1. : Complements Vectorisation NEC Lussor

17/04/02 : Adresse support développeur  

Une nouvelle adresse:  elsA-infodev@onera.fr 
vient d'être créée pour les développeurs d'elsA souhaitant adresser à 
l'équipe elsA des courriers concernant le développement. 

Cette adresse s'ajoute à l'adresse de la liste de discussion 
elsA-dev@onera.fr permettant les échanges directs entre développeurs. 

Tout comme l'adresse elsa-info@onera.fr a pour mission de répondre 
aux questions des utilisateurs et aux difficultés qu'ils rencontrent, 
l'adresse elsa-infodev@onera.fr accueillera toute question 
et toute information concernant le développement d'elsA. 

Plus précisement, cette adresse devra être utilisée: 

 -1- pour les demandes d'intégration et échanges à ce propos, 
 -2- pour tous les échanges de courrier concernant les 
     revues d'intégration et en particulier les remarques de relecture, 
 -3- pour toute question concernant le support développeur : 
     - compréhension et apprentissage du code, 
     - méthode de développement, 
     - tout problème rencontré au cours d'un développement. 

Cette adresse facilitera ainsi le partage de l'information 
à l'interieur de l'équipe elsA et permettra une meilleure continuité 
de l'assistance aux développeurs. 
  

Pour l'instant cette adresse est renvoyée sur tous les membres 
permanents de l'unité DSNA/ELSA.

16/04/02 : Mise à jour de la base de validation  

De la part de Sebastien Heib

La base de validation a ete mise a jour et un tag lui a ete appose.
Ce tag est le V2210.


15/04/02 : Production de la version V2.2.11  

From Michel Gazaix.
Une nouvelle production elsA est disponible : ELSAVERSION = v2.2.11 cvs rtag : I2211f Date : 15 Avril 2002 Cette production a ete compilee avec : ELSAPROD = sgi, sgi_dbg, nec_r8, linux, hp_r8, hp_r8_mpi, hp_r8_dbg ---------------------------- NOUVELLES FONCTIONNALITES ---------------------------- 1. F. Cayre : Reduction encombrement memoire (calcul des gradients) (D02.08.02) ------------ rtag I2211a ------------- 2. C. Benoit, X. Juvigny, G. Jeanfaivre: Chimere + multigrille et diverses ameliorations chimeres (D02.10.01) ------------ rtag I2211b ------------- 3. A. Jolles: Reduction encombrement memoire [constitution Rhs], RI-NI-02-04-02 (D01-44-04) ------------ rtag I2211d ------------- 4. S. Plot, A. Roure: Complements extraction multi-etage instationnaire (D02-13-19) 5. M. Costes: Couplage faible HOST (D02-15-01) 6. R. Houdeville: Calcul des epaisseurs integrales dans les coins (D02-15-03) ------------ rtag I2211e ------------- 7. M. Costes: Couplage faible HOST (D02-15-01) ------------ rtag I2211f ------------- ---------------------------- CORRECTIONS BUGS, PORTAGE ---------------------------- 1. S. Beaugrand: Correction bug extraction liste de points (apparu dans la v2.2.08)(D02.14.04) 2. M. Gazaix: 'wallslip': Assoc. FxcUpwind + BndEuler a partir de l'interface( 3. M.C. Le Pape: Correction bug Chimere (D02-11-04, D02.14.01) 4. C. Benoit: Correction bug ALE (Tbx/Function/TbxsoftblademeshDeformationF.for) (D02.14.03) 5. M. Montagnac: bug LussorSca et LussorMat (diffusiveSpectralRadiusBorder) (D02.14.02) 6. J. Delbove: JoinMatchParP.C (D02.11.02) 7. M. Gazaix: Assert GeoConnect, FxcUpwind, BndNS.C (mode Debug) (D02-11.03- D02-13-11) 8. J. Peter: van Albada (D02.11.09) 9. A. Roure: Raccords partellement coincidants (D02-11-05, D02-11-06, D02-11-07) 10. M. Gazaix: Corrections linux. HP. 11. S. Heib: Changement codage flux de Roe (E_MIN_SURFACE) (D02-14-05). 12. A. Roure: Multigrille + nearmatch (D02-15-06) 13. R. Houdeville: Vieux bug FactSio.C (coeffmutinit, utiliser BndNS au lieu de BndSupOut) 14. A. Jolles: Lussor + periodicite (D02-13-23) ---------------------------- MODIFICATIONS INTERFACES ---------------------------- 1. 'extrap=-1' : Force l'appel a BndFirstOrderExt si 'wallslip' et schema decentre 2. 'winnow_type' : Remplace 'typevanne' (egalement: "winnow_ref_pres", "winnow_ref_mflow", "winnow_relax") 3. 'full_cooling' : 'porosity_file' au lieu de 'zone_file' (egalement: "active_control' au lieu de 'activecontrol", ...) (D02-11-13) 4. 'jtopo', 'periodic', 'noperiodic' : raccord partiellement coincidant ('ptype' devient obsolete) 5. 'specrad_s_mean' : valeur par defaut modifiee (moyenne arithmetique au lieu de moyenne de Roe) --> gain CPU ---------------------------- OPTIMISATION CPU ---------------------------- 1. M. Montagnac : Vectorisation NEC Lussor (D02-13-14)


01/03/02 : Production de la version V2.2.10  

From Michel Gazaix et Anne Jollès.
Une nouvelle production elsA est disponible : ELSAVERSION = v2.2.10 cvs rtag : I2210j Date : 1 mars 2002 Cette production a ete compilee avec : ELSAPROD = sgi, sgi_r8, sgi_dbg, nec_r8, linux, linux_mpi, hp_r8 Reference de la base de non-regression ElsaNoReg disponible pour les modes de production sgi : Ref2210 et sgi_r8 : Ref2210_r8 ---------------------------- NOUVELLES FONCTIONNALITES ---------------------------- 1. L. Castillon / R. Houdeville: Lois de paroi + reperes mobiles (D01.50.05) 2. R. Houdeville: Corrections Spalart (I2210a) 3. J.C. Boniface: Complements Actuator Disk (D01.42.05) 4. J.C. Boniface: Definition axes helico en ALE et chimere (D02.04.01) 5. J.C. Boniface: Tourbillon prescrit dans une grille(D01.25.08) (I2210b) 6. M.C. Le Pape: Chimere : Correction 2 erreurs + extraction champ de points masques (D02.04.08) 7. V. Brunet: Conditions de porosite type Darcy (D01.50.06) 8. X. Juvigny: Graphes des communications entre processeurs (I2210c) 9. R. Houdeville: Blindages SST (D02.02.01) 10. S. Heib : Condition aval de debit global (D01.50.01) (I2210d) 11. F. Cayre (SNECMA) : Vectorisation Turbulent Michel (delta3v_) 12. M. Montagnac / J. Delbove : Divers HMR, Baldwin, MPI + Corrections Bugs (D02.04.02) (I2210e) 13. R. Houdeville : Vectorisation loi de paroi (bndinvwlv_) (D01.04.03) (I2210f) 14. R. Houdeville : Corrections modele K-omega AVTAC (variante Kok) (D02.08.03) (I2210g) 15. S. Plot / A. Roure : Interface instationnaire avec reduction du nombre de canaux (COMET3) (D01.25.10) (I2210g) 16. S. Beaugrand : Extraction Flux en instationnaire (D01.48.03) (I2210i) 17. A. Jolles: Association LUSSOR / Periodicite (D02.06.02) (I2210j)


18/01/02 : Production de la version V2.2.09  

From Michel Gazaix et Anne Jollès.
Une nouvelle production elsA est disponible : ELSAVERSION = v2.2.09 cvs rtag : I2209h Date : 18 janvier 2002 Cette production a ete compilee avec : ELSAPROD = sgi, sgi_dbg, nec_r8, linux, linux_mpi, hp_r8 ---------------------------- NOUVELLES FONCTIONNALITES ---------------------------- 1. M. Costes : Extraction corde de section de pale (D50.02) 2. JC Jouhaud : Complements AMR (D50.04) 3. JC Boniface : Complements Actuator Disk (D42.05) 4. MC Le Pape : Initialisation et Extraction Chimere (D12.05) -------------- CORRECTIONS -------------- 1. AJ : Association pas de temps dual + lois de paroi 2. AJ : Bug instationnaire visqueux + mouvement corps rigide dans la classe BndNSRelFrame (cf mail a elsa-dev: [elsa-dev] Bug v2.2.08 : instationnaire visqueux + mouvement corps rigide) 3. RH : Correction de la c.l. de pression subsonique en basse vitesse (BndWbPresLowspeedF.for) 4. AR : Bug OperGradGradVelocity (dimension wb1) 5. AR : Bug raccord "NoMatch" en mode DEBUG (JoinNoMatchLineP.C, TmoPrepareNoMatchLine.C) 6. SB : Bug "pinf+frottement" (objets partages) 7. SB : Indices fichier tecplot pour extractions sur interfaces 8. SB : Bug iterations sur grandeurs partiellement integrees 9. SB : PR-2000-07-08 : bug sur certaines configurations en parallele pour des extractions utilisant les raccords (node et cellfict) 10.SB : initialisation du SioZoneElem dans SioComputeBase ----------------------- AMELIORATION QUALITE ----------------------- 1. AJ: Amelioration codage (Geo) ------------------------------------------- MODIFICATIONS DE L'INTERFACE UTILISATEUR ------------------------------------------- I - Nouvelles cles Pythons Utilisateur -1- Developpement JCB : D42.05 "Complements actuator-disc". La prise en compte de l'actuator-disc peut actuellement se faire de deux facons differentes: -a- en utilisant le nouveau developpement. Dans ce cas le disque d'action est considere comme un raccord coincident specifique. Un nouvel attribut a ete rajoute aux c.l. de type raccord: 'actuator_type'. L'attribut 'reference_state' n'est pas utilise. Pour la frontiere "Lower": bnd.set('type','join') ---> type de la c.l. bnd.set('jtype','match') bnd.set('blkrac':'Blk0003', 'wndrac':'Bnd0022_Blk0003') bnd.set('actuator_type','actuatorhelico') ---> nouvelle cle bnd.set('file','Force_Actu_1.dat') bnd.set('format','fmt_tp') Et pour la frontiere "Upper": bnd.set('type','join') ---> type de la c.l. bnd.set('jtype','match') bnd.set('blkrac':'Blk0001', 'wndrac':'Bnd0005_Blk0001') -b- toujours comme avant (c'est a dire jusqu'a la v2.2.08 comprise). Dans ce cas la condition de disque d'action est traitee comme une condition limite physique : Les cles Python utilisateur n'ont pas change (MU p 72,80,85). Pour la frontiere "Lower": bnd.set('type','actuatorhelico') ---> type de la c.l. bnd.set('jtype','match') bnd.set('blkrac':'Blk0003', 'wndrac':'Bnd0022_Blk0003') bnd.set('reference_state','lower') bnd.set('file','Force_Actu_1.dat') bnd.set('format','fmt_tp') Et pour la frontiere "Upper": bnd.set('type','actuatorhelico') ---> type de la c.l. bnd.set('jtype','match') bnd.set('reference_state','upper') bnd.set('blkrac':'Blk0001', 'wndrac':'Bnd0005_Blk0001') La nouvelle methode -a- est conseillee, etant plus stable numeriquement. Elle possede l'avantage tres interessant d'autoriser l'association avec le preconditionnement aux basses vitesses, ce que ne permet pas la methode -b-. La methode -b- a pour l'instant ete conservee a des fins de non regression. -2- Developpement MCLP (Chimere): "Reglage des parametres de la grille cartesienne d'interpolation". -2.1- Nouvelle methode d'information display() accessible a l'utilisateur pour les objets block, mask, mesh. Ex. d'utilisation: blk.display() -2.2- Nouvelles cles: 'interpol_tool' : preconditioning for interpolation cell search. valeurs possibles : 'icg' ou 'adt' ('icg' par defaut). 'icg_refn_depth' : refinement depth for the 'icg' grid. valeur entiere positive, 5 par defaut. (0 signifie absence de raffinement). 'icg_refn_cut' : refinement cut-off for the 'icg' grid. valeur entiere strictement positive, 30 par defaut. La valeur par defaut valait 20 dans la v2.2.06, ou cet attribut etait alors impose par le noyau. 'icg_dim_factor' : factor for basic icg number of cells. valeur reelle strictement positive, 1. par defaut. Remarque de MCLP: Dans la v2.2.06, l'attribut 'icg_dim_factor' n'existait pas sous cette forme. Les dimensions topologiques de l'icg etaient les memes dans les 3 directions et etaient chacune egale a environ 1.8 du nombre total de points du domaine asssocie. Pour se rapprocher (avec la v2.2.09) de ces valeurs, il faut s'arranger pour que le nombre de cellules de l'icg soit presque les memes pour les deux versions. II - Modification des cles Python utilisateur Developpement MCLP (Chimere): "Reglage des parametres de la grille cartesienne d'interpolation". La cle 'type' n'a plus de valeur par defaut. III - Modification des cles Python API (uniquement pour des cles n'existant pas encore en Python Utilisateur) Developpement RH "Criteres de transition": "transcriteria" a ete modifiee en "trans_crit" "transi_exit" a ete modifiee en "trans_exit" "transcriteria_file" a ete modifiee en "trans_crit_file"


16/01/02 : Bilan Apps de la version V2.2.08  

From Anne Jollès.
La base des tests de developpement Apps a ete mise a jour avec la version v2.2.08. (cvs tag: I2208) Ci-joint 2 listings de l'execution de la base: apps2208_sgi.gz (mode sgi) apps2208_sgi_dbg.gz (mode sgi_dbg) Remarques: (1) Les 3 tests de Descp ont ete remis a jour (depuis le bilan v2.2.05). (2) Le test TmoSquaredNozzle3D.C ne fonctionne plus. (3) Le tests SouDts.C ne fonctionne toujours pas. (4) Le test OperGradMultiDom.C ne fonctionne etrangement qu'en mode sgi_dbg. (5) Le test GlobBorderIsConnexTest.C est deja corrige pour la v2.2.09 (il ne fonctionne pas en mode sgi_dbg pour la v2.2.08) (6) De nombreux tests de Lhs ne marchent pas. Ils ont ete corriges pour fonctionner dans la v2.2.09.


14-18/12/01 : Production de la version V2.2.08  

From Anne Jollès.
Une nouvelle production elsA est disponible : ELSAVERSION = v2.2.08 cvs rtag : I2208l Cette production a ete compilee avec : ELSAPROD = sgi, sgi_r8, nec_r8 Les references de non-regression (ElsaNoReg) correspondant a cette production sont: Ref2208 pour le mode sgi Ref2208_r8 pour le mode sgi_r8. Cette production comprend les developpements suivants: ------------------------------------------------------- 1. J.Delbove : cvs rtag I2208a D14.05 : Complements Ale + SSOR : Modifications apportees a l'implicite Lussor (scalaire et matriciel) permettant de prendre en compte la vitesse de grille (entrainement et deformation). Introduction du nouveau cas de non-regression: nacaProfile/naca64_unsteady_motion10_bw_lussormat_ns_dt.py .................................................................... 2. M.Montagnac : cvs rtag I2208b D41.01 : Complements "patch grid" : Extension multigrille et parallhle des raccords adjacents non cooncidents, Codage algorithme de "clipping"pour le calcul d'intersection de polygones. .................................................................... 3. S.Beaugrand : cvs rtag I2208c 1) D00.31.03 : C.L. avec cartes de valeurs sur grilles grossieres. 2) D25.02 : Extraction fonction d'intermittence a la paroi (utile pour verifier le fichier d'intermittence lorsque la transition est imposee, ou pour connaitre le resultat du critere de transition lorsqu'on calcule la transition). 3) D40.02 : Generalisation de la methode FactSio::initFld(). .................................................................... 4. A.Jolles : cvs rtag I2208d D04.12 : Complements et correction d'un bug pour reduction du nombre de cellules minimum dans les grilles grossieres. .................................................................... 5. H.Pascal-Jenny : cvs rtag I2208e D30.04 : Implicitation LU-SSOR du Spalart-Allmaras. Resolution de l'equation de transport de la turbulence proposee par Spalart-Allmaras a l'aide d'un solveur Backward-Euler + LU-SSOR scalaire. Introduction d'un schema centre avec dissipation artificielle pour traiter les flux convectifs de l'equation de transport (uniquement pour Spalart-Allmaras). Une note de specification-conception NI-01079 est disponible. "Complement au Manuel Theorique d'elsA: Implicitation du solveur de l'equation de transport de Spalart-Allmaras". Introduction de 2 nouveaux cas de non-regression: flatPlateSA/fpSA_ssor.py et RAE2822/rae2822_ssor.py .................................................................... 6. A.Roure : cvs rtag I2208g D42.01 : Probleme "reprises" condition de vorticite. Ce developpement resout le probleme identifie par Frederic Moens (DAAP/ACI). Introduction du nouveau cas de non-regression: OAT15A/OAT15A-Vrt-reprise.py .................................................................... 7. R.Houdeville : cvs rtag I2208h 1) D22.02 : Adaptation modele de Spalart-Allmaras (APAM) : a- Fonction correctrice du terme de production de nu_tilde pour Spalart pour prendre en compte les effets de courbure. APAM b- Optimisation place memoire du calcul de grad( grad(U) ) pour Spalart SARC (evite une augmetation de 40% de la place mmoire). Utilise aussi pour le terme bas Reynolds equation pour epsilon (modele k-epsilon : gain 25% - 30%) 2) D38.02 : Modeles k-omega AVTAC : a- modele k-omega Wilcox : adjonction terme de diffusion croise grad(k).grad(omega) pour supprimer la sensibilite a omega_infini (AVTAC, J.C. Kok) b- modele k-omega de J.C. Kok - condition de paroi soit en 1/y**2 au voisinage de la paroi soit omega_paroi = 18*nu/y**2 - terme de diffusion croise grad(k).grad(omega) pour supprimer la sensibilite a omega_infini - limitation de Prod(k) a 20*epsilon c- modele k-tau de J.C. Kok. Identique au k-omega de J.C. Kok mais l'equation pour omega est resolue en tau=1/omega. Pour l'instant, ne donne pas de bons resultats avec transition imposee. Deux notes internes sont disponibles en version provisoire : NI-01062 : Modeles k-omega et k-tau, NI-01063 : Effets de courbure et de rotation dans le modele de Spalart-Allmaras. Correction bugs: ------------------- J.Delbove: correction d'un bug dans le calcul du pas de temps pour tous les cas avec mouvement (mouvement de corps rigide et/ou deformation). J-C Boniface: correction d'un bug pour la c.l. pression subsonique avec preconditionnement basse vitesse. Amelioration du source: ----------------------- A.Jolles: - blindage c.l. glissement et preconditionnement - simplifications methodes Fld M.Montagnac: - nettoyage methode Tmo cvs rtag I2208j -------------------------- Premiere mise à disposition de la production 2.2.08, le 14 décembre 2001. cvs rtag I2208l -------------------------- Deuxième mise à disposition de la production 2.2.08, le 18 décembre 2001, comprenant des modifications mineures supplémentaires concernant les fichiers : cfg/prods/Make_nec.mk et api/Py/EpAttrDefs.py . Modification cles Python API: ---------------------------- Attention les cles suivantes (concernant la reduction de viscosite artificielle) n'existent plus. << 'av_lower' 'av_lower_exp' 'av_lower_var' 'av_lower_varinf' >> Elles ont ete remplacees par les cles : << 'av_lsr_dir' 'av_lsr_exp' 'av_lsr_var' 'av_lsr_val' >> Consequence au niveau Python utilisateur: ---------------------------------------- Il faut utiliser la-aussi les nouvelles cles 'av_lsr_dir', 'av_lsr_exp', 'av_lsr_var', 'av_lsr_val' . On peut cependant continuer d'utiliser les cles preconisees par le Manuel Utilisateur c'est a dire: 'av_lsr_red' 'av_lsr_red_exp' 'av_lsr_red_var' 'av_lsr_red_val' MAIS en utilisant la commande: elsa -f script.epy --allow_obsolete.


26/11/01 : Production de la version V2.2.07  

From Michel Gazaix.
Une nouvelle production elsA est disponible : ELSAVERSION = v2.2.07 cvs rtag : I2207a Cette production a ete compilee avec : ELSAPROD = sgi, nec_r8 Cette production n'apporte pas de nouvelles fonctionnalites. ------------------- Corrections : ------------------- Bug Turbulent Michel : < IF (test(m)) STOP 'Tur/Alg/TurAlgDeltaT-3F.for : bad test(m)' --- > IF (test(m)) EXIT ---------------------------------------------------------- Bug Constructeur RhsEquation < _rhs(grid.getNbCell(), eqnNumber), --- > _rhs(TYPE_CELL, grid.getNbCell(), eqnNumber), (Avait comme consequence l'arret du calcul en mode DEBUG lorqu'on demandait residual_type != 0) ---------------------------------------------------------- bug DefError.C: < if (mesg) cerr << mesg << endl; --- > if (mesg) cerr << *mesg << endl; ---------------------------------------------------------- Quelques modifications du traitement des cutoffs sur les valeurs conservatives ---------------------------------------------------------- Suppression de pas mal de divisions par zero potentielles dans les FORTRAN de Bnd Exemple : dans BndSymF.for on remplace < invsfn = ONE/snorm(ling) par : > invsfn = ONE/MAX(snorm(ling),E_MIN_SURFACE) Par souci d'uniformisation, j'ai remplace les cutoff dans Fxc par E_MIN_SURFACE Exemple dans FxcMatNumDissDirF.for : < nosurf = MAX(E_CUTOFF,surfno(indint)) --- > nosurf = MAX(E_MIN_SURFACE,surfno(indint)) ---------------------------------------------------------- Pour portage IBM, remplacement d'une constante entiere par un enum dans OperBase.h : < static const E_Int E_MAX_OPERTOCOMP = 7; --- > enum { E_MAX_OPERTOCOMP = 7 }; ---------------------------------------------------------- Correction de nommage : < magic(restrict_visc_fields) --- > magic(restrict_vfields) ---------------------------------------------------------- Suppression de pas mal de valeurs par defaut "en dur" (regroupement dans EpKernelDefVal.py) ----------------------------------------------------------


20-21/11/01 : Production de la version V2.2.06  

From Michel Gazaix.
Une nouvelle production elsA est disponible : ELSAVERSION = v2.2.06 cvs rtag : I2206k Cette production a ete compilee avec : ELSAPROD = sgi, sgi_dbg, sgi_r8, sgi_mpi, nec_r8, nec, hp, hp_r8_mpi, linux, linux_mpi, dec, dec_opt, dec_opt_mpi, dec_opt_r8_mpi Cette production apporte de nouvelles fonctionnalites importantes, en particulier dans le cadre du DTP CAIMEN : .................................................................... 1. RH : cvs rtag I2206a Transition (PRC COMET) Concerne les criteres de transition non locaux (Note NI-01026) avec en plus: - effet de la turbulence exterieure sur la couche limite laminaire (aubes de turbines) - taux de turbulence exterieur variable (aubes, profils) - application de la fonction d'intermittence uniquement dans les couches limites (application aux geometries complexes, aux multi-corps) .................................................................... 2. SB : cvs rtag I2206b 2.1- Diverses corrections de problemes En particulier : a) D39.05 : extraction loc=node dans le cas d'un raccord et en parallele, b) correction portant sur la prise en compte de pinf pour l'integration des flux 2.2- Ecritures periodiques de fichiers pour extractions periodiques de type "group" qui n'avaient pas ete envisagees (writingmode=3); (voir SioWriteCommand::execute) 2.3- D39.04: Possibilite de faire des extractions de type flux en parallele separees par processeur (reference : PR-2001-09-23) (voir FactSio::buildSioFilters). 2.4- D35.02: Extraction node pour contraintes (voir SioFilterFlux::writeLocalValues ) (sans application des raccords) 2.5- D39.02: Fichier commun pression et contraintes - voir SioFilterFlux::writeLocalValues - sans application des CL - comme prevu cette nouvelle possibilite a un peu degrade la qualite du code, de plus un warning apparait. Notamment, lorsque "loc=1" pour un DesExtractGroup on passe outre les verifications sur les noms de variables. 2.6- D39.06 Extraction des normales (voir SioFilterFlux::writeLocalValues) 2.7- D40.05 Extraction du temps pour l'instationnaire ( voir SioFilterBase::whenReady SioFilterFlux::writeGlobalValues SioWritePoint::Update ) .................................................................... 3. AJ : D41.02 : cvs rtag I2206c Ameliorations pour implicitation terme source en turbulent (dont correction Bug formule 3.200 MT) Remarque : il semble que la correction de ce bug entraine la divergence du calcul de P. Viscat (2D KL, avec preconditionement et implicitation terme source) : a suivre ... .................................................................... 4. L. Castillon (& AR) : D36.01 : cvs rtag I2206d Condition limite "vanne" .................................................................... 5. AJ : D04.12 : cvs rtag I2206e Reduction du nombre minimum de cellules dans les grilles grossieres (de 3 a 2) Ce developpement est limite a l'adaptation de l'operateur FxcScaNumDiss dans sa 1ere variante (d3Type = FxcScaNumDiss::rspec_d3), avec raccords coincidants. Il faut donc utiliser un schema centre avec dissipation artificielle scalaire de type 4eme ordre: 'rspec_d3', et des raccords coincidents pour resoudre: - un systeme a 5 equations (Euler, Navier Stokes laminaire ou turbulent avec uniquement modele algebrique) sur des grilles pouvant n'avoir que 2 cellules (grilles fines ou grilles grossieres en multigrille) - le systeme decouple des equations de Navier Stokes turbulentes avec equations de transport en multigrille, les grilles grossieres pouvant n'avoir que 2 cellules. Par contre, toute grille fine doit posseder comme avant au moins 3 cellules dans chaque direction. La note NI-01053 presente la methode retenue et les choix effectues. .................................................................... 6. AJ : D42.04 : Basses vitesses et Vorticite : cvs rtag I2206f .................................................................... 7. CB/GJ : D44.02 : Reduction CPU chimere : cvs rtag I2206f .................................................................... 8. JP : CAIMEN D30.05 : cvs rtag I2206g * multigrille ------------- - restriction semi-synchrone : dans cette variante de restriction, on ne recalcule que le flux convectif (qu'on ajoute aux autres flux - 'asynchrones'- calcules precedemment) avant la restriction. [Cette variante est parfois tres fructueuse avec les schemas Backward-Euler] - appel a la phase implicite a la fin de la prolongation. Au lieu ajouter l'increment prolonge de grille fine de la methode standard (dwst), on appelle la phase implicite avec dwst comme membre de droite --> il en ressort un dwi qui est ajoute a la place de dwst. [ Cette variante existait dans canari, elle n'ameliore pas la convergence des calculs qui fonctionnent sans y recourir. Elle permet, parait-il, de demarrer et de d'obtenir cetains calculs tres complexes ] * Implicites LURELAX -------------------- Suite du developpement des implicites LURELAX pour des problemes decouples (5+2) equations avec notamment des linearisations sur cinq points du flux visqueux explicite. Desormais, six variantes sont disponibles : RELUS3 (qui existait deja) RELUM3 (qui existait seulement pour 5 equations sans terme visqueux), RELUMS3, RELUM5, RELUM5, RELUS5 (nouvelles) [ Bonne efficacite des phases implictes RELUM5, RELUMS5 RELUM3 pour les calculs a viscosite matricielle, bonne efficacite de RELUS5 pour les calculs a viscosite scalaire ] .................................................................... 9. D. Jouvenot : RVSLIB : cvs rtag I2206h Possibilite de realiser directement des films sur cronos, sans passer par un traitement de fichiers (Tecplot ou autre). .................................................................... 10. MG : cvs rtag I2206i 10.1 : D43.01 Ameliorations Schemas decentres (vectorisation, introduction 2 nouveaux flux : coquel_d, coquel_i) 10.2 : D42.02 Cutoffs et Limitation increments (A terminer, non utilisable actuellement en turbulent) 10.3 : D42.03 Tecplot binaire (en ecriture) Disponible sur sgi (oneroa13, akira), Compaq (feder), HP (octopus), linux (oneroa36) Non disponible sur NEC (cronos) 10.4 : D40.03 Test Metrique et message volumes negatifs (methode DesBlock::show()) Concerne uniquement les grilles fines. 10.5 : D42.06 2D plan ou axi avec un seul plan ------------------- Corrections BUG : ------------------- F. Cayre (Snecma), AR+GB : Bug modele de Michel MG : Bug basse vitesse si interface de surface nulle) (calcul de P. Viscat) OR/XJ : Bug rayon spectral turbulent K-Epsilon SB : Correction 'period' en cas de reprise SB : extraction sur liste de points en double precision -------------------------------------- Ultime correction : cvs rtag I2206k -------------------------------------- MG : correction d'une erreur dans la preparation de la version // MPI de la v2.2.06. Il faut changer, dans Sio/Api/SioApi.C, un '||' par un '&&'. --------------- Divers --------------- Passage sur COMPAQ (Serveur DAAP (Lille) "feder") en mode parallele (MPI) : sur un cas 2D, 40 blocs, 3,7 Millions de points, turbulent KL multigrille (2 grilles grossieresa), en simple precision : dec_opt : 199 secondes par iteration dec_opt_mpi : 132 secondes par iteration (2 processeurs, reparttion "au pif", probablement ameliorable) nec : 37 secondes par iteration (800 MegaFloips) Pour ceux que cela interesse, vous pouvez consulter le profiler correspondant au calcul sur un processeur.


05/11/01 : Production de la version V2.1-p4  

From Anne-Marie Vuillot.
La production v2.1.15-p4 correspondant au patch V2.1-p4 est disponible sur NEC . Un executable est disponible : elsA.x (sans profiler).


18/10/01 : Conception / modèle objet  

From Anne Jollès.
Une action generale d'evaluation et d'amelioration de la conception dans elsA est actuellement en cours. M. Depitre, consultant de la societe Softeam, travaille en etroite collaboration avec nous afin de realiser un audit general de la conception, et dans un premier temps de nous aider a mettre a jour le modele objet. Ce nouveau modele sera realise a l'aide de l'outil Rational Rose, a partir d'une operation de reverse-engineering sur le code C++ de la version v2.2.05 d'elsA.


19/09/01 : Correction de la version V2.2.05  

From Michel Gazaix.
Philippe Girodroux-Lavigne a cru bon d'essayer de produire elsA v2.2.05 sur une station HP, en mode MPI (ELSAPROD=hp_mpi_r8) Il a mis en evidence que la v2.2.05 ne compilait plus, pour plusieurs raisons: 1. un 'explicit' malheureux devant le constructeur de FldSlice. C'est le compilo HP qui a raison. Par consequent, sgi, g++, sxc++, cxx avaient tort : nous sommes en bonne compagnie (enfin ca depend des gouts) (cela dit, pas de panique : le code genere est identique). 2. des 'const_iterator' que le compilo HP n'accepte pas: il prefere 'iterator'. Pour l'instant, on regle ce pb en utilisant la macro E_CONST_ITERATOR (qui s'expand en 'iterator' sur HP, 'const_iterator' ailleurs). Je ne sais pas si le compilateur HP a encore raison sur ce point (jusqu'a maintenant, nous avons observe que le bougre se trompe rarement...). Si oui, il vaudrait mieux, a terme, remplacer les 'const_iterator' concernes par 'iterator'. 3. l'utilisation de la fonction puissance 'pow' demande a etre precisee sur HP, on a introduit les cast qui vont bien. 4. le parallele MPI ne compilait plus (Bien sur, on aurait du essayer avant d'annoncer la production v2.2.05, meme si je n'avais pas explicitement annonce que sgi_mpi etait disponible...). Les fichiers modifies sont : Tmo/SetOfSolver/TmoSetOfSolver.C Sio/Compute/SioComputeNode.C D'autre part, trois bugs ont ete trouves depuis hier : createBndSubInOut (intervient dans certains calculs actuator-disk) FxcRoeCorr.C: intervient seulement dans les calculs chimeres turbulents. FxcUpwind.C : ne compile pas si _E_FORTRAN_LOOPS_ n'est pas defini DesWindow.C : tous les cas paralleles --> core dump J'ai hesite entre creer une production v2.2.06 incorporant ces corrections, ou produire une version v2.2.05 avec un nouveau tag : I2205h. Les deux solutions ont leurs avantages et inconvenients... Je choisis pour cette fois la deuxieme solution, ce qui signifie concretement: 1. que le repository CVS est modifie; si vous vous etiez synchronise des hier, vous avez probablement une re-synchro mineure (et donc indolore) a effectuer. 2. ceux qui utilisent les header et les librairies installes sur oneroa13 dans /v2.2.05/Dist/bin verront automatiquement les nouvelles librairies (je m'occupe de recompiler et de re-installer). En promotion, j'en profite pour ajouter la production sgi_mpi. 3. sur NEC, je remplace les 2 executables (elsA.x et elsA_prof.x) a priori, de toute facon, il ne doit y avoir aucune difference (a part les calculs chimeres turbulents, qui etaient reserves de toute facon a un duo d'intrepides calculateurs). PS : pour les curieux, je joins les differences entre le tag CVS I2205g (correspondant a l'annonce d'hier), et le tag CVS I2205h (correspondant au repository aujourd'hui).


18/09/01 : Bilan Apps de la version V2.2.05  

From Anne Jollès.
La base des tests de developpement Apps a ete mise a jour avec la version v2.2.05. (cvs tag: I2205) Un listing de l'execution de la base en mode de production sgi est donne. Remarques: (1) les 3 tests de Descp ne marchent plus. Leur remise a jour est prevue dans les jours qui viennent. (2) Les tests SouDts.C et OperGradMultiDom.C ne fonctionnent egalement plus. L'aide des auteurs de ces tests serait grandement utile a leur remise en etat! (3) les tests qui ne marchaient plus lors du bilan Apps de la version v2.2.02 fonctionnent a nouveau.


18/09/01 : Production de la version V2.2.05  

From Michel Gazaix.
Une nouvelle production elsA est disponible : ELSAVERSION = v2.2.05 ELSAPROD = sgi, sgi_dbg, sgi_r8, nec_r8 cvs rtag : I2205g Cette production apporte de nouvelles fonctionnalites importantes : .................................................................... 1. JP : CAIMEN C1 (cvs rtag : I2205a) (developpement decrit avec plus de precision dans un message de JP du 7 Septembre a elsA-dev) a - Reduction de la dissipation artificielle par produit du flux dissipatif par un facteur sense detecter la couche limite (cette reduction n'est faite que dans la/les direction(s) orthogonale(s) aux parois). JP a code cette fonctionnalite avec un facteur sur la vitesse. SP a code une formule sur la pression d'arret. Ces fonctionnalites ont ete testees avec succes par VB (As28G Nacelle Lara) b - D12.04 : Codage d'un flux visqueux avec gradient corrige aux interfaces (methode MSD Borrel). D20.01 : XJ a ecrit le complement raccord necessaire, VB a teste cette fonctionnalite. c - D04-06 : Separation de FxdTur en deux classes : FxdTurMeanFlow et FxdTurTurbVar --> remplacement de TurBase::compDiffFluxDensity par deux methodes: compDiffFluxDens_gradCen compDiffFluxDens_gradCenCorr 6 nouveaux cas test d'integration : rae2822_VRed.py rae2822_mfc.py rae2822_mftvc.py rae2822_tvc.py sqnz_ns_adia_mfc.py sqnz_ns_adia_paRed.py Nouvelles clefs python utilisateur A PREVOIR (et documenter dans le MU) .................................................................... 2. CB : D14.04 : Chimere points orphelins + Preconditionnement de la recherche de la cellule d'interpolation pour une grille cylindrique 3 nouveaux cas tests d'integration : rae-lam.py, rae-ko.py, rae-ko-orph.py (cvs rtag I2205b) .................................................................... 3. XJ : D10.01 : Developpement "super-frontiere": Adressage super-frontiere et gestion de clusters de processeurs. cvs rtag I2205c .................................................................... 4. MG : D03.02 : Vectorisation Flux Decentres cvs rtag I2205d Sur un cas 3D (256 X 49 X 33) multigrille (1 grille grossiere), flux de Roe, implicite relaxsca, limiteur van Albada, en multigrille (1 grille grossiere, nitercoarse = 2) on avoisine 1900 MFlops (2,3 micro-seconde par pt et par iter) en monogrille : 2 GFlops, (1.2 micro-seconde par pt et par it) Vous pouvez consulter les resultats multigrille et monogrille du profiler sur NEC . 15 nouveaux cas test d'integrations (10 shockTube, flux roe/vleer, limiteurs null/minmod/vanalbada/vanleer/superbe, 5 as28g : flux de roe avec els 5 limiteurs)) .................................................................... 5. AJ : D30.02 : Implicitation IRS terme source turbulent (Lhs) cvs rtag I2205f 4 nouveaux cas tests d'integration (plaque plane en KL, KO, KEps, SA) Plus de precisions dans la fiche de revue RI-NI-01-09-05. .................................................................... 6. MM : D25.05 : Patch grid (raccords non coincidents) - raccords sequentiels non coincident total - correction d'un bug sur les super-frontieres - optimisation pour les raccords non coincidents par ligne 2 nouveaux cas test d'integration (naca9 et tuyere 8 domaines) ............................................................................ Les extractions ont legerement changees (SB) ............................................................................ L'attribut 'mode' n'est plus pris en compte par le code, le mode (periodique ou non) des extractions se regle tout seul. D'autre part, Pour des extractions de type flux, avant on mettait : "w1+w2-w3-w4". Maintenant on peut mettre : "w1 w2 w3 w4", le sens des fenetres sera donne par le sens de la normale exterieure au fluide. Mais si on veut quand meme fixer le sens (par exemple en dehors d'une paroi, par exemple pour un debit), il faut faire attention a ne pas oublier le signe pour la premiere fenetre : "+w1+w2-w3-w4". .................................................................... Correction bugs principaux : .................................................................... CB : GeoCfdField::removeCfdField Bug logique : on ne rendait jamais la memoire utilisee. Impact : probablement assez faible. VB/MG : TurKO::computeWallForWallLaw bug de calcul d'addresse Attention aux calculs KO avec loi de paroi realises avec les anciennes productions ... AR : RhsJoinNearMatchFine attribut _isComputedByCoarse non initialise Ce bug a ete trouve apres BEAUCOUP d'efforts, en utilisant purify (sur SUN de DRIS). Peut expliquer les anomalies du cas farceur naca-LMG.py Cependant, il semble qu'il subsiste des anomalies. (message de Barrera du 18 Septembre a elsA-info) MG : Si on associe une window a un block invalide, ou si on cree une boundary associee a une window ou a un block invalide, on a mainteneant un message d'erreur.


14/09/01 : Production de la version V2.1-p3  

From Anne Jollès.
La production v2.1.15-p3 correspondant au patch V2.1-p3 est disponible sur NEC . Deux executables sont disponibles : elsA.x (sans profiler) et elsA_prof.x (avec profiler).


30/08/01 : Production de la version V2.2.04  

From Michel Gazaix.
Une nouvelle production elsA est disponible : ELSAVERSION = v2.2.04 ELSAPROD = sgi, sgi_dbg, nec cvs rtag : I2204e Cette production comprend : RH : D99.22.04 : Implicitation IRS du terme source en turbulent (Systeme Tur) RH : D25-03 : Modele K-omega AVTAC (debut) RH : D25-01 : Orthogonalistion des distances 4 nouveaux cas test d'integration : flatPlateKL/flatPlateKL_incli_gridline.py flatPlateKL/flatPlateKL_incli_gridline_ortho.py flatPlateKL/flatPlateKL_incli_mininterf.py flatPlateKL/flatPlateKL_incli_mininterf_ortho.py VB : D12.06 : Loi de paroi en multigrille 3 nouveaux cas de non regression : AS28G/as28g_wlmg.py flatPlateTurb/flatPlateKL_wall1.py flatPlateTurb/flatPlateKL_wall2.py AR : D25.04 (partiel) : raccord multietages avec eqs transport SB : D12.03 : Extractions RANS2 (lot 1) SB : D25.09 : Extractions RANS2 (lot 2) SB : D22.03 (partiel) : Extraction node pour mut


24/08/01 : Production de la version V2.1-p2  

From Anne Jollès.
La production v2.1.15-p2 correspondant au patch V2.1-p2 est disponible sur NEC . Deux executables sont disponibles : elsA.x (sans profiler) et elsA_prof.x (avec profiler).


23/08/01 : Production de la version V2.2.03  

From Michel Gazaix.
Une nouvelle production elsA est disponible: ELSAVERSION = 2.2.03 cvs rtag : I2203a Cette version n'est pas actuellement disponible sur NEC. Cette version n'apporte pas de nouvelles fonctionnalites; elle corrige plusieurs bugs, et introduit les modifs du patch 2.1.15-p2 dans la branche principale. Principaux bugs corriges: 1. AR : raccords non coincidents a lignes coincidentes (cf son message a elsa-user) 2. AJ : Bug Condition limite supersonique avec fichier (createBndSupInOut) Bug mis en evidence par A. Chatelain (CORIA). 3. MG : Bug dans l'implicite ADI. Se manifeste par une division par zero (Floating point Exception). Ce bug apparait lorqu'on utilise l'otion padding, sur les machines qui refusent de diviser par zero (avec les options presentes dans les makefiles actuels, c'est par exemple le cas sur DEC, et sur SGI en mode sgi_dbg). Ce bug n'affecte pas les resultats (c'est une division qui intervient sur les cellules fictives introduites par le padding). Donc pas depanique, sur NEC les resultats sont corrects ! 4. Bug Makefile pour fujitsu (complements apportes par Ph. Chanez, SNECMA) Pour la prise en compte du patch v2.1.5-p2 : 1. J'ai pris en compte le patch de Marc Montagnac (Pcm, FactTurb.C) J'ai, sur sa suggestion, introduit un "static_cast(nbCell)" sur un argument (entier) de PcmReduce::apply, utilise dans SioFilterResidual.C (je relaie son mail a elsa-dev pour plus de precisions) 2, Je crois avoir pris en compte les modifs de Marc Lazareff concernant l'interface (Kernel/api/Py) pour 2 fichiers : EpAttrDefs.py et EpAttrDict.py. En son absence, je n'ai pas pu faire plus, il faudra penser a verifier a son retour.


06/07/01 : Bilan Apps de la version V2.2.02  

From Anne Jollès.
Apps a ete mis a jour avec la version v2.2.02. cvs tag: I2202 Un listing de l'execution de la base en mode de production sgi est donne ci-joint. On remarque que les test: BndWallLaw1.C, TurMKFLC2_wl.C, TurWallLaw1.C, TurMichelTyp1B.C, TurMichelTyp2.C, TurMichel.C, TurTwoLayer.C, TurTransiLoc.C, TurTransiNoLoc1.C, TurTransiNoLoc3.C, TurTransiNoLoc4.C ne fonctionnent plus. Ils sont en cours de correction.


06/07/01 : Production de la version V2.2.02  

From Michel Gazaix.
Une nouvelle production elsA est disponible: ELSAVERSION=v2.2.02 cvs rtag = I2202e ELSAPRODUCTION: sgi sgi_dbg sgi_r8 sgi_r8_mpi linux nec_r8 (avec et sans profiler) Cette production comprend: 1- Marc Montagnac : Nouveaux complements Dtw (D22.05) 2. Xavier Juvigny : Introduction de la notion de groupe ("MPI_Comm") en parallele (D00.26.01) 3- Anne Jolles : Complements actuator-disk en centre et extension en visqueux et turbulent (D15.02) 4. MG : Une nouvelle gestion des valeurs par defauts il faut maintenant utiliser les methodes (static) FactDataBase::getDef[IFS](...) Les valeurs par defauts sont centralisees dans le module python EpKernelDefVal.py : l'utilisateur averti peut donc modifier dynamiquement ces valeurs, ou simplement les controler. 5. MG : Les fichiers FactDataBase.dat et FactCfdDataBase.dat sont maintenant inutiles, et remplaces par le module python EpKernelClassName.py. En python-API, on n'a plus besoin de ELSAHOME (ou ELSAROOT). 6. Corrections de bug : 1) AgCoord (E_FUNC::fEqual aulieu de fEqual) (MM) 2) Transi/TurCompLam2F.for : correction de signe (RH) 3) BndWbNRARturF.for : epso --> eps0(0:intb-1) (MM) 4) BndNoReflUnst.C (Julien Delbove) : 5 bOut(lint, eqInd) = _state(lint,eqInd); au lieu de wbOut(lint, eqInd) = _state[lint,eqInd]; 6) Correction Fact*.dat (bien sur reporte dans EpKernelClassName.py : peux-tu verifier, RH ?) 7) "covariant return" non supporte sur DEC, FUJITSU (en plus de CRAY) (Luis Barrera), dans BlkMotionBlock.[hC]. 8) ML : Contournement d'un bug (trouve par LBa) en python-utilisateur 9) ML : ajout de quelques elements manquants en python-utilisateur 10) Bug option compilation g++ sur linux (MG) 7. Introduction d'enum (cf DesNumerics.h): enum FrameFormulation { AbsoluteFrame, RelativeFrame }; enum Deform { Undeformable, Deformable }; enum VelocityFormulation { AbsoluteVel, RelativeVel }; enum Motion { Fixed, Mobile }; Ces enum sont deja utilises dans Fact; Il faudra prevoir une action pour les utiliser systematiquement. 8. Simplification / Rationalisation des nommages dans Fact 9. Remplacement de pas mal de cles definies dans DesDef1.h : Je rappelle que les clefs doivent commencer par le prefixe 'KEY_' (et non pas 'E_'). Ce sont les valeurs associees aux clefs qui commencent par 'E_'. Exemple : 'E_TRANSL_SPEED_X' devient 'KEY_TRANSL_SPEED_X'


22/06/01 : Production V2.2.01 / Branche M2115  

From Michel Gazaix.
Une nouvelle production est disponible : ELSAVERSION=v2.2.01 cvs rtag : I2201a Cette production integre le developpement "Multigrille Local" effectuee par Antoine Roure. (D99.40.2 ). Par ailleurs, afin de preparer la livraison de la version finale 2.1, une branche (au sens CVS) a ete creee. L'identificateur CVS de la branche est : M2115. A priori, la creation de cette branche n'a pas d'impact sur les developpeurs: il leur est toujours conseille de se synchroniser sur la production la plus recente du tronc commun.


29/06/01 : Bilan Apps de la version V2.1.15  

From Anne Jollès.
La mise a jour des tests de developpement Apps/test a ete effectuee sur la version de production v2.1.15. Le Repository de Apps a ete modifie: identificateur CVS: I2115 Le bilan des tests de Apps executes avec la production v2.1.15 dans le mode de production sgi est maintenant disponible. Quelques remarques : Les 2 tests LhsAdi_17.C et LhsAdi_18.C qui ne marchent pas seront supprimes lors de la prochaine mise a jour d'Apps avec la v2.2.01 (ceci etant du a l'integration du developpement Multigrille local de A. Roure dans la v2.2.01). Les 3 tests: DtwFldWall.C, DtwFldWall1.C, DtwFldWall2.C seront remis a jour lors de la prochaine integration de M.Montagnac .


29/06/01 : Production de la version V2.1.15  

From Anne Jollès.
La production v2.1.15 correspondant a la nouvelle livraison v2.1 decrite dans le Proces verbal de livraison /ELSA/PVL-01038 (pdf), est disponible sur NEC et SGI. ELSAVERSION=v2.1.15 cvs rtag : M2115-p1 (Cette production a ete faite a partir de la branche CVS cree pour cette livraison : M2115) Sur SGI, les modes de production disponibles sont les suivants: sgi, sgi_dbg, sgi_r8, sgi_mpi Sur NEC, 2 executables sont disponibles : elsA.x (sans profiler) et elsA_prof.x (avec profiler).


18/06/01 : Production de la version V2.1.14  

From Michel Gazaix.
Une nouvelle production elsA est disponible : cvs rtag : I2114e ELSAVERSION=v2.1.14 ELSAPROD=nec_r8, sgi, sgi_r8, sgi_mpi, sgi_r8_mpi Cette version n'apporte pas de fonctionalites supplementaires. Les principales corrections par rapport a la v2.1.13 sont: 1. AR : Correction d'un bug pour les raccords non coincidents a lignes coincidentes (ouf). Cas test concerne : turbine VEGA 2. L'introduction d'une clef, "vectplan_i" permettant de choisir de parcourir certaines boucles par ligne I, ou par plan J-K. Cette clef est "experimentale", et n'est pas decrite, volontairement, dans le Manuel Utilisateur. 3. VB : Une modification de la prise en compte des composantes du vecteur vitesse pour le calcul des flux diffusifs en axisymetrique. Cette modif ameliore le cas "tuyere MARCO". 4. Modifs pour pouvoir passer en mode "dbg" Beaucoup de ces modifs tournent autour de l'usage des methodes infernales GeoCfdField::getCfdField et GeoCfdField::createCfdField. L'interface de GeoCfdField est probablement a clarifier, avec peut-etre seulement une seule methode public pour acceder aux champs geres par GeoCfdField. 5. MG : Des modifications relatives au parallele: a) Pour les raccords coincidents, remplacement de MPI_BUFFER_Send par MPI_Sendrecv_replace; Ca semble eliminer des deadlocks sur Fujitsu et HP, et continuer a marcher sur SGI... Au cours du debogage. J 'ai du ecrire quelques cas tests unitaires pour cerner le probleme. Il faudra valider completement cette modif. La validation sur le gros cas fuselage est en cours sur la SGI d'EUROCOPTER. Pour ceux que ca interesse, Une justification du remplacement se trouve dans les commentaires de la methode PcmStreamOutImplMPI::initBuffer() b) La clef "node" n'est plus indispensable pour les objets DesMesh (on se sert de la clef du DesBlock auquel est "attache" l'objet mesh). Si la clef 'node' est absente, il faut eviter de faire un 'submit' explicite (donc, pas de "myMesh.submit()") c) Un script parallele, donc comportant des clefs "node", est maintenant utilisable sans modif par un executable non parallele. Ceci devrait simplifier la gestion des cas tests. 6. MM : Vectorisation Lussor 7. MM : Correction fuite memoire (en turbulent BL parallele)


07/06/01 : Bilan Apps de la version V2.1.13  

From Anne Jollès.
La mise a jour des tests de developpement Apps/test a ete effectuee sur la version de production v2.1.13. Le Repository de Apps a ete modifie: identificateur CVS: I2113 Le bilan des tests de Apps executes avec la production v2.1.13 dans le mode de production sgi est maintenant disponible. On constate que 2 cas qui marchaient avec la v2.1.10: LhsAdi_17.C et LhsAdi_18.C ne marchent plus. Ce probleme est en cours de traitement par l'auteur de ces tests. On constate egalement que les 3 tests: DtwFldWall.C, DtwFldWall1.C, DtwFldWall2.C ne marchent pas. Leur mise a jour sera publiee lors de la prochaine integration de M.Montagnac . Enfin, la mise a jour du test FldField.C a egalement ete promise par son auteur.


01/06/01 : Production de la version V2.1.13  

From Michel Gazaix.
Une nouvelle production est disponible : ELSAHOME = v2.1.13 cvs rtag = I2113b ELSAPROD : sgi, nec_r8 Cette production corrige des bugs mineurs. On a egaleemnt supprime pratiquement tous les "cast" non qualifies.


23/05/01 : Production de la version V2.1.12  

From Michel Gazaix.
Une nouvelle production est disponible: ELSAVERSION = v2.1.12 cvs rtag = I2112e Cette production est disponible: sur oneroa13 (sgi, sgi_r8, sgi_mpi, sgi_r8_mpi, sgi_dbg) sur cronos (nec_r8, nec_r8_mpi) sur oneroa36 (linux, linux_mpi, linux_r8) On a egalement verifie la production d'executable: sur HP (hp_r8, hp_r8_mpi) sur SUN (sun, sun_gcc) Une validation a ete effectuee sur cronos (A. Jolles) Cette production comprend: 1. D42.3 D3.2 D12.2 (S. Beaugrand) 2. D4.1, D5.2 (R. Houdeville) 3. Corrections de bug: 3.1 Erreur en dissipation matricielle (FxcMatNumDiss), sur NEC (specifique v2.1.11) (AJ) 3.2 Erreur algorithmique dans l'ordre des messages envoyes en parallele (Bug trouve par Marc Montagnac, suite a l'anomalie detectee par EUROCOPTER) 3.3 Bug Parallele dans l'appel de MPI_Buffer_detach (MG) 3.4 Bug en turbulent avec mouvement (bug trouve par AJ, specifique a la v2.1.11) 3.5 Bug ecrasement memoire dans createFxc par utilisation inconsideree de sprintf (trouve par MP) 3.6 Bug GeoConnect (JP) 3.7 Bug raccords non coincidents (XJ, AR): cependant, le cas caimen-08 (?) ne fonctionne toujours pas... 3.8 CFL variable en python utilisateur 3.9 Nombreux complements en python-utilisateur (ML, AJ)


25/04/01 : Production de la version V2.1.11  

From Michel Gazaix.
Une nouvelle production : ELSAVERSION = v2.1.11 est disponible. Sur oneroa13 : sgi, sgi_dbg, sgi_r8 Sur cronos : nec_r8 Sur oneroa36 : linux, linux_mpi, linux_so ======================================================== Cette production comprend : 1. JC Boniface (D10.6) ALE + decentre, ALE +dismat, pas de temps dual + maillage entraine,... 2. AJ et MG : Nettoyage / Simplification de GeoCfdField 3. G. Bastin (D13.1)Condition au limite "Full Cooling" et paroi isoth mouvement 4. XJ : Condition debit global (remis a jour) 5. VB : D4.10 Debit local 6. VB : D2.2 Modele bi-couche 7. Steve Champagneux (CERFACS) : Uniformisation des directives de compilation 8. ML D10.7 : Interface relax LU, dissipation, mouvement,... 9. MG / AJ : D13.5, D13.6 Optimisation CPU NEC Pour le cas Aile M6 60 blocs (sans mouvement), sur NEC, on passe de: 2.8 secondes par iter a 2.05 microsecondes L'executable NEC ne genere plus de fichier mon.out : on ne peut donc plus obtenir de profiler. J'ai fait ce changement pour gagner un peu de temps CPU. Si on veut faire un profiler, il faut utiliser elsA_prof.x (se trouve dans le meme directory que elsA.x)


23/04/01 : Bilan Apps de la version V2.1.10  

From Anne Jollès.
La mise a jour des tests de developpement Apps/test a ete effectuee sur la version de production v2.1.10. Le Repository de Apps a ete modifie: identificateur CVS: I2110a Le bilan des tests de Apps executes avec la production v2.1.10 dans le mode de production sgi est maintenant disponible. Cette mise a jour a demande un travail important a plusieurs d'entre nous, ce qui explique sa publication bien tardive. J'espere que la gene que cela a pu vous occasionner n'a pas ete trop importante. Attention: de facon a accelerer la mise a jour de Apps, je commence des maintenant a modifier a nouveau son Repository pour une mise a jour sur la version v2.1.11 en preparation.


12/03/01 : Production de la version V2.1.10  

From Michel Gazaix.
Une nouvelle production: cvs rtag : I2110a ELSAPROD=sgi (ELSAPROD, et non plus ELSAPRODUCTION !!!) est disponible sur SGI. Cette production comprend: ............................................ 1. Marc Poinot : D00.20.5 ............................................ Cette nouvelle structure de Makefile vise a plus de generalite. MP enverra un mail pour donner tous details utiles. Vous pouvez sans difficulte continuer temporairement a utiliser les anciens Makefile, si vous le souhaitez. ............................................ 2. Marc Montagnac: ............................................ 2.1 Diverses corrections pour le Lussor (implicitation correcte des C.L., acces en python-utilisateur) 2.2 Modif dans le calcul des rayons spectraux. Voici ses explications : Je rappelle rapidement la situation : Dans les flux de dissipation scalaires, le calcul du rayon spectral aux interfaces fait intervenir les memes methodes que ce soit pour les raccords ou pour des CL physiques. (FxcScaNumDiss::computeField, case SPECRADIUS_INT :) or getCellValues renvoie une valeur au centre d'une cellule pour les raccords alors que compBoundaryValues renvoie une valeur a l'interface pour une CL. J'ai donc differencie les traitements entre raccords et CL dans cette partie. fichiers touches : Eos/Base/EosIdealGas.C Eos/Base/EosIdealGas.h Eos/Base/EosIdealGasP.h Fxc/Centered/FxcScaNumDiss.C Fxc/Centered/FxcScaNumDiss.h Fxc/Centered/FxcScaNumDissGridMotion.h Fxc/Centered/FxcScaNumDissGridMotion.C fichiers nouveaux : Eos/Base/EosRadiusIntBorBndF.for Eos/Base/EosRadiusIntBorBndF.h Eos/Base/EosRadiusIntBorBndMotionOrAleF.for Eos/Base/EosRadiusIntBorBndMotionAndAleF.for Tous ces developpements peuvent se resoudre plus simplement si les conditions limites ont la possibilite de renvoyer des valeurs dans une cellule fictive mais il n'est sans doute pas preferable de melanger les genres. Quelques remarques et details annexes au fil du developpement : a. Dans FxcScaNumDissGridMotion::compRadiusIntBor, une branche conditionnelle peut etre supprimee. (mais peut-etre est-elle presente pour bien differencier les aspects vitesses entrainement et deformation. J'ai donc repris le meme modele) b. EosRadiusIntBorMotionF.for et EosRadiusIntBorAleF.for se ressemblent fortement. Il ne change que par le nom de la vitesse. (mais peut-etre par souci de differenciation). Pour ma part, j'ai factorise mon developpement. c. Les portions de codes utilisant une moyenne arithmetique n'ont pas besoin de la notion de sens et peuvent certainement se simplifier. Pour finir, je precise que ce developpement n'a pas ete teste ni valide et que ces deux etapes risquent fort de prendre un certain temps vu l'impact numerique.


12/03/01 : Production de la version V2.1.09  

From Michel Gazaix.
Une nouvelle production du logiciel elsA est disponible: cvs rtag : I2109k ELSAVERSION=v2.1.09 Sur oneroa13: ELSAPRODUCTION=sgi_cc_64 Sur NEC: nec_sxcc8 Cette production contient: ........................................................ 1. Sebastien Beaugrand ........................................................ 1.1 : D48.2 : Epaisseurs integrales extractions de "bl_quantities_2d" et "bl_quantities_3d" (loc=4) 2d correspond a : DELTA1, THETA11, HI, DELTA, RUNIT 3d correspond a : DELTA2, THETA12, THETA22, BETA0 1.2 : D4.5 : Chimere extraction fichier loc=CELL de la nature de la cellule: nature de cellule : - 1 : classique (i.e calculee par le schema, par CL autres que Chimere) ; - 0 : a interpoler (au sens Chimere) ; - -1 : masquee. 1.3 : Verification qu'on est bien en turbulent lors d'extractions de ce type. Sinon message (avant ca plantait cf. PR-2001-01-21). 1.4 : Titre dans un fichier d' "ExtractGroup" (PR-2001-01-20) ........................................................ ........................................................ 2. Sylvie Plot : ........................................................ D6.4 : Correction de bugs (trouves par BC et JCB) concernant le calcul de la vitesse d'entrainement du repere relatif. Attention au changement de signe pour les valeurs asociees aux clefs : transl_speed_[xyz] (Ex: script naca64_rigidMotion.py) ........................................................ 3. Jacques Peter : ........................................................ Ce developpement est decrit en detail dans la fiche de revue RI-NI-01-01-03. 3.1: D6.5 : transformation LDU --> LURelax, et linearisation des dissipations d'ordre 4 (LuRelax scalaire --> D4 scalaire LuRelax matriciel --> D4 scalaire ou matriciel ) + correction d'un bug dans LhsAdiVL7 + renommage .. 3.2 : D00.14.2 : vectorisation CL paroi Euler avec extrapolation multidimensionnelle de pression 3.3.1 : mise a disposition des trois versions de la dissipation d'ordre 4, a- rspec d3 b- d1 rspec d2, c- d2 rspec d1 dans la classe FxcScaNumDiss, et suppression de la classe FxcNumDiss. Suite a ce developpement la dissipation champ moyen devient une fonctionnalite a quatre variantes independantes (hors parametre scalaire) - la formule du senseur - la formule du d3 - l'utilisation de la correction de Martinelli - l'utilisation de la reduction de dissipation 3.3.2 : diminution de la dissipation la ou la vitesse est faible; 2.3.3 : supression du fichier fortran "SGI" FxcScaNumDissDirF.for Pour rendre ce developpement accessible en Python-utilisateur, Anne Jolles et Marc Lazareff ont modifie le fichier EpAttrDefs.py. L'utilisation de 'disnum' est donc aujourd'hui obsolete. L'ancienne ligne (Python utilisateur): << num.set('artviscosity','disnum') >> est a remplacer par les 2 lignes: << num.set('artviscosity','dissca') num.set('av_type','d1_rspec_d2') >> L'utilisation de l'option de Barrera devient aujourd'hui: << num.set('artviscosity','dissca') num.set('av_type','d2_rspec_d1') >> Pour resumer, la cle 'av_type' peut prendre 3 valeurs: 'rspec_d3' (ou 0 en API, valeur par defaut), 'd1_rspec_d2' (ou 1 en API), 'd2_rspec_d1' (ou 2 en API). Enfin les nouvelles cles rendant compte de la fonctionnalite "lowering of dissipation " sont les cles: av_lower, av_lower_exp, av_lower_velinf et s'utilisent de la facon suivante: (je ne sais pas si c'est un exemple intelligent) num.set('artviscosity','dissca') num.set('av_type','rspec_d3') num.set('av_lower','lower_i') num.set('av_lower_exp',2) num.set('av_lower_velinf',1.) En tout etat de cause : une validation complete est inevitable, avec bien sur introduction de nouveaux cas tests. (cvs rtag I2109a) ........................................................ 4- Anne Jolles et Vincent Gleize ........................................................ D23.3 : Complement Basse Vitesse Ce developpement est decrit en detail dans la fiche de revue RI-NI-01-01-02. (cvs rtag I2109b) ........................................................ ........................................................ 5. Christophe Benoit et Gaelle Jeanfaivre ........................................................ D43.2 : Chimere instationnaire Partie 2 (cvs rtag I2109c) ........................................................ ........................................................ 6. Robert Houdeville ........................................................ D38.1 : Initialisation des modele de turbulence a partir de: 6.1- mut = coeffmutinit * mu : comme avant, sans regression ------------------------------------------------------------- 6.2- meme chose, mais a l'iteration "n" --------------------------------------- Evite de faire un premier calcul en laminaire, ou pseudo laminaire, c'est-a-dire avec un modele de turbulence qui ne developpe pas la turbulence. En .py, il faut: timeInteg.setI("iter_ini_tur", 10) 6.3- mut donne par fichier et hypothese de bradshaw --------------------------------------------------- Permet de passer d'un modele algebrique a un modele a equations de transport: ini1_blk1 = DesWindow("ini1_blk1","blk1") ini1_blk1.setS("all","yes") ini1 = DesInit('ini1','ini1_blk1') # initialisation champ aero ini1.setS("file",INIT_FILE_NAME) ini1.setS("format",INIT_FORMAT) ini1.setS("var",INIT_VARIABLES) # FROM_MU=1, FROM_MUT=2, FROM_KL=3, FROM_KO=4, FROM_SA=5 ini1.setI("initialmodel",2) # fichier de mut ini2 = DesInit("ini2","ini1_blk1") ini2.setS("file","mutKL_45x50.v3d") ini2.setS("format","fmt_v3d") ini2.setS("var","viscturb") 6.4- initialisation d'un modele par un autre -------------------------------------------- a- modeles k-eps ou k-omega a partir de k-l mut connu (fichier) : pas d'hypothese distance connue (fichier) : pas d'hypothese ni mut, ni distance : hypothese fmu = 1 (effets bas reynolds neglige pour calcul de mut depuis k et l. b- modele k-eps a partir de k-omega pas d'hypothese c- modele Spalart Allmaras mut connu (fichier) : pas d'hypothese k-omega : pas d'hypothese k-l : non prevu, passer par fichier mut (on peut faire sans mais avec hypothese si demande) On peut toujours imaginer d'autres combinaisons. On verra a l'usage de l'interet. 6.5- Initialisation a partir d'une solution sur un maillage grossier -------------------------------------------------------------------- Permet de repartir d'une solution obtenue un point sur 2 dans 1,2 ou 3 directions: FactSio::fillSolFromCoarse Ce n'est pas finalise, mais ca marche. Il faudra modifier pour eviter de reecrire la convention de numerotation des cellules par rapport a i,j,k. Le but est d'avoir rapidement un champ approche pour mieux initialise un maillage fin. Pour le RAE822, je suis parti d'un maillage 129x33. Apres 3000 ite, c'est converge en 10'. La reprise en maillage fin donne a peu pres les memes residus que le calcul des 3000 ite maillage fin, 4 fois plus long. Pour plus d'info: ELSA/NI-00079 NE PAS OUBLIER : Plusieurs cas de validation seront a prendre en compte. 6.6- Dans FactSio::initSolutionFromDesc --------------------------------------- Le calcul de grad(V) pour Bradshaw qui utilisait un prolongement du champ de vitesse sur toutes les frontieres est maintenant fait en utilisant les parois pour la condition de non glissement. Ca devrait accelerer un peu le developpement de la turbulence (pour toutes les autres frontieres, on continue a prolonger la vitesse). Cette modif provoque une regression pour tous les cas qui utilisent "coeffmutinit". 6.7- Dans Tur -------------- Avec Anne Jolles (cf point numero 9 ci-dessous), J'ai commence a enlever quelques champs inutiles (sqrt(k) et son gradient). (cvs rtag I2109d) ........................................................ ........................................................ 7. Michel Gazaix ........................................................ 7.1 : Introduction de nouveaux operateurs dans FldFieldF. Correlativement, suppression de certaines methodes (add, sub, setNewArray, multFract, multNum, multDen,...) Une description / justification est en fichier attache. 7.2 : Introduction (experimentale ...) de deux namespace: E_FUNC: contient les fonctions utilitaires definies dans DefFunction.h E_FACT: contient les celebres createXXX, qui n'appartiennent a aucune classe. 7.3 : Remplacement de certains EosBase en EosIdealGas: on evite ainsi pas mal de downcast disgracieux, et surtout completement inutiles aujourd'hui, puisqu'elsA ne sait pas encore traiter autre chose que des gaz parfaits. 7.4 : Deplacement de l'attribut (protected) _grid de BndBase vers un attribut private dans les classes Bnd concretes. 7.5 : Une fonction fEqual(a,b) est maintenant disponible: inline E_Boolean fEqual(E_Float lhs, E_Float rhs, E_Float precision = E_CUTOFF ) { return (E_abs(lhs - rhs) < precision) ? E_True : E_False; } Cette fonction devrait etre utilisee systematiquement lorsqu'on doit comparer des nombres reels (et donc en particulier dans les tests unitaires). On ne devrait donc plus voir une orgie de "E_CUTOFF", E_TOL, "precision", "epsilon", ... dissemines dans tout le code. (cvs rtag I2109e) ........................................................ ........................................................ 8. Bernard Cantaloube ........................................................ D6.9 : ALE pale souple Attention : la clef "gcl_law" devient "vol_comp_meth". Penser a introduire un (ou plusieurs) cas de validation. (cvs rtag I2109g) ........................................................ ........................................................ 9. Anne Jolles et Robert Houdeville ........................................................ D10.2 : Encombrement memoire : Tur et Fxd Economie de plusieurs champs inutiles dans GeoCfdField. D'autre part, le choix a ete fait de ne plus stocker la vitesse (identificateur dans GeoCfdField: VELO_CELL) ainsi que la temperature (TEMP_CELL). Ceci a evidemment pour consequence des recalculs supplementaires. Si une validation complete montre une degradation jugee trop grande des performances, il faudra revenir sur ce choix. (cvs rtag I2109h) ........................................................ ........................................................ 10. Anne Jolles ........................................................ D6.8 : Impression en cas de pression ou temperature negative Au cours du calcul, des impressions ont ete rajoutees, informant l'utilisateur des que le calcul de la pression ou celui de la temperature donne une valeur negative. Les informations suivantes permettant de localiser dans le temps et dans l'espace ces valeurs negatives, sont alors donnees: - le niveau de grille, - le numero du bloc, - le numero du noeud (en parallele), - le numero de l'iteration, - le numero du sous-pas, - le numero du systeme, - la fonctionnalite correspondant a la premiere observation d'une pression ou temperature negative au cours du calcul, - la liste des cellules ou interfaces (liste limitee a 10 elements) ou la pression ou la temperature calculee est negative avec leur numero global ainsi que leur position sur le maillage. Ce developpement est pour l'instant "experimental". Il pourra par la suite etre modifie apres retour ou reaction des utilisateurs/developpeurs. (cvs rtag I2109j) ........................................................ Je signale que la production v2.1.09 : 1- corrige un bug relatif au choix du senseur en python-utilisateur (developpement Girodroux) 2- doit permettre d'ecrire des scripts "multi-etage" en python-utilisateur 3- doit fonctionner en MUSE-elsA (pour les fonctionnalites disponibles dans ce langage, biensur) Les points 2 et 3 n'ont ete que partiellement testes.


02/02/01 : Production de la version V2.1.08  

From Michel Gazaix.
Une nouvelle production : v2.1.08 est disponible. cvs rtag : I2108c Sur oneroa13: ELSAPRODUCTION=sgi_cc_64 ELSAPRODUCTION=sgi_cc_64_i8r8 ELSAPRODUCTION=sgi_cc_64_DEBUG ELSAPRODUCTION=sgi_cc_64_double_MPI Sur NEC: nec_sxcc8 Sur oneroa36: linux_egcs Cette production comprend: 1. Modele algebrique generalisee (Marc Montagnac, CERFACS). 2. Complements LUSSOR (Marc Montagnac) L'integration des developpemnts 1 et 2 correspond au tag CVS I2108a. 3. Preconditionnement basse vitesse (A. Jolles, V. Gleize, O. Rouzaud) (tag cvs : I2108b). 4. Simplification des classes FldField (M. Gazaix) Les classes FldInt[FI], FldCell[FI], FldNode[FI] sont supprimees. On aura donc moins de code a maintenir, ou a assimiler pour les nouveaux developpeurs. On supprime egalement tous les "dynamic_cast, ce qui allege le code. Pour les remplacer, on a definit des "typedef" (dans FldField.h): typedef FldFieldF FldNodeF; typedef FldFieldF FldCellF; typedef FldFieldF FldIntF; Il est fortement conseille de continuer a choisir le type adapte lorsqu'on construit un objet Field --> meilleure comprehension, documentation "gratuite". Attention: la suppression des classes specialisees (Node, Cell, Int) impose l'initialisation correcte de l'attribut _geoLocation, lorsque celui-ci est utilise (periodicite, par exemple). --------------------------------------------------- Remarques: 1. La perte d'efficacite CPU sur NEC specifique a la version 2.1.07 (resultant de l'option "-ftrace" passee au compilateur sxc++) a ete corrigee. 2. Toujours sur NEC, la version du compilateur sxc++ a changee. Je n'ai pas trouve d'anomalies pour l'instant... 3. J'ai re-mesure le temps par point et par iteration pour le cas "Rotor-7A" (multigrille, 2 grilles grossieres, modele de michel, moyenne arithmetique, avec padding) On obtient 4.81 micro-secondes par point et par iteration. 4. J'ai passe pas mal de temps a corriger des bugs (que j'avais introduit moi-meme) dans le traitement des conditions aux limites Multi-Etages. Ce qui a rendu la recherche de bug tres lente, c'est qu'on ne dispose pas (sauf erreur de ma part) de cas test de taille moyenne correspondant a la fonctionnalite multi-etage. Il faut environ 15 minutes pour faire 10 iterations, Heureusement que les autres cas de la base de non regression - environ une centaine - sont plus courts ! Nous sommes TRES demandeur d'un (ou de plusieurs) cas test de non-regression. Avis aux amateurs!!!


01/02/01 : Bilan Apps de la version V2.1.07  

From Anne Jollès.
La mise a jour des tests de developpement Apps a ete effectuee (production v2.1.07). Le Repository de Apps a ete modifie: identificateur CVS: I2107b Le bilan des tests de Apps executes avec la production v2.1.07 dans le mode de production sgi_cc_64 est maintenant disponible.


26/01/01 : Production de la version V2.1.07  

From Michel Gazaix.
Une nouvelle production : v2.1.07 est disponible. Sur oneroa13: ELSAPRODUCTION=sgi_cc_64 Sur NEC: nec_sxcc8 Sur oneroa36: linux_egcs Cette production comprend: 1. Complements modeles de turbulence (R. Houdeville) : Enrichissement du modele k-omega: Modeles k-omega Menter BSL, k-eps grand Reynolds avec possibilite de correction SST, finalisation criteres de transition locaux. Ce developpement est decrit en detail dans la fiche de revue RI-NI-01-01-01. 2. Complements vectorisation transferts multigrilles (J. Peter). 3. Complements Mouvement (JC Boniface). 4. Complements Actuator Disk, correction bug condition au limite Froude en multiblock, coefficients de forces et moments integres sur des sections de pales (S. Beaugrand). 5. Correction de bugs (introduits apres la derniere livraison) relatifs aux raccords multi-etages (A. Roure, M. Gazaix). 6. Corrections et modifications de l'interface (M. Lazareff) : - le limiteur de Zheng a pour cle' zhenglim (et non plus limzhen) ; - il est de nouveau possible de creer des structures de controle en Python-utilisateur, mais pour l'instant les constructeurs doivent dans ce contexte renseigner l'argument optionnel 'name', ex : i=3 if i == 2: print "doing nothing" elif i ==3: m1=model(name='m1') else: pass - les noms de fichiers de script peuvent contenir des points (ex: eu8dom.log.epy) ; - les variantes de macro-attributs utilisent un codage de leur longueur sur 2 chiffres (utile pour les harmoniques de la classe function).


11/01/01 : Production de la version V2.0-p2  

From Anne Jollès.
Pour la version patchee v2.0-p2 incluant les 2 patchs p1 et p2 apportes a la version v2.0, les productions suivantes sont disponibles : sgi_cc_64, sgi_cc_64_MPI, sgi_cc_64_i8r8, nec_sxcc8. La version v2.0-p2 a ete taguee "M2017_P3". Tests effectues: --------------- (1) Les tests de developpement (Apps) ont ete passes avec succes. (2) La non-regression (ElsaNoReg) egalement. (3) Sur SGI: quelques tests ont ete passes sur quelques iterations pour verifier les corrections du patch. Le cas de validation parallele SQNZ-MDPAR-01 a ete passe. (4) Sur NEC: les cas ROTOR37-MICHEL et VEGA-2HOH ont ete passes. La base de validation n'a pas ete executee.


22/12/00 : Production de la version V2.1.06  

From Michel Gazaix.
La derniere production du millenaire : v2.1.06 est disponible. cvs rtag : I2106d ELSAPRODUCTION: sgi_cc_64, nec_sxcc8 Cette production comprend: 1. La possibilite d'utiliser un senseur pression-entropie (au lieu de pression-vitesse) pour la viscosite artificielle (Ph. Girodroux-Lavigne /DDSS). Ce senseur semble indispensable pour calculer des profiles "aigues" (?), type PFSU (cette conclusion provient d'une comparaison avec des calculs COLIBRI effectues par DDSS). Ce developpement est accessible en Python-Utilisateur, pas en MUSE. 2. Des complements relatifs a l'actuator disk (J. Brezillon et S. Beaugrand). 3. L'introduction du modele ASM (V. Brunet), ainsi qu'un renommage dans les systemes Tur et Fxd (TurBoussinesq --> TurBase, FxdTurBoussionesq --> FxdTur)). Bonnes fetes a tous


15/12/00 : Production de la version V2.1.05  

From Michel Gazaix.
Une nouvelle production : v2.1.05 est disponible. cvs rtag : I2105e ELSAPRODUCTION: sgi_cc_64, sgi_cc_64_i8r8 sgi_cc_64_DEBUG, nec_sxcc8 Cette production comprend: 1. Possibilite de reduire le pas de temps pres des bords (Action D31.2) : A. Roure 2. Limiteur du mut variable en fonction des iterations (Action D45.1) : A. Roure 3. CFL par bloc (Action D45.2): A. Roure 4. Une optimisation CPU, consistant a calculer et a stocker la numerotation des interfaces associees a une condition au limite (au lieu de recalculer le tableau d'indirection a chaque iteration) : M. Gazaix. 5. Une action de "qualite logicielle", portant sur les "cast", comportant 2 volets a) De nombreux cast inutiles ont ete supprimes Pour cela, on a ete amene a corriger (ajouter ou enlever) pas mal de "const". b) Parmi les cast restant, un grand nombre ont ete precises: l'idee est de se forcer a ecrire precisement pourquoi on a besoin de caster. Cela permet de clarifier, et de retrouver plus facilement les occurences de cast. Pour cela on a utilise les 4 types de cast du C++ (static_cast, dynamic_cast, const_cast, reinterpret_cast) Pour le cas particulier du "dynamic_cast", on a introduit la macro E_DYNAMIC_CAST, afin d'eviter, si on le souhaite, tout surcout eventuel (temps CPU et memoire) introduit par la verification du type a l'execution. Dans la version v2.0.17, il y a (approximativement) 570 casts. Dans la nouvelle version, il subsiste 183 cast non qualifies (il faudra finir le menage ...). La presence de 232 (!) dynamic_cast suggere qu'on a aussi a faire une re-conception de certaines classes, ou en tout cas des conditions dans lesquelles nous les employons. 6. Une action de qualite logicielle dans le systeme Bnd (simplification de signatures, suppression de l'argument sizewb1 pour certaines routines Fortran).


15/12/00 : Mode de production pour V2.0-p1  

From Marie-Claire Le Pape.
Pour le patch v2.0-p1, decrit dans la page News/Project et deja mentionne precedemment, une nouvelle production est disponible : sgi_cc_64_double_MPI


08/12/00 : Production de la version V2.1.04  

From Michel Gazaix.
Une nouvelle production: v2.1.04 est disponible ELSAPRODUCTION: sgi_cc_64, nec_sxcc9 cvs rtag : I2104b Cette production comprend: 1. Codage "parabolique" de la correction Harten (cf Manuel theorique): JC Boniface Les cas test utilisant les flux de Roe sont a re-valider. 2. Complement de conditions aux limites (decentre, mouvement,...) Des cas tests specifiques sont a prendre encompte dans la base de Validation (cylindre, as28g) Pour 1 et 2, les tests Apps seront a verifier (JCB s'en occupe). 3. Codage Gear : J Peter 4. Amelioration / factorisation LU : J Peter Simplification des tests Lhs dans Apps


05/12/00 : Production de la version V2.1.03  

From Michel Gazaix.
Une nouvelle production: v2.1.03 est disponible. ELSAPRODUCTION=sgi_cc_64, nec_sxcc8 cvs rtag : I2103 Cette production contient les modifs necessaires au portage d'elsA sur IBM.


04/12/00 : Production de la version V2.0-p1  

From Marie-Claire Le Pape.
Pour le patch v2.0-p1, decrit dans la page News/Project, les productions suivantes sont disponibles : sgi_cc_64 sgi_cc_64_i8r8 sgi_cc_64_MPI nec_sxcc8 La version v2.0-p1 a ete taguee "M2017_P2". La base de validation n'a pas ete executee. Pour verification de la version patchee, les test-unitaires et la non-regression ont ete executes avec succes (par comparaison avec les versions v2.0.16 et v2.0.17) ; quelques cas-tests ont ete executes sur quelques iterations pour verifier certaines corrections de bugs. Le cas-test SQNZ-MDPAR-01 a ete execute sur SGI (cas-test parallele) avec succes.


01/12/00 : Production de la version V2.1.02  

From Michel Gazaix.
Une nouvelle production est disponible: ELSAVERSION = v2.1.02 ELSAPRODUCTION: sgi_cc_64, sgi_cc_64_DEBUG, sgi_cc_64_i8r8 nec_sxcc8 cvs rtag: I2102 Cette production comprend: 1. Prise en compte de criteres de transition (RH) 2. Correction de bugs (en gros les memes que ceux corriges dans le patch v2-P1)


23/11/00 : Production de la version V2.1.01  

From Michel Gazaix.
Une nouvelle production est disponible: ELSADIST=v2.1.01 ELSAPRODUCTION=sgi_cc_64, nec_sxcc8 Identificateur CVS: I2101b Cette production n'apporte que tres peu de modifications. La "raison d'etre" de cette production reside esentiellement dans le portage d'elsA sur HP (portage realise en collaboration avec P. Girodroux / DDSS). Le compilateur C++ de HP a permis de detecter quelques bugs ou maladresses, On espere par la meme occasion que le portage sur le nouveau serveur de calcul IBM, acquis par SNECMA, sera facilite. Cette production N'EST PAS destinee a corriger les "imperfections" de la derniere version livree (v2.0): une version "patchee" specifique est en cours de realisation, et sera mise a disposition des utilisateurs prochainement.


20/11/00 : Mode de production pour V2.0.17  

From Anne Jollès.
Pour la production v2.0.17, decrite dans le Proces verbal de livraison /ELSA/PVL-00071 (pdf) et deja mentionnee precedemment, la production dec_cxx_optimize a ete reinstallee, une production de python 1.5 ayant ete refaite sur oneroa12. ELSAROOT=/data/Export/v2.0.17


15/11/00 : Bilan Apps de la version V2.0.17  

From Anne Jollès.
Le Repository de Apps a ete modifie. L'identification CVS est: I2017 Les bilans des tests de Apps executes avec la production v2.0.17 sont disponibles dans les modes de production : sgi_cc_64et sgi_cc_64_DEBUG. On remarque des resultats differents pour certains tests entre ces 2 modes de production. En particulier: 1) OperGradMultiDom.C marche en mode sgi_cc_64_DEBUG mais pas en mode sgi_cc_64 2) BlkOverlappingData.C , LhsAdiVL71_3D, LhsAdi21.C, LhsAdi22.C, LhsAdi23.C, LhsAdi24.C, LhsAdi25.C, LhsAdi26.C, LhsLussor2D2Sca1.C, LhsLussor2D7Sca1.C, LhsLussor3D2Sca1.C marchent eux en mode sgi_cc_64 mais pas en mode sgi_cc_64_DEBUG.


15/11/00 : Production de la version V2.0.17  

From Michel Gazaix.
Une nouvelle production v2.0.17 est disponible. Cette production correspond a la version elsA 2.0, qui a ete livree le 14 Novembre 2000 et qui est decrite dans le Proces verbal de livraison /ELSA/PVL-00071 (pdf). L'identification CVS est: I2017 Les librairies sont disponibles pour les modes de production: sgi_cc_64 sgi_cc_64_DEBUG sgi_cc_64_i8r8 sgi_cc_64_double_MPI nec_sxcc8


27/10/00 : Production de la version V2.0.14  

From Michel Gazaix.
Une nouvelle production v2.0.14 est disponible. (cvs rtag: I2014c) Cette version est disponible en: sgi_cc_64 sgi_cc_64_DEBUG nec_sxcc8 Cette version correspond "pratiquement" a la prochaine livraison.


10/10/00 : Bilan Apps de la version V2.0.11  

From Michel Gazaix.
ATTENTION =========== Le Repository de Apps a ete modifie. (cvs rtag: I2011a) (Penser a faire un "cvs update") Introduction ------------ Ces modifications ont permis de remettre a jour plusieurs tests qui: -soit ne marchaient plus depuis la v2.0.09 - soit "semblaient " marcher lorsque l'on utilisait la commande "make testall" //' mais en realite ne marchaient pas lorsqu'on les lancait "a la main" par exemple par la //' commande: Obj/Src/Test.x Contributions de VB, MG, AR, XJ, AJ. La correction de "make testall" a ete effectuee par MP. Il reste a l'heure actuelle encore plusieurs tests qui //' ne marchent pas, dont certains commentes dans les Makefile. Il est evidemment important de ne PAS commenter dans un Makefile un test qui ne marcherait plus, notre but etant au contraire d'arriver a faire marcher //' TOUS les tests presents dans les Makefile, y compris ceux qui sont aujourd'hui encore commentes car //' necessitant une reconsideration relativement importante. Ce travail de remise en etat et augmentation du nombre des cas tests de developpement etant reconnu par tous comme indispensable dans une finalite de non-regression du source, j'espere que vous ferez bon acceuil a de //' nouvelles sollicitations ! Lancement des tests de Apps (provisoire) ---------------------------------------- Jusqu'a l' annonce de la prochaine production, une installation particuliere de cfg est a effectuer dans Apps afin de beneficier de la correction de MP et de pouvoir s'"appuyer" sur la commande: //' "make testall". Celle-ci est la suivante: cd .../DevWork/Apps/test 1- detruire le lien cfg s'il existe //' sinon ne PAS effectuer la commande "source.me" creant ce lien 2- faire une extraction de cfg par: cvs -d$CVSROOTKER co -P cfg 3- Ne surtout PAS faire de "cvs add" du repertoire cfg. Bilan des tests de Apps ----------------------- Modes de production : sgi_cc_64 et sgi_cc_64_DEBUG version: v2.0.11 Systeme Bnd =========== Test BndInj2Sub.C ---> FAILED ----------------------------- *** Test > ***Run : Obj/sgi_cc_64_DEBUG/Src/BndInj2Sub.x *** Test > : -- ################################# FAILED Process 7026785 (BndInj2Sub.x) stopped on signal SIGFPE: Floating point exception (default) at [BndInj2Sub::run(void):119 +0x60,0x10011bd0] 119 if (fabs(wb1(i,eqn)-vcons[eqn-1])>tol) Systeme Dtw =========== Test DtwThick.C ---> FAILED ---------------------------- *** Test > : *** Test > ***Run : Obj/sgi_cc_64_DEBUG/Src/DtwThick.x *** Test > : -- ################################# FAILED !!!DtwThick::computeOneWall: gradVelo must be calculated to compute the boundary layer thicknesses!!! Assertion failed: _data != E_NULLPTR, file /v2.0.11/elsA/inc/Fld/Base/FldArray.h, line 233, pid v26575 Abort(coredump) Systeme Lhs : (14 tests ne fonctionnent pas) =========== Test LhsAdiVL71_3D.C ---> FAILED -------------------------------- *** Test > : *** Test > ***Run : Obj/sgi_cc_64_DEBUG/Src/LhsAdiVL71_3D.x *** Test > : -- ################################# FAILED Process 6900502 (LhsAdiVL71_3D.x) stopped on signal SIGFPE: Floating point exception (default) at [coeffivl:169 +0x10,0x400047f588] Test LhsAdiVL72_3D.C ---> FAILED -------------------------------- *** Test > : *** Test > ***Run : Obj/sgi_cc_64_DEBUG/Src/LhsAdiVL72_3D.x *** Test > : -- ################################# FAILED probleme lhs. indice ij k ch 1 1 1 1 dq calcule -11.3126 dq prevu -6.16 ...................... ...................... probleme lhs. indice ij k ch 4 3 2 7 dq calcule -17.5594 dq prevu -17.56 nombre d'erreur(s) phase implicite = 167 //' Test LhsAdiVL73_3D.C ---> FAILED -------------------------------- *** Test > : *** Test > ***Run : Obj/sgi_cc_64_DEBUG/Src/LhsAdiVL73_3D.x *** Test > : -- ################################# FAILED Process 6820041 (LhsAdiVL73_3D.x) stopped on signal SIGFPE: Floating point exception (default) at [EosIdealGas::EosIdealGas (float,float,float,float,float,float,int):164+0x8,0x40000d4b14] Test LhsAdiVL74_3D.C ---> FAILED -------------------------------- *** Test > : *** Test > ***Run : Obj/sgi_cc_64_DEBUG/Src/LhsAdiVL74_3D.x *** Test > : -- ################################# FAILED probleme lhs. indice ij k ch 1 1 1 1 dq calcule -11.2407 dq prevu -6.16 .............. .............. probleme lhs. indice ij k ch 4 3 2 7 dq calcule -17.5594 dq prevu -17.56 nombre d'erreur(s) phase implicite = 166 //' Test LhsAdi20.C ---> FAILED ---------------------------- *** Test > : *** Test > ***Run : Obj/sgi_cc_64_DEBUG/Src/LhsAdi20.x *** Test > : -- ################################# FAILED Process 4550210 (LhsAdi20.x) stopped on signal SIGFPE: Floating point exception (default) at [solve2i:116 +0xbc,0x40004c53b8] Test LhsAdi21.C ---> FAILED ---------------------------- *** Test > : *** Test > ***Run : Obj/sgi_cc_64_DEBUG/Src/LhsAdi21.x *** Test > : -- ################################# FAILED Process 6581184 (LhsAdi21.x) stopped on signal SIGFPE: Floating point exception (default) at [solve5i:121 +0xbc,0x40004ca628] Test LhsAdi22.C ---> FAILED ---------------------------- *** Test > : *** Test > ***Run : Obj/sgi_cc_64_DEBUG/Src/LhsAdi22.x *** Test > : -- ################################# FAILED Process 7084777 (LhsAdi22.x) stopped on signal SIGFPE: Floating point exception (default) at [solve2i:72 +0xbc,0x40004c4708] Test LhsAdi23.C ---> FAILED ---------------------------- *** Test > : *** Test > ***Run : Obj/sgi_cc_64_DEBUG/Src/LhsAdi23.x *** Test > : -- ################################# FAILED Process 6614471 (LhsAdi23.x) stopped on signal SIGFPE: Floating point exception (default) at [solve5i:68 +0xbc,0x40004c8c5c] Test LhsAdi24.C ---> FAILED ---------------------------- *** Test > : *** Test > ***Run : Obj/sgi_cc_64_DEBUG/Src/LhsAdi24.x *** Test > : -- ################################# FAILED Process 7077849 (LhsAdi24.x) stopped on signal SIGFPE: Floating point exception (default) at [solve5i:68 +0xbc,0x40004c8c5c] Test LhsAdi25.C ---> FAILED ---------------------------- *** Test > : *** Test > ***Run : Obj/sgi_cc_64_DEBUG/Src/LhsAdi25.x *** Test > : -- ################################# FAILED Process 7087607 (LhsAdi25.x) stopped on signal SIGFPE: Floating point exception (default) at [solve2i:72 +0xbc,0x40004c4708] Test LhsAdi26.C ---> FAILED ---------------------------- *** Test > : *** Test > ***Run : Obj/sgi_cc_64_DEBUG/Src/LhsAdi26.x *** Test > : -- ################################# FAILED [LhsAdi::compCoefficient (float*,float*,float*,float*,float*,float*,float*,float*,float*,float*, int,int,int,int,int,int):925 +0x90,0x40004b5284] Test LhsLussor2D2Sca1.C ---> FAILED ----------------------------------- *** Test > : *** Test > ***Run : Obj/sgi_cc_64_DEBUG/Src/LhsLussor2D2Sca1.x *** Test > : -- ################################# FAILED Assertion failed: fld <= _nfldLoc, file /v2.0.11/elsA/inc/Fld/Base/FldArray.h, line 261, pid v84350 Process 7084350 (LhsLussor2D2Sca1.x) stopped on signal SIGABRT: Abort (see abort(3c)) (default) at [_prctl:15 +0x8,0xda4b0f8] Test LhsLussor2D7Sca1.C ---> FAILED ----------------------------------- *** Test > : *** Test > ***Run : Obj/sgi_cc_64_DEBUG/Src/LhsLussor2D7Sca1.x *** Test > : -- ################################# FAILED Process 7029510 (LhsLussor2D7Sca1.x) stopped on signal SIGFPE: Floating point exception (default) at [lussorsca7:387 +0x8,0x4000875930] Test LhsLussor3D2Sca1.C ---> FAILED ----------------------------------- *** Test > : *** Test > ***Run : Obj/sgi_cc_64_DEBUG/Src/LhsLussor3D2Sca1.x *** Test > : -- ################################# FAILED Process 7087143 (LhsLussor3D2Sca1.x) stopped on signal SIGSEGV: Segmentation violation (default) at [lhsupfnp:59 +0x21c,0x40004a800c] Systeme Sou =========== Test SouDts.C ---> FAILED ------------------------- *** Test > ***Run : Obj/sgi_cc_64_DEBUG/Src/SouDts.x *** Test > : -- ################################# FAILED Process 7082294 (SouDts.x) stopped on signal SIGSEGV: Segmentation violation (default) at [GeoCfdField::getDtsSolution(int):256 +0x20,0x40008dc370] 256 return *(_dtsSolution[order-1]); Systeme Tur =========== Test TurBlx1.C ---> FAILED --------------------------- *** Test > ***Run : Obj/sgi_cc_64_DEBUG/Src/TurBlx1.x *** Test > : -- ################################# FAILED Avec l'option verbose=2, on obtient: //' ===>nbErrors =8 wanted tolerance =0.001 obtained = 0.998919 Errors on mut : 8 -------------------------------------------------------- -------------------------------------------------------- **** TurBlx1 failed **** Test TurMichel.C ---> FAILED ---------------------------- *** Test > : *** Test > ***Run : Obj/sgi_cc_64_DEBUG/Src/TurMichel.x *** Test > : -- ################################# FAILED Avec l'option verbose=2, on obtient: //' ===>NbErrors = 11 wanted tolerance = 5e-07 obtained = 1 Errors on mut : 11


06/10/00 : Production de la version V2.0.11  

From Anne Jollès.
Une nouvelle production v2.0.11 est disponible. ------------------------------------------ sur oneroa13 : ELSAPRODUCTION=sgi_cc_64 sgi_cc_64_i8r8 sgi_cc_64_DEBUG sgi_cc_64_DEBUG_NEC nec_sxcc8 nec_sxcc4 ------------------------------------------- sur NEC : ELSAPRODUCTION=nec_sxcc8 nec_sxcc4 ------------------------------------------- Nouvelle version scx++ utilisee (Rev.035) cvs rtag : I2011 (Apps a ete tagge avec le meme tag) Cette production apporte des complements pour l'interface graphique PyGelsA, ameliore la vectorisation du multigrille, donne l'acces au limiteur de Zheng (en Python-API), corrige quelques (!) bugs, et modifie (legerement) les classes FldArray, FldField et FldMemory. 1. Bug dans BndSym::compWb1 (trouve par AJ) ------------------------ J'ai introduit ce bug a partir de la v2.0.01. Ce bug n'apparaissait que (c'est deja pas mal) sur NEC. Une fois corige, l'aile F4 en turbulent K-L, qui comporte une condition limite "sym" re-marche sur NEC. Attention donc si vous aviez effectue des calculs sur NEC avce condition de symetrie, avec les versions de production v2.0.01 a v2.0.10. 2. Bug dans FldArrayF::operator (trouve par SB) ---------------------------- Meme genre de bug, meme coupable, a priori lui-aussi nefaste seulement (?!!!?) sur NEC. Consequence(s) de ce bug: indeterminees... 3. Bug sur NEC (fxccen1.f) (trouve par AR et MG) ----------------------- Une ligne de trop ..., bug particulierement recalcitrant. Ce bug illustre s'il en etait encore besoin notre cruel manque de cas turbulent 3D, en non-regression et en validation. Impact: calcul Spalart 3D. 4. Modification de la destruction des objets de description -------------------------------------------------------- XJ et AJ avaient trouve un bug dans DesBase::unDesRegister: on detruisait un objet local au lieu de detruire le bon objet. Apres correction (dans la version 2.0.10), on avait un comportement "bizarre" sur NEC (dans certains cas(?!), le code ne s'arrete plus...). On revient pour l'instant a un laissez-aller bien reposant (plus de destruction, les objets mourront comme ils voudront). 5. Fusion en un seul fichier de Geo/Grid/GeoAdrCell2F.h et Geo/Grid/GeoIndCell2F.h ------------------------------------------------------------------------------- Ca permet de compiler sur NEC et linux (gcc). 6. Modification des impressions de FldMemory::statistics() ------------------------------------------------------- En fait, au lieu des 3 valeurs (flottants, entiers, booleens) la valeur reellement interessante est la valeur maxima globale (qui est en general legerement inferieure a la somme des trois precedentes). Une nouvelle methode a aussi ete ajoutee: FldMemory::max_memory() A mon avis, c'est la seule qui devrait etre utilisee par les utilisateurs, et donc decrite dans le Manuel Utilisateur; cette methode ne donne QUE la valeur max: le nombre de new et autres details ne les regardent pas, ca nous sert seulement pour voir les evolutions entre productions. 7. Changement de l'interface de FldArray et FldField ------------------------------------------------- Comme cette interface est particulierement importante, je donne ci-dessous quelques precisions. a) setAllValues(E_Float) devient setAllValuesAt b) Les constructeurs a deux arguments entiers dont le dernier est optionnel (exemple: FldArrayF(E_Int size, E_Int nfld=1); ) sont precedes par le keyword "explicit" --> evite des conversions derriere le dos du programmeur. Exemple: class A { method1(float x); method2(FldArrayF f); }; ... A a; / / faute d'ortographe / / on voulait method1 a.method2(0.); Le compilateur ne se demonte pas: il convertit le flottant 0. en l'entier 0, puis l'entier 0 en l'objet FldArray(0,1): c'est pas forcement ce qu'on voulait. Le keyword "explicit" interdit la conversion de 0 en FldArray.


22/09/00 : Production de la version V2.0.10  

From Michel Gazaix.
Une nouvelle production v2.0.10 est disponible. ------------------------------------------ sur oneroa13 : ELSAPRODUCTION=sgi_cc_64 sgi_cc_64_DEBUG sgi_cc_64_MPI sgi_cc_64_i8r8 ------------------------------------------- sur NEC : ELSAPRODUCTION=nec_sxcc8 nec_sxcc4 On rappelle qu'il faut definir pour lancer un calcul sur le NEC les 4 variables: ELSAROOT, ELSAPRODUCTION, PYTHONPATH, et PYTHONHOME. ------------------------------------------- cvs rtag : I2010e La non-regression a ete effectuee en sgi_cc_64 et sgi_cc_64_DEBUG. La validation de cette production sera effectuee la semaine prochaine. Les tests C++ (Apps) n'ont pas ete verifies (par manque de temps). Il est evident qu'il faudra de toute urgence mettre a jour les tests Apps. Je rappelle qu'il faut debugger le "make testall" (difference entre lancement des tests en interactif ou par make testall). Les principales differences % a la v2.0.09: 1. D35.4 : Enrichissement interface (M. Lazareff) Relecture non effectuee, donc a prevoir. Prise en compte complete des descriptions de mouvement (harmoniques,...) Prise en compte des raccords multi-etages, ... Cas tests (.epy et .muse) a prevoir. 2. D26.3 : Adaptation pas de temps dual (J. Raddatz) Ce travail n'avait pas pu etre integre avant le depart de JR. AJ a realise la mise a jour de l'espace de developpement (resolution des conflits, passage des test de non-regression, test de validation). Question : est-ce-que c'etait forcement a un membre de l'unite elsA de finaliser ce travail? --> Ajout d'un cas test de non-regression --> Cas test de validation a prendre en compte (il faudra resoudre un probleme pour une extraction specifique a ce cas). 3. D30.1 : Complements multi-echelle (V. Brunet) --> Modification du test de non-regression et du cas de validation Ajout/Mise a jour des tests unitaires C++ 4. D37.2 : Complements helicopteres (B. Cantaloube, JC. Boniface, A. Jolles) Ce developement permet d'introduire des mouvements plus generaux, diminue le temps CPU sur NEC (traduction C++ --> Fortran), et corrige des bugs. Plusieurs cas test de validation a prevoir. 5. D31.1 : Modifications schemas decentres (J. Peter) Le developpement a ete integre, des "ameliorations" (delicat euphemisme) devront etre apportees pour prendre en compte les remarques des relecteurs. Correction de bug (flux de Roe a 5 et 7 equations) Correction de bug calcul des residus (option ???) Une modif dans FactOper et/ou dans les descriptions sera necessaire pour pouvoir brancher correctement les aiguillages couple/non couple, flux de roe 5 ou 7 equations... --> un cas test de non-regression a ete ajoute. --> les cas test de validation avec flux de Roe seront a examiner avec attention. --> Il faut ajouter un cas de validation. 6. Complements chimere (C. Benoit, G. Jeanfaivre) Amelioration CPU, diminution drastique du nombre de new (drastique = facteur ~ 100). Nouvelle option pour l'interpolation : ADT (Alternative Digital Tree). Attention : un certain nombre de methodes ont ete integrees (en particulier dans Tmo). alors qu'elles ne sont pas utilisees dans elsA aujourd'hui. Ces methodes seront utilisees lorsque l'espace de "recherche" de CB et GJ sera integre. L'integration "en avance de phase" devrait faciliter les mises a jour ulterieures: ce sera a verifier le moment venu. En principe, un commentaire a ete insere pour justifier l'existence de ces methodes. Les nettoyeurs a l'affut des interfaces trop grasses sont donc invites a agir avec circonspection. Un renommage/rationalisation des classes chimere semble souhaitable a terme (module Blk). 7. D20.3 : Suppression TbxDico (M. Poinot) Apres cette simplification, sauf erreur, il ne reste plus de template C++ dans elsA en dehors de la STL. 8. D37.3 : Modification passage des arguments pour le calcul des rayons spectraux : Permet de pouvoir utiliser l'option "check_bound" (verification des debordements de tableaux) en Fortran. 9. D37.1 : Verification des declarations des conditions aux limites (MG): on verifie que toutes les facettes frontieres sont definies exactement une fois. (Cas particulier: condition "overlap", reflexion a faire sur la pertinence du choix effectue) 10. D28.1 : Modification Makefiles pour simplification parallele (M. Poinot) La gymnastique avec les 2 fichiers Pcm/Make_obj.mk devient inutile. Il faudra modifier la procedure decrivant les productions. 11. Modification du "make install" (M. Poinot) Il faut desormais definir la variable ELSADIST, ce qui permet de recopier librariries, executables et headers C++ dans un espace completement distinct. Par defaut on ne genere pas la grammaire ("make installmuse" genere la grammaire, "make installdoc" genere la doc HTML). Les headers sont maintenant installes (dans $ELSADIST/Dist/include). Le repository CVS Dist devient inutile. Les fichiers FactCfdDataBase.Dat et FactDataBase.dat sont maintenant dans elsA/src/Fact. Pour la livraison, il faudra penser a livrer la doc HTML avec les docs ps/pdf (auparavant, la doc HTML etait geree en CVS dans le repository Dist). 12. Corrections de Bugs Spalart + injection, paroi isotherme + intermittence, raccord avec facettes degenerees (nacelle LARA)


28/08/00 : Production de la version V2.0.09  

From Anne Jollès.
Une nouvelle production v2.0.09 est disponible. ------------------------------------------ sur oneroa13 : ELSAPRODUCTION=sgi_cc_64 ------------------------------------------- sur NEC : ELSAPRODUCTION=nec_sxcc8 PYTHONPATH =$ELSAROOT/Dist/lib/py. ------------------------------------------- cvs rtag : I2009c (pour elsA et Apps) Cette version contient: ====================== 1- "Correction du bug de l'integration de l'ALE dans la production v2.0.07 et creation de routines Fortran". (BC) 2- "Corrections pour le parallele et correction bug Sio ". (MG) 3- "Cartes de valeurs". (SB) Bilan des tests de Apps pour la production sgi_cc_64: ===================================================== Fxc : FxcNumMat3.x : OK mais avec le message: runSingle[29]: 15429772 Killed Lhs : LhsLussor3D2Sca1.x : OK mais avec le message: runSingle[29]: 14921097 Memory fault(coredump) Mask : MaskCartesianElements.x : FAILED Tur : TurMKFLC2_C.x : OK mais avec le message: runSingle[29]: 14454825 Killed Bilan de la non-regresion (execution py) ======================================== Seuls les 2 tests ALE presentent des differences dans les residus (naca64_ale.py et pf1_iso_ale.py). Des calculs complets restent maintenant a faire concernant l'ALE afin de valider la v2.0.09 . Pour la plupart des autres tests: legere diminution du nombre de creations de FldArray. Pas de validation effectuee. ============================


25/07/00 : Production de la version V2.0.08  

From Marie-Claire Le Pape.
Une nouvelle production est disponible: ELSAPRODUCTION=sgi_cc_64 cvs rtag : I2008d (pour elsA et Apps) Cette production comprend : ------------------------- Nb | Revue | Theme/Developpeur | -------------------------------------------------------------- 1 | RI-NI-00-06-03 | "Epaisseurs de deplacement | | | Sio" | | | Sebastien Beaugrand | -------------------------------------------------------------- 2 | RI-NI-00-06-04 | "Complements pour actuator | | | disk" | | | Sebastien Beaugrand | -------------------------------------------------------------- 3 | RI-NI-00-07-02 | "Reprises avec condition de | | | Froude et condition de | | | vorticite en multi-bloc" | | | Sebastien Beaugrand | -------------------------------------------------------------- 4 | RI-NI-00-07-03 |"Choix entre 3 formules de residus, | | | destruction de l'increment implicite, | | | prefixage des enumeres de Eos par eos" | | | Jacques Peter | -------------------------------------------------------------- 5 | RI-NI-00-06-07 | "Complements LU-SSOR | | | CERFACS" | | | Jacques Peter | -------------------------------------------------------------- 6 | | Introduction de la combinaison | | | FxcUpwind+insup | Bilan des tests : --------------- a. APPS : Compilations/editions de liens : OK Executions : FAILED : LhsLdu3D7S5.x MaskCartesianElements.x killed : TurMKFLC2_C.x autres messages : don't know how to make #Src/FxdTest.o (bu42). don't know how to make #Src/HstMemory.o (bu42). don't know how to make #Src/SouDts.x (bu42). b. Non-regression : Execution PY : OK, sauf quelques ecarts dans les residus : rae_lam_lss_mat1.py rae_lam_lss_mat2.py rae_lam_lss_sca1.py (faibles ecarts) oat15a_vrt.py (faibles ecarts - precedente integration) naca_eu_lss_sca1.py (faibles ecarts) naca_eu_lss_sca2.py (faibles ecarts) naca_eu_lss_mat.py flatPlateMKFLC2.py (faibles ecarts) differences lu-ssor matriciel acceptees par le CERFACS. gain en place memoire d'environ 3 - 4% pour certains cas. legere augmentation du nombre de creations de FldArra. c. Validation : non effectuee.


18/07/00 : Validation de la version V2.0.06  

From Marie-Claire Le Pape.
La version v2.0.06 d´elsA a ete validee sur SGI (mode de production sgi_cc_64_double) et sur NEC (mode de production nec_sxcc8). Les resultats obtenus sont dans l´ensemble corrects. Pour connaitre les cas-tests qui ont pose des difficultes, consulter les dossiers de validation correspondants. Les dossiers de validation correspondants n´ont ete tires qu´en un seul exemplaire chacun ; il sont disponibles en salle de reunion elsA : SGI sgi_cc_64_double : PTST-00051 NEC nec_sxcc8 : PTST-00052.


12/07/00 : Production de la version V2.0.07  

From Michel Gazaix.
Une nouvelle production est disponible: ELSAPRODUCTION=sgi_cc_64, nec_sxcc8 cvs rtag : I2007f Cette production comprend: 1. Simplification et ajouts dans Sio (RI-NI-00-02-04) 1.1 Possibilite de donner un format d'ecriture (en "fmt-v3d" ou "fmt_tp") 1.2 Traitement du format "fmt_v3d" tres peu modifie 1.3 Changement complet du traitement du format "fmt_tp" (utilisation de la librairie Fortran, a la place des iostream C++) 1.3.1 L'utilisation de la cle "overwrite" sur NEC est maintenant inutile. 1.3.2 Le nom de la fenetre associee a l'objet DesExtract apparait maintenant dans le nom de la zone Tecplot, en deuxieme ligne du fichier Tecplot. 1.3.3 Si on ne precise pas de titre, le nom de l'objet DesExtract est utilise (au lieu de "elsA" precedemment) 1.3.4 Possibilite (conservee momentanement par securite, en cas de difficulte constatee avec le nouveau traitement) de recourir a l'ancien traitement avec l'option "fmt_tp_block" 1.3.5 Mise a jour des messages d'erreur 2. Modele multi-ichelles MKFLC2 (RI-NI-00-05-03) 2 tests unitaires, un cas non regression 3. Complements ALE 4. Correction tenseur visqueux en axi 5. Correction d'un bug relatif au multi-domaine (bug qui pouvait conduire a des problemes memoire pour certains calculs multi-domaines et dans le cadre du choix de certaines options de compilation) 6. Les mots-cles correspondant a "farfield" ont ete modifies, de maniere compatible avec le Manuel Utilisateur. Dans les scripts Python (Python API ET Python utilisateur), il faut utiliser les mots-cles : "farfieldtype" "froudemod" "froudelength" Un exemple de syntaxe sur ce sujet peut etre consulte dans : /data/Wksp/jolles/Val/Validation/script/ROTOR-Michel-E2


04/07/00 : Production de la version V2.0.06  

From Michel Gazaix.
Une nouvelle production est disponible: ELSAPRODUCTION=sgi_cc_64, sgi_cc_64_double, nec_sxcc8 cvs rtag : I2006f La base de non-regression a ete passee sur SGI. La base de validation n'a pas ete passee. Cette production comprend: 1. Complements ALE (RI-NI-00-05-06): BC Pas de cas test unitaire, Cas test de non-regression a clarifier. 2. Portage sur CRAY J90 (RI-NI-00-06-05): MP elsA fonctionne sur J90, il faut probablement quelques optimisations (par exemple, directives de vectorisation FORTRAN). Une validation, meme partielle, semble indispensable. 3. Conditions de vorticite (RI-NI-00-06-02): RH Ajout d'un cas test unitaire, et d'un cas de non-regression. 4. Resolution de systemes couples - Phase 1 (RI-NI-00-06-01) : JP Un cas test de non-regression ajoute. Quelques corrections de bugs: - variable utilisee avant initialisation dans FxcRoe7F.for - Bug specifique au NEC pour Spalart-Almaras - Amelioration des messages d'erreurs lorsque le nombre de points est insuffisant Remarque 1 ------------------------------------- Il faut maintenant inclure "Def/Error/DefError.h" explicitement, lorsqu'on utilise un objet DefError. Auparavant, "DefInclude.h" incluait DefError.h. Par contre, dans tout le noyau (ou presque...) on peut (on doit) se passer de "Tbx/Libraries/TbxString.h", qui est maintenant inclus par "DefError.h". Remarque 2 -------------------------------- Les cas test unitaires: Mask/Src/MaskCartesianElements.C Fxc/Src/FxcNumMat3.C ne fonctionnent plus, je ne sais pas pourquoi.


26/05/00 : Production de la version V2.0.05  

From Michel Gazaix.
Une nouvelle production est disponible: ELSAPRODUCTION=sgi_cc_64, nec_sxcc8, sgi_cc_64_i8r8 cvs rtag : I2005a La base de non-regression a ete passee sur SGI (sgi_cc_64). Les tests unitaires Apps ont ete passes (sgi_cc_64, sgi_cc_64_i8r8). Un cas test unitaire ne fonctionne pas en sgi_cc_64_i8r8 (DesMask). La base de validation n'a pas ete passee. Cette version contient: 1. Les premiers elements correspondant a la fonctionnalite "Actuator Disk" (S. Beaugrand & J. Brezillon (DAAP)). Il n'existe aucun cas test unitaire (C++) ou d'integration. Ce developpement sera complete au retour de vacances de Brezillon. 2. Des complements au developpement "reduction consommation memoire". Ces complements ont ete realises par S. Beaugrand, suite a la revue effectuee par J. Peter. 3. Prise en compte de la cinematique d'un rotor et introduction d'une condition limite champ lointain instationnaire (J.C. Boniface, A. Jolles). 4. K-OMEGA avec lois de parois (R. Houdeville (DMAE)). Il faudra penser a ajouter un cas test de validation. 5. Un renommage de la classe DesCfdSys en EosSysEq (M. Poinot). Les fichiers associes ont ete deplaces dans Eos/Base. Cette modification doit permettre de mieux separer ce qui releve du noyau CFD (kernel) de ce qui releve de la description d'un calcul. Cette modification, certes mineure -) pourrait peut-etre entrainer certaines difficultes aux developpeurs intrepides qui voudraient se synchroniser "la fleur au fusil". En effet, plusieurs centaines de fichiers sont modifies. 3 possibilites existent: a) utiliser Sniff pour effectuer le remplacement systematique, puis cvs. Incidemment, l'operation pourrait constituer un premier exemple d'utilisation de Sniff. b) utiliser un script (appelant sed, ou awk, ou python) pour faire les substitutions, puis cvs. c) faire confiance a CVS, pour resoudre les conflits "automatiquement". J'ai fait un essai sur l'espace de B. Cantaloube: la solution c) a marche du tonnerre, mais il s'etait "pre-integre" sur la version qui precedait immediatement le renommage de DesCfdSys. Son cas n'est donc pas forcement le cas general. En tout cas, je conseille de faire d'abord un essai avec c), sur un seul fichier. Si ca marche, continuer... Si ca ne marche pas convenablement, n'hesitez pas a contacter Marc Poinot. 6. Amelioration des temps d'interpretaion en python-utilisateur. + commandes relatives au chimere + manuel relatif aux commandes python-utilisateur (M. Lazareff). 7. Optimisation CPU des routines d'inversion ADI sur NEC (G. Bastin) Cette optimisation est decrite dans la "Note Technique elsA NI-00025". L'optimisation principale consiste a vectoriser par plan, au lieu de vectoriser par ligne. Sur NEC, malheureusement, cette vectorisation produit parfois des resultats TRES decevants, par suite d'une grande sensibilite de cet ordinateur aux problemes de "conflits de bancs memoire". Pour detecter ce probleme, la seule solution est de regarder attentivement le resultat de la commande "ntrace". Si le nombre de MegaFlops pour la routine solve5i est de l'ordre de 300, alors il y a comme un defaut (il apparait aussi un nombre anormalement eleve dans la colonne "BANK CONFLICT"). Ce probleme est malheureusement assez frequent, en particulier pour les maillages multigrilles (voir note NI-00025 pour l'explication). On propose une solution (partielle) a ce probleme, qui consiste a supprimer les conflits d'acces memoire en rangeant les donnees utilisees dans la phase implicite de maniere NON CONSECUTIVES. (on consomme donc un peu plus de memoire). Cettepossibilite d'optimisation sur NEC est encore EXPERIMENTALE: il est fortement conseille de verifier que le resultat est identique, et que les performances sont vraiment ameliorees. N'hesitez pas a nous communiquer les resultats de vos experiences... Concretement, il suffit de modifier le fichier de donnee en ajoutant(ici, desMesh est suppose etre le nom d'un objet maillage de type DesMesh): desMesh.setI("padding_i", 1) # python-elsA_API desMesh.set("padding_i",1) # python-elsA Cette optimisation a necessite une reorganisation (partielle) de Lhs (M. Gazaix): - factorisation de code commun dans LhsBaseImp - uniformisation du passage des donnees de Eos (gamma, prandtl,...) dans les constructeurs - transformation (encore imcomplete) d'adressage de type "lhs1,lhs2, ,lhs5" en lhs(ncell, neq). Cet adressage permettra dans l'avenir une optimisation differente sur machine vectoriel et machine RISC. - Transformation d'enum globaux en enum de classe. Cela evite de polluer le namespace global, et identification claire du lieu de definition. Ex: LhsAdi::Euler au lieu de LhsAdiEuler - introduction d'un nouveau Fortran dans Tmo (tmocpsoln_) appele par TmoSolverBaseP::advanceStep Il se peut que le multigrille local ne fonctionne plus; je suis pret a faire l'effort de le faire remarcher, si on me fournit les indications necessaires (en particulier, des cas tests). Notons encore que les implicites non ADI (Lussor, Ldu) ne sont pas du tout optimises sur NEC. 8. Suppression des deux ALLOCATE/DEALLOCATE (delta3_ et delta4_). Je n'ai malheureusement pas pu tester la modif, n'ayant pas trouve de cas tests. Quelqu'un en aurait-il dans sa musette? Plus generalement, nous manquons cruellement d'outils de couverture de test.


22/05/00 : Validation de la version V2.0.03  

From Marie-Claire Le Pape.
La version v2.0.03 d'elsA a ete validée avec succès sur SGI, mode de production sgi_cc_64. Le dossier de validation correspondant n´a été tiré qu´en un seul exemplaire, disponible en salle de réunion elsA (5f07) : document relié /ELSA/PTST-00034, v1.0 .


05/05/00 : Production de la version V2.0.03  

From Marc Poinot.
Une nouvelle production est disponible : ELSAPRODUCTION=sgi_cc_64 cvs rtag : I2003 Les nouveautés de cette production ne résident que dans les fichiers Python : - Chimère en Python utilisateur, - corrections de bugs. Attention, il y a toujours un bug dans la conversion des groupes Python API en Python utilisateur !


02/05/00 : Production de la version V2.0.02  

From Michel Gazaix.
Une nouvelle production est disponible: ELSAPRODUCTION=sgi_cc_64 cvs rtag: I2002 La principale nouveaute de cette production concerne l'implicite LDU-SSOR, introduit par Marc Montagnac, du CERFACS. (L'adaptation a ete realisee par J. Peter, qui a aussi ajoute des elements relatifs aux implicites visqueux). Des complements chimeres (description) ont ete egalement ajoutes (MCLP). Les cas tests unitaires de Apps ont ete augmentes (par 8 tests relatifs a l'implicite LU-SSOR) et verifies. Une non-regression a ete effectuee sur la base de non-regression de l'integration. Par contre, il n'a pas ete effectue de validation complete.


20/04/00 : Production de la version V2.0.01  

From Michel Gazaix.
Une nouvelle production est disponible: ELSAPRODUCTION=sgi_cc_c4, sgi_cc_64_DEBUG, nec_sxcc8 cvs rtag: I2001 La principale nouveaute de cette production concerne le chimere (MCLP/CB) Egalement, introduction d'une nouvelle configuration (tmic8) pour le modele de Michel (JCB), complement pour le calcul des epaisseurs integrales (R. Houdeville), vectorisation du calcul des distances (updist1) (Houdeville) condition limite "farfield" avec equation de transport (AJ). Correction de plusieurs bugs: - calcul du rayon spectral avec moyenne arithmetique - initialisation du modele de Michel Les cas tests unitaires de Apps ont ete augmentes et verifies. Une non-regression a ete effectuee: sur la base de non-regression pour l'integration, il n'y a aucune difference entre v1.2.07 et v2.0.01. Par contre, aucune validation complete n'a ete effectuee.


19/04/00 : Production de la version V1.2.07  

From Michel Gazaix.
La version v1.2.07 est disponible pour les developpeurs: -- sur SGI (oneroa13) ELSAPRODUCTION=sgi_cc_64, sgi_cc_64_MPI, sgi_cc_64_i8r8 -- sur NEC: nec_sxcc8 Dans la mesure ou cette version a ete livree, et a donc ete decrite de maniere relativement exhaustive dans le document accompagnant la livraison, nous ne donnons pas d'informations supplementaires sur le contenu. Nous profitons de cette nouvelle production pour annoncer que les tests unitaires ont ete reorganises. Tous les tests unitaires sont rassembles dans une arborescence UNIX separee (Apps), et sont geres par CVS dans un repository distinct (CVSROOT=$CVSROOTAPP) On a verifie que tous les tests explicitement declares dans les fichiers Make_tst.mk de Apps fonctionnent correctement avec les librairies v1.2.07. A l'occasion de chaque integration, nous verifierons que l'ensemble de ces cas tests continuent a fonctionner. Bien sur, beaucoup d'integrations ajouteront de nouveaux cas tests unitaires, qui seront ajoutes a notre "banque" de cas test. En outre, un certain nombre de cas tests, qui demeurent pertinents, n'ont pas pu, faute de temps, etre pris en compte: en general, ils ne compilaient plus, ou ils renvoyaient des resultats faux ou ambigus a l'execution. Ces cas tests sont pour la plupart desactives, dans les fichiers Make_tst.mk. Il est evident qu'il est indispensable de recuperer ces cas tests, qui permettent de completer la couverture de test d'elsA. Anne Jolles a la responsabilite de suivre la mise a jour de ces cas tests: elle demandera a certains d'entre vous (peut-etre l'a-t-elle deja fait dans certains cas) une "assistance" pour la mise a jour. En general, la personne qui a ecrit le test est la mieux placee pour le mettre a jour.


14/04/00 : Developer's Guide  

From Anne-Marie Vuillot.
La première version du guide du développeur est disponible en version PostScript (5.4 Mo). Elle sera diffusée à partir du lundi 17/04/00 associée dans un classeur avec les planches du modèle objet pour former le vade-mecum du developpeur elsA.


08/03/00 : Production de la version V1.2.04  

From Michel Gazaix.
Une nouvelle production est disponible sur SGI (ELSAPRODUCTION=sgi_cc_64, tag CVS = I1204). Pour la non-regression, il faut positionner la variable OLDREF a Ref1204f (dans le Makefile de non-regression resultant de la commande "make install ..."). Cette version contient probablement un bug pour le pas de temps dual. Cette version n'a pas ete completement validee, elle est donc "instable". Sa mise a disposition "anticipee" vise a permettre une synchronisation aussi rapide que possible, notamment pour les developpements attendus pour la prochaine livraison : chimere, gain memoire, calcul residus, optimisation CPU en multigrille...).


17/02/00 : Production de la version V1.2.02  

From Michel Gazaix.
Une nouvelle production est disponible : v1.2.02 Par manque de temps, elle n'a pu etre produite qu'en sgi_cc_64 (simple precison SGI). Cette version contient: 1. Complement ALE 2. Complement Pas de temps dual 3. Nouvelle condition au limite BndInj2, BndSubInj2,... 4. Des elements d'optimisation CPU - calcul des rayons spectraux - vectorisation Cond. aux limites repere mobile 5. Complement lois de paroi et Spallart 6. Corrections bugs: - initialisation en turbulent - calcul axi visqueux 7. Complement MUSE-elsA


02/02/00 : Production de la version V1.2.01  

From Michel Gazaix.
Une nouvelle production : v1.2.01 (CVS: rtag=I1201) est disponible. ------------------------------------------------------- ELSAPRODUCTION=sgi_cc_64 ou sgi_cc_64_double --------------------------------------------------------- ELSAPRODUCTION=nec_sxcc8 et nec_sxcc4, --------------------------------------------------------- ELSAPRODUCTION=dec_cxx_optimize, ------------------------------------------------- ! ATTENTION - NEC ! : ------------------------------------------------- Nous attirons l'attention des heureux utilisateurs d'elsA sur NEC sur le point TRES IMPORTANT suivant: La nouvelle production v1.2.01 a ete produite avec le nouveau compilateur C++ NEC. En consequence, on peut craindre certaines anomalies: prudence donc. Si cela s'averait necessaire, on pourra toujours utiliser l'ancien compilateur. Pour l'instant, Il est fortement conseille d'utiliser le mode python-API. ------------------------------------------------- NOUVEAUTES ------------------------------------------------- 1) Premiere integration des lois de parois ------------------------------------------ Developpement non acheve 2) Developpement de la condition limite "champ lointain" ("farfield") --------------------------------------------------------------------- Possibilite de prendre en compte les vitesses de Froude dans le cas : repere relatif, formulation vitesses absolue, cas stationnaire , vitesse de translation nulle. Pas de restriction sur l'utilisation de la condition "champ lointain" sans Froude. L'association multigrille et c.l. champ lointain avec (et sans) Froude fonctionne. 3) Action d'amelioration de la qualite dans le module Fxd --------------------------------------------------------- - Suppression de compIntOfBorder dans FxdFlux et FxdTurBoussinesq - Modification des signes dans les formules de calcul de flux dans un souci de coherence avec les notations du Manuel Theorique - Amelioration notable du respect de la nomenclature - Suppression de FxdDiffIntF.for (a confirmer) - Verification et amelioration des commentaires - Ameliorations mineures de codage (par ex. , suppression de scalaire temporaire inutile) 4) Modification du calcul de muTur en debut de calcul ou en cas de reprise -------------------------------------------------------------------------- Dans les versions precedentes, muTur etait initialise a 0 en debut de calcul (methode solutionNEqualSolutionP) au lieu d'etre calcule. Maintenant muTur est calcule, avec un traitement encore insatisfaisant dans le cas des modeles algebriques. En effet dans ce cas, le calcul du gradient de la vitesse necessaire au calcul du muTur est effectue en ne tenant pas compte des conditions de raccord (mais seulement des conditions limites physiques). Ceci sera bien sur a ameliorer lorsque la nouvelle version des raccords sera livree. 5)! Fichiers de donnees modeles de turbulence algebriques a MODIFIER ! ---------------------------------------------------------------------- Un bug ayant ete corrige, il faut desormais indiquer explicitement le type de distance choisi pour effectuer un calcul en modele algebrique (cle walldistcompute). 6) Possibilite d'associer Flux decentre (Roe, van Leer) ------------------------------------------------------- et calcul turbulent avec equation de transport. ----------------------------------------------- 7) Elements de la methode Chimere relatifs aux transferts entre grille ---------------------------------------------------------------------- 8) Extensions des sorties ------------------------- Calcul des moments,possibilites d'ajouter les flux ... Acces MUSE-elsA pour ces extensions. 9) Developpements permettant de diminuer le temps d'interpretation des ---------------------------------------------------------------------- commandes elsA. --------------- 10) Correction de bug divers et varies --------------------------------------


17/01/00 : Production de la version V1.1.12  

From Michel Gazaix.
Une nouvelle production : v1.1.12 (rtag=I1112) est disponible. Cette production fera l'objet d'une LIVRAISON, avec Proces Verbal de Livraison, qui decrira de maniere plus complete le contenu de la livraison. ------------------------------------------------------- ELSAPRODUCTION=sgi_cc_64 ou sgi_cc_64_double ou sgi_cc_64_MPI --------------------------------------------------------- ELSAPRODUCTION=NEC_edg_Opt8 et NEC_edg_Opt4, --------------------------------------------------------- ELSAPRODUCTION=dec_cxx_optimize, --------------------------------------------------------- Par rapport a la production precedente, cette version corrige des bugs subtils concernant le multi-grille, et ameliore notablement la consommation memoire. A titre indicatif, entre la presente version et la derniere version livree a l'exterieur de DSNA (c'est-a-dire la v1.0.07) la memoire max. est divisee par un facteur 1.8 (cas fuselage helico, 10 domaines, turbulent k-L, multigrille). Liste des modifications (sauf oubli ...) =========================================================================== Contribution JP/OR--> V1112 a - suppression de la copie du champ conservatif dans OperFlux::compute b - correction du bug multigrille/validation c - calcul du pas de temps argument de advanceIteration, calcul du pas de temps dans le cycle multigrille. d - suppresion des calculs de metrique par restriction pour les grilles grossieres. Creation dynamique des maillages grossiers dans BlkElemBlock. e - raccourcissement de l'interface publique de FldTrans et du source appelant f - Modification GeoCfdField des raccords (un champ aero et non plus deux ) g - suppresion de la possibilite de multigrille sans forcing function (interet tres limite pour la complexite induite) ------------------------------------------------------------------ --> le travail fait dans BlkElemBlock appelera sans doute un travail equivalent dans BlkMotionBlock que je n'ai pas fait (me demander le cas echeant). --> la modification de Join est a reporter dans les nouvelles versions ------------------------------------------------------------------- --> le gain de memoire induit par a/ et f/ induit - malgre e/ en multigrille - est de l'ordre de 5 a 10% sur l'ensemble des calculs de la base de validation des integrations. -------------------------------------------------------------------- --> il faudrait des specifications tres precises sur les besoins futurs des uns et des autres (ale, chimere,turbulence) pour savoir si on peut detruire le maillage des niveaux grossiers apres construction des grilles. Qui a besoin, sur les niveaux grossier, de calculer des centres, de detruire la metrique et de la recalculer a chaques iteration .... ---------------------------------------------------------------------- ======================================================================== Contribution AJ/BC Supression de la double creation des grilles (gain en place memoire) ========================================================================= Contribution AJ -- Branchement du calcul en Fortran du pas de temps turbulent sur NEC (gain en CPU) -- Limitation de l'invalidation de muTur conforme a CANARI: muTur est invalide apres resolution des 4 sous-pas de Runge-Kutta du systeme turbulent (donc une seule fois par iteration) -- Dimensionnement du vecteur aero limite a 5 champs sur les grilles grossieres (en cas de multigrille) au lieu de 7 (en turbulent) (gain en place memoire) ============================================================ Contribution SB -- Suppresion de la duplication des variables conservatives dans Blk. -- Mise a jour de Sio ============================================================= Attention: modification d'une cle utile aux calculs visqueux en reperes tournants!! La cle "etapm" est maintenant obsolete et doit etre remplacee par "mobile_coef" Acces (en python-API uniquement) a la possibilite de faire des sous-iterations sur les systemes champs moyens et champs turbulents. Mise a jour de la documentation HTML des sources (doc++). La documentation HTML parait relativement correcte, et merite sans doute d'etre plus consultee.


20/12/99 : Production de la version V1.1.11  

From Michel Gazaix.
Une nouvelle production : v1.1.11 (rtag=I1111) est disponible: ------------------------------------------------------- ELSAPRODUCTION=sgi_cc_64 ou sgi_cc_64_double --------------------------------------------------------- ELSAPRODUCTION=NEC_edg_Opt8 et NEC_edg_Opt4 --------------------------------------------------------- ELSAPRODUCTION=dec_cxx_optimize --------------------------------------------------------- Cette version contient: Avec revue: =================================================================== RI-NI-99-12-04 : Implicite van Leer RI-NI-99-12-03 : Residus, portance pour Froude, reduction memoire Sans revue: =================================================================== Diverses corrections de bug : ---------------------- 1. Algebrique Michel pour rotor 2. Correction FxcRoeCorr et FxcHarCorr si maillage degenere 3. Bug sur NEC pour l'application du limiteur en turbulent decouple 4. Bug(s) intervenant sur la double precision: il en reste surement encore, mais on a un "compromis" qui semble marcher (sur la base de non-regression d'integration en tout cas) sur SGI (sgi_cc_64_double). Sur SGI: Les entiers sont sur 4 octets, les reels sur 8 Une validation complete sur SGI en sgi_cc_64 et sgi_cc_64_double semble s'imposer.


07/12/99 : Production de la version V1.1.09  

From Michel Gazaix.
Une nouvelle production : v1.1.09 (rtag=I1109a) est disponible: ------------------------------------------------------- ELSAPRODUCTION=sgi_cc_64 ou sgi_cc_64_double --------------------------------------------------------- ELSAPRODUCTION=NEC_edg_Opt8 et NEC_edg_Opt4 --------------------------------------------------------- ELSAPRODUCTION=dec_cxx_optimize --------------------------------------------------------- Cette version contient: Avec revue: =================================================================== RI-NI-99-12-01 : Pas de temps dual pour l'instationnaire (Il faudra finaliser les remarques faites au cours de la revue.) Sans revue: =================================================================== 1. Modif "Qualite" : --------------- - BndBase remplace desormais BndOperator (le nom "BndOperator" etait source de confusion, puisque les Bnd ne sont PAS des opeateurs, au sens d'elsA). - menage proprement dit (voir liste ci-dessous) - correction include fortran (copmilation DEC) - correction du bug dans FxcHarrCorrF.for pour NEC - descente de wb1 dans OperFlux::computeAllBorders, (disparition de RhsTerm) disparition de wb1 des classes gradients (sauf dans OperGradGradVelocity qui a un mecanisme particulier et ou j'ai du faire une tambouille un peu malheureuse). wb1 apparait encore sous sa forme ectoplasmique dans TurBousinesq. Il faudrait continuer a le faire disparaitre en descendant jusqu'au niveau ou on le dimensionne avant de l'utiliser effectivement. 2. Pour l'axi: 2.1/ prise en compte des grandeurs turbulentes dans la c.l. BndAxiSym 2.2/ association LDU et Axisym rendue possible 3. Pour les reperes relatifs: 3.1/correction calcul de la vitesse d'entrainement : 1 seule fois 3.2/ implementation du pas de temps convectif en Fortran 4. Correction bug compMutInModel dans TurKO.C et TurKEps.C. (devrait permettre un traitement correct de l'intermittence en KO et K-EPS) 5. Variante du modele de Michel pour les helico 6. Affichage plus complet en cas de probleme memoire: On affiche maintenant le nombre de bytes demande lorsque l'allocation a echoue. ainsi que le maximum de memoire precedent (Ceci repond au mail Snecma (F. Bastin) du 3 Decembre). 7. "Le rab de Join" (SIC) decele par J. Peter est corrige (mail JP du 6 Decembre) On obtient un gain en temps et memoire (Ex. sur SGI: helico 10 blocks NS (1 200 000 / 8 ) points t(CPU) : avant 308s (pour 10 iter), apres : 273s memoire max avant : 259 MegaOctets, apres : 226 MegaOctets



(updated $Date: 2016/04/28 13:06:51 $)Home Site map Contacts