41       SUBROUTINE GMRES(nFaces, dof, nNo, nnz, mynNo, commu, cS, face,   argument
46       INTEGER, INTENT(IN) :: nFaces, dof, nNo, nnz, mynNo  local
52       REAL*8, INTENT(IN) :: Val(dof*dof,nnz), R(dof,nNo)
53       REAL*8, INTENT(OUT) :: X(dof,nNo)
61       ALLOCATE(h(ls%sD+1,ls%sD), u(dof,nNo,ls%sD+1), unCondU(dof,nNo), 
73             CALL SPARMULVV(dof, nNo, nnz, commu, cS, rowPtr, colPtr, 
76             CALL ADDBCMUL(BCOP_TYPE_ADD, nFaces, dof, nNo, mynNo, commu,
83             CALL ADDBCMUL(BCOP_TYPE_PRE, nFaces, dof, nNo, mynNo, commu,
87          err(1)   = NORMV(dof, mynNo, commu, u(:,:,1))
103             CALL SPARMULVV(dof, nNo, nnz, commu, cS, rowPtr, colPtr, 
106             CALL ADDBCMUL(BCOP_TYPE_ADD, nFaces, dof, nNo, mynNo, commu,
111                CALL ADDBCMUL(BCOP_TYPE_PRE, nFaces, dof, nNo, mynNo, 
115                h(j,i) = DOTV(dof, mynNo, commu, u(:,:,i+1), u(:,:,j))
118             h(i+1,i)   = NORMV(dof, mynNo, commu, u(:,:,i+1))
163       SUBROUTINE GMRESV(nFaces, dof, nNo, nnz, mynNo, commu, cS, face,   argument
168       INTEGER, INTENT(IN) :: nFaces, dof, nNo, nnz, mynNo  local
174       REAL*8, INTENT(IN) :: Val(dof*dof,nnz)
175       REAL*8, INTENT(INOUT) :: R(dof,nNo)
183       ALLOCATE(h(ls%sD+1,ls%sD), u(dof,nNo,ls%sD+1), X(dof,nNo), 
188       eps       = NORMV(dof, mynNo, commu, R)
203          CALL SPARMULVV(dof, nNo, nnz, commu, cS, rowPtr, colPtr, 
205          CALL ADDBCMUL(BCOP_TYPE_ADD, nFaces, dof, nNo, mynNo, commu,
209          err(1)   = NORMV(dof, mynNo, commu, u(:,:,1))
213             CALL SPARMULVV(dof, nNo, nnz, commu, cS, rowPtr, colPtr, 
215             CALL ADDBCMUL(BCOP_TYPE_ADD, nFaces, dof, nNo, mynNo, commu,
219                h(j,i) = DOTV(dof, mynNo, commu, u(:,:,i+1), u(:,:,j))
222             h(i+1,i)   = NORMV(dof, mynNo, commu, u(:,:,i+1))