Nos partenaires et nous-mêmes utilisons différentes technologies, telles que les cookies, pour personnaliser les contenus et les publicités, proposer des fonctionnalités sur les réseaux sociaux et analyser le trafic. Merci de cliquer sur le bouton ci-dessous pour donner votre accord. Vous pouvez changer d’avis et modifier vos choix à tout moment. Informations RGPD
#!/usr/bin/env python2 # -*- coding: utf-8 -*- """ Created on Sun Apr 9 14:11:37 2017 @author: john """ import os import shutil import time as t import numpy as n ti=t.time() ########################################## filE='2D.dgibi' #parametres para1='P' para2='NO' Para=[para1,para2] Obs=['MXSSVM','DAX','ATAN'] X1=list(n.linspace(0.2,0.74,8)) #X1=[1,10,100,1000] X2=list(n.linspace(0,2500,8)) ########################################## cwdo=os.getcwd() os.chdir(cwdo) #creation du fichier stock if 'stock' not in os.listdir(cwdo): os.mkdir(cwdo+'/stock') print('stock n existait pas') os.chdir(cwdo+'/stock') cwd=os.getcwd() #creation du fichier experience try: shutil.rmtree(cwd+'/stock '+para1+'-'+para2) print('/stock '+para1+'-'+para2+'/ para remove') except: pass print('/stock '+para1+'-'+para2+'/ n existait pas') os.mkdir(cwd+'/stock '+para1+'-'+para2) def modif_original(Para,Obs): os.chdir(cwdo) file_handle = open(filE, 'r') N = file_handle.read().split('\n') file_handle.close() M=[] for k in range(len(N)): M+=[N[k]] if ('FIN' in N[k].upper()) and ('*' not in N[k]): del M[k] file_handle.close() M+=["i=0;COMM para;"] M+=["chsort = CHAINE 'sort' i;"] M+=["OPTION SORT chsort;"] M+=["TAB=TABLE;"] mp='' for k in Para: mp+=(k + ' \' ' +k +' \' ') mo='' for k in Obs: mo+=(k + ' \' ' +k +' \' ') M+=["TAB.X1=CHAINE " +mp+";"] M+=["SORT 'CHAI' (TAB.X1);"] M+=["TAB.X2=CHAINE " +mo+";"] M+=["SORT 'CHAI' (TAB.X2);"] M+=["fin;"] os.chdir(cwd+'/stock '+para1+'-'+para2+'/') file_handle = open('2D_base.dgibi', 'w') file_handle.write('\n'.join(M)) file_handle.close() def remplpara(Para, New, fileindex): filE='2D_base.dgibi' file_handle = open(filE, 'r') M = file_handle.read().split('\n') file_handle.close() for p in range(len(Para)): i=0 c=True while c: if (Para[p] in M[i]) and ('COMM para' in M[i]): M[i]=Para[p] + ' = ' + str(New[p]) +'; COMM para;' c=False i+=1 file_handle = open('2D'+str(fileindex)+'.dgibi', 'w') file_handle.write('\n'.join(M)) file_handle.close() modif_original(Para,Obs) print('modif_original : ok') os.chdir(cwd+'/stock '+para1+'-'+para2) L=[] i=0 tr=[] tc=[] for k in X1: for l in X2: i+=1 remplpara(Para+['i'],[k,l,i],i) L+=[''] L[-1] = os.popen('castem16 '+'2D'+str(i)+'.dgibi') print('remplpara') tf=t.time() print(tf-ti)