1 subroutine mv(m
,nc
,irp
,ja
,val
,x
,ldx
,y
,ldy
,acc
)
5 integer, parameter :: ipk_
= int32
6 integer, parameter :: spk_
= real32
7 complex(spk_
), parameter :: czero
=(0.0_spk_
,0.0_spk_
)
9 integer(ipk_
), intent(in
) :: m
,ldx
,ldy
,nc
,irp(*),ja(*)
10 complex(spk_
), intent(in
) :: x(ldx
,*),val(*)
11 complex(spk_
), intent(inout
) :: y(ldy
,*)
12 complex(spk_
), intent(inout
) :: acc(*)
13 integer(ipk_
) :: i
,j
,k
, ir
, jc
17 do j
=irp(i
), irp(i
+1)-1
18 acc(1:nc
) = acc(1:nc
) + val(j
) * x(ja(j
),1:nc
)
20 y(i
,1:nc
) = -acc(1:nc
)