control2d.py

À propos

Type de fichier

Fichier PY de 3 Ko (text/x-python)

Confidentialité

Fichier public, envoyé le 21 avril 2017 à 17:55, depuis l'adresse IP 162.38.x.x (France)

Sécurité

Ne contient aucun Virus ou Malware connus - Dernière vérification: 3 jours

Statistiques

La présente page de téléchargement a été vue 214 fois depuis l'envoi du fichier

Page de téléchargement

https://www.petit-fichier.fr/2017/04/21/control2d/   Copier

Aperçu du fichier

#!/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)

Partager le fichier

Télécharger control2d.py