11 integer, parameter :: fp
= selected_real_kind(13)
12 complex(kind
=fp
) :: c(:,:)
13 where( abs( aimag( c
) ) < 1.e
-10_fp
) &
14 & c
= cmplx( real( c
, fp
) , 0._fp
, fp
)
15 where( abs( real( c
, fp
) ) < 1.e
-10_fp
) &
16 & c
= cmplx( 0._fp
, aimag( c
) , fp
)