2 ! PR 60834 - this used to ICE.
10 subroutine solve(this
, x
)
11 class(t
), intent(in
) :: this
12 real, intent(in
) :: x(:)
14 integer, parameter :: n(1:5)=[(i
,i
=1, 5)]
16 associate( nu
=>this
%diffusion
)
17 associate( exponential
=>exp(-(x(i
)-n
) ))