223031012
/
DAMASK_EICMD
mirror of https://github.com/achalhp/DAMASK_EICMD.git
1
0
Fork 0
Code Issues Projects Releases Wiki Activity

Merge branch 'development' of magit1.mpie.de:damask/DAMASK into development

This commit is contained in:
Martin Diehl 2018-02-03 09:02:20 +01:00
parent d80e15bd76 911fb84e81
commit ba6fd70654
5 changed files with 78 additions and 49 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

2
PRIVATE

@ -1 +1 @@
Subproject commit bdbf2da71cd9e0825d17f673ec2fbabc2c8027f8 Subproject commit 5fc3188c86ea1f4159db87529ac3e3169fb56e5d

2
VERSION
View File

@ -1 +1 @@
v2.0.1-1029-ge304ce3 v2.0.1-1035-gd80a255

53
lib/damask/test/test.py
View File

@ -322,8 +322,10 @@ class Test():
cur1Name = self.fileInCurrent(cur1) cur1Name = self.fileInCurrent(cur1)
return self.compare_Array(cur0Name,cur1Name) return self.compare_Array(cur0Name,cur1Name)
def compare_Table(self,headings0,file0,headings1,file1,normHeadings='',normType=None, def compare_Table(self,headings0,file0,
absoluteTolerance=False,perLine=False,skipLines=[]): headings1,file1,
normHeadings='',normType=None,
absoluteTolerance=False,perLine=False,skipLines=[]):
import numpy as np import numpy as np
logging.info('\n '.join(['comparing ASCII Tables',file0,file1])) logging.info('\n '.join(['comparing ASCII Tables',file0,file1]))
@ -337,7 +339,7 @@ class Test():
data = [[] for i in range(dataLength)] data = [[] for i in range(dataLength)]
maxError = [0.0 for i in range(dataLength)] maxError = [0.0 for i in range(dataLength)]
absTol = [absoluteTolerance for i in range(dataLength)] absTol = [absoluteTolerance for i in range(dataLength)]
column = [[1 for i in range(dataLength)] for j in range(2)] column = [[1 for i in range(dataLength)] for j in range(2)]
norm = [[] for i in range(dataLength)] norm = [[] for i in range(dataLength)]
normLength = [1 for i in range(dataLength)] normLength = [1 for i in range(dataLength)]
@ -368,11 +370,11 @@ class Test():
key1 = ('1_' if length[i]>1 else '') + headings1[i]['label'] key1 = ('1_' if length[i]>1 else '') + headings1[i]['label']
normKey = ('1_' if normLength[i]>1 else '') + normHeadings[i]['label'] normKey = ('1_' if normLength[i]>1 else '') + normHeadings[i]['label']
if key0 not in table0.labels(raw = True): if key0 not in table0.labels(raw = True):
raise Exception('column {} not found in 1. table...\n'.format(key0)) raise Exception('column "{}" not found in first table...\n'.format(key0))
elif key1 not in table1.labels(raw = True): elif key1 not in table1.labels(raw = True):
raise Exception('column {} not found in 2. table...\n'.format(key1)) raise Exception('column "{}" not found in second table...\n'.format(key1))
elif normKey not in table0.labels(raw = True): elif normKey not in table0.labels(raw = True):
raise Exception('column {} not found in 1. table...\n'.format(normKey)) raise Exception('column "{}" not found in first table...\n'.format(normKey))
else: else:
column[0][i] = table0.label_index(key0) column[0][i] = table0.label_index(key0)
column[1][i] = table1.label_index(key1) column[1][i] = table1.label_index(key1)
@ -400,9 +402,9 @@ class Test():
norm[i] = [1.0 for j in range(line0-len(skipLines))] norm[i] = [1.0 for j in range(line0-len(skipLines))]
absTol[i] = True absTol[i] = True
if perLine: if perLine:
logging.warning('At least one norm of {} in 1. table is 0.0, using absolute tolerance'.format(headings0[i]['label'])) logging.warning('At least one norm of "{}" in first table is 0.0, using absolute tolerance'.format(headings0[i]['label']))
else: else:
logging.warning('Maximum norm of {} in 1. table is 0.0, using absolute tolerance'.format(headings0[i]['label'])) logging.warning('Maximum norm of "{}" in first table is 0.0, using absolute tolerance'.format(headings0[i]['label']))
line1 = 0 line1 = 0
while table1.data_read(): # read next data line of ASCII table while table1.data_read(): # read next data line of ASCII table
@ -414,7 +416,7 @@ class Test():
norm[i][line1-len(skipLines)]) norm[i][line1-len(skipLines)])
line1 +=1 line1 +=1
if (line0 != line1): raise Exception('found {} lines in 1. table but {} in 2. table'.format(line0,line1)) if (line0 != line1): raise Exception('found {} lines in first table but {} in second table'.format(line0,line1))
logging.info(' ********') logging.info(' ********')
for i in range(dataLength): for i in range(dataLength):
@ -561,25 +563,28 @@ class Test():
return allclose return allclose
def compare_TableRefCur(self,headingsRef,ref,headingsCur='',cur='',normHeadings='',normType=None,\ def compare_TableRefCur(self,headingsRef,ref,headingsCur='',cur='',
absoluteTolerance=False,perLine=False,skipLines=[]): normHeadings='',normType=None,
absoluteTolerance=False,perLine=False,skipLines=[]):
if cur == '': cur = ref return self.compare_Table(headingsRef,
if headingsCur == '': headingsCur = headingsRef self.fileInReference(ref),
refName = self.fileInReference(ref) headingsRef if headingsCur == '' else headingsCur,
curName = self.fileInCurrent(cur) self.fileInCurrent(ref if cur == '' else cur),
return self.compare_Table(headingsRef,refName,headingsCur,curName,normHeadings,normType, normHeadings,normType,
absoluteTolerance,perLine,skipLines) absoluteTolerance,perLine,skipLines)
def compare_TableCurCur(self,headingsCur0,Cur0,Cur1,headingsCur1='',normHeadings='',normType=None,\ def compare_TableCurCur(self,headingsCur0,Cur0,Cur1,
absoluteTolerance=False,perLine=False,skipLines=[]): headingsCur1='',
normHeadings='',normType=None,
absoluteTolerance=False,perLine=False,skipLines=[]):
if headingsCur1 == '': headingsCur1 = headingsCur0 return self.compare_Table(headingsCur0,
cur0Name = self.fileInCurrent(Cur0) self.fileInCurrent(Cur0),
cur1Name = self.fileInCurrent(Cur1) headingsCur0 if headingsCur1 == '' else headingsCur1,
return self.compare_Table(headingsCur0,cur0Name,headingsCur1,cur1Name,normHeadings,normType, self.fileInCurrent(Cur1),
absoluteTolerance,perLine,skipLines) normHeadings,normType,absoluteTolerance,perLine,skipLines)
def report_Success(self,culprit): def report_Success(self,culprit):

6
src/IO.f90
View File

@ -1268,7 +1268,7 @@ integer(pInt) function IO_countNumericalDataLines(fileUnit)
line = IO_read(fileUnit) line = IO_read(fileUnit)
chunkPos = IO_stringPos(line) chunkPos = IO_stringPos(line)
tmp = IO_lc(IO_stringValue(line,chunkPos,1_pInt)) tmp = IO_lc(IO_stringValue(line,chunkPos,1_pInt))
if (scan(tmp,"abcdefghijklmnopqrstuvwxyz")/=0) then ! found keyword if (verify(trim(tmp) ,"0123456789")/=0) then ! found keyword
line = IO_read(fileUnit, .true.) ! reset IO_read line = IO_read(fileUnit, .true.) ! reset IO_read
exit exit
else else
@ -1839,12 +1839,12 @@ integer(pInt) function IO_verifyIntValue (string,validChars,myName)
invalidWhere = verify(string,validChars) invalidWhere = verify(string,validChars)
if (invalidWhere == 0_pInt) then if (invalidWhere == 0_pInt) then
read(UNIT=string,iostat=readStatus,FMT=*) IO_verifyIntValue ! no offending chars found read(UNIT=string,iostat=readStatus,FMT=*) IO_verifyIntValue ! no offending chars found
if (readStatus /= 0_pInt) & ! error during string to float conversion if (readStatus /= 0_pInt) & ! error during string to integer conversion
call IO_warning(203_pInt,ext_msg=myName//'"'//string//'"') call IO_warning(203_pInt,ext_msg=myName//'"'//string//'"')
else else
call IO_warning(202_pInt,ext_msg=myName//'"'//string//'"') ! complain about offending characters call IO_warning(202_pInt,ext_msg=myName//'"'//string//'"') ! complain about offending characters
read(UNIT=string(1_pInt:invalidWhere-1_pInt),iostat=readStatus,FMT=*) IO_verifyIntValue ! interpret remaining string read(UNIT=string(1_pInt:invalidWhere-1_pInt),iostat=readStatus,FMT=*) IO_verifyIntValue ! interpret remaining string
if (readStatus /= 0_pInt) & ! error during string to float conversion if (readStatus /= 0_pInt) & ! error during string to integer conversion
call IO_warning(203_pInt,ext_msg=myName//'"'//string(1_pInt:invalidWhere-1_pInt)//'"') call IO_warning(203_pInt,ext_msg=myName//'"'//string(1_pInt:invalidWhere-1_pInt)//'"')
endif endif

64
src/mesh.f90
View File

@ -1642,11 +1642,12 @@ subroutine mesh_marc_get_matNumber(fileUnit)
if (len(trim(line))/=0_pInt) then if (len(trim(line))/=0_pInt) then
chunkPos = IO_stringPos(line) chunkPos = IO_stringPos(line)
data_blocks = IO_intValue(line,chunkPos,1_pInt) data_blocks = IO_intValue(line,chunkPos,1_pInt)
endif endif
allocate(Marc_matNumber(data_blocks))
do i=1_pInt,data_blocks ! read all data blocks do i=1_pInt,data_blocks ! read all data blocks
read (fileUnit,610,END=620) line read (fileUnit,610,END=620) line
chunkPos = IO_stringPos(line) chunkPos = IO_stringPos(line)
Marc_matNumber = (/Marc_matNumber, IO_intValue(line,chunkPos,1_pInt)/) Marc_matNumber(i) = IO_intValue(line,chunkPos,1_pInt)
do j=1_pint,2_pInt + hypoelasticTableStyle ! read 2 or 3 remaining lines of data block do j=1_pint,2_pInt + hypoelasticTableStyle ! read 2 or 3 remaining lines of data block
read (fileUnit,610,END=620) line read (fileUnit,610,END=620) line
enddo enddo
@ -1809,12 +1810,11 @@ subroutine mesh_marc_count_cpElements(fileUnit)
do i=1_pInt,3_pInt+hypoelasticTableStyle ! Skip 3 or 4 lines do i=1_pInt,3_pInt+hypoelasticTableStyle ! Skip 3 or 4 lines
read (fileUnit,610,END=620) line read (fileUnit,610,END=620) line
enddo enddo
mesh_NcpElems = mesh_NcpElems + IO_countContinuousIntValues(fileUnit) ! why not simply mesh_NcpElems = IO_countContinuousIntValues(fileUnit)? keyword hypoelastic might appear several times so the next line probably should not be there mesh_NcpElems = mesh_NcpElems + IO_countContinuousIntValues(fileUnit) ! why not simply mesh_NcpElems = IO_countContinuousIntValues(fileUnit)? not fully correct as hypoelastic can have multiple data fields, needs update
exit exit
endif endif
enddo enddo
else ! Marc2017 and later else ! Marc2017 and later
call IO_error(error_ID=701_pInt)
do do
read (fileUnit,610,END=620) line read (fileUnit,610,END=620) line
chunkPos = IO_stringPos(line) chunkPos = IO_stringPos(line)
@ -1839,6 +1839,7 @@ subroutine mesh_marc_map_elements(fileUnit)
use math, only: math_qsort use math, only: math_qsort
use IO, only: IO_lc, & use IO, only: IO_lc, &
IO_intValue, &
IO_stringValue, & IO_stringValue, &
IO_stringPos, & IO_stringPos, &
IO_continuousIntValues IO_continuousIntValues
@ -1847,7 +1848,8 @@ subroutine mesh_marc_map_elements(fileUnit)
integer(pInt), intent(in) :: fileUnit integer(pInt), intent(in) :: fileUnit
integer(pInt), allocatable, dimension(:) :: chunkPos integer(pInt), allocatable, dimension(:) :: chunkPos
character(len=300) :: line character(len=300) :: line, &
tmp
integer(pInt), dimension (1_pInt+mesh_NcpElems) :: contInts integer(pInt), dimension (1_pInt+mesh_NcpElems) :: contInts
integer(pInt) :: i,cpElem = 0_pInt integer(pInt) :: i,cpElem = 0_pInt
@ -1856,25 +1858,47 @@ subroutine mesh_marc_map_elements(fileUnit)
610 FORMAT(A300) 610 FORMAT(A300)
contInts = 0_pInt
rewind(fileUnit) rewind(fileUnit)
do do
read (fileUnit,610,END=660) line read (fileUnit,610,END=660) line
chunkPos = IO_stringPos(line) chunkPos = IO_stringPos(line)
if( IO_lc(IO_stringValue(line,chunkPos,1_pInt)) == 'hypoelastic' ) then if (MarcVersion < 13) then ! Marc 2016 or earlier
do i=1_pInt,3_pInt+hypoelasticTableStyle ! skip three (or four if new table style!) lines if( IO_lc(IO_stringValue(line,chunkPos,1_pInt)) == 'hypoelastic' ) then
read (fileUnit,610,END=660) line do i=1_pInt,3_pInt+hypoelasticTableStyle ! skip three (or four if new table style!) lines
enddo read (fileUnit,610,END=660) line
contInts = IO_continuousIntValues(fileUnit,mesh_NcpElems,mesh_nameElemSet,& enddo
contInts = IO_continuousIntValues(fileUnit,mesh_NcpElems,mesh_nameElemSet,&
mesh_mapElemSet,mesh_NelemSets) mesh_mapElemSet,mesh_NelemSets)
do i = 1_pInt,contInts(1) exit
cpElem = cpElem+1_pInt endif
mesh_mapFEtoCPelem(1,cpElem) = contInts(1_pInt+i) else ! Marc2017 and later
mesh_mapFEtoCPelem(2,cpElem) = cpElem if ( IO_lc(IO_stringValue(line,chunkPos,1_pInt)) == 'connectivity') then
enddo read (fileUnit,610,END=660) line
endif chunkPos = IO_stringPos(line)
enddo if(any(Marc_matNumber==IO_intValue(line,chunkPos,6_pInt))) then
do
660 call math_qsort(mesh_mapFEtoCPelem,1_pInt,int(size(mesh_mapFEtoCPelem,2_pInt),pInt)) ! should be mesh_NcpElems read (fileUnit,610,END=660) line
chunkPos = IO_stringPos(line)
tmp = IO_lc(IO_stringValue(line,chunkPos,1_pInt))
if (verify(trim(tmp),"0123456789")/=0) then ! found keyword
exit
else
contInts(1) = contInts(1) + 1_pInt
read (tmp,*) contInts(contInts(1)+1)
endif
enddo
endif
endif
endif
enddo
660 do i = 1_pInt,contInts(1)
cpElem = cpElem+1_pInt
mesh_mapFEtoCPelem(1,cpElem) = contInts(1_pInt+i)
mesh_mapFEtoCPelem(2,cpElem) = cpElem
enddo
call math_qsort(mesh_mapFEtoCPelem,1_pInt,int(size(mesh_mapFEtoCPelem,2_pInt),pInt)) ! should be mesh_NcpElems
end subroutine mesh_marc_map_elements end subroutine mesh_marc_map_elements