From 056f154a13cd677a181b9a0b3159b42928376da4 Mon Sep 17 00:00:00 2001
From: Salvatore Filippone <salvatore.filippone@cranfield.ac.uk>
Date: Tue, 10 Apr 2007 14:23:46 +0000
Subject: [PATCH] *** empty log message ***

---
 base/serial/csr/dcsrmm.f | 80 ++++++++++++++++++++--------------------
 1 file changed, 40 insertions(+), 40 deletions(-)

diff --git a/base/serial/csr/dcsrmm.f b/base/serial/csr/dcsrmm.f
index 76df9493..1e3d36f8 100644
--- a/base/serial/csr/dcsrmm.f
+++ b/base/serial/csr/dcsrmm.f
@@ -35,7 +35,7 @@ c
 *     
 *     
       SUBROUTINE DCSRMM(TRANSA,M,K,N,ALPHA,DESCRA,AR,
-     *     JA,IA,B,LDB,BETA,C,LDC,WORK,LWORK,IERROR)
+     *  JA,IA,B,LDB,BETA,C,LDC,WORK,LWORK,IERROR)
 C     
 C     
 C     .. Scalar Arguments ..
@@ -66,23 +66,23 @@ C
       IF (DESCRA(1:1).EQ.'S' .AND. DESCRA(2:2).EQ.'U') TRANS = 'U'
       IF (DESCRA(1:1).EQ.'S' .AND. DESCRA(2:2).EQ.'L') TRANS = 'L'
 c     
-c     Does DSRMV manage this case too? <????????????????????????????
+c     Does DSRMV manage this case too? 
 c     
       IF (DESCRA(1:1).EQ.'D') THEN
-         IF (DESCRA(3:3).EQ.'U') THEN
-            DO 40 I = 1, K
-               DO 20 J = 1, M
-                  C(J,I) = BETA*C(J,I) + ALPHA*B(J,I)
- 20            CONTINUE
- 40         CONTINUE
-         ELSE
-            DO 80 I = 1, K
-               DO 60 J = 1, M
-                  C(J,I) = BETA*C(J,I) + ALPHA*AR(J)*B(J,I)
- 60            CONTINUE
- 80         CONTINUE
-         END IF
-         RETURN
+        IF (DESCRA(3:3).EQ.'U') THEN
+          DO 40 I = 1, K
+            DO 20 J = 1, M
+              C(J,I) = BETA*C(J,I) + ALPHA*B(J,I)
+ 20         CONTINUE
+ 40       CONTINUE
+        ELSE
+          DO 80 I = 1, K
+            DO 60 J = 1, M
+              C(J,I) = BETA*C(J,I) + ALPHA*AR(J)*B(J,I)
+ 60         CONTINUE
+ 80       CONTINUE
+        END IF
+        RETURN
       END IF
 C     
       IF (DESCRA(3:3).EQ.'N') DIAG = 'N'
@@ -95,38 +95,38 @@ C
       K4 = MOD(K,NB)
       SELECT CASE(K4)
       CASE(1)
-         IF (K==1) THEN
-            CALL DCSRMV(TRANS,DIAG,M,N,ALPHA,AR,JA,IA,B(1,1),
-     +           BETA,C(1,1),WORK,LWORK,IERROR)
-         ELSE 
-            CALL DCSRMV2(TRANS,DIAG,M,N,ALPHA,AR,JA,IA,
-     +           B(1,1),LDB,BETA,C(1,1),LDC,WORK,LWORK,IERROR)
-            CALL DCSRMV3(TRANS,DIAG,M,N,ALPHA,AR,JA,IA,
-     +           B(1,3),LDB,BETA,C(1,3),LDC,WORK,LWORK,IERROR)
-            K4 = K4 + 4 
-         ENDIF
+        IF (K==1) THEN
+          CALL DCSRMV(TRANS,DIAG,M,N,ALPHA,AR,JA,IA,B(1,1),
+     +      BETA,C(1,1),WORK,LWORK,IERROR)
+        ELSE 
+          CALL DCSRMV2(TRANS,DIAG,M,N,ALPHA,AR,JA,IA,
+     +      B(1,1),LDB,BETA,C(1,1),LDC,WORK,LWORK,IERROR)
+          CALL DCSRMV3(TRANS,DIAG,M,N,ALPHA,AR,JA,IA,
+     +      B(1,3),LDB,BETA,C(1,3),LDC,WORK,LWORK,IERROR)
+          K4 = K4 + 4 
+        ENDIF
       CASE(2)
-         CALL DCSRMV2(TRANS,DIAG,M,N,ALPHA,AR,JA,IA,
-     +        B(1,1),LDB,BETA,C(1,1),LDC,WORK,LWORK,IERROR)
+        CALL DCSRMV2(TRANS,DIAG,M,N,ALPHA,AR,JA,IA,
+     +    B(1,1),LDB,BETA,C(1,1),LDC,WORK,LWORK,IERROR)
       CASE(3)
-         CALL DCSRMV3(TRANS,DIAG,M,N,ALPHA,AR,JA,IA,
-     +        B(1,1),LDB,BETA,C(1,1),LDC,WORK,LWORK,IERROR)
+        CALL DCSRMV3(TRANS,DIAG,M,N,ALPHA,AR,JA,IA,
+     +    B(1,1),LDB,BETA,C(1,1),LDC,WORK,LWORK,IERROR)
       END SELECT
 
       IF(IERROR.NE.0) THEN
-         INT_VAL(1)=IERROR
-         CALL FCPSB_ERRPUSH(4012,NAME,INT_VAL)
-         GOTO 9999
+        INT_VAL(1)=IERROR
+        CALL FCPSB_ERRPUSH(4012,NAME,INT_VAL)
+        GOTO 9999
       END IF
 
       DO I=K4+1,K,NB
-         CALL DCSRMV4(TRANS,DIAG,M,N,ALPHA,AR,JA,IA,
-     +        B(1,I),LDB,BETA,C(1,I),LDC,WORK,LWORK,IERROR)
+        CALL DCSRMV4(TRANS,DIAG,M,N,ALPHA,AR,JA,IA,
+     +    B(1,I),LDB,BETA,C(1,I),LDC,WORK,LWORK,IERROR)
       ENDDO
       IF(IERROR.NE.0) THEN
-         INT_VAL(1)=IERROR
-         CALL FCPSB_ERRPUSH(4012,NAME,INT_VAL)
-         GOTO 9999
+        INT_VAL(1)=IERROR
+        CALL FCPSB_ERRPUSH(4012,NAME,INT_VAL)
+        GOTO 9999
       END IF
 
       CALL FCPSB_ERRACTIONRESTORE(ERR_ACT)
@@ -136,8 +136,8 @@ C
       CALL FCPSB_ERRACTIONRESTORE(ERR_ACT)
 
       IF ( ERR_ACT .NE. 0 ) THEN 
-         CALL FCPSB_SERROR()
-         RETURN
+        CALL FCPSB_SERROR()
+        RETURN
       ENDIF
 
       RETURN