Skip to content

Commit

Permalink
FDS Source: bit of clean up for hypre
Browse files Browse the repository at this point in the history
  • Loading branch information
rmcdermo committed Oct 2, 2024
1 parent 51576ef commit 32bc865
Showing 1 changed file with 1 addition and 11 deletions.
12 changes: 1 addition & 11 deletions Source/pres.f90
Original file line number Diff line number Diff line change
Expand Up @@ -1712,18 +1712,8 @@ SUBROUTINE ULMAT_SOLVE_ZONE(NM,IPZ)
#endif

#ifdef WITH_HYPRE /*---------------------------------------------------------------------------------------*/

! Set values
! ZM%X_H = 0._EB
CALL HYPRE_IJVECTORSETVALUES(ZM%HYPRE_ZM%F_H, ZM%NUNKH, ZM%HYPRE_ZM%INDICES, ZM%F_H, HYPRE_IERR)
!CALL HYPRE_IJVECTORSETVALUES(ZM%HYPRE_ZM%X_H, ZM%NUNKH, ZM%HYPRE_ZM%INDICES, ZM%X_H, HYPRE_IERR)
! Assemble vectors
CALL HYPRE_IJVECTORASSEMBLE(ZM%HYPRE_ZM%F_H, HYPRE_IERR)
!CALL HYPRE_IJVECTORASSEMBLE(ZM%HYPRE_ZM%X_H, HYPRE_IERR)
! Get rhs and soln objects
!CALL HYPRE_IJVECTORGETOBJECT(ZM%HYPRE_ZM%F_H, ZM%HYPRE_ZM%PAR_F_H, HYPRE_IERR)
!CALL HYPRE_IJVECTORGETOBJECT(ZM%HYPRE_ZM%X_H, ZM%HYPRE_ZM%PAR_X_H, HYPRE_IERR)
! Solve
CALL HYPRE_PARCSRPCGSOLVE(ZM%HYPRE_ZM%SOLVER, ZM%HYPRE_ZM%PARCSR_A_H, ZM%HYPRE_ZM%PAR_F_H, ZM%HYPRE_ZM%PAR_X_H, HYPRE_IERR)
IF (CHECK_POISSON .AND. HYPRE_SOLVER_SETPRINTLEVEL>0) THEN
CALL HYPRE_PARCSRPCGGETNUMITERATIONS(ZM%HYPRE_ZM%SOLVER, ZM%HYPRE_ZM%NUM_ITERATIONS, HYPRE_IERR)
Expand Down Expand Up @@ -2796,7 +2786,7 @@ SUBROUTINE ULMAT_H_MATRIX_SOLVER_SETUP(NM,IPZ)
! Allocate Solution and RHS vectors:
IF (ALLOCATED(ZM%X_H)) DEALLOCATE(ZM%X_H)
IF (ALLOCATED(ZM%F_H)) DEALLOCATE(ZM%F_H)
ALLOCATE( ZM%X_H(ZM%NUNKH) , ZM%F_H(ZM%NUNKH) ); ZM%F_H(:) = 0._EB; ZM%X_H(:) = 0.
ALLOCATE( ZM%X_H(ZM%NUNKH) , ZM%F_H(ZM%NUNKH) ); ZM%F_H(:) = 0._EB; ZM%X_H(:) = 0._EB

#ifdef WITH_HYPRE /* --------------------------------------------------------------------------------------------------*/

Expand Down

0 comments on commit 32bc865

Please sign in to comment.