# -*- coding: UTF-8 no BOM -*- import re import os class Section(): def __init__(self,data = {'__order__':[]},part = ''): classes = { 'homogenization':Homogenization, 'microstructure':Microstructure, 'crystallite':Crystallite, 'phase':Phase, 'texture':Texture, } self.parameters = {} for key in data: if type(data[key]) is not list: self.parameters[key] = [data[key]] #this is creating a dictionary. So this is a way to create elements in a dictionary else: self.parameters[key] = data[key] if '__order__' not in self.parameters: self.parameters['__order__'] = list(self.parameters.keys()) if part.lower() in classes: self.__class__ = classes[part.lower()] self.__init__(data) def add_multiKey(self,key,data): multiKey = '(%s)'%key if multiKey not in self.parameters: self.parameters[multiKey] = [] if multiKey not in self.parameters['__order__']: self.parameters['__order__'] += [multiKey] if type(data) == list: self.parameters[multiKey] += [[item] for item in data] else: self.parameters[multiKey] += [[data]] def data(self): return self.parameters class Homogenization(Section): def __init__(self,data = {'__order__':[]}): Section.__init__(self,data) class Crystallite(Section): def __init__(self,data = {'__order__':[]}): Section.__init__(self,data) class Phase(Section): def __init__(self,data = {'__order__':[]}): Section.__init__(self,data) class Microstructure(Section): def __init__(self,data = {'__order__':[]}): Section.__init__(self,data) class Texture(Section): def __init__(self,data = {'__order__':[]}): Section.__init__(self,data) def add_component(self,theType,properties): if 'scatter' not in list(map(str.lower,list(properties.keys()))): scatter = 0.0 else: scatter = properties['scatter'] if 'fraction' not in list(map(str.lower,list(properties.keys()))): fraction = 1.0 else: fraction = properties['fraction'] try: multiKey = theType.lower() except AttributeError: pass if multiKey == 'gauss': self.add_multiKey(multiKey,'phi1 %g\tPhi %g\tphi2 %g\tscatter %g\tfraction %g'%( properties['eulers'][0], properties['eulers'][1], properties['eulers'][2], scatter, fraction, ) ) if multiKey == 'fiber': self.add_multiKey(multiKey,'alpha1 %g\talpha2 %g\tbeta1 %g\tbeta2 %g\tscatter %g\tfraction %g'%( properties['eulers'][0], properties['eulers'][1], properties['eulers'][2], properties['eulers'][3], scatter, fraction, ) ) class Material(): """Reads, manipulates and writes material.config files""" def __init__(self,verbose=True): """Generates ordered list of parts""" self.parts = [ 'homogenization', 'microstructure', 'crystallite', 'phase', 'texture', ] self.data = {\ 'homogenization': {'__order__': []}, 'microstructure': {'__order__': []}, 'crystallite': {'__order__': []}, 'phase': {'__order__': []}, 'texture': {'__order__': []}, } self.verbose = verbose def __repr__(self): """Returns current configuration to be used as material.config""" me = [] for part in self.parts: if self.verbose: print('doing '+part) me += ['','#-----------------------------#','<%s>'%part,'#-----------------------------#',] for section in self.data[part]['__order__']: me += ['','[%s] %s'%(section,'-'*max(0,27-len(section))),'',] for key in self.data[part][section]['__order__']: if key.startswith('(') and key.endswith(')'): # multiple (key) me += ['%s\t%s'%(key,' '.join(values)) for values in self.data[part][section][key]] else: # plain key me += ['%s\t%s'%(key,' '.join(map(str,self.data[part][section][key])))] return '\n'.join(me) def parse_data(self, part=None, sections=[], content=None): re_part = re.compile(r'^<(.+)>$') # pattern for part re_sec = re.compile(r'^\[(.+)\]$') # pattern for section re_include = re.compile(r'^{(.+)}$') # pattern for include name_section = '' active = False for line in content: line = line.split('#')[0].strip() # kill comments and extra whitespace line = line.split('/echo/')[0].strip() # remove '/echo/' tags line = line.lower() # be case insensitive if line: # content survives... match_part = re_part.match(line.split()[0]) if match_part: # found separator active = (match_part.group(1) == part) # only active in continue if active: match_sec = re_sec.match(line.split()[0]) if match_sec: # found [section] name_section = match_sec.group(1) # remember name ... if '__order__' not in self.data[part]: self.data[part]['__order__'] = [] self.data[part]['__order__'].append(name_section) # ... and position self.data[part][name_section] = {'__order__':[]} continue # This part accesses the data in the files mentioned in the curly brackets. For now the file should be present in the working folder #trying to generalize the file and its location. if line: match_include = re_include.match(line.split()[0]) if match_include: sub_file = match_include.group(1).split('/')[1] #get the file name which has to be opened to get the data #code to search for file when it is located anywhere path_file = os.environ['PATH'] path_file = path_file[:(path_file.find("DAMASK") + len("DAMASK"))] for root, dirs, files in os.walk(path_file): if sub_file in files: target_path = os.path.join(root, sub_file) #this would give the path for the sub_file with open(target_path,'r') as f1: #using the target_path to open the file subdata = f1.readlines() subdata = [a for a in subdata if a != '\n'] #removes '\n' elements from the list if subdata: for entry in subdata: if sections == [] or name_section in sections: # respect subset items = entry.split() if items[0] not in self.data[part][name_section]: # first encounter of key? self.data[part][name_section][items[0]] = [] # create item self.data[part][name_section]['__order__'].append(items[0]) if items[0].startswith('(') and items[0].endswith(')'): # multiple "(key)" self.data[part][name_section][items[0]].append(items[1:]) else: # plain key self.data[part][name_section][items[0]] = items[1:] continue if sections == [] or name_section in sections: # respect subset items = line.split() if items[0] not in self.data[part][name_section]: # first encounter of key? self.data[part][name_section][items[0]] = [] # create item self.data[part][name_section]['__order__'].append(items[0]) if items[0].startswith('(') and items[0].endswith(')'): # multiple "(key)" self.data[part][name_section][items[0]].append(items[1:]) else: # plain key self.data[part][name_section][items[0]] = items[1:] def read(self,file=None): re_include = re.compile(r'^{(.+)}$') # pattern for include c = [] with open(file,'r') as f: for entry in f: match_include = re_include.match(entry.split()[0]) if match_include: sub_file = match_include.group(1).split('/')[-1] #get the file name which has to be opened to get the data, #code to search for file when it is located anywhere path_file = os.environ['PATH'] path_file = path_file[:(path_file.find("DAMASK") + len("DAMASK"))] for root, dirs, files in os.walk(path_file): if sub_file in files: target_path = os.path.join(root, sub_file) #this would give the path for the sub_file with open(target_path,'r') as f1: #using the target_path to open the file for data in f1: c = c.append(data) else: c = c.append(entry) for p in self.parts: self.parse_data(part=p, content=c) def write(self,file='material.config', overwrite=False): import os i = 0 saveFile = file while not overwrite and os.path.exists(saveFile): i += 1 saveFile = file+'_%i'%i if self.verbose: print('Writing material data to file %s'%saveFile) f=open(saveFile,'w') f.write(str(self)+'\n') #newline at end f.close() return saveFile def add_section(self, part=None, section=None, initialData=None, merge = False): """adding/updating""" part = part.lower() section = section.lower() if part not in self.parts: raise Exception('invalid part %s'%part) if type(initialData) is not dict: initialData = initialData.data() if section not in self.data[part]: self.data[part]['__order__'] += [section] if section in self.data[part] and merge: for existing in self.data[part][section]['__order__']: # replace existing if existing in initialData['__order__']: if existing.startswith('(') and existing.endswith(')'): # multiple (key) self.data[part][section][existing] += initialData[existing] # add new multiple entries to existing ones else: # regular key self.data[part][section][existing] = initialData[existing] # plain replice for new in initialData['__order__']: # merge new content if new not in self.data[part][section]['__order__']: self.data[part][section][new] = initialData[new] self.data[part][section]['__order__'] += [new] else: self.data[part][section] = initialData def add_microstructure(self, section='', components={}, # dict of phase,texture, and fraction lists ): """Experimental! Needs expansion to multi-constituent microstructures...""" microstructure = Microstructure() # make keys lower case (http://stackoverflow.com/questions/764235/dictionary-to-lowercase-in-python) components=dict((k.lower(), v) for k,v in components.items()) for key in ['phase','texture','fraction','crystallite']: if type(components[key]) is not list: try: components[key] = [components[key].lower()] except AttributeError: components[key] = [components[key]] else: for i, x in enumerate(components[key]): try: components[key][i] = x.lower() except AttributeError: pass for (phase,texture,fraction,crystallite) in zip(components['phase'],components['texture'], components['fraction'],components['crystallite']): microstructure.add_multiKey('constituent','phase %i\ttexture %i\tfraction %g\ncrystallite %i'%( self.data['phase']['__order__'].index(phase)+1, self.data['texture']['__order__'].index(texture)+1, fraction, self.data['crystallite']['__order__'].index(crystallite)+1)) self.add_section('microstructure',section,microstructure) def change_value(self, part=None, section=None, key=None, value=None): if not isinstance(value,list): if not isinstance(value,str): value = '%s'%value value = [value] newlen = len(value) oldval = self.data[part.lower()][section.lower()][key.lower()] oldlen = len(oldval) print('changing %s:%s:%s from %s to %s '%(part.lower(),section.lower(),key.lower(),oldval,value)) self.data[part.lower()][section.lower()][key.lower()] = value if newlen is not oldlen: print('Length of value was changed from %i to %i!'%(oldlen,newlen))