SetDefaults.F

#include "parmdefines.h"
c**********************************************************************
#include "author.inc"
c*    $Id: SetDefaults.F,v 1.16 1996/10/09 02:13:03 turner Exp $
c*
c*    Sets reasonable default values for iparm and rparm.
c*
c*    <PARAMETER LIST>
c*
c*     Output:
c*      iparm - array of integer parameters
c*      rparm - array of floating point parameters
c*      status - return status
c*
c*    <SUBROUTINES REQUIRED>
c*
c*     JT_FillVectorFloat
c*     JT_FillVectorInt
c*
#include "parmdesc.inc"
c*
#include "copyright.inc"
c**********************************************************************
      subroutine JT_SetDefaults (iparm, rparm, status)
      implicit none
c
c ... Parameters.
#include "numbers.inc"
c
c ... Output:
      integer status
      integer iparm(_JT_no_of_iparms_)
      real rparm(_JT_no_of_rparms_)
c
      call JT_FillVectorInt (_JT_no_of_iparms_, 0, iparm, status)
      iparm(_JT_out_) = _JT_out_warnings_  ! by default output errors and warnings
      iparm(_JT_luout_) = _JT_stdout_   ! logical unit for output
#ifdef hpux
      iparm(_JT_luerr_) = _JT_stderr_hp_   ! logical unit for errors
#else
      iparm(_JT_luerr_) = _JT_stderr_   ! logical unit for errors
#endif
      iparm(_JT_norm_) = _JT_norm_2_  ! use L-2 norm
      iparm(_JT_steps_) = _JT_steps_default_   ! steps in m-step preconditioners
      iparm(_JT_nold_) = _JT_nold_default_   ! vectors to use in truncated methods
c
      call JT_FillVectorFloat (_JT_no_of_rparms_, zero, rparm, status)
      rparm(_JT_eps_) = _JT_eps_default_
      rparm(_JT_epspre_) = _JT_epspre_default_
      rparm(_JT_omega_) = _JT_omega_default_
      rparm(_JT_tiny_) = _JT_tiny_default_
c
      status = 0
      return
      end