Qprod Subroutine

  1. util.f90
  2. Qprod

subroutine Qprod(n, nmax, Q, x, b)

forms b=M.x where Q is nxn, stored by columns, with stride nmax

Arguments

Type IntentOptional AttributesName
integer :: n
integer :: nmax
real :: Q
real :: x
real :: b

Called by

proc~~qprod~~CalledByGraph proc~qprod Qprod proc~eptsol eptsol proc~eptsol->proc~qprod proc~updatese updateSE proc~updatese->proc~eptsol

Contents

Source Code

Qprod

Source Code

      subroutine Qprod(n,nmax,Q,x,b)
      implicit double precision (a-h,o-z)
      dimension Q(*),x(*),b(*)
      do i=1,n
        b(i)=0.D0
      end do
      i1=1
      do i=1,n
        call mysaxpy(x(i),Q(i1),b,n)
        i1=i1+nmax
      end do
      end subroutine Qprod