5       subroutine mtx_colm(a_colm,a_rowp,  argument
9       integer,intent(in) :: a_colm(m+1),b_colm(l+1)  local
11       integer,intent(in) :: a_rowp(a_colm(m+1)-1)
21          do k=a_colm(i),a_colm(i+1)-1
36       subroutine mtx_rowp(a_colm,a_rowp,b_colm,b_rowp,  argument
39       integer,intent(in) :: a_colm(m+1),b_colm(l+1)  local
42       integer,intent(in) :: a_rowp(a_colm(m+1)-1)
57          do k=a_colm(i),a_colm(i+1)-1
72       subroutine mtx_mtx(a_colm,a_rowp,a_mtx,  argument
77       integer,intent(in) :: a_colm(m+1),b_colm(l+1)  local
80       integer,intent(in) :: a_rowp(a_colm(m+1)-1)
83       real(kind=8),intent(in) :: a_mtx(a_colm(m+1)-1)
96          do k=a_colm(i),a_colm(i+1)-1