Use array setter in place of loop for particle coordinate lattice mapping

This commit is contained in:
Will Boyd 2015-05-26 23:28:53 -04:00
parent c93b8a666f
commit c8e614346d

View file

@ -266,10 +266,7 @@ contains
if (.not. associated(p % coord % mapping)) then
allocate(p % coord % mapping(n_maps))
end if
do j = 1, n_maps
p % coord % mapping(j) = lat % offset(j, i_xyz(1), i_xyz(2), i_xyz(3))
end do
p % coord % mapping(:) = lat % offset(:, i_xyz(1), i_xyz(2), i_xyz(3))
end if
call find_cell(p, found)