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