1 #if !defined(__pcbddc_h) 2 #define __pcbddc_h 3 4 #include <../src/ksp/pc/impls/is/pcis.h> 5 #include <../src/ksp/pc/impls/bddc/bddcstructs.h> 6 7 #if !defined(PETSC_PCBDDC_MAXLEVELS) 8 #define PETSC_PCBDDC_MAXLEVELS 8 9 #endif 10 11 PETSC_EXTERN PetscLogEvent PC_BDDC_Topology[PETSC_PCBDDC_MAXLEVELS]; 12 PETSC_EXTERN PetscLogEvent PC_BDDC_LocalSolvers[PETSC_PCBDDC_MAXLEVELS]; 13 PETSC_EXTERN PetscLogEvent PC_BDDC_LocalWork[PETSC_PCBDDC_MAXLEVELS]; 14 PETSC_EXTERN PetscLogEvent PC_BDDC_CorrectionSetUp[PETSC_PCBDDC_MAXLEVELS]; 15 PETSC_EXTERN PetscLogEvent PC_BDDC_CoarseSetUp[PETSC_PCBDDC_MAXLEVELS]; 16 PETSC_EXTERN PetscLogEvent PC_BDDC_ApproxSetUp[PETSC_PCBDDC_MAXLEVELS]; 17 PETSC_EXTERN PetscLogEvent PC_BDDC_ApproxApply[PETSC_PCBDDC_MAXLEVELS]; 18 PETSC_EXTERN PetscLogEvent PC_BDDC_CoarseSolver[PETSC_PCBDDC_MAXLEVELS]; 19 PETSC_EXTERN PetscLogEvent PC_BDDC_AdaptiveSetUp[PETSC_PCBDDC_MAXLEVELS]; 20 PETSC_EXTERN PetscLogEvent PC_BDDC_Scaling[PETSC_PCBDDC_MAXLEVELS]; 21 PETSC_EXTERN PetscLogEvent PC_BDDC_Schurs[PETSC_PCBDDC_MAXLEVELS]; 22 23 /* Private context (data structure) for the BDDC preconditioner. */ 24 typedef struct { 25 /* First MUST come the folowing line, for the stuff that is common to FETI and Neumann-Neumann. */ 26 PC_IS pcis; 27 /* Coarse stuffs needed by BDDC application in KSP */ 28 Vec coarse_vec; 29 KSP coarse_ksp; 30 Mat coarse_phi_B; 31 Mat coarse_phi_D; 32 Mat coarse_psi_B; 33 Mat coarse_psi_D; 34 PetscInt local_primal_size; 35 PetscInt coarse_size; 36 PetscInt* global_primal_indices; 37 VecScatter coarse_loc_to_glob; 38 /* Local stuffs needed by BDDC application in KSP */ 39 Vec vec1_P; 40 Vec vec1_C; 41 Mat local_auxmat1; 42 Mat local_auxmat2; 43 Vec vec1_R; 44 Vec vec2_R; 45 IS is_R_local; 46 VecScatter R_to_B; 47 VecScatter R_to_D; 48 KSP ksp_R; 49 KSP ksp_D; 50 /* Quantities defining constraining details (local) of the preconditioner */ 51 /* These quantities define the preconditioner itself */ 52 PetscInt n_vertices; 53 Mat ConstraintMatrix; 54 PetscBool new_primal_space; 55 PetscBool new_primal_space_local; 56 PetscInt *primal_indices_local_idxs; 57 PetscInt local_primal_size_cc; 58 PetscInt *local_primal_ref_node; 59 PetscInt *local_primal_ref_mult; 60 PetscBool use_change_of_basis; 61 PetscBool use_change_on_faces; 62 PetscBool fake_change; 63 Mat ChangeOfBasisMatrix; 64 Mat user_ChangeOfBasisMatrix; 65 PetscBool change_interior; 66 Mat switch_static_change; 67 Vec work_change; 68 Vec original_rhs; 69 Vec temp_solution; 70 Mat local_mat; 71 PetscBool use_exact_dirichlet_trick; 72 PetscBool exact_dirichlet_trick_app; 73 PetscBool ksp_guess_nonzero; 74 PetscBool rhs_change; 75 PetscBool temp_solution_used; 76 /* benign subspace trick */ 77 PetscBool benign_saddle_point; 78 PetscBool benign_have_null; 79 PetscBool benign_skip_correction; 80 PetscBool benign_compute_correction; 81 Mat benign_change; 82 Mat benign_original_mat; 83 IS *benign_zerodiag_subs; 84 Vec benign_vec; 85 Mat benign_B0; 86 PetscSF benign_sf; 87 PetscScalar *benign_p0; 88 PetscInt benign_n; 89 PetscInt *benign_p0_lidx; 90 PetscInt *benign_p0_gidx; 91 PetscBool benign_null; 92 PetscBool benign_change_explicit; 93 PetscBool benign_apply_coarse_only; 94 95 /* Some defaults on selecting vertices and constraints*/ 96 PetscBool use_local_adj; 97 PetscBool use_vertices; 98 PetscBool use_faces; 99 PetscBool use_edges; 100 101 /* Some customization is possible */ 102 PetscBool corner_selection; 103 PetscBool corner_selected; 104 PetscBool recompute_topography; 105 PetscBool graphanalyzed; 106 PCBDDCGraph mat_graph; 107 PetscInt graphmaxcount; 108 MatNullSpace onearnullspace; 109 PetscObjectState *onearnullvecs_state; 110 PetscBool NullSpace_corr[4]; 111 IS user_primal_vertices; 112 IS user_primal_vertices_local; 113 PetscBool use_nnsp; 114 PetscBool use_nnsp_true; 115 PetscBool use_qr_single; 116 PetscBool user_provided_isfordofs; 117 PetscInt n_ISForDofs; 118 PetscInt n_ISForDofsLocal; 119 IS *ISForDofs; 120 IS *ISForDofsLocal; 121 IS NeumannBoundaries; 122 IS NeumannBoundariesLocal; 123 IS DirichletBoundaries; 124 IS DirichletBoundariesLocal; 125 PetscBool eliminate_dirdofs; 126 PetscBool switch_static; 127 PetscInt coarsening_ratio; 128 PetscInt coarse_adj_red; 129 PetscInt current_level; 130 PetscInt max_levels; 131 PetscInt coarse_eqs_per_proc; 132 PetscInt coarse_eqs_limit; 133 IS coarse_subassembling; 134 PetscBool use_coarse_estimates; 135 PetscBool symmetric_primal; 136 PetscInt vertex_size; 137 PCBDDCInterfaceExtType interface_extension; 138 139 /* no-net-flux */ 140 PetscBool compute_nonetflux; 141 Mat divudotp; 142 PetscBool divudotp_trans; 143 IS divudotp_vl2l; 144 145 /* nedelec */ 146 Mat discretegradient; 147 PetscInt nedorder; 148 PetscBool conforming; 149 PetscInt nedfield; 150 PetscBool nedglobal; 151 Mat nedcG; 152 IS nedclocal; 153 154 /* local disconnected subdomains */ 155 PetscBool detect_disconnected; 156 PetscBool detect_disconnected_filter; 157 PetscInt n_local_subs; 158 IS *local_subs; 159 160 /* scaling */ 161 Vec work_scaling; 162 PetscBool use_deluxe_scaling; 163 PCBDDCDeluxeScaling deluxe_ctx; 164 PetscBool deluxe_zerorows; 165 PetscBool deluxe_singlemat; 166 167 /* schur complements on interface's subsets */ 168 PCBDDCSubSchurs sub_schurs; 169 PetscBool sub_schurs_rebuild; 170 PetscBool sub_schurs_exact_schur; 171 PetscInt sub_schurs_layers; 172 PetscBool sub_schurs_use_useradj; 173 PetscBool computed_rowadj; 174 175 /* adaptive selection of constraints */ 176 PetscBool adaptive_selection; 177 PetscBool adaptive_userdefined; 178 PetscReal adaptive_threshold[2]; 179 PetscInt adaptive_nmin; 180 PetscInt adaptive_nmax; 181 PetscInt* adaptive_constraints_n; 182 PetscInt* adaptive_constraints_idxs; 183 PetscInt* adaptive_constraints_idxs_ptr; 184 PetscScalar* adaptive_constraints_data; 185 PetscInt* adaptive_constraints_data_ptr; 186 187 /* For verbose output of some bddc data structures */ 188 PetscInt dbg_flag; 189 PetscViewer dbg_viewer; 190 } PC_BDDC; 191 192 #endif /* __pcbddc_h */ 193