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: 15 18 15 15 9 5*c4762a1bSJed Brown 1-cells: 22 24 22 22 12 6*c4762a1bSJed Brown 2-cells: 8 8 8 8 4 7*c4762a1bSJed BrownLabels: 8*c4762a1bSJed Brown depth: 3 strata with value/size (0 (15), 1 (22), 2 (8)) 9*c4762a1bSJed Brown celltype: 3 strata with value/size (3 (8), 0 (15), 1 (22)) 10*c4762a1bSJed Brown marker: 2 strata with value/size (1 (9), 4 (5)) 11*c4762a1bSJed Brown Face Sets: 2 strata with value/size (1 (6), 4 (3)) 12*c4762a1bSJed Brown markerBottom: 1 strata with value/size (1 (9)) 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 (5)) 16*c4762a1bSJed BrownIntegral of pressure: -0.0387881 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: shell 41*c4762a1bSJed Brown FETI-DP multipliers 42*c4762a1bSJed Brown Dirichlet preconditioner (just from rank 0) 43*c4762a1bSJed Brown Mat Object: 1 MPI processes 44*c4762a1bSJed Brown type: schurcomplement 45*c4762a1bSJed Brown rows=24, cols=24 46*c4762a1bSJed Brown Schur complement A11 - A10 inv(A00) A01 47*c4762a1bSJed Brown A11 48*c4762a1bSJed Brown Mat Object: 1 MPI processes 49*c4762a1bSJed Brown type: seqaij 50*c4762a1bSJed Brown rows=24, cols=24 51*c4762a1bSJed Brown total: nonzeros=200, allocated nonzeros=200 52*c4762a1bSJed Brown total number of mallocs used during MatSetValues calls=0 53*c4762a1bSJed Brown using I-node routines: found 16 nodes, limit used is 5 54*c4762a1bSJed Brown A10 55*c4762a1bSJed Brown Mat Object: 1 MPI processes 56*c4762a1bSJed Brown type: seqaij 57*c4762a1bSJed Brown rows=24, cols=76 58*c4762a1bSJed Brown total: nonzeros=470, allocated nonzeros=470 59*c4762a1bSJed Brown total number of mallocs used during MatSetValues calls=0 60*c4762a1bSJed Brown not using I-node routines 61*c4762a1bSJed Brown KSP of A00 62*c4762a1bSJed Brown KSP Object: (fetidp_bddc_pc_bddc_dirichlet_) 1 MPI processes 63*c4762a1bSJed Brown type: preonly 64*c4762a1bSJed Brown maximum iterations=10000, initial guess is zero 65*c4762a1bSJed Brown tolerances: relative=1e-05, absolute=1e-50, divergence=10000. 66*c4762a1bSJed Brown left preconditioning 67*c4762a1bSJed Brown using NONE norm type for convergence test 68*c4762a1bSJed Brown PC Object: (fetidp_bddc_pc_bddc_dirichlet_) 1 MPI processes 69*c4762a1bSJed Brown type: svd 70*c4762a1bSJed Brown All singular values smaller than 1e-12 treated as zero 71*c4762a1bSJed Brown Provided essential rank of the matrix 0 (all other eigenvalues are zeroed) 72*c4762a1bSJed Brown linear system matrix = precond matrix: 73*c4762a1bSJed Brown Mat Object: (fetidp_bddc_pc_bddc_dirichlet_) 1 MPI processes 74*c4762a1bSJed Brown type: seqaij 75*c4762a1bSJed Brown rows=76, cols=76 76*c4762a1bSJed Brown total: nonzeros=1576, allocated nonzeros=1576 77*c4762a1bSJed Brown total number of mallocs used during MatSetValues calls=0 78*c4762a1bSJed Brown not using I-node routines 79*c4762a1bSJed Brown A01 80*c4762a1bSJed Brown Mat Object: 1 MPI processes 81*c4762a1bSJed Brown type: seqaij 82*c4762a1bSJed Brown rows=76, cols=24 83*c4762a1bSJed Brown total: nonzeros=470, allocated nonzeros=470 84*c4762a1bSJed Brown total number of mallocs used during MatSetValues calls=0 85*c4762a1bSJed Brown using I-node routines: found 48 nodes, limit used is 5 86*c4762a1bSJed Brown linear system matrix = precond matrix: 87*c4762a1bSJed Brown Mat Object: F (fetidp_) 5 MPI processes 88*c4762a1bSJed Brown type: shell 89*c4762a1bSJed Brown rows=60, cols=60 90*c4762a1bSJed Brown Inner BDDC solver details 91*c4762a1bSJed Brown PC Object: (fetidp_bddc_) 5 MPI processes 92*c4762a1bSJed Brown type: bddc 93*c4762a1bSJed Brown Use verbose output: 0 94*c4762a1bSJed Brown Use user-defined CSR: 0 95*c4762a1bSJed Brown Use local mat graph: 1 96*c4762a1bSJed Brown Connectivity graph topological dimension: 2 97*c4762a1bSJed Brown Graph max count: 2 98*c4762a1bSJed Brown Use vertices: 1 (vertex size 3) 99*c4762a1bSJed Brown Use edges: 1 100*c4762a1bSJed Brown Use faces: 1 101*c4762a1bSJed Brown Use true near null space: 0 102*c4762a1bSJed Brown Use QR for single constraints on cc: 0 103*c4762a1bSJed Brown Use change of basis on local edge nodes: 0 104*c4762a1bSJed Brown Use change of basis on local face nodes: 0 105*c4762a1bSJed Brown User defined change of basis matrix: 0 106*c4762a1bSJed Brown Has change of basis matrix: 1 107*c4762a1bSJed Brown Eliminate dirichlet boundary dofs: 1 108*c4762a1bSJed Brown Switch on static condensation ops around the interface preconditioner: 0 109*c4762a1bSJed Brown Use exact dirichlet trick: 1 110*c4762a1bSJed Brown Interface extension: DIRICHLET 111*c4762a1bSJed Brown Multilevel max levels: 0 112*c4762a1bSJed Brown Multilevel coarsening ratio: 8 113*c4762a1bSJed Brown Use estimated eigs for coarse problem: 0 114*c4762a1bSJed Brown Use deluxe scaling: 0 115*c4762a1bSJed Brown Use deluxe zerorows: 0 116*c4762a1bSJed Brown Use deluxe singlemat: 0 117*c4762a1bSJed Brown Rebuild interface graph for Schur principal minors: 0 118*c4762a1bSJed Brown Number of dofs' layers for the computation of principal minors: -1 119*c4762a1bSJed Brown Use user CSR graph to compute successive layers: 0 120*c4762a1bSJed Brown Adaptive constraint selection threshold (active 0, userdefined 0): 0. 121*c4762a1bSJed Brown Min constraints / connected component: 0 122*c4762a1bSJed Brown Max constraints / connected component: 0 123*c4762a1bSJed Brown Invert exact Schur complement for adaptive selection: 0 124*c4762a1bSJed Brown Symmetric computation of primal basis functions: 1 125*c4762a1bSJed Brown Num. Procs. to map coarse adjacency list: 0 126*c4762a1bSJed Brown Coarse eqs per proc (significant at the coarsest level): 1 127*c4762a1bSJed Brown Detect disconnected: 1 (filter 0) 128*c4762a1bSJed Brown Benign subspace trick: 1 (change explicit 1) 129*c4762a1bSJed Brown Benign subspace trick is active: 1 130*c4762a1bSJed Brown Algebraic computation of no-net-flux: 1 131*c4762a1bSJed Brown ********************************** STATISTICS AT LEVEL 0 ********************************** 132*c4762a1bSJed Brown Global dofs sizes: all 394 interface 68 coarse 23 133*c4762a1bSJed Brown Coarsening ratios: all/coarse 17 interface/coarse 2 134*c4762a1bSJed Brown Active processes : 5 135*c4762a1bSJed Brown Total subdomains : 6 136*c4762a1bSJed Brown Benign subs : 6 137*c4762a1bSJed Brown Dofs type : MIN MAX MEAN 138*c4762a1bSJed Brown Interior dofs : 36 76 65 139*c4762a1bSJed Brown Interface dofs : 16 40 28 140*c4762a1bSJed Brown Primal dofs : 5 14 9 141*c4762a1bSJed Brown Local dofs : 52 109 94 142*c4762a1bSJed Brown Local subs : 1 2 143*c4762a1bSJed Brown --- Interior solver (rank 0) 144*c4762a1bSJed Brown KSP Object: (fetidp_bddc_pc_bddc_dirichlet_) 1 MPI processes 145*c4762a1bSJed Brown type: preonly 146*c4762a1bSJed Brown maximum iterations=10000, initial guess is zero 147*c4762a1bSJed Brown tolerances: relative=1e-05, absolute=1e-50, divergence=10000. 148*c4762a1bSJed Brown left preconditioning 149*c4762a1bSJed Brown using NONE norm type for convergence test 150*c4762a1bSJed Brown PC Object: (fetidp_bddc_pc_bddc_dirichlet_) 1 MPI processes 151*c4762a1bSJed Brown type: svd 152*c4762a1bSJed Brown All singular values smaller than 1e-12 treated as zero 153*c4762a1bSJed Brown Provided essential rank of the matrix 0 (all other eigenvalues are zeroed) 154*c4762a1bSJed Brown linear system matrix = precond matrix: 155*c4762a1bSJed Brown Mat Object: (fetidp_bddc_pc_bddc_dirichlet_) 1 MPI processes 156*c4762a1bSJed Brown type: seqaij 157*c4762a1bSJed Brown rows=76, cols=76 158*c4762a1bSJed Brown total: nonzeros=1576, allocated nonzeros=1576 159*c4762a1bSJed Brown total number of mallocs used during MatSetValues calls=0 160*c4762a1bSJed Brown not using I-node routines 161*c4762a1bSJed Brown --- Correction solver (rank 0) 162*c4762a1bSJed Brown KSP Object: (fetidp_bddc_pc_bddc_neumann_) 1 MPI processes 163*c4762a1bSJed Brown type: preonly 164*c4762a1bSJed Brown maximum iterations=1000, initial guess is zero 165*c4762a1bSJed Brown tolerances: relative=1e-10, absolute=1e-10, divergence=10000. 166*c4762a1bSJed Brown left preconditioning 167*c4762a1bSJed Brown using NONE norm type for convergence test 168*c4762a1bSJed Brown PC Object: (fetidp_bddc_pc_bddc_neumann_) 1 MPI processes 169*c4762a1bSJed Brown type: svd 170*c4762a1bSJed Brown All singular values smaller than 1e-12 treated as zero 171*c4762a1bSJed Brown Provided essential rank of the matrix 0 (all other eigenvalues are zeroed) 172*c4762a1bSJed Brown linear system matrix = precond matrix: 173*c4762a1bSJed Brown Mat Object: (fetidp_bddc_pc_bddc_neumann_) 1 MPI processes 174*c4762a1bSJed Brown type: seqaij 175*c4762a1bSJed Brown rows=96, cols=96 176*c4762a1bSJed Brown total: nonzeros=2424, allocated nonzeros=2424 177*c4762a1bSJed Brown total number of mallocs used during MatSetValues calls=0 178*c4762a1bSJed Brown not using I-node routines 179*c4762a1bSJed Brown --- Coarse solver 180*c4762a1bSJed Brown KSP Object: (fetidp_bddc_pc_bddc_coarse_) 5 MPI processes 181*c4762a1bSJed Brown type: preonly 182*c4762a1bSJed Brown maximum iterations=1000, initial guess is zero 183*c4762a1bSJed Brown tolerances: relative=1e-10, absolute=1e-10, divergence=10000. 184*c4762a1bSJed Brown left preconditioning 185*c4762a1bSJed Brown using NONE norm type for convergence test 186*c4762a1bSJed Brown PC Object: (fetidp_bddc_pc_bddc_coarse_) 5 MPI processes 187*c4762a1bSJed Brown type: redundant 188*c4762a1bSJed Brown First (color=0) of 5 PCs follows 189*c4762a1bSJed Brown linear system matrix = precond matrix: 190*c4762a1bSJed Brown Mat Object: (fetidp_bddc_pc_bddc_coarse_) 5 MPI processes 191*c4762a1bSJed Brown type: mpiaij 192*c4762a1bSJed Brown rows=23, cols=23 193*c4762a1bSJed Brown total: nonzeros=365, allocated nonzeros=414 194*c4762a1bSJed Brown total number of mallocs used during MatSetValues calls=0 195*c4762a1bSJed Brown has attached null space 196*c4762a1bSJed Brown using I-node (on process 0) routines: found 3 nodes, limit used is 5 197*c4762a1bSJed Brown KSP Object: (fetidp_bddc_pc_bddc_coarse_redundant_) 1 MPI processes 198*c4762a1bSJed Brown type: preonly 199*c4762a1bSJed Brown maximum iterations=10000, initial guess is zero 200*c4762a1bSJed Brown tolerances: relative=1e-05, absolute=1e-50, divergence=10000. 201*c4762a1bSJed Brown left preconditioning 202*c4762a1bSJed Brown using NONE norm type for convergence test 203*c4762a1bSJed Brown PC Object: (fetidp_bddc_pc_bddc_coarse_redundant_) 1 MPI processes 204*c4762a1bSJed Brown type: svd 205*c4762a1bSJed Brown All singular values smaller than 1e-12 treated as zero 206*c4762a1bSJed Brown Provided essential rank of the matrix 0 (all other eigenvalues are zeroed) 207*c4762a1bSJed Brown linear system matrix = precond matrix: 208*c4762a1bSJed Brown Mat Object: 1 MPI processes 209*c4762a1bSJed Brown type: seqaij 210*c4762a1bSJed Brown rows=23, cols=23 211*c4762a1bSJed Brown total: nonzeros=365, allocated nonzeros=365 212*c4762a1bSJed Brown total number of mallocs used during MatSetValues calls=0 213*c4762a1bSJed Brown using I-node routines: found 16 nodes, limit used is 5 214*c4762a1bSJed Brown linear system matrix = precond matrix: 215*c4762a1bSJed Brown Mat Object: 5 MPI processes 216*c4762a1bSJed Brown type: is 217*c4762a1bSJed Brown rows=394, cols=394 218*c4762a1bSJed Brown total: nonzeros=11830, allocated nonzeros=11830 219*c4762a1bSJed Brown total number of mallocs used during MatSetValues calls=0 220*c4762a1bSJed Brown has attached near null space 221*c4762a1bSJed Brown maximum iterations=10000, initial guess is zero 222*c4762a1bSJed Brown tolerances: relative=1e-12, absolute=1e-10, divergence=10000. 223*c4762a1bSJed Brown left preconditioning 224*c4762a1bSJed Brown using NONE norm type for convergence test 225*c4762a1bSJed Brown PC Object: 5 MPI processes 226*c4762a1bSJed Brown type: none 227*c4762a1bSJed Brown linear system matrix followed by preconditioner matrix: 228*c4762a1bSJed Brown Mat Object: 5 MPI processes 229*c4762a1bSJed Brown type: is 230*c4762a1bSJed Brown rows=394, cols=394 231*c4762a1bSJed Brown total: nonzeros=11830, allocated nonzeros=14730 232*c4762a1bSJed Brown total number of mallocs used during MatSetValues calls=0 233*c4762a1bSJed Brown has attached null space 234*c4762a1bSJed Brown Mat Object: 5 MPI processes 235*c4762a1bSJed Brown type: is 236*c4762a1bSJed Brown rows=394, cols=394 237*c4762a1bSJed Brown total: nonzeros=11830, allocated nonzeros=14730 238*c4762a1bSJed Brown total number of mallocs used during MatSetValues calls=0 239*c4762a1bSJed BrownL_2 Error: 0.00418 [1.2e-05, 0.00418] 240*c4762a1bSJed BrownIntegral of pressure: 0. 241