DumpSystem.F

#include "iadefines.h"
c**********************************************************************
#include "author.inc"
c*    $Id: DumpSystem.F,v 1.3 1996/07/08 00:47:36 turner Exp $
c*
c*    Subroutine to dump an entire linear system.
c*
c*    <PARAMETER LIST>
c*
c*     Input:
c*      lu - logical unit to which to write
c*      a - coefficient matrix
c*      ia - integer vector containing info about how "a" is stored
c*        NOTE: see description of ia below
c*      ja - map array
c*      x - vector of unknowns
c*      b - vector containing right-hand side
c*
c*     Output:
c*      status - return status
c*        -1  ==>  invalid argument(s)
c*         0  ==>  success
c*
#include "iadesc.inc"
c*
c*    <SUBROUTINES REQUIRED>
c*
c*     JT_DumpCoeff
c*
#include "copyright.inc"
c**********************************************************************
      subroutine JT_DumpSystem (lu, a, ia, ja, x, b, status)
      implicit none
c
c ... Input:
      integer lu
      integer ia(_JT_no_of_storage_parameters_)
      integer ja(*)
      real a(*), x(*), b(*)
c
c ... Output:
      integer status
c
c ... Local:
      integer i
c
      call JT_DumpCoeff (lu, a, ia, ja, status)
      do i=1,ia(2)
       write(lu,100) x(i), b(i)
      enddo
c
      return
  100 format (1p, 2(1x, g15.8))
      end