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

replaces Nast with Nreps in hybridIA 

formerly, if the number of orientations requested was (far) less than the number of nonzero bins in the texture representation, you only got a subset of low phi1 values..! (that is the top corner in Euler space -- Hi Bing ;-) )
This commit is contained in:
Denny Tjahjanto 2009-08-13 13:21:22 +00:00
parent 0ef4ec27b3
commit cafc39b183
1 changed files with 4 additions and 7 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

11
code/IO.f90
View File

@ -201,7 +201,7 @@
enddo enddo
enddo enddo
enddo enddo
dV_V = dV_V/sum_dV_V ! normalize to 1 dV_V = dV_V/sum_dV_V ! normalize to 1
!--- now fix bounds --- !--- now fix bounds ---
@ -229,7 +229,7 @@
exit exit
endif endif
enddo enddo
allocate(binSet(Nreps)) allocate(binSet(Nreps))
bin = 0 ! bin counter bin = 0 ! bin counter
i = 1 ! set counter i = 1 ! set counter
@ -247,7 +247,7 @@
do i=1,Nast do i=1,Nast
if (i < Nast) then if (i < Nast) then
call random_number(rnd) call random_number(rnd)
j = nint(rnd*(Nast-i)+i+0.5_pReal,pInt) j = nint(rnd*(Nreps-i)+i+0.5_pReal,pInt)
else else
j = i j = i
endif endif
@ -858,7 +858,7 @@ endfunction
case (265) case (265)
msg = 'Limit for crystallite loop too small' msg = 'Limit for crystallite loop too small'
case (266) case (266)
msg = 'Limit for crystallite state loop too small' msg = 'Limit for state loop too small'
case (267) case (267)
msg = 'Limit for stress loop too small' msg = 'Limit for stress loop too small'
case (268) case (268)
@ -884,9 +884,6 @@ endfunction
case (277) case (277)
msg = 'Non-positive relevant mismatch in RGC' msg = 'Non-positive relevant mismatch in RGC'
case (279)
msg = 'Limit for materialpoint state loop too small'
case (300) case (300)
msg = 'This material can only be used with elements with three direct stress components' msg = 'This material can only be used with elements with three direct stress components'
case (500) case (500)