FDS Source: add HYPRE_PARCSRMATRIXPRINT to HYPRE interface

Author: rmcdermo

Source/cons.f90

@@ -194,6 +194,7 @@ MODULE GLOBAL_CONSTANTS
 LOGICAL :: ACCUMULATE_WATER=.FALSE.         !< Indicates that integrated liquid outputs are specified
 LOGICAL :: WRITE_XYZ=.FALSE.                !< Indicates that a Plot3D geometry file is specified by user
 LOGICAL :: CHECK_POISSON=.FALSE.            !< Check the accuracy of the Poisson solver
+LOGICAL :: WRITE_PARCSRPCG_MATRIX=.FALSE.   !< If true, write out matrix for UGLMAT HYPRE solver
 LOGICAL :: TWO_D=.FALSE.                    !< Perform a 2-D simulation
 LOGICAL :: SETUP_ONLY=.FALSE.               !< Indicates that the calculation should be stopped before time-stepping
 LOGICAL :: CHECK_MESH_ALIGNMENT=.FALSE.     !< Indicates that the user wants to check the mesh alignment and then stop
@@ -561,7 +562,7 @@ MODULE GLOBAL_CONSTANTS
 ! Logical units and output file names
 
 INTEGER                              :: LU_ERR=ERROR_UNIT,LU_END=2,LU_GIT=3,LU_SMV=4,LU_INPUT=5,LU_OUTPUT=6,LU_STOP=7,LU_CPU=8,&
-                                        LU_CATF=9,LU_RDIR=10,LU_GDIR=11,LU_SETCC=12,LU_BINGEOM=13
+                                        LU_CATF=9,LU_RDIR=10,LU_GDIR=11,LU_SETCC=12,LU_BINGEOM=13,LU_PARCSRPCG_MATRIX=14
 INTEGER                              :: LU_MASS,LU_HRR,LU_STEPS,LU_NOTREADY,LU_VELOCITY_ERROR,LU_CFL,LU_LINE=-1,LU_CUTCELL, &
                                         LU_CVODE_SUBSTEPS
 INTEGER                              :: LU_HISTOGRAM,LU_HVAC

Source/geom.f90

@@ -22175,7 +22175,7 @@ SUBROUTINE READ_GEOM
 
 CHARACTER(LABEL_LENGTH) :: ID,MATL_ID,TEXTURE_MAPPING, &
                            DEVC_ID,CTRL_ID,SURF_IDS(3),SURF_ID6(6),MOVE_ID
-CHARACTER(MESSAGE_LENGTH) :: BUFFER,FN_BINGEOM,BINARY_FILE
+CHARACTER(FN_LENGTH) :: BUFFER,FN_BINGEOM,BINARY_FILE
 CHARACTER(LABEL_LENGTH),  ALLOCATABLE, DIMENSION(:) :: SURF_ID
 REAL(EB), ALLOCATABLE, DIMENSION(:) :: ZVALS,TFACES
 REAL(EB), ALLOCATABLE, TARGET, DIMENSION(:) :: VERTS,VERTS_AUX

Source/imkl.f90

@@ -463,6 +463,13 @@ SUBROUTINE HYPRE_PARCSRPCGGETFINALRELATIVE(SOLVER, FINAL_RES_NORM, IERR)
          INTEGER,         INTENT(OUT)    :: IERR
       END SUBROUTINE HYPRE_PARCSRPCGGETFINALRELATIVE
 
+      SUBROUTINE HYPRE_PARCSRMATRIXPRINT(PARCSR_A, FILE_NAME, FILE_NAME_SIZE, IERR)
+         INTEGER(KIND=8), INTENT(IN)     :: PARCSR_A
+         CHARACTER(*),    INTENT(IN)     :: FILE_NAME
+         INTEGER,         INTENT(IN)     :: FILE_NAME_SIZE
+         INTEGER,         INTENT(OUT)    :: IERR
+      END SUBROUTINE HYPRE_PARCSRMATRIXPRINT
+
    END INTERFACE
 
    ! HYPRE info for ZONE_MESH defined in HYPRE_ZM, ZONE_MESH_TYPE, type.f90.
@@ -530,7 +537,8 @@ END SUBROUTINE HYPRE_PARCSRPCGGETFINALRELATIVE
              HYPRE_PARCSRPCGSETUP,                &  ! subroutine in HYPRE library
              HYPRE_PARCSRPCGSOLVE,                &  ! subroutine in HYPRE library
              HYPRE_PARCSRPCGGETNUMITERATIONS,     &  ! subroutine in HYPRE library
-             HYPRE_PARCSRPCGGETFINALRELATIVE         ! subroutine in HYPRE library
+             HYPRE_PARCSRPCGGETFINALRELATIVE,     &  ! subroutine in HYPRE library
+             HYPRE_PARCSRMATRIXPRINT                 ! subroutine in HYPRE library
 
 END MODULE HYPRE_INTERFACE
 #endif /* WITH_HYPRE */

Source/pres.f90

@@ -3933,6 +3933,7 @@ SUBROUTINE GET_H_MATRIX_LUDCMP
 #ifdef WITH_HYPRE
 INTEGER :: ONEV(1), END_ROW
 REAL(EB) :: ZEROV(1)
+CHARACTER(FN_LENGTH) :: FN_PARCSRPCG_MATRIX
 #endif
 !.. All other variables
 INTEGER :: INNZ, IROW, JCOL, IERR
@@ -4188,6 +4189,10 @@ SUBROUTINE GET_H_MATRIX_LUDCMP
    CALL HYPRE_BOOMERAMGSETMAXITER(ZSL%HYPRE_ZSL%PRECOND, HYPRE_PRECOND_MAXITER, HYPRE_IERR)
    CALL HYPRE_PARCSRPCGSETPRECOND(ZSL%HYPRE_ZSL%SOLVER, HYPRE_PRECOND_ID, ZSL%HYPRE_ZSL%PRECOND, HYPRE_IERR)
    ! Solver setup
+   IF(WRITE_PARCSRPCG_MATRIX) THEN
+      WRITE(FN_PARCSRPCG_MATRIX,'(A,A,A)') './',TRIM(CHID),'_hypre_matrix.txt'
+      CALL HYPRE_PARCSRMATRIXPRINT(ZSL%HYPRE_ZSL%PARCSR_A_H, TRIM(FN_PARCSRPCG_MATRIX), LEN(TRIM(FN_PARCSRPCG_MATRIX)), HYPRE_IERR)
+   ENDIF
    CALL HYPRE_PARCSRPCGSETUP(ZSL%HYPRE_ZSL%SOLVER, ZSL%HYPRE_ZSL%PARCSR_A_H,&
                              ZSL%HYPRE_ZSL%PAR_F_H, ZSL%HYPRE_ZSL%PAR_X_H, HYPRE_IERR)
 #endif

Source/read.f90

@@ -9826,7 +9826,8 @@ SUBROUTINE READ_PRES
 
 NAMELIST /PRES/ BAROCLINIC,CHECK_POISSON,FISHPAK_BC,ITERATION_SUSPEND_FACTOR, &
                 MAX_PRESSURE_ITERATIONS,PRESSURE_RELAX_TIME,PRESSURE_TOLERANCE, &
-                RELAXATION_FACTOR,SOLVER,SUSPEND_PRESSURE_ITERATIONS,TUNNEL_PRECONDITIONER,VELOCITY_TOLERANCE
+                RELAXATION_FACTOR,SOLVER,SUSPEND_PRESSURE_ITERATIONS,TUNNEL_PRECONDITIONER,VELOCITY_TOLERANCE, &
+                WRITE_PARCSRPCG_MATRIX
 
 ! Read the single PRES line