Loading src/LinearSolverTypes.f90 +4 −0 Original line number Diff line number Diff line Loading @@ -577,7 +577,11 @@ SUBROUTINE init_LinearSolverType_Base(solver,Params,A) IF(TPLType==PETSC) THEN #ifdef FUTILITY_HAVE_PETSC !create and initialize KSP #ifdef HAVE_MPI CALL KSPCreate(solver%MPIparallelEnv%comm,solver%ksp,iperr) #else CALL KSPCreate(PETSC_COMM_SELF,solver%ksp,iperr) #endif !set iterative solver type SELECTCASE(solverMethod) Loading src/LinearSolverTypes_Multigrid.f90 +3 −2 Original line number Diff line number Diff line Loading @@ -35,6 +35,7 @@ MODULE LinearSolverTypes_Multigrid USE IntrType USE ParameterLists USE ParallelEnv USE VectorTypes USE MatrixTypes USE LinearSolverTypes Loading Loading @@ -177,7 +178,7 @@ SUBROUTINE init_LinearSolverType_Multigrid(solver,Params,A) !Pull LS data from the parameter list TPLType=-1 MPI_Comm_ID=-1 MPI_Comm_ID=PE_COMM_NULL matType=-1 matEngine=-1 numberOMP=-1 Loading @@ -202,7 +203,7 @@ SUBROUTINE init_LinearSolverType_Multigrid(solver,Params,A) CALL validParams%clear() !Initialize parallel environments based on input IF(MPI_Comm_ID /= -1) CALL solver%MPIparallelEnv%init(MPI_Comm_ID) IF(MPI_Comm_ID /= PE_COMM_NULL) CALL solver%MPIparallelEnv%init(MPI_Comm_ID) IF(numberOMP > 0) CALL solver%OMPparallelEnv%init(numberOMP) IF(TPLType /= PETSC) THEN Loading src/MatrixTypes.f90 +1 −1 Original line number Diff line number Diff line Loading @@ -882,7 +882,7 @@ SUBROUTINE pop_recv(acc,valBuf,thisMat,ctBuf,idx,req,f) !> Array of buffer sizes INTEGER(SIK), INTENT(INOUT) :: ctBuf(MATVEC_SLOTS) !> The buffer slot popped INTEGER(SIK), INTENT(OUT) :: idx INTEGER(SIK), INTENT(INOUT) :: idx !> List of MPI Requests INTEGER(SIK), INTENT(INOUT) :: req(MATVEC_SLOTS) !> Optional argument to flush all buffers Loading src/MatrixTypes_PETSc.f90 +4 −0 Original line number Diff line number Diff line Loading @@ -121,7 +121,11 @@ SUBROUTINE init_PETScMatrixParam(matrix,Params) IF(validParams%has('MatrixType->isSym')) & CALL validParams%get('MatrixType->isSym',isSym) CALL validParams%get('MatrixType->matType',matType) #ifdef HAVE_MPI CALL validParams%get('MatrixType->MPI_COMM_ID',MPI_COMM_ID) #else MPI_COMM_ID = PETSC_COMM_SELF #endif CALL validParams%get('MatrixType->nlocal',nlocal) m=n Loading src/ParameterLists.f90 +2 −2 Original line number Diff line number Diff line Loading @@ -172,8 +172,8 @@ INTEGER(SIK),PARAMETER :: VALIDTYPE_VERIFYTEST=1 INTEGER(SIK),PARAMETER :: VALIDTYPE_VERIFYLIST=2 !> Hashing parameters INTEGER(SLK),PARAMETER,PRIVATE :: PL_HASH_POLYNOMIAL_BASE=97_SIK INTEGER(SLK),PARAMETER,PRIVATE :: PL_HASH_MODULUS=1000000009_SIK INTEGER(SLK),PARAMETER,PRIVATE :: PL_HASH_POLYNOMIAL_BASE=97_SLK INTEGER(SLK),PARAMETER,PRIVATE :: PL_HASH_MODULUS=1000000009_SLK !> Exception handler for the module TYPE(ExceptionHandlerType),SAVE :: eParams Loading Loading
src/LinearSolverTypes.f90 +4 −0 Original line number Diff line number Diff line Loading @@ -577,7 +577,11 @@ SUBROUTINE init_LinearSolverType_Base(solver,Params,A) IF(TPLType==PETSC) THEN #ifdef FUTILITY_HAVE_PETSC !create and initialize KSP #ifdef HAVE_MPI CALL KSPCreate(solver%MPIparallelEnv%comm,solver%ksp,iperr) #else CALL KSPCreate(PETSC_COMM_SELF,solver%ksp,iperr) #endif !set iterative solver type SELECTCASE(solverMethod) Loading
src/LinearSolverTypes_Multigrid.f90 +3 −2 Original line number Diff line number Diff line Loading @@ -35,6 +35,7 @@ MODULE LinearSolverTypes_Multigrid USE IntrType USE ParameterLists USE ParallelEnv USE VectorTypes USE MatrixTypes USE LinearSolverTypes Loading Loading @@ -177,7 +178,7 @@ SUBROUTINE init_LinearSolverType_Multigrid(solver,Params,A) !Pull LS data from the parameter list TPLType=-1 MPI_Comm_ID=-1 MPI_Comm_ID=PE_COMM_NULL matType=-1 matEngine=-1 numberOMP=-1 Loading @@ -202,7 +203,7 @@ SUBROUTINE init_LinearSolverType_Multigrid(solver,Params,A) CALL validParams%clear() !Initialize parallel environments based on input IF(MPI_Comm_ID /= -1) CALL solver%MPIparallelEnv%init(MPI_Comm_ID) IF(MPI_Comm_ID /= PE_COMM_NULL) CALL solver%MPIparallelEnv%init(MPI_Comm_ID) IF(numberOMP > 0) CALL solver%OMPparallelEnv%init(numberOMP) IF(TPLType /= PETSC) THEN Loading
src/MatrixTypes.f90 +1 −1 Original line number Diff line number Diff line Loading @@ -882,7 +882,7 @@ SUBROUTINE pop_recv(acc,valBuf,thisMat,ctBuf,idx,req,f) !> Array of buffer sizes INTEGER(SIK), INTENT(INOUT) :: ctBuf(MATVEC_SLOTS) !> The buffer slot popped INTEGER(SIK), INTENT(OUT) :: idx INTEGER(SIK), INTENT(INOUT) :: idx !> List of MPI Requests INTEGER(SIK), INTENT(INOUT) :: req(MATVEC_SLOTS) !> Optional argument to flush all buffers Loading
src/MatrixTypes_PETSc.f90 +4 −0 Original line number Diff line number Diff line Loading @@ -121,7 +121,11 @@ SUBROUTINE init_PETScMatrixParam(matrix,Params) IF(validParams%has('MatrixType->isSym')) & CALL validParams%get('MatrixType->isSym',isSym) CALL validParams%get('MatrixType->matType',matType) #ifdef HAVE_MPI CALL validParams%get('MatrixType->MPI_COMM_ID',MPI_COMM_ID) #else MPI_COMM_ID = PETSC_COMM_SELF #endif CALL validParams%get('MatrixType->nlocal',nlocal) m=n Loading
src/ParameterLists.f90 +2 −2 Original line number Diff line number Diff line Loading @@ -172,8 +172,8 @@ INTEGER(SIK),PARAMETER :: VALIDTYPE_VERIFYTEST=1 INTEGER(SIK),PARAMETER :: VALIDTYPE_VERIFYLIST=2 !> Hashing parameters INTEGER(SLK),PARAMETER,PRIVATE :: PL_HASH_POLYNOMIAL_BASE=97_SIK INTEGER(SLK),PARAMETER,PRIVATE :: PL_HASH_MODULUS=1000000009_SIK INTEGER(SLK),PARAMETER,PRIVATE :: PL_HASH_POLYNOMIAL_BASE=97_SLK INTEGER(SLK),PARAMETER,PRIVATE :: PL_HASH_MODULUS=1000000009_SLK !> Exception handler for the module TYPE(ExceptionHandlerType),SAVE :: eParams Loading