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