forms b=M'.x where Q is nxn, stored by columns, with stride nmax
Type | Intent | Optional | Attributes | Name | ||
---|---|---|---|---|---|---|
integer | :: | n | ||||
integer | :: | nmax | ||||
real | :: | Q | ||||
real | :: | x | ||||
real | :: | b |
subroutine Qtprod(n,nmax,Q,x,b)
implicit double precision (a-h,o-z)
dimension Q(*),x(*),b(*)
i1=1
do i=1,n
b(i)=scpr(0.D0,Q(i1),x,n)
i1=i1+nmax
end do
end subroutine Qtprod