Replace -pc_air_max_dd_ratio with new CF splitting type

Makefile

@@ -105,7 +105,10 @@ ifeq ($(PETSC_HAVE_KOKKOS),1)
 export OBJS := $(OBJS) $(SRCDIR)/PETSc_Helperk.o \
 							  $(SRCDIR)/Grid_Transferk.o \
 							  $(SRCDIR)/VecISCopyLocalk.o \
-							  $(SRCDIR)/PMISR_DDCk.o \
+							  $(SRCDIR)/Device_Datak.o \
+							  $(SRCDIR)/MatDiagDomk.o \
+							  $(SRCDIR)/PMISR_Modulek.o \
+							  $(SRCDIR)/DDC_Modulek.o \
 							  $(SRCDIR)/Gmres_Polyk.o
 endif	
 
@@ -116,9 +119,13 @@ OBJS := $(OBJS) $(SRCDIR)/PETSc_Helper.o \
 		  $(SRCDIR)/AIR_MG_Stats.o \
 		  $(SRCDIR)/SAI_Z.o \
 		  $(SRCDIR)/Constrain_Z_or_W.o \
-		  $(SRCDIR)/PMISR_DDC.o \
+		  $(SRCDIR)/MatDiagDom.o \
+		  $(SRCDIR)/SAbs.o \
+		  $(SRCDIR)/DDC_Module.o \
+		  $(SRCDIR)/PMISR_Module.o \
 		  $(SRCDIR)/Aggregation.o \
 		  $(SRCDIR)/CF_Splitting.o \
+		  $(SRCDIR)/MatDiagDomSubmatrix.o \
 		  $(SRCDIR)/Repartition.o \
 		  $(SRCDIR)/Timers.o \
 		  $(SRCDIR)/Weighted_Jacobi.o \

Makefile.deps

@@ -7,15 +7,18 @@ src/AIR_Operators_Setup.o : src/AIR_Operators_Setup.F90 src/PETSc_Helper.o src/S
 src/Approx_Inverse_Setup.o : src/Approx_Inverse_Setup.F90 src/Matshell_Data_Type.o src/PETSc_Helper.o src/Repartition.o src/SAI_Z.o src/Weighted_Jacobi.o src/Neumann_Poly.o src/Gmres_Poly_Newton.o src/Gmres_Poly.o src/Pflare_Parameters.o src/TSQR.o 
 src/Binary_Tree.o : src/Binary_Tree.F90 
 src/C_Fortran_Bindings.o : src/C_Fortran_Bindings.F90 src/AIR_Data_Type_Routines.o src/CF_Splitting.o src/Approx_Inverse_Setup.o src/PCAIR_Shell.o src/PCAIR_Data_Type.o 
-src/CF_Splitting.o : src/CF_Splitting.F90 src/PETSc_Helper.o src/Aggregation.o src/C_PETSc_Interfaces.o src/PMISR_DDC.o src/Pflare_Parameters.o 
+src/CF_Splitting.o : src/CF_Splitting.F90 src/PETSc_Helper.o src/Aggregation.o src/C_PETSc_Interfaces.o src/SAbs.o src/DDC_Module.o src/PMISR_Module.o src/Pflare_Parameters.o 
 src/Constrain_Z_or_W.o : src/Constrain_Z_or_W.F90 src/PETSc_Helper.o src/C_PETSc_Interfaces.o 
 src/C_PETSc_Interfaces.o : src/C_PETSc_Interfaces.F90 
+src/DDC_Module.o : src/DDC_Module.F90 src/MatDiagDom.o src/Pflare_Parameters.o src/PMISR_Module.o src/SAbs.o src/C_PETSc_Interfaces.o src/PETSc_Helper.o 
 src/FC_Smooth.o : src/FC_Smooth.F90 src/Matshell_Data_Type.o src/PETSc_Helper.o src/AIR_Data_Type.o src/C_PETSc_Interfaces.o 
 src/Gmres_Poly_Data_Type.o : src/Gmres_Poly_Data_Type.F90 src/TSQR.o 
 src/Gmres_Poly.o : src/Gmres_Poly.F90 src/PETSc_Helper.o src/Gmres_Poly_Data_Type.o src/TSQR.o src/Matshell_Data_Type.o src/Pflare_Parameters.o src/C_PETSc_Interfaces.o src/Sorting.o 
 src/Gmres_Poly_Newton.o : src/Gmres_Poly_Newton.F90 src/C_PETSc_Interfaces.o src/Gmres_Poly.o 
 src/Grid_Transfer.o : src/Grid_Transfer.F90 src/PETSc_Helper.o src/C_PETSc_Interfaces.o 
 src/Grid_Transfer_Improve.o : src/Grid_Transfer_Improve.F90 src/PETSc_Helper.o src/Pflare_Parameters.o src/Timers.o 
+src/MatDiagDom.o : src/MatDiagDom.F90 src/Pflare_Parameters.o src/C_PETSc_Interfaces.o src/PETSc_Helper.o 
+src/MatDiagDomSubmatrix.o : src/MatDiagDomSubmatrix.F90 src/Pflare_Parameters.o src/PETSc_Helper.o src/CF_Splitting.o 
 src/Matshell_Data_Type.o : src/Matshell_Data_Type.F90 src/Pflare_Parameters.o src/AIR_Data_Type.o 
 src/Neumann_Poly.o : src/Neumann_Poly.F90 src/Pflare_Parameters.o src/TSQR.o src/Matshell_Data_Type.o src/PETSc_Helper.o src/Gmres_Poly.o 
 src/PCAIR_C_Fortran_Bindings.o : src/PCAIR_C_Fortran_Bindings.F90 src/PCAIR_Interfaces.o 
@@ -26,8 +29,9 @@ src/PCPFLAREINV_Interfaces.o : src/PCPFLAREINV_Interfaces.F90
 src/PETSc_Helper.o : src/PETSc_Helper.F90 src/C_PETSc_Interfaces.o 
 src/PFLARE.o : src/PFLARE.F90 src/Pflare_Parameters.o src/PCAIR_Interfaces.o src/PCPFLAREINV_Interfaces.o 
 src/Pflare_Parameters.o : src/Pflare_Parameters.F90 
-src/PMISR_DDC.o : src/PMISR_DDC.F90 src/Pflare_Parameters.o src/C_PETSc_Interfaces.o src/PETSc_Helper.o 
+src/PMISR_Module.o : src/PMISR_Module.F90 src/Pflare_Parameters.o src/C_PETSc_Interfaces.o src/PETSc_Helper.o 
 src/Repartition.o : src/Repartition.F90 src/PETSc_Helper.o src/C_PETSc_Interfaces.o 
+src/SAbs.o : src/SAbs.F90 src/PETSc_Helper.o 
 src/SAI_Z.o : src/SAI_Z.F90 src/Pflare_Parameters.o src/PETSc_Helper.o src/C_PETSc_Interfaces.o src/Sorting.o src/Binary_Tree.o 
 src/Sorting.o : src/Sorting.F90 src/Binary_Tree.o 
 src/Timers.o : src/Timers.F90 src/Pflare_Parameters.o 

README.md

@@ -25,6 +25,7 @@ PFLARE adds new methods to PETSc, including:
 1) Polynomial approximate inverses, e.g., GMRES and Neumann polynomials
 2) Reduction multigrids, e.g., AIRG, nAIR and lAIR
 3) CF splittings, e.g., PMISR DDC
+4) Extracting diagonally dominant submatrices
 
 ## Quick start
 
@@ -46,12 +47,13 @@ For details about PFLARE, please see:
 | [docs/gpus.md](docs/gpus.md) | Using GPUs with PFLARE |
 | [docs/reuse.md](docs/reuse.md) | Re-using components of PFLARE |
 | [docs/options.md](docs/options.md) | List of the options available in PFLARE |   
+| [docs/faq.md](docs/faq.md) | Frequently asked questions and help! |   
 
 and the Jupyter notebooks:
 
 | Path | Contents |
 |---|---|
-| [notebooks/01_getting_started.ipynb](notebooks/01_getting_started.ipynb) | Introduce PFLARE |
+| [notebooks/01_getting_started.ipynb](notebooks/01_getting_started.ipynb) | Introduction to PFLARE |
 | [notebooks/02_pcpflareinv.ipynb](notebooks/02_pcpflareinv.ipynb) | Examine some of the approximate inverses found in PCPFLAREINV |
 | [notebooks/03_cf_splitting.ipynb](notebooks/03_cf_splitting.ipynb) | Visualise the C/F splitting and explore the PMISR-DDC algorithm |
 | [notebooks/04_pcair.ipynb](notebooks/04_pcair.ipynb) | Introduction to PCAIR and the AIRG method |