* PACKAGE LGLU1  !" 線型解析 : LU 分解による連立方程式の解
*
*"  [HIS] 93/07/10(takepiro) 「行列計算ソフトウェア」使用(GLU1)
*
**********************************************************************
      SUBROUTINE LGLU
     M            ( AMAT   , 
     M              BVEC   ,
     O              IERR   , 
     I              MATDIM , OLU   )
*
#ifdef SYS_IBMS
      INCLUDE   (ZCDIM)                      !" 
      INCLUDE   (ZLDIM)                      !" 格子点数, 波数
#else
#include        "zcdim.F"                    !" 
#include        "zldim.F"                    !" 係数行列の大きさ
#endif
*
*   [MODIFY] 
      REAL      AMAT  ( MATDMX, MATDMX )     !" 係数行列
      REAL      BVEC  ( MATDMX )             !" 右辺ベクトル / 解
*
*   [OUTPUT] 
      INTEGER   IERR                         !" サブルーチンエラーコード
*
*   [INPUT] 
      INTEGER   MATDIM                       !" 計算する次元
      LOGICAL   OLU                          !" LU 分解するか否か
*
*   [INTERNAL WORK] 
      INTEGER   IP    ( MATDMX )             !" ピボット情報
      REAL      WORK  ( MATDMX )             !" 作業配列
      REAL      EPS                          !" 特異性の判定基準
      DATA      EPS  / 3.52D-15 /
*
      INTEGER   IUNIT
      LOGICAL   OFIRST 
      DATA      OFIRST  / .TRUE. /
*
      SAVE      EPS, OFIRST
*
      IF ( OFIRST )THEN
         OFIRST = .FALSE.
         CALL GLPGET( 'MSGUNIT', IUNIT )
         WRITE ( IUNIT, * ) ' GAUSSIAN ELIMINATION METHOD : 93/07/10'
      ENDIF
*
      IERR = 0
*
      IF ( OLU ) THEN
         CALL GLU1
     M          ( AMAT   , 
     I            MATDIM , MATDMX , 
     O            EPS    , 
     W            WORK   , 
     O            IP     , IERR     )
      ENDIF
*
      IF ( IERR .EQ. 0 ) THEN
         CALL GSLV1
     I          ( AMAT   , 
     I            MATDIM , MATDMX , 
     M            BVEC   , 
     I            IP                )
      ELSE
         CALL MSGDMP
     &         ( 'W', 'LGLU', 'SOL. NOT CALCULATED BECAUSE OF ERRORS' )
      END IF
*
      RETURN
      END