anchor.diff:        write(iDump, "(I0)") idx
       end do
       ! item
  +    write(iDump,"(A)") '% item(1:nnonzero)'
       do i = 1, n
         ! Lower
         do j = hecMAT%indexL(i-1)+1,hecMAT%indexL(i)
  
positive.code:     ! header
      write(iDump,"(A)") '%%Block-CSR matrix real general'
      nnz = n + hecMAT%indexL(n) + hecMAT%indexU(n)
      write(iDump,"(I0,' ',I0,' ',I0,' ',I0)") n, np, nnz, ndof
      ! index
      idx = 0
      write(iDump, "(I0)") idx
      do i = 1, n
  
negative.code:     ! header
      write(iDump,"(A)") '%%CSR matrix real general'
      nnz = ndof2 * (n + hecMAT%indexL(n) + hecMAT%indexU(n))
      write(iDump,"(A)") '% nrow ncol nnonzero'
      write(iDump,"(I0,' ',I0,' ',I0)") n*ndof, np*ndof, nnz
      ! index
      write(iDump,"(A)") '% index(0:nrow)'
      idx = 0
      write(iDump, "(I0)") idx
      do i = 1, n
  