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

exception in function math_expand in module math.f90

This commit is contained in:
Sharan Roongta 2018-05-17 19:46:38 +02:00
parent 8d2c0ec857
commit e107a5e004
1 changed files with 3 additions and 1 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

4
src/math.f90
View File

@ -82,6 +82,7 @@ module math
math_crossproduct, & math_crossproduct, &
math_tensorproduct33, & math_tensorproduct33, &
math_mul3x3, & math_mul3x3, &
math_expand, &
math_mul6x6, & math_mul6x6, &
math_mul33xx33, & math_mul33xx33, &
math_mul3333xx33, & math_mul3333xx33, &
@ -384,7 +385,8 @@ pure function math_expand(what,how)
integer(pInt), dimension(:), intent(in) :: how integer(pInt), dimension(:), intent(in) :: how
real(pReal), dimension(sum(how)) :: math_expand real(pReal), dimension(sum(how)) :: math_expand
integer(pInt) :: i integer(pInt) :: i
if(sum(how)==0) &
return
do i = 1_pInt, size(how) do i = 1_pInt, size(how)
math_expand(sum(how(1:i-1))+1:sum(how(1:i))) = what(mod(i-1_pInt,size(what))+1_pInt) math_expand(sum(how(1:i-1))+1:sum(how(1:i))) = what(mod(i-1_pInt,size(what))+1_pInt)
enddo enddo