remove trailing white spaces

Author: pmocz

mtx/public/mtx_bcyclic_dble_decsol.dek

@@ -1,4 +1,4 @@
-   
+
       subroutine bcyclic_dble_decsolblk( &
             iop,caller_id,nvar,nz,lblk,dblk,ublk,brhs,ipiv,lrd,rpar_decsol,lid,ipar_decsol,ierr)
          use mtx_def, only: lapack
@@ -8,7 +8,7 @@
          real(dp), dimension(:), pointer, intent(inout) :: lblk, dblk, ublk ! =(nvar,nvar,nz)
          ! row(i) of mtx has lblk(:,:,i), dblk(:,:,i), ublk(:,:,i)
          ! lblk(:,:,1) is not used; ublk(:,:,nz) is not used.
-         real(dp), pointer, intent(inout)  :: brhs(:) ! =(nvar,nz)     
+         real(dp), pointer, intent(inout)  :: brhs(:) ! =(nvar,nz)
          real(dp), pointer, intent(inout) :: rpar_decsol(:) ! (lrd)
          integer, pointer, intent(inout) :: ipar_decsol(:) ! (lid)
          integer, intent(out) :: ierr
@@ -35,7 +35,7 @@
          end if
       end subroutine bcyclic_dble_decsolblk
 
-      
+
       subroutine bcyclic_dble_work_sizes(nvar,nz,lrd,lid)
          integer, intent(in) :: nvar,nz
          integer, intent(out) :: lrd,lid

mtx/public/mtx_dble_decsol.dek

@@ -1,4 +1,4 @@
-   
+
       subroutine lapack_decsol(iop,n,ndim,a,ml,mu,b,ip,lrd,rpar_decsol,lid,ipar_decsol,ier)
          use const_def, only: dp
          integer, intent(in) :: iop, n, ndim, lrd, lid
@@ -37,7 +37,7 @@
          end if
       end subroutine lapack_decsol
 
-      
+
       subroutine lapack_work_sizes(n,lrd,lid)
          integer, intent(in) :: n
          integer, intent(out) :: lrd,lid
@@ -46,10 +46,10 @@
       end subroutine lapack_work_sizes
 
 
-      
+
       ! complex(dp) routines
 
-      
+
       subroutine lapack_zwork_sizes(n,mu,ml,lcd,lrd,lid)
          integer, intent(in) :: n,mu,ml
          integer, intent(out) :: lcd,lrd,lid
@@ -63,7 +63,7 @@
          lrd = 0
          lcd = lda*(n+1) + n
       end subroutine lapack_zwork_sizes
-      
+
 
 
       subroutine lapack_decsolc(iop,n,ndim,ar,ai,ml,mu,br,bi,ip,lcd,cpar_decsol,lrd,rpar_decsol,lid,ipar_decsol,ier)
@@ -112,20 +112,20 @@
          real(dp), intent(inout), pointer :: rpar_decsol(:) ! (lrd)
          integer, intent(inout), pointer :: ipar_decsol(:) ! (lid)
          integer, intent(out) :: ier
-         
+
          ! lcd must be >= n*ndim
-         
+
          complex(dp), pointer :: a(:,:)
          integer :: i, j
-         
+
          i=1
          a(1:ndim,1:n) => cpar_decsol(i:i+n*ndim-1)
          i=i+n*ndim
          if (i-1 > lcd) then
             ier = -1
             return
          end if
-         
+
          do j=1,n
             do i=1,ndim
                a(i,j) = cmplx(ar(i,j),ai(i,j), kind=dp)
@@ -138,11 +138,11 @@
                ai(i,j) = AIMAG(a(i,j))
             end do
          end do
-         
+
       end subroutine lapack_decc
-         
-      
-      
+
+
+
       subroutine lapack_decbc(n,ndim,ar,ai,ml,mu,ip,lcd,cpar_decsol,lrd,rpar_decsol,lid,ipar_decsol,ier)
          use const_def, only: dp
          integer, intent(in) :: n, ndim, ml, mu, lcd, lrd, lid
@@ -152,20 +152,20 @@
          real(dp), intent(inout), pointer :: rpar_decsol(:) ! (lrd)
          integer, intent(inout), pointer :: ipar_decsol(:) ! (lid)
          integer, intent(out) :: ier
-         
+
          ! ldc must be >= n*ndim
-         
+
          complex(dp), pointer :: a(:,:)
          integer :: i, j
-         
+
          i=1
          a(1:ndim,1:n) => cpar_decsol(i:i+n*ndim-1)
          i=i+n*ndim
          if (i-1 > lcd) then
             ier = -1
             return
          end if
-         
+
          do j=1,n
             do i=1,ndim
                a(i,j) = cmplx(ar(i,j),ai(i,j), kind=dp)
@@ -178,10 +178,10 @@
                ai(i,j) = AIMAG(a(i,j))
             end do
          end do
-         
+
       end subroutine lapack_decbc
-      
-      
+
+
       subroutine lapack_solc(n,ndim,ar,ai,br,bi,ip,lcd,cpar_decsol,lrd,rpar_decsol,lid,ipar_decsol,ier)
          integer, intent(in) :: n, ndim, lcd, lrd, lid
          integer, intent(in) :: ip(n)
@@ -191,12 +191,12 @@
          real(dp), intent(inout), pointer :: rpar_decsol(:) ! (lrd)
          integer, intent(inout), pointer :: ipar_decsol(:) ! (lid)
          integer, intent(out) :: ier
-         
+
          ! ldc must be >= n*(ndim+1)
-         
+
          complex(dp), pointer :: a(:,:), b(:)
          integer :: i, j
-         
+
          i = 1
          a(1:ndim,1:n) => cpar_decsol(i:i+n*ndim-1)
          i=i+n*ndim
@@ -205,7 +205,7 @@
             ier = -1
             return
          end if
-         
+
          do j=1,n
             do i=1,ndim
                a(i,j) = cmplx(ar(i,j),ai(i,j), kind=dp)
@@ -232,12 +232,12 @@
          real(dp), intent(inout), pointer :: rpar_decsol(:) ! (lrd)
          integer, intent(inout), pointer :: ipar_decsol(:) ! (lid)
          integer, intent(out) :: ier
-         
+
          ! ldc must be >= n*(ndim+1)
-         
+
          complex(dp), pointer :: a(:,:), b(:)
          integer :: i, j
-         
+
          i = 1
          a(1:ndim,1:n) => cpar_decsol(i:i+n*ndim-1)
          i=i+n*ndim
@@ -246,7 +246,7 @@
             ier = -1
             return
          end if
-         
+
          do j=1,n
             do i=1,ndim
                a(i,j) = cmplx(ar(i,j),ai(i,j), kind=dp)
@@ -260,8 +260,8 @@
                bi(j) = AIMAG(b(j))
             end do
          end do
-         
-         
+
+
       end subroutine lapack_solbc
-      
-      
+
+

mtx/public/mtx_debug_decsol.dek

@@ -1,4 +1,4 @@
-   
+
 
       subroutine debug_decsol(iop,n,ndim,a,ml,mu,b,ip,lrd,rpar,lid,ipar,ier)
          integer, intent(in) :: iop, n, ndim, lrd, lid, ml, mu
@@ -27,7 +27,7 @@
             end if
          end if
       end subroutine debug_decsol
-      
+
 
       subroutine debug_decsolc(iop,n,ndim,ar,ai,ml,mu,br,bi,ip,lcd,cpar_decsol,lrd,rpar_decsol,lid,ipar_decsol,ier)
          integer, intent(in) :: iop, n, ndim, lcd, lrd, lid, ml, mu

mtx/public/mtx_decsolblk_dble.dek

@@ -4,7 +4,7 @@
          integer, intent(in) :: caller_id, nvar, nz, iop, lrd, lid
          integer, pointer, intent(inout) :: ipiv(:) ! =(nvar,nz)
          real(dp), dimension(:), pointer, intent(inout) :: lblk, dblk, ublk ! =(nvar,nvar,nz)
-         real(dp), pointer, intent(inout)  :: brhs(:) ! =(nvar,nz)     
+         real(dp), pointer, intent(inout)  :: brhs(:) ! =(nvar,nz)
          real(dp), pointer, intent(inout) :: rpar_decsol(:) ! (lrd)
          integer, pointer, intent(inout) :: ipar_decsol(:) ! (lid)
          integer, intent(out) :: ierr