rtsol Subroutine

  1. util.f90
  2. rtsol

subroutine rtsol(n, nn, nmax, R, b)

solves Rt.x=b with same conventions as above nn is not required on entry but is set on exit

Arguments

Type IntentOptional AttributesName
integer :: n
integer :: nn
integer :: nmax
real :: R
real :: b

Called by

proc~~rtsol~~CalledByGraph proc~rtsol rtsol proc~eptsol eptsol proc~eptsol->proc~rtsol proc~checkt checkT proc~checkt->proc~rtsol proc~updatese updateSE proc~updatese->proc~eptsol

Contents

Source Code

rtsol

Source Code

      subroutine rtsol(n,nn,nmax,R,b)
      implicit double precision (a-h,o-z)
      dimension R(*),b(*)
      n2=nmax+2
      nn=1
      b(1)=b(1)/R(1)
      do i=2,n
        i1=i-1
        call mysaxpy(-b(i1),R(nn+1),b(i),n-i1)
        nn=nn+n2-i
        b(i)=b(i)/R(nn)
      end do
      end subroutine rtsol