JSOC/proj/vfisv/apps/svbksb.f90 - view

Return to svbksb.f90 CVS log

Up to [Development] / JSOC / proj / vfisv / apps

File: [Development] / JSOC / proj / vfisv / apps / svbksb.f90 (download)
Revision: 1.5, Tue Apr 10 21:17:28 2012 UTC (11 years, 1 month ago) by keiji
Branch: MAIN
CVS Tags: Ver_LATEST, Ver_9-5, Ver_9-41, Ver_9-4, Ver_9-3, Ver_9-2, Ver_9-1, Ver_9-0, Ver_8-8, Ver_8-7, Ver_8-6, Ver_8-5, Ver_8-4, Ver_8-3, Ver_8-2, Ver_8-12, Ver_8-11, Ver_8-10, Ver_8-1, Ver_8-0, Ver_7-1, Ver_7-0, Ver_6-4, Ver_6-3, Ver_6-2, HEAD
Changes since 1.4: +0 -0 lines

*** empty log message ***

      SUBROUTINE SVBKSB(U,W,V,M,N,MP,NP,B,X)

      USE CONS_PARAM
      Real(dp), Dimension(:), Allocatable :: TMP
      Integer :: Status
      real(dp) U(MP,NP),W(NP),V(NP,NP),B(MP),X(NP)
      
      Allocate (TMP(N), Stat=Status)
      If (Status .ne. 0) then
        Print *,'Unable to allocate memory for TMP in svbksb.f90'
        Stop
      End if

      DO 12 J=1,N

        S=0.

        IF(W(J).NE.0.)THEN

          DO 11 I=1,M

            S=S+U(I,J)*B(I)

11        CONTINUE

          S=S/W(J)

        ENDIF

        TMP(J)=S

12    CONTINUE

      DO 14 J=1,N

        S=0.

        DO 13 JJ=1,N

          S=S+V(J,JJ)*TMP(JJ)

13      CONTINUE

        X(J)=S

14    CONTINUE

      Deallocate (TMP)

      RETURN

      END

!CVSVERSIONINFO "$Id: svbksb.f90,v 1.5 2012/04/10 22:17:28 keiji Exp $"

Karen Tian