1DM Object: Distributed Mesh 5 MPI processes 2 type: plex 3Distributed Mesh in 2 dimensions: 4 0-cells: 29 29 30 15 15 5 1-cells: 60 60 60 30 30 6 2-cells: 32 32 32 16 16 7Labels: 8 celltype: 3 strata with value/size (1 (60), 3 (32), 0 (29)) 9 depth: 3 strata with value/size (0 (29), 1 (60), 2 (32)) 10 marker: 2 strata with value/size (1 (8), 4 (17)) 11 Face Sets: 2 strata with value/size (1 (7), 4 (14)) 12 markerBottom: 1 strata with value/size (1 (8)) 13 markerRight: 0 strata with value/size () 14 markerTop: 0 strata with value/size () 15 markerLeft: 1 strata with value/size (1 (17)) 16Integral of pressure: 0. 17SNES Object: 5 MPI processes 18 type: newtonls 19 tolerances: relative=1e-08, absolute=1e-50, solution=1e-08 20 norm schedule ALWAYS 21 SNESLineSearch Object: 5 MPI processes 22 type: bt 23 interpolation: cubic 24 alpha=1.000000e-04 25 maxstep=1.000000e+08, minlambda=1.000000e-12 26 tolerances: relative=1.000000e-08, absolute=1.000000e-15, lambda=1.000000e-08 27 maximum iterations=40 28 KSP Object: 5 MPI processes 29 type: fetidp 30 fully redundant: 0 31 saddle point: 1 32 Inner KSP solver details 33 KSP Object: (fetidp_) 5 MPI processes 34 type: cg 35 maximum iterations=10000, initial guess is zero 36 tolerances: relative=1e-12, absolute=1e-10, divergence=10000. 37 left preconditioning 38 using NATURAL norm type for convergence test 39 PC Object: (fetidp_) 5 MPI processes 40 type: fieldsplit 41 FieldSplit with Schur preconditioner, factorization DIAG 42 Preconditioner for the Schur complement formed from user provided matrix 43 Split info: 44 Split number 0 Defined by IS 45 Split number 1 Defined by IS 46 KSP solver for A00 block 47 KSP Object: (fetidp_fieldsplit_lag_) 5 MPI processes 48 type: preonly 49 maximum iterations=10000, initial guess is zero 50 tolerances: relative=1e-05, absolute=1e-50, divergence=10000. 51 left preconditioning 52 using NONE norm type for convergence test 53 PC Object: (fetidp_fieldsplit_lag_) 5 MPI processes 54 type: shell 55 FETI-DP multipliers 56 Dirichlet preconditioner (just from rank 0) 57 Mat Object: 1 MPI processes 58 type: schurcomplement 59 rows=59, cols=59 60 Schur complement A11 - A10 inv(A00) A01 61 A11 62 Mat Object: 1 MPI processes 63 type: seqaij 64 rows=59, cols=59 65 total: nonzeros=651, allocated nonzeros=651 66 total number of mallocs used during MatSetValues calls=0 67 using I-node routines: found 25 nodes, limit used is 5 68 A10 69 Mat Object: 1 MPI processes 70 type: seqaij 71 rows=59, cols=121 72 total: nonzeros=551, allocated nonzeros=551 73 total number of mallocs used during MatSetValues calls=0 74 using I-node routines: found 25 nodes, limit used is 5 75 KSP of A00 76 KSP Object: (fetidp_bddc_pc_bddc_dirichlet_) 1 MPI processes 77 type: preonly 78 maximum iterations=10000, initial guess is zero 79 tolerances: relative=1e-05, absolute=1e-50, divergence=10000. 80 left preconditioning 81 using NONE norm type for convergence test 82 PC Object: (fetidp_bddc_pc_bddc_dirichlet_) 1 MPI processes 83 type: svd 84 All singular values smaller than 1e-12 treated as zero 85 Provided essential rank of the matrix 0 (all other eigenvalues are zeroed) 86 linear system matrix = precond matrix: 87 Mat Object: (fetidp_bddc_pc_bddc_dirichlet_) 1 MPI processes 88 type: seqaij 89 rows=121, cols=121 90 total: nonzeros=2225, allocated nonzeros=2225 91 total number of mallocs used during MatSetValues calls=0 92 using I-node routines: found 64 nodes, limit used is 5 93 A01 94 Mat Object: 1 MPI processes 95 type: seqaij 96 rows=121, cols=59 97 total: nonzeros=551, allocated nonzeros=551 98 total number of mallocs used during MatSetValues calls=0 99 using I-node routines: found 53 nodes, limit used is 5 100 linear system matrix = precond matrix: 101 Mat Object: (fetidp_fieldsplit_lag_) 5 MPI processes 102 type: submatrix 103 rows=112, cols=112 104 KSP solver for S = A11 - A10 inv(A00) A01 105 KSP Object: (fetidp_fieldsplit_p_) 5 MPI processes 106 type: preonly 107 maximum iterations=10000, initial guess is zero 108 tolerances: relative=1e-05, absolute=1e-50, divergence=10000. 109 left preconditioning 110 using NONE norm type for convergence test 111 PC Object: (fetidp_fieldsplit_p_) 5 MPI processes 112 type: jacobi 113 linear system matrix followed by preconditioner matrix: 114 Mat Object: (fetidp_fieldsplit_p_) 5 MPI processes 115 type: schurcomplement 116 rows=31, cols=31 117 Schur complement A11 - A10 inv(A00) A01 118 A11 119 Mat Object: (fetidp_fieldsplit_p_) 5 MPI processes 120 type: submatrix 121 rows=31, cols=31 122 A10 123 Mat Object: 5 MPI processes 124 type: submatrix 125 rows=31, cols=112 126 KSP of A00 127 KSP Object: (fetidp_fieldsplit_lag_) 5 MPI processes 128 type: preonly 129 maximum iterations=10000, initial guess is zero 130 tolerances: relative=1e-05, absolute=1e-50, divergence=10000. 131 left preconditioning 132 using NONE norm type for convergence test 133 PC Object: (fetidp_fieldsplit_lag_) 5 MPI processes 134 type: shell 135 FETI-DP multipliers 136 Dirichlet preconditioner (just from rank 0) 137 Mat Object: 1 MPI processes 138 type: schurcomplement 139 rows=59, cols=59 140 Schur complement A11 - A10 inv(A00) A01 141 A11 142 Mat Object: 1 MPI processes 143 type: seqaij 144 rows=59, cols=59 145 total: nonzeros=651, allocated nonzeros=651 146 total number of mallocs used during MatSetValues calls=0 147 using I-node routines: found 25 nodes, limit used is 5 148 A10 149 Mat Object: 1 MPI processes 150 type: seqaij 151 rows=59, cols=121 152 total: nonzeros=551, allocated nonzeros=551 153 total number of mallocs used during MatSetValues calls=0 154 using I-node routines: found 25 nodes, limit used is 5 155 KSP of A00 156 KSP Object: (fetidp_bddc_pc_bddc_dirichlet_) 1 MPI processes 157 type: preonly 158 maximum iterations=10000, initial guess is zero 159 tolerances: relative=1e-05, absolute=1e-50, divergence=10000. 160 left preconditioning 161 using NONE norm type for convergence test 162 PC Object: (fetidp_bddc_pc_bddc_dirichlet_) 1 MPI processes 163 type: svd 164 All singular values smaller than 1e-12 treated as zero 165 Provided essential rank of the matrix 0 (all other eigenvalues are zeroed) 166 linear system matrix = precond matrix: 167 Mat Object: (fetidp_bddc_pc_bddc_dirichlet_) 1 MPI processes 168 type: seqaij 169 rows=121, cols=121 170 total: nonzeros=2225, allocated nonzeros=2225 171 total number of mallocs used during MatSetValues calls=0 172 using I-node routines: found 64 nodes, limit used is 5 173 A01 174 Mat Object: 1 MPI processes 175 type: seqaij 176 rows=121, cols=59 177 total: nonzeros=551, allocated nonzeros=551 178 total number of mallocs used during MatSetValues calls=0 179 using I-node routines: found 53 nodes, limit used is 5 180 linear system matrix = precond matrix: 181 Mat Object: (fetidp_fieldsplit_lag_) 5 MPI processes 182 type: submatrix 183 rows=112, cols=112 184 A01 185 Mat Object: 5 MPI processes 186 type: submatrix 187 rows=112, cols=31 188 Mat Object: 5 MPI processes 189 type: is 190 rows=31, cols=31 191 total: nonzeros=216, allocated nonzeros=216 192 total number of mallocs used during MatSetValues calls=0 193 linear system matrix = precond matrix: 194 Mat Object: G (fetidp_) 5 MPI processes 195 type: shell 196 rows=143, cols=143 197 Inner BDDC solver details 198 PC Object: (fetidp_bddc_) 5 MPI processes 199 type: bddc 200 Use verbose output: 0 201 Use user-defined CSR: 0 202 Use local mat graph: 1 203 Connectivity graph topological dimension: 2 204 Graph max count: 2 205 Use vertices: 1 (vertex size 3) 206 Use edges: 1 207 Use faces: 0 208 Use true near null space: 0 209 Use QR for single constraints on cc: 0 210 Use change of basis on local edge nodes: 0 211 Use change of basis on local face nodes: 0 212 User defined change of basis matrix: 0 213 Has change of basis matrix: 0 214 Eliminate dirichlet boundary dofs: 1 215 Switch on static condensation ops around the interface preconditioner: 1 216 Use exact dirichlet trick: 0 217 Interface extension: DIRICHLET 218 Multilevel max levels: 0 219 Multilevel coarsening ratio: 8 220 Use estimated eigs for coarse problem: 0 221 Use deluxe scaling: 0 222 Use deluxe zerorows: 0 223 Use deluxe singlemat: 0 224 Rebuild interface graph for Schur principal minors: 0 225 Number of dofs' layers for the computation of principal minors: -1 226 Use user CSR graph to compute successive layers: 0 227 Adaptive constraint selection threshold (active 0, userdefined 0): 0. 228 Min constraints / connected component: 0 229 Max constraints / connected component: 0 230 Invert exact Schur complement for adaptive selection: 0 231 Symmetric computation of primal basis functions: 1 232 Num. Procs. to map coarse adjacency list: 0 233 Coarse eqs per proc (significant at the coarsest level): 1 234 Detect disconnected: 1 (filter 0) 235 Benign subspace trick: 0 (change explicit 0) 236 Benign subspace trick is active: 0 237 Algebraic computation of no-net-flux: 0 238 ********************************** STATISTICS AT LEVEL 0 ********************************** 239 Global dofs sizes: all 591 interface 149 coarse 6 240 Coarsening ratios: all/coarse 98 interface/coarse 24 241 Active processes : 5 242 Total subdomains : 8 243 Dofs type : MIN MAX MEAN 244 Interior dofs : 45 121 88 245 Interface dofs : 40 95 62 246 Primal dofs : 3 5 4 247 Local dofs : 95 195 151 248 Local subs : 1 2 249 --- Interior solver (rank 0) 250 KSP Object: (fetidp_bddc_pc_bddc_dirichlet_) 1 MPI processes 251 type: preonly 252 maximum iterations=10000, initial guess is zero 253 tolerances: relative=1e-05, absolute=1e-50, divergence=10000. 254 left preconditioning 255 using NONE norm type for convergence test 256 PC Object: (fetidp_bddc_pc_bddc_dirichlet_) 1 MPI processes 257 type: svd 258 All singular values smaller than 1e-12 treated as zero 259 Provided essential rank of the matrix 0 (all other eigenvalues are zeroed) 260 linear system matrix = precond matrix: 261 Mat Object: (fetidp_bddc_pc_bddc_dirichlet_) 1 MPI processes 262 type: seqaij 263 rows=121, cols=121 264 total: nonzeros=2225, allocated nonzeros=2225 265 total number of mallocs used during MatSetValues calls=0 266 using I-node routines: found 64 nodes, limit used is 5 267 --- Correction solver (rank 0) 268 KSP Object: (fetidp_bddc_pc_bddc_neumann_) 1 MPI processes 269 type: preonly 270 maximum iterations=1000, initial guess is zero 271 tolerances: relative=1e-10, absolute=1e-10, divergence=10000. 272 left preconditioning 273 using NONE norm type for convergence test 274 PC Object: (fetidp_bddc_pc_bddc_neumann_) 1 MPI processes 275 type: svd 276 All singular values smaller than 1e-12 treated as zero 277 Provided essential rank of the matrix 0 (all other eigenvalues are zeroed) 278 linear system matrix = precond matrix: 279 Mat Object: (fetidp_bddc_pc_bddc_neumann_) 1 MPI processes 280 type: seqaij 281 rows=176, cols=176 282 total: nonzeros=3856, allocated nonzeros=3856 283 total number of mallocs used during MatSetValues calls=0 284 using I-node routines: found 85 nodes, limit used is 5 285 --- Coarse solver 286 KSP Object: (fetidp_bddc_pc_bddc_coarse_) 5 MPI processes 287 type: preonly 288 maximum iterations=1000, initial guess is zero 289 tolerances: relative=1e-10, absolute=1e-10, divergence=10000. 290 left preconditioning 291 using NONE norm type for convergence test 292 PC Object: (fetidp_bddc_pc_bddc_coarse_) 5 MPI processes 293 type: redundant 294 First (color=0) of 5 PCs follows 295 linear system matrix = precond matrix: 296 Mat Object: (fetidp_bddc_pc_bddc_coarse_) 5 MPI processes 297 type: mpiaij 298 rows=6, cols=6 299 total: nonzeros=34, allocated nonzeros=36 300 total number of mallocs used during MatSetValues calls=0 301 using I-node (on process 0) routines: found 1 nodes, limit used is 5 302 KSP Object: (fetidp_bddc_pc_bddc_coarse_redundant_) 1 MPI processes 303 type: preonly 304 maximum iterations=10000, initial guess is zero 305 tolerances: relative=1e-05, absolute=1e-50, divergence=10000. 306 left preconditioning 307 using NONE norm type for convergence test 308 PC Object: (fetidp_bddc_pc_bddc_coarse_redundant_) 1 MPI processes 309 type: cholesky 310 out-of-place factorization 311 tolerance for zero pivot 2.22045e-14 312 matrix ordering: nd 313 factor fill ratio given 5., needed 1. 314 Factored matrix follows: 315 Mat Object: 1 MPI processes 316 type: seqsbaij 317 rows=6, cols=6 318 package used to perform factorization: petsc 319 total: nonzeros=20, allocated nonzeros=20 320 block size is 1 321 linear system matrix = precond matrix: 322 Mat Object: 1 MPI processes 323 type: seqaij 324 rows=6, cols=6 325 total: nonzeros=34, allocated nonzeros=34 326 total number of mallocs used during MatSetValues calls=0 327 using I-node routines: found 3 nodes, limit used is 5 328 linear system matrix = precond matrix: 329 Mat Object: 5 MPI processes 330 type: is 331 rows=591, cols=591 332 total: nonzeros=17047, allocated nonzeros=17047 333 total number of mallocs used during MatSetValues calls=0 334 has attached near null space 335 maximum iterations=10000, initial guess is zero 336 tolerances: relative=1e-12, absolute=1e-10, divergence=10000. 337 left preconditioning 338 using NONE norm type for convergence test 339 PC Object: 5 MPI processes 340 type: none 341 linear system matrix followed by preconditioner matrix: 342 Mat Object: 5 MPI processes 343 type: is 344 rows=591, cols=591 345 total: nonzeros=17047, allocated nonzeros=20789 346 total number of mallocs used during MatSetValues calls=0 347 has attached null space 348 Mat Object: 5 MPI processes 349 type: is 350 rows=591, cols=591 351 total: nonzeros=17047, allocated nonzeros=20789 352 total number of mallocs used during MatSetValues calls=0 353L_2 Error: 0.00151 [9.32e-06, 0.00151] 354Integral of pressure: 0. 355