(t_initf) Read in prof_inparam namelist from: input PERF_SETOPTS: PAPI library not linked in. Request to enable PAPI ignored. (t_initf) Using profile_disable= F profile_timer= 2 (t_initf) profile_depth_limit= 99999 profile_detail_limit= 1 (t_initf) profile_barrier= F profile_outpe_num= 1 (t_initf) profile_outpe_stride= 1 profile_single_file= F (t_initf) profile_global_stats= T profile_papi_enable= F call petsc_init call adios_init use adios reader v2 # of OMP threads = 1 Total simulation processors : 32 setting up......... 0 Read input file sml_param ptl_param eq_param col_param diag_param smooth_param mon_param Interplane-Major ordering for sml_comm communicator plane mpi communication inter-plane mpi communication adios mpi communication? end of commnunicator setup reading...... eq end_flag -1 ptl_num for each CPU = 100000 memory allocation init_interpolation libhugetlbfs [nid00594:27475]: ERROR: RTLD_NEXT used in code not dynamically loadedlibhugetlbfs [nid00595:7974]: ERROR: RTLD_NEXT used in code not dynamically loadedlibhugetlbfs [nid00594:27461]: ERROR: RTLD_NEXT used in code not dynamically loadedlibhugetlbfs [nid00594:27474]: ERROR: RTLD_NEXT used in code not dynamically loadedlibhugetlbfs [nid00594:27464]: ERROR: RTLD_NEXT used in code not dynamically loadedlibhugetlbfs [nid00594:27470]: ERROR: RTLD_NEXT used in code not dynamically loadedlibhugetlbfs [nid00594:27472]: ERROR: RTLD_NEXT used in code not dynamically loadedlibhugetlbfs [nid00594:27468]: ERROR: RTLD_NEXT used in code not dynamically loadedlibhugetlbfs [nid00594:27477]: ERROR: RTLD_NEXT used in code not dynamically loadedlibhugetlbfs [nid00594:27462]: ERROR: RTLD_NEXT used in code not dynamically loadedlibhugetlbfs [nid00594:27465]: ERROR: RTLD_NEXT used in code not dynamically loadedlibhugetlbfs [nid00594:27480]: ERROR: RTLD_NEXT used in code not dynamically loadedlibhugetlbfs [nid00594:27471]: ERROR: RTLD_NEXT used in code not dynamically loadedlibhugetlbfs [nid00594:27463]: ERROR: RTLD_NEXT used in code not dynamically loadedlibhugetlbfs [nid00594:27479]: ERROR: RTLD_NEXT used in code not dynamically loadedlibhugetlbfs [nid00594:27481]: ERROR: RTLD_NEXT used in code not dynamically loadedlibhugetlbfs [nid00594:27469]: ERROR: RTLD_NEXT used in code not dynamically loadedlibhugetlbfs [nid00594:27473]: ERROR: RTLD_NEXT used in code not dynamically loadedlibhugetlbfs [nid00594:27466]: ERROR: RTLD_NEXT used in code not dynamically loadedlibhugetlbfs [nid00594:27478]: ERROR: RTLD_NEXT used in code not dynamically loadedlibhugetlbfs [nid00594:27460]: ERROR: RTLD_NEXT used in code not dynamically loadedlibhugetlbfs [nid00594:27467]: ERROR: RTLD_NEXT used in code not dynamically loadedlibhugetlbfs [nid00594:27483]: ERROR: RTLD_NEXT used in code not dynamically loadedlibhugetlbfs [nid00594:27482]: ERROR: RTLD_NEXT used in code not dynamically loadedlibhugetlbfs [nid00594:27476]: ERROR: RTLD_NEXT used in code not dynamically loadedlibhugetlbfs [nid00595:7980]: ERROR: RTLD_NEXT used in code not dynamically loadedlibhugetlbfs [nid00595:7978]: ERROR: RTLD_NEXT used in code not dynamically loadedlibhugetlbfs [nid00595:7975]: ERROR: RTLD_NEXT used in code not dynamically loadedlibhugetlbfs [nid00595:7979]: ERROR: RTLD_NEXT used in code not dynamically loadedlibhugetlbfs [nid00595:7976]: ERROR: RTLD_NEXT used in code not dynamically loadedlibhugetlbfs [nid00595:7973]: ERROR: RTLD_NEXT used in code not dynamically loadedlibhugetlbfs [nid00595:7977]: ERROR: RTLD_NEXT used in code not dynamically loaded third difinition part eq_x_slope= 0.000000000000000E+000 eq_x2_slope= 0.000000000000000E+000 R range : 1.03053998549000 2.36946001451000 Z range : -0.669460014510200 0.669460014510200 bounce setup get_mid_r setup diag_efld setup setting up grid system..... 0 Guess min 1.03053998549000 -0.669460014510200 Guess max 2.36946001451000 0.669460014510200 init_grid open node file (*.node) : ../XGC-1_inputs/cyclon16mm_32.1.node open connectivity file (*.ele) : ../XGC-1_inputs/cyclon16mm_32.1.ele grid%nnode = 3705 sum(tr_vol) 1.71167703906022 sum(node_area) 1.28200418926821 sum(node_vol) 1.71167703906022 init_guess_table init bfollow, init_ff diagnosis_1d init Monte-Carlo volume calculation... smoothing initializing load 0 nonzero n-mode, potential boundary extended bd nseg 2 iseg( 1 )= 1 1 iseg( 2 )= 2873 3705 nonzero n-mode, charge boundary extended bd nseg 2 iseg( 1 )= 1 29 iseg( 2 )= 2706 3705 n=0-mode, potential boundary extended bd nseg 1 iseg( 1 )= 2873 3705 n=0-mode, charge boundary extended bd nseg 2 iseg( 1 )= 1 1 iseg( 2 )= 2706 3705 loading end 0 Init poisson Initialize 00 solver - 00 solver memory create FSA SOLVER: 20 flux surfaces, (doit= T ) make A (A_00) mat. FSA: Add to A mat. with 7309 triangles. FSA: create B mat. M,N= 20 3705 FSA: make B. low= 0 , high= 927 Initialize turb solver - turb solver memory create start init_1field_solver prepare gyro-averaging matrix setup poloidal decomposition initial diagnosis dump_bfield background elec charge **************** solve_poisson_private static variable initialized ************ **** End of setup main loop started partition_opt: (using Alg. 2) 0.25028E+00 partition_opt: (alt. Alg. 1) 0.25195E+00 0.24495E+00 10 - do loop start step,ratio,# of ion 1 1.007620 3200000 petsc_solve: setting up Schur complement Linear fsa_fieldsplit_lambda_inner_ solve converged due to CONVERGED_RTOL iterations 12 Linear fsa_fieldsplit_lambda_inner_ solve converged due to CONVERGED_RTOL iterations 13 Linear fsa_fieldsplit_lambda_inner_ solve converged due to CONVERGED_RTOL iterations 13 Linear fsa_fieldsplit_lambda_inner_ solve converged due to CONVERGED_RTOL iterations 14 Linear fsa_fieldsplit_lambda_inner_ solve converged due to CONVERGED_RTOL iterations 14 Linear fsa_fieldsplit_lambda_inner_ solve converged due to CONVERGED_RTOL iterations 14 Linear fsa_fieldsplit_lambda_inner_ solve converged due to CONVERGED_RTOL iterations 15 Linear fsa_fieldsplit_lambda_inner_ solve converged due to CONVERGED_RTOL iterations 15 Linear fsa_fieldsplit_lambda_inner_ solve converged due to CONVERGED_RTOL iterations 15 Linear fsa_fieldsplit_lambda_inner_ solve converged due to CONVERGED_RTOL iterations 15 Linear fsa_fieldsplit_lambda_inner_ solve converged due to CONVERGED_RTOL iterations 15 Linear fsa_fieldsplit_lambda_inner_ solve converged due to CONVERGED_RTOL iterations 15 Linear fsa_fieldsplit_lambda_inner_ solve converged due to CONVERGED_RTOL iterations 15 Linear fsa_fieldsplit_lambda_inner_ solve converged due to CONVERGED_RTOL iterations 15 Linear fsa_fieldsplit_lambda_inner_ solve converged due to CONVERGED_RTOL iterations 15 Linear fsa_fieldsplit_lambda_inner_ solve converged due to CONVERGED_RTOL iterations 15 Linear fsa_fieldsplit_lambda_inner_ solve converged due to CONVERGED_RTOL iterations 15 Linear fsa_fieldsplit_lambda_inner_ solve converged due to CONVERGED_RTOL iterations 15 Linear fsa_fieldsplit_lambda_inner_ solve converged due to CONVERGED_RTOL iterations 15 Linear fsa_fieldsplit_lambda_inner_ solve converged due to CONVERGED_RTOL iterations 15 0 SNES Function norm 2.787917906837e-11 Residual norms for fsa_ solve. 0 KSP Residual norm 2.787917906837e-11 Linear fsa_fieldsplit_phi_ solve converged due to CONVERGED_RTOL iterations 19 Residual norms for fsa_fieldsplit_lambda_ solve. 0 KSP Residual norm 1.137395651284e-03 Linear fsa_fieldsplit_lambda_inner_ solve converged due to CONVERGED_RTOL iterations 14 1 KSP Residual norm 1.128554060335e-03 Linear fsa_fieldsplit_lambda_inner_ solve converged due to CONVERGED_RTOL iterations 15 2 KSP Residual norm 1.052344323087e-03 Linear fsa_fieldsplit_lambda_inner_ solve converged due to CONVERGED_RTOL iterations 15 3 KSP Residual norm 1.002723605629e-03 Linear fsa_fieldsplit_lambda_inner_ solve converged due to CONVERGED_RTOL iterations 15 4 KSP Residual norm 9.742118406180e-04 Linear fsa_fieldsplit_lambda_inner_ solve converged due to CONVERGED_RTOL iterations 15 5 KSP Residual norm 8.805098337403e-04 Linear fsa_fieldsplit_lambda_inner_ solve converged due to CONVERGED_RTOL iterations 15 6 KSP Residual norm 7.126652505928e-04 Linear fsa_fieldsplit_lambda_inner_ solve converged due to CONVERGED_RTOL iterations 15 7 KSP Residual norm 6.422221670734e-04 Linear fsa_fieldsplit_lambda_inner_ solve converged due to CONVERGED_RTOL iterations 15 8 KSP Residual norm 6.096647312050e-04 Linear fsa_fieldsplit_lambda_inner_ solve converged due to CONVERGED_RTOL iterations 15 9 KSP Residual norm 5.448833672044e-04 Linear fsa_fieldsplit_lambda_inner_ solve converged due to CONVERGED_RTOL iterations 15 10 KSP Residual norm 5.085599658599e-04 Linear fsa_fieldsplit_lambda_inner_ solve converged due to CONVERGED_RTOL iterations 15 11 KSP Residual norm 4.500226838175e-04 Linear fsa_fieldsplit_lambda_inner_ solve converged due to CONVERGED_RTOL iterations 15 12 KSP Residual norm 3.890160785864e-04 Linear fsa_fieldsplit_lambda_inner_ solve converged due to CONVERGED_RTOL iterations 15 13 KSP Residual norm 3.421662052661e-04 Linear fsa_fieldsplit_lambda_inner_ solve converged due to CONVERGED_RTOL iterations 15 14 KSP Residual norm 2.795271537693e-04 Linear fsa_fieldsplit_lambda_inner_ solve converged due to CONVERGED_RTOL iterations 15 15 KSP Residual norm 2.518425576162e-04 Linear fsa_fieldsplit_lambda_inner_ solve converged due to CONVERGED_RTOL iterations 15 16 KSP Residual norm 2.042190098742e-04 Linear fsa_fieldsplit_lambda_inner_ solve converged due to CONVERGED_RTOL iterations 15 17 KSP Residual norm 1.775002211922e-04 Linear fsa_fieldsplit_lambda_inner_ solve converged due to CONVERGED_RTOL iterations 15 18 KSP Residual norm 5.014002458711e-05 Linear fsa_fieldsplit_lambda_inner_ solve converged due to CONVERGED_RTOL iterations 15 19 KSP Residual norm 1.027247491223e-05 Linear fsa_fieldsplit_lambda_inner_ solve converged due to CONVERGED_RTOL iterations 15 20 KSP Residual norm 7.612495277727e-17 Linear fsa_fieldsplit_lambda_ solve converged due to CONVERGED_RTOL iterations 20 Linear fsa_fieldsplit_phi_ solve converged due to CONVERGED_RTOL iterations 20 1 KSP Residual norm 1.414274792264e-21 Linear fsa_fieldsplit_phi_ solve converged due to CONVERGED_RTOL iterations 22 Residual norms for fsa_fieldsplit_lambda_ solve. 0 KSP Residual norm 1.235140093850e-04 Linear fsa_fieldsplit_lambda_inner_ solve converged due to CONVERGED_RTOL iterations 15 1 KSP Residual norm 5.916065076726e-05 Linear fsa_fieldsplit_lambda_inner_ solve converged due to CONVERGED_RTOL iterations 15 2 KSP Residual norm 4.591862058126e-05 Linear fsa_fieldsplit_lambda_inner_ solve converged due to CONVERGED_RTOL iterations 15 3 KSP Residual norm 3.548051010066e-05 Linear fsa_fieldsplit_lambda_inner_ solve converged due to CONVERGED_RTOL iterations 14 4 KSP Residual norm 2.936599475239e-05 Linear fsa_fieldsplit_lambda_inner_ solve converged due to CONVERGED_RTOL iterations 15 5 KSP Residual norm 2.563081567874e-05 Linear fsa_fieldsplit_lambda_inner_ solve converged due to CONVERGED_RTOL iterations 15 6 KSP Residual norm 2.415591000692e-05 Linear fsa_fieldsplit_lambda_inner_ solve converged due to CONVERGED_RTOL iterations 14 7 KSP Residual norm 2.325948998846e-05 Linear fsa_fieldsplit_lambda_inner_ solve converged due to CONVERGED_RTOL iterations 15 8 KSP Residual norm 2.266243522248e-05 Linear fsa_fieldsplit_lambda_inner_ solve converged due to CONVERGED_RTOL iterations 15 9 KSP Residual norm 2.128465835692e-05 Linear fsa_fieldsplit_lambda_inner_ solve converged due to CONVERGED_RTOL iterations 15 10 KSP Residual norm 1.995419543224e-05 Linear fsa_fieldsplit_lambda_inner_ solve converged due to CONVERGED_RTOL iterations 15 11 KSP Residual norm 1.816619936019e-05 Linear fsa_fieldsplit_lambda_inner_ solve converged due to CONVERGED_RTOL iterations 15 12 KSP Residual norm 1.614264873374e-05 Linear fsa_fieldsplit_lambda_inner_ solve converged due to CONVERGED_RTOL iterations 15 13 KSP Residual norm 1.439168777342e-05 Linear fsa_fieldsplit_lambda_inner_ solve converged due to CONVERGED_RTOL iterations 15 14 KSP Residual norm 1.201523848017e-05 Linear fsa_fieldsplit_lambda_inner_ solve converged due to CONVERGED_RTOL iterations 15 15 KSP Residual norm 1.036617539918e-05 Linear fsa_fieldsplit_lambda_inner_ solve converged due to CONVERGED_RTOL iterations 14 16 KSP Residual norm 8.574558177929e-06 Linear fsa_fieldsplit_lambda_inner_ solve converged due to CONVERGED_RTOL iterations 15 17 KSP Residual norm 7.660619967582e-06 Linear fsa_fieldsplit_lambda_inner_ solve converged due to CONVERGED_RTOL iterations 15 18 KSP Residual norm 7.352150700734e-06 Linear fsa_fieldsplit_lambda_inner_ solve converged due to CONVERGED_RTOL iterations 15 19 KSP Residual norm 2.142008990262e-07 Linear fsa_fieldsplit_lambda_inner_ solve converged due to CONVERGED_RTOL iterations 15 20 KSP Residual norm 1.368965502943e-18 Linear fsa_fieldsplit_lambda_ solve converged due to CONVERGED_RTOL iterations 20 Linear fsa_fieldsplit_phi_ solve converged due to CONVERGED_RTOL iterations 21 2 KSP Residual norm 5.260079183868e-27 Linear fsa_ solve converged due to CONVERGED_RTOL iterations 2 KSP Object:(fsa_) 4 MPI processes type: fgmres GMRES: restart=50, using Classical (unmodified) Gram-Schmidt Orthogonalization with no iterative refinement GMRES: happy breakdown tolerance 1e-30 maximum iterations=3, initial guess is zero tolerances: relative=1e-12, absolute=1e-50, divergence=10000 right preconditioning using UNPRECONDITIONED norm type for convergence test PC Object:(fsa_) 4 MPI processes type: fieldsplit FieldSplit with Schur preconditioner, factorization FULL Preconditioner for the Schur complement formed from A11 Split info: Split number 0 Defined by IS Split number 1 Defined by IS KSP solver for A00 block KSP Object: (fsa_fieldsplit_phi_) 4 MPI processes type: cg maximum iterations=10000, initial guess is zero tolerances: relative=1e-12, absolute=1e-50, divergence=10000 left preconditioning using PRECONDITIONED norm type for convergence test PC Object: (fsa_fieldsplit_phi_) 4 MPI processes type: gamg MG: type is MULTIPLICATIVE, levels=2 cycles=v Cycles per PCApply=1 Using Galerkin computed coarse grid matrices Coarse grid solver -- level ------------------------------- KSP Object: (fsa_fieldsplit_phi_mg_coarse_) 4 MPI processes type: preonly maximum iterations=1, initial guess is zero tolerances: relative=1e-05, absolute=1e-50, divergence=10000 left preconditioning using NONE norm type for convergence test PC Object: (fsa_fieldsplit_phi_mg_coarse_) 4 MPI processes type: bjacobi block Jacobi: number of blocks = 4 Local solve is same for all blocks, in the following KSP and PC objects: KSP Object: (fsa_fieldsplit_phi_mg_coarse_sub_) 1 MPI processes type: preonly maximum iterations=1, initial guess is zero tolerances: relative=1e-05, absolute=1e-50, divergence=10000 left preconditioning using NONE norm type for convergence test PC Object: (fsa_fieldsplit_phi_mg_coarse_sub_) 1 MPI processes type: lu LU: out-of-place factorization tolerance for zero pivot 2.22045e-14 using diagonal shift on blocks to prevent zero pivot [INBLOCKS] matrix ordering: nd factor fill ratio given 5, needed 4.18681 Factored matrix follows: Mat Object: 1 MPI processes type: seqaij rows=298, cols=298 package used to perform factorization: petsc total: nonzeros=13582, allocated nonzeros=13582 total number of mallocs used during MatSetValues calls =0 not using I-node routines linear system matrix = precond matrix: Mat Object: 1 MPI processes type: seqaij rows=298, cols=298 total: nonzeros=3244, allocated nonzeros=3244 total number of mallocs used during MatSetValues calls =0 not using I-node routines linear system matrix = precond matrix: Mat Object: 4 MPI processes type: mpiaij rows=298, cols=298 total: nonzeros=3244, allocated nonzeros=3244 total number of mallocs used during MatSetValues calls =0 not using I-node (on process 0) routines Down solver (pre-smoother) on level 1 ------------------------------- KSP Object: (fsa_fieldsplit_phi_mg_levels_1_) 4 MPI processes type: richardson Richardson: damping factor=1 maximum iterations=2 tolerances: relative=1e-05, absolute=1e-50, divergence=10000 left preconditioning using nonzero initial guess using NONE norm type for convergence test PC Object: (fsa_fieldsplit_phi_mg_levels_1_) 4 MPI processes type: sor SOR: type = local_symmetric, iterations = 1, local iterations = 1, omega = 1 linear system matrix = precond matrix: Mat Object: (fsa_fieldsplit_phi_) 4 MPI processes type: mpiaij rows=3705, cols=3705 total: nonzeros=20597, allocated nonzeros=555750 total number of mallocs used during MatSetValues calls =0 not using I-node (on process 0) routines Up solver (post-smoother) same as down solver (pre-smoother) linear system matrix = precond matrix: Mat Object: (fsa_fieldsplit_phi_) 4 MPI processes type: mpiaij rows=3705, cols=3705 total: nonzeros=20597, allocated nonzeros=555750 total number of mallocs used during MatSetValues calls =0 not using I-node (on process 0) routines KSP solver for S = A11 - A10 inv(A00) A01 KSP Object: (fsa_fieldsplit_lambda_) 4 MPI processes type: gmres GMRES: restart=30, using Classical (unmodified) Gram-Schmidt Orthogonalization with no iterative refinement GMRES: happy breakdown tolerance 1e-30 maximum iterations=10000, initial guess is zero tolerances: relative=1e-12, absolute=1e-50, divergence=10000 left preconditioning using PRECONDITIONED norm type for convergence test PC Object: (fsa_fieldsplit_lambda_) 4 MPI processes type: bjacobi block Jacobi: number of blocks = 4 Local solve is same for all blocks, in the following KSP and PC objects: KSP Object: (fsa_fieldsplit_lambda_sub_) 1 MPI processes type: preonly maximum iterations=10000, initial guess is zero tolerances: relative=1e-05, absolute=1e-50, divergence=10000 left preconditioning using NONE norm type for convergence test PC Object: (fsa_fieldsplit_lambda_sub_) 1 MPI processes type: lu LU: out-of-place factorization tolerance for zero pivot 2.22045e-14 matrix ordering: nd factor fill ratio given 5, needed 1 Factored matrix follows: Mat Object: 1 MPI processes type: seqaij rows=20, cols=20 package used to perform factorization: petsc total: nonzeros=20, allocated nonzeros=20 total number of mallocs used during MatSetValues calls =0 not using I-node routines linear system matrix = precond matrix: Mat Object: (fsa_fieldsplit_lambda_) 1 MPI processes type: seqaij rows=20, cols=20 total: nonzeros=20, allocated nonzeros=20 total number of mallocs used during MatSetValues calls =0 not using I-node routines linear system matrix followed by preconditioner matrix: Mat Object: (fsa_fieldsplit_lambda_) 4 MPI processes type: schurcomplement rows=20, cols=20 Schur complement A11 - A10 inv(A00) A01 A11 Mat Object: (fsa_fieldsplit_lambda_) 4 MPI processes type: mpiaij rows=20, cols=20 total: nonzeros=20, allocated nonzeros=40 total number of mallocs used during MatSetValues calls =0 not using I-node (on process 0) routines A10 Mat Object: 4 MPI processes type: mpiaij rows=20, cols=3705 total: nonzeros=9501, allocated nonzeros=19002 total number of mallocs used during MatSetValues calls =0 using I-node (on process 0) routines: found 12 nodes, limit used is 5 KSP of A00 KSP Object: (fsa_fieldsplit_lambda_inner_) 4 MPI processes type: cg maximum iterations=10000, initial guess is zero tolerances: relative=1e-12, absolute=1e-50, divergence=10000 left preconditioning using PRECONDITIONED norm type for convergence test PC Object: (fsa_fieldsplit_lambda_inner_) 4 MPI processes type: gamg MG: type is MULTIPLICATIVE, levels=2 cycles=v Cycles per PCApply=1 Using Galerkin computed coarse grid matrices Coarse grid solver -- level ------------------------------- KSP Object: (fsa_fieldsplit_lambda_inner_mg_coarse_) 4 MPI processes type: preonly maximum iterations=1, initial guess is zero tolerances: relative=1e-05, absolute=1e-50, divergence=10000 left preconditioning using NONE norm type for convergence test PC Object: (fsa_fieldsplit_lambda_inner_mg_coarse_) 4 MPI processes type: bjacobi block Jacobi: number of blocks = 4 Local solve is same for all blocks, in the following KSP and PC objects: KSP Object: (fsa_fieldsplit_lambda_inner_mg_coarse_sub_) 1 MPI processes type: preonly maximum iterations=1, initial guess is zero tolerances: relative=1e-05, absolute=1e-50, divergence=10000 left preconditioning using NONE norm type for convergence test PC Object: (fsa_fieldsplit_lambda_inner_mg_coarse_sub_) 1 MPI processes type: lu LU: out-of-place factorization tolerance for zero pivot 2.22045e-14 using diagonal shift on blocks to prevent zero pivot [INBLOCKS] matrix ordering: nd factor fill ratio given 5, needed 4.18681 Factored matrix follows: Mat Object: 1 MPI processes type: seqaij rows=298, cols=298 package used to perform factorization: petsc total: nonzeros=13582, allocated nonzeros=13582 total number of mallocs used during MatSetValues calls =0 not using I-node routines linear system matrix = precond matrix: Mat Object: 1 MPI processes type: seqaij rows=298, cols=298 total: nonzeros=3244, allocated nonzeros=3244 total number of mallocs used during MatSetValues calls =0 not using I-node routines linear system matrix = precond matrix: Mat Object: 4 MPI processes type: mpiaij rows=298, cols=298 total: nonzeros=3244, allocated nonzeros=3244 total number of mallocs used during MatSetValues calls =0 not using I-node (on process 0) routines Down solver (pre-smoother) on level 1 ------------------------------- KSP Object: (fsa_fieldsplit_lambda_inner_mg_levels_1_) 4 MPI processes type: chebyshev Chebyshev: eigenvalue estimates: min = 0.147014, max = 3.0873 maximum iterations=2 tolerances: relative=1e-05, absolute=1e-50, divergence=10000 left preconditioning using nonzero initial guess using NONE norm type for convergence test PC Object: (fsa_fieldsplit_lambda_inner_mg_levels_1_) 4 MPI processes type: sor SOR: type = local_symmetric, iterations = 1, local iterations = 1, omega = 1 linear system matrix = precond matrix: Mat Object: (fsa_fieldsplit_phi_) 4 MPI processes type: mpiaij rows=3705, cols=3705 total: nonzeros=20597, allocated nonzeros=555750 total number of mallocs used during MatSetValues calls =0 not using I-node (on process 0) routines Up solver (post-smoother) same as down solver (pre-smoother) linear system matrix = precond matrix: Mat Object: (fsa_fieldsplit_phi_) 4 MPI processes type: mpiaij rows=3705, cols=3705 total: nonzeros=20597, allocated nonzeros=555750 total number of mallocs used during MatSetValues calls =0 not using I-node (on process 0) routines A01 Mat Object: 4 MPI processes type: mpiaij rows=3705, cols=20 total: nonzeros=9303, allocated nonzeros=9303 total number of mallocs used during MatSetValues calls =0 using I-node (on process 0) routines: found 190 nodes, limit used is 5 Mat Object: (fsa_fieldsplit_lambda_) 4 MPI processes type: mpiaij rows=20, cols=20 total: nonzeros=20, allocated nonzeros=40 total number of mallocs used during MatSetValues calls =0 not using I-node (on process 0) routines linear system matrix = precond matrix: Mat Object: 4 MPI processes type: nest rows=3725, cols=3725 Matrix object: type=nest, rows=2, cols=2 MatNest structure: (0,0) : prefix="fsa_fieldsplit_phi_", type=mpiaij, rows=3705, cols=3705 (0,1) : type=mpiaij, rows=3705, cols=20 (1,0) : type=mpiaij, rows=20, cols=3705 (1,1) : prefix="fsa_fieldsplit_lambda_", type=mpiaij, rows=20, cols=20 1 SNES Function norm 4.441813144645e-24 Nonlinear fsa_ solve converged due to CONVERGED_ITS iterations 1 |C phi - | / || = 3.8475E-16, ||= 4.1991E-12, ||/|phi|= 9.8870E-01 2d particle and heat flux is not ready yet for deltaf-f method 0 SNES Function norm 2.963555088136e+02 Residual norms for fsa_ solve. 0 KSP Residual norm 2.963555088136e+02 Linear fsa_fieldsplit_phi_ solve converged due to CONVERGED_RTOL iterations 19 Residual norms for fsa_fieldsplit_lambda_ solve. 0 KSP Residual norm 4.995891851558e-06 Linear fsa_fieldsplit_lambda_inner_ solve converged due to CONVERGED_RTOL iterations 15 1 KSP Residual norm 2.715409504032e-06 Linear fsa_fieldsplit_lambda_inner_ solve converged due to CONVERGED_RTOL iterations 15 2 KSP Residual norm 1.929214177358e-06 Linear fsa_fieldsplit_lambda_inner_ solve converged due to CONVERGED_RTOL iterations 15 3 KSP Residual norm 1.609858508337e-06 Linear fsa_fieldsplit_lambda_inner_ solve converged due to CONVERGED_RTOL iterations 15 4 KSP Residual norm 1.372457921000e-06 Linear fsa_fieldsplit_lambda_inner_ solve converged due to CONVERGED_RTOL iterations 15 5 KSP Residual norm 1.281510967656e-06 Linear fsa_fieldsplit_lambda_inner_ solve converged due to CONVERGED_RTOL iterations 15 6 KSP Residual norm 1.091376459209e-06 Linear fsa_fieldsplit_lambda_inner_ solve converged due to CONVERGED_RTOL iterations 15 7 KSP Residual norm 1.071024388980e-06 Linear fsa_fieldsplit_lambda_inner_ solve converged due to CONVERGED_RTOL iterations 15 8 KSP Residual norm 1.039932347970e-06 Linear fsa_fieldsplit_lambda_inner_ solve converged due to CONVERGED_RTOL iterations 15 9 KSP Residual norm 1.017009602837e-06 Linear fsa_fieldsplit_lambda_inner_ solve converged due to CONVERGED_RTOL iterations 15 10 KSP Residual norm 9.832222040058e-07 Linear fsa_fieldsplit_lambda_inner_ solve converged due to CONVERGED_RTOL iterations 15 11 KSP Residual norm 9.323902437711e-07 Linear fsa_fieldsplit_lambda_inner_ solve converged due to CONVERGED_RTOL iterations 15 12 KSP Residual norm 8.949854068652e-07 Linear fsa_fieldsplit_lambda_inner_ solve converged due to CONVERGED_RTOL iterations 15 13 KSP Residual norm 8.119936645246e-07 Linear fsa_fieldsplit_lambda_inner_ solve converged due to CONVERGED_RTOL iterations 15 14 KSP Residual norm 6.649600103835e-07 Linear fsa_fieldsplit_lambda_inner_ solve converged due to CONVERGED_RTOL iterations 15 15 KSP Residual norm 4.498846283438e-07 Linear fsa_fieldsplit_lambda_inner_ solve converged due to CONVERGED_RTOL iterations 15 16 KSP Residual norm 3.160271742987e-07 Linear fsa_fieldsplit_lambda_inner_ solve converged due to CONVERGED_RTOL iterations 15 17 KSP Residual norm 1.563923650135e-07 Linear fsa_fieldsplit_lambda_inner_ solve converged due to CONVERGED_RTOL iterations 15 18 KSP Residual norm 1.093472275581e-07 Linear fsa_fieldsplit_lambda_inner_ solve converged due to CONVERGED_RTOL iterations 15 19 KSP Residual norm 1.045080205752e-07 Linear fsa_fieldsplit_lambda_inner_ solve converged due to CONVERGED_RTOL iterations 15 20 KSP Residual norm 7.524210239844e-20 Linear fsa_fieldsplit_lambda_ solve converged due to CONVERGED_RTOL iterations 20 Linear fsa_fieldsplit_phi_ solve converged due to CONVERGED_RTOL iterations 19 1 KSP Residual norm 1.005146975949e-10 Linear fsa_ solve converged due to CONVERGED_RTOL iterations 1 KSP Object:(fsa_) 4 MPI processes type: fgmres GMRES: restart=50, using Classical (unmodified) Gram-Schmidt Orthogonalization with no iterative refinement GMRES: happy breakdown tolerance 1e-30 maximum iterations=3, initial guess is zero tolerances: relative=1e-12, absolute=1e-50, divergence=10000 right preconditioning using UNPRECONDITIONED norm type for convergence test PC Object:(fsa_) 4 MPI processes type: fieldsplit FieldSplit with Schur preconditioner, factorization FULL Preconditioner for the Schur complement formed from A11 Split info: Split number 0 Defined by IS Split number 1 Defined by IS KSP solver for A00 block KSP Object: (fsa_fieldsplit_phi_) 4 MPI processes type: cg maximum iterations=10000, initial guess is zero tolerances: relative=1e-12, absolute=1e-50, divergence=10000 left preconditioning using PRECONDITIONED norm type for convergence test PC Object: (fsa_fieldsplit_phi_) 4 MPI processes type: gamg MG: type is MULTIPLICATIVE, levels=2 cycles=v Cycles per PCApply=1 Using Galerkin computed coarse grid matrices Coarse grid solver -- level ------------------------------- KSP Object: (fsa_fieldsplit_phi_mg_coarse_) 4 MPI processes type: preonly maximum iterations=1, initial guess is zero tolerances: relative=1e-05, absolute=1e-50, divergence=10000 left preconditioning using NONE norm type for convergence test PC Object: (fsa_fieldsplit_phi_mg_coarse_) 4 MPI processes type: bjacobi block Jacobi: number of blocks = 4 Local solve is same for all blocks, in the following KSP and PC objects: KSP Object: (fsa_fieldsplit_phi_mg_coarse_sub_) 1 MPI processes type: preonly maximum iterations=1, initial guess is zero tolerances: relative=1e-05, absolute=1e-50, divergence=10000 left preconditioning using NONE norm type for convergence test PC Object: (fsa_fieldsplit_phi_mg_coarse_sub_) 1 MPI processes type: lu LU: out-of-place factorization tolerance for zero pivot 2.22045e-14 using diagonal shift on blocks to prevent zero pivot [INBLOCKS] matrix ordering: nd factor fill ratio given 5, needed 4.18681 Factored matrix follows: Mat Object: 1 MPI processes type: seqaij rows=298, cols=298 package used to perform factorization: petsc total: nonzeros=13582, allocated nonzeros=13582 total number of mallocs used during MatSetValues calls =0 not using I-node routines linear system matrix = precond matrix: Mat Object: 1 MPI processes type: seqaij rows=298, cols=298 total: nonzeros=3244, allocated nonzeros=3244 total number of mallocs used during MatSetValues calls =0 not using I-node routines linear system matrix = precond matrix: Mat Object: 4 MPI processes type: mpiaij rows=298, cols=298 total: nonzeros=3244, allocated nonzeros=3244 total number of mallocs used during MatSetValues calls =0 not using I-node (on process 0) routines Down solver (pre-smoother) on level 1 ------------------------------- KSP Object: (fsa_fieldsplit_phi_mg_levels_1_) 4 MPI processes type: richardson Richardson: damping factor=1 maximum iterations=2 tolerances: relative=1e-05, absolute=1e-50, divergence=10000 left preconditioning using nonzero initial guess using NONE norm type for convergence test PC Object: (fsa_fieldsplit_phi_mg_levels_1_) 4 MPI processes type: sor SOR: type = local_symmetric, iterations = 1, local iterations = 1, omega = 1 linear system matrix = precond matrix: Mat Object: (fsa_fieldsplit_phi_) 4 MPI processes type: mpiaij rows=3705, cols=3705 total: nonzeros=20597, allocated nonzeros=555750 total number of mallocs used during MatSetValues calls =0 not using I-node (on process 0) routines Up solver (post-smoother) same as down solver (pre-smoother) linear system matrix = precond matrix: Mat Object: (fsa_fieldsplit_phi_) 4 MPI processes type: mpiaij rows=3705, cols=3705 total: nonzeros=20597, allocated nonzeros=555750 total number of mallocs used during MatSetValues calls =0 not using I-node (on process 0) routines KSP solver for S = A11 - A10 inv(A00) A01 KSP Object: (fsa_fieldsplit_lambda_) 4 MPI processes type: gmres GMRES: restart=30, using Classical (unmodified) Gram-Schmidt Orthogonalization with no iterative refinement GMRES: happy breakdown tolerance 1e-30 maximum iterations=10000, initial guess is zero tolerances: relative=1e-12, absolute=1e-50, divergence=10000 left preconditioning using PRECONDITIONED norm type for convergence test PC Object: (fsa_fieldsplit_lambda_) 4 MPI processes type: bjacobi block Jacobi: number of blocks = 4 Local solve is same for all blocks, in the following KSP and PC objects: KSP Object: (fsa_fieldsplit_lambda_sub_) 1 MPI processes type: preonly maximum iterations=10000, initial guess is zero tolerances: relative=1e-05, absolute=1e-50, divergence=10000 left preconditioning using NONE norm type for convergence test PC Object: (fsa_fieldsplit_lambda_sub_) 1 MPI processes type: lu LU: out-of-place factorization tolerance for zero pivot 2.22045e-14 matrix ordering: nd factor fill ratio given 5, needed 1 Factored matrix follows: Mat Object: 1 MPI processes type: seqaij rows=20, cols=20 package used to perform factorization: petsc total: nonzeros=20, allocated nonzeros=20 total number of mallocs used during MatSetValues calls =0 not using I-node routines linear system matrix = precond matrix: Mat Object: (fsa_fieldsplit_lambda_) 1 MPI processes type: seqaij rows=20, cols=20 total: nonzeros=20, allocated nonzeros=20 total number of mallocs used during MatSetValues calls =0 not using I-node routines linear system matrix followed by preconditioner matrix: Mat Object: (fsa_fieldsplit_lambda_) 4 MPI processes type: schurcomplement rows=20, cols=20 Schur complement A11 - A10 inv(A00) A01 A11 Mat Object: (fsa_fieldsplit_lambda_) 4 MPI processes type: mpiaij rows=20, cols=20 total: nonzeros=20, allocated nonzeros=40 total number of mallocs used during MatSetValues calls =0 not using I-node (on process 0) routines A10 Mat Object: 4 MPI processes type: mpiaij rows=20, cols=3705 total: nonzeros=9501, allocated nonzeros=19002 total number of mallocs used during MatSetValues calls =0 using I-node (on process 0) routines: found 12 nodes, limit used is 5 KSP of A00 KSP Object: (fsa_fieldsplit_lambda_inner_) 4 MPI processes type: cg maximum iterations=10000, initial guess is zero tolerances: relative=1e-12, absolute=1e-50, divergence=10000 left preconditioning using PRECONDITIONED norm type for convergence test PC Object: (fsa_fieldsplit_lambda_inner_) 4 MPI processes type: gamg MG: type is MULTIPLICATIVE, levels=2 cycles=v Cycles per PCApply=1 Using Galerkin computed coarse grid matrices Coarse grid solver -- level ------------------------------- KSP Object: (fsa_fieldsplit_lambda_inner_mg_coarse_) 4 MPI processes type: preonly maximum iterations=1, initial guess is zero tolerances: relative=1e-05, absolute=1e-50, divergence=10000 left preconditioning using NONE norm type for convergence test PC Object: (fsa_fieldsplit_lambda_inner_mg_coarse_) 4 MPI processes type: bjacobi block Jacobi: number of blocks = 4 Local solve is same for all blocks, in the following KSP and PC objects: KSP Object: (fsa_fieldsplit_lambda_inner_mg_coarse_sub_) 1 MPI processes type: preonly maximum iterations=1, initial guess is zero tolerances: relative=1e-05, absolute=1e-50, divergence=10000 left preconditioning using NONE norm type for convergence test PC Object: (fsa_fieldsplit_lambda_inner_mg_coarse_sub_) 1 MPI processes type: lu LU: out-of-place factorization tolerance for zero pivot 2.22045e-14 using diagonal shift on blocks to prevent zero pivot [INBLOCKS] matrix ordering: nd factor fill ratio given 5, needed 4.18681 Factored matrix follows: Mat Object: 1 MPI processes type: seqaij rows=298, cols=298 package used to perform factorization: petsc total: nonzeros=13582, allocated nonzeros=13582 total number of mallocs used during MatSetValues calls =0 not using I-node routines linear system matrix = precond matrix: Mat Object: 1 MPI processes type: seqaij rows=298, cols=298 total: nonzeros=3244, allocated nonzeros=3244 total number of mallocs used during MatSetValues calls =0 not using I-node routines linear system matrix = precond matrix: Mat Object: 4 MPI processes type: mpiaij rows=298, cols=298 total: nonzeros=3244, allocated nonzeros=3244 total number of mallocs used during MatSetValues calls =0 not using I-node (on process 0) routines Down solver (pre-smoother) on level 1 ------------------------------- KSP Object: (fsa_fieldsplit_lambda_inner_mg_levels_1_) 4 MPI processes type: chebyshev Chebyshev: eigenvalue estimates: min = 0.147014, max = 3.0873 maximum iterations=2 tolerances: relative=1e-05, absolute=1e-50, divergence=10000 left preconditioning using nonzero initial guess using NONE norm type for convergence test PC Object: (fsa_fieldsplit_lambda_inner_mg_levels_1_) 4 MPI processes type: sor SOR: type = local_symmetric, iterations = 1, local iterations = 1, omega = 1 linear system matrix = precond matrix: Mat Object: (fsa_fieldsplit_phi_) 4 MPI processes type: mpiaij rows=3705, cols=3705 total: nonzeros=20597, allocated nonzeros=555750 total number of mallocs used during MatSetValues calls =0 not using I-node (on process 0) routines Up solver (post-smoother) same as down solver (pre-smoother) linear system matrix = precond matrix: Mat Object: (fsa_fieldsplit_phi_) 4 MPI processes type: mpiaij rows=3705, cols=3705 total: nonzeros=20597, allocated nonzeros=555750 total number of mallocs used during MatSetValues calls =0 not using I-node (on process 0) routines A01 Mat Object: 4 MPI processes type: mpiaij rows=3705, cols=20 total: nonzeros=9303, allocated nonzeros=9303 total number of mallocs used during MatSetValues calls =0 using I-node (on process 0) routines: found 190 nodes, limit used is 5 Mat Object: (fsa_fieldsplit_lambda_) 4 MPI processes type: mpiaij rows=20, cols=20 total: nonzeros=20, allocated nonzeros=40 total number of mallocs used during MatSetValues calls =0 not using I-node (on process 0) routines linear system matrix = precond matrix: Mat Object: 4 MPI processes type: nest rows=3725, cols=3725 Matrix object: type=nest, rows=2, cols=2 MatNest structure: (0,0) : prefix="fsa_fieldsplit_phi_", type=mpiaij, rows=3705, cols=3705 (0,1) : type=mpiaij, rows=3705, cols=20 (1,0) : type=mpiaij, rows=20, cols=3705 (1,1) : prefix="fsa_fieldsplit_lambda_", type=mpiaij, rows=20, cols=20 1 SNES Function norm 1.005118954668e-10 Nonlinear fsa_ solve converged due to CONVERGED_ITS iterations 1 |C phi - | / || = 5.3109E-14, ||= 5.0171E-02, ||/|phi|= 7.8371E-02 step,ratio,# of ion 2 1.007470 3200000 0 SNES Function norm 5.925636187376e+02 Residual norms for fsa_ solve. 0 KSP Residual norm 5.925636187376e+02 Linear fsa_fieldsplit_phi_ solve converged due to CONVERGED_RTOL iterations 19 Residual norms for fsa_fieldsplit_lambda_ solve. 0 KSP Residual norm 5.790987941564e-06 Linear fsa_fieldsplit_lambda_inner_ solve converged due to CONVERGED_RTOL iterations 15 1 KSP Residual norm 2.827834278506e-06 Linear fsa_fieldsplit_lambda_inner_ solve converged due to CONVERGED_RTOL iterations 15 2 KSP Residual norm 1.767762837351e-06 Linear fsa_fieldsplit_lambda_inner_ solve converged due to CONVERGED_RTOL iterations 15 3 KSP Residual norm 1.338411768727e-06 Linear fsa_fieldsplit_lambda_inner_ solve converged due to CONVERGED_RTOL iterations 15 4 KSP Residual norm 1.041579294454e-06 Linear fsa_fieldsplit_lambda_inner_ solve converged due to CONVERGED_RTOL iterations 15 5 KSP Residual norm 9.812686352687e-07 Linear fsa_fieldsplit_lambda_inner_ solve converged due to CONVERGED_RTOL iterations 15 6 KSP Residual norm 7.014959547015e-07 Linear fsa_fieldsplit_lambda_inner_ solve converged due to CONVERGED_RTOL iterations 15 7 KSP Residual norm 6.259175841640e-07 Linear fsa_fieldsplit_lambda_inner_ solve converged due to CONVERGED_RTOL iterations 15 8 KSP Residual norm 5.779742850396e-07 Linear fsa_fieldsplit_lambda_inner_ solve converged due to CONVERGED_RTOL iterations 15 9 KSP Residual norm 5.591462868660e-07 Linear fsa_fieldsplit_lambda_inner_ solve converged due to CONVERGED_RTOL iterations 15 10 KSP Residual norm 5.479546540146e-07 Linear fsa_fieldsplit_lambda_inner_ solve converged due to CONVERGED_RTOL iterations 15 11 KSP Residual norm 5.344499562971e-07 Linear fsa_fieldsplit_lambda_inner_ solve converged due to CONVERGED_RTOL iterations 15 12 KSP Residual norm 5.239386806233e-07 Linear fsa_fieldsplit_lambda_inner_ solve converged due to CONVERGED_RTOL iterations 15 13 KSP Residual norm 4.959956080119e-07 Linear fsa_fieldsplit_lambda_inner_ solve converged due to CONVERGED_RTOL iterations 15 14 KSP Residual norm 4.686268010470e-07 Linear fsa_fieldsplit_lambda_inner_ solve converged due to CONVERGED_RTOL iterations 15 15 KSP Residual norm 4.459275043164e-07 Linear fsa_fieldsplit_lambda_inner_ solve converged due to CONVERGED_RTOL iterations 15 16 KSP Residual norm 4.144536528804e-07 Linear fsa_fieldsplit_lambda_inner_ solve converged due to CONVERGED_RTOL iterations 15 17 KSP Residual norm 3.779721432521e-07 Linear fsa_fieldsplit_lambda_inner_ solve converged due to CONVERGED_RTOL iterations 15 18 KSP Residual norm 3.416289975754e-07 Linear fsa_fieldsplit_lambda_inner_ solve converged due to CONVERGED_RTOL iterations 15 19 KSP Residual norm 3.322368903701e-07 Linear fsa_fieldsplit_lambda_inner_ solve converged due to CONVERGED_RTOL iterations 15 20 KSP Residual norm 5.332676678708e-20 Linear fsa_fieldsplit_lambda_ solve converged due to CONVERGED_RTOL iterations 20 Linear fsa_fieldsplit_phi_ solve converged due to CONVERGED_RTOL iterations 19 1 KSP Residual norm 1.492596808614e-10 Linear fsa_ solve converged due to CONVERGED_RTOL iterations 1 KSP Object:(fsa_) 4 MPI processes type: fgmres GMRES: restart=50, using Classical (unmodified) Gram-Schmidt Orthogonalization with no iterative refinement GMRES: happy breakdown tolerance 1e-30 maximum iterations=3, initial guess is zero tolerances: relative=1e-12, absolute=1e-50, divergence=10000 right preconditioning using UNPRECONDITIONED norm type for convergence test PC Object:(fsa_) 4 MPI processes type: fieldsplit FieldSplit with Schur preconditioner, factorization FULL Preconditioner for the Schur complement formed from A11 Split info: Split number 0 Defined by IS Split number 1 Defined by IS KSP solver for A00 block KSP Object: (fsa_fieldsplit_phi_) 4 MPI processes type: cg maximum iterations=10000, initial guess is zero tolerances: relative=1e-12, absolute=1e-50, divergence=10000 left preconditioning using PRECONDITIONED norm type for convergence test PC Object: (fsa_fieldsplit_phi_) 4 MPI processes type: gamg MG: type is MULTIPLICATIVE, levels=2 cycles=v Cycles per PCApply=1 Using Galerkin computed coarse grid matrices Coarse grid solver -- level ------------------------------- KSP Object: (fsa_fieldsplit_phi_mg_coarse_) 4 MPI processes type: preonly maximum iterations=1, initial guess is zero tolerances: relative=1e-05, absolute=1e-50, divergence=10000 left preconditioning using NONE norm type for convergence test PC Object: (fsa_fieldsplit_phi_mg_coarse_) 4 MPI processes type: bjacobi block Jacobi: number of blocks = 4 Local solve is same for all blocks, in the following KSP and PC objects: KSP Object: (fsa_fieldsplit_phi_mg_coarse_sub_) 1 MPI processes type: preonly maximum iterations=1, initial guess is zero tolerances: relative=1e-05, absolute=1e-50, divergence=10000 left preconditioning using NONE norm type for convergence test PC Object: (fsa_fieldsplit_phi_mg_coarse_sub_) 1 MPI processes type: lu LU: out-of-place factorization tolerance for zero pivot 2.22045e-14 using diagonal shift on blocks to prevent zero pivot [INBLOCKS] matrix ordering: nd factor fill ratio given 5, needed 4.18681 Factored matrix follows: Mat Object: 1 MPI processes type: seqaij rows=298, cols=298 package used to perform factorization: petsc total: nonzeros=13582, allocated nonzeros=13582 total number of mallocs used during MatSetValues calls =0 not using I-node routines linear system matrix = precond matrix: Mat Object: 1 MPI processes type: seqaij rows=298, cols=298 total: nonzeros=3244, allocated nonzeros=3244 total number of mallocs used during MatSetValues calls =0 not using I-node routines linear system matrix = precond matrix: Mat Object: 4 MPI processes type: mpiaij rows=298, cols=298 total: nonzeros=3244, allocated nonzeros=3244 total number of mallocs used during MatSetValues calls =0 not using I-node (on process 0) routines Down solver (pre-smoother) on level 1 ------------------------------- KSP Object: (fsa_fieldsplit_phi_mg_levels_1_) 4 MPI processes type: richardson Richardson: damping factor=1 maximum iterations=2 tolerances: relative=1e-05, absolute=1e-50, divergence=10000 left preconditioning using nonzero initial guess using NONE norm type for convergence test PC Object: (fsa_fieldsplit_phi_mg_levels_1_) 4 MPI processes type: sor SOR: type = local_symmetric, iterations = 1, local iterations = 1, omega = 1 linear system matrix = precond matrix: Mat Object: (fsa_fieldsplit_phi_) 4 MPI processes type: mpiaij rows=3705, cols=3705 total: nonzeros=20597, allocated nonzeros=555750 total number of mallocs used during MatSetValues calls =0 not using I-node (on process 0) routines Up solver (post-smoother) same as down solver (pre-smoother) linear system matrix = precond matrix: Mat Object: (fsa_fieldsplit_phi_) 4 MPI processes type: mpiaij rows=3705, cols=3705 total: nonzeros=20597, allocated nonzeros=555750 total number of mallocs used during MatSetValues calls =0 not using I-node (on process 0) routines KSP solver for S = A11 - A10 inv(A00) A01 KSP Object: (fsa_fieldsplit_lambda_) 4 MPI processes type: gmres GMRES: restart=30, using Classical (unmodified) Gram-Schmidt Orthogonalization with no iterative refinement GMRES: happy breakdown tolerance 1e-30 maximum iterations=10000, initial guess is zero tolerances: relative=1e-12, absolute=1e-50, divergence=10000 left preconditioning using PRECONDITIONED norm type for convergence test PC Object: (fsa_fieldsplit_lambda_) 4 MPI processes type: bjacobi block Jacobi: number of blocks = 4 Local solve is same for all blocks, in the following KSP and PC objects: KSP Object: (fsa_fieldsplit_lambda_sub_) 1 MPI processes type: preonly maximum iterations=10000, initial guess is zero tolerances: relative=1e-05, absolute=1e-50, divergence=10000 left preconditioning using NONE norm type for convergence test PC Object: (fsa_fieldsplit_lambda_sub_) 1 MPI processes type: lu LU: out-of-place factorization tolerance for zero pivot 2.22045e-14 matrix ordering: nd factor fill ratio given 5, needed 1 Factored matrix follows: Mat Object: 1 MPI processes type: seqaij rows=20, cols=20 package used to perform factorization: petsc total: nonzeros=20, allocated nonzeros=20 total number of mallocs used during MatSetValues calls =0 not using I-node routines linear system matrix = precond matrix: Mat Object: (fsa_fieldsplit_lambda_) 1 MPI processes type: seqaij rows=20, cols=20 total: nonzeros=20, allocated nonzeros=20 total number of mallocs used during MatSetValues calls =0 not using I-node routines linear system matrix followed by preconditioner matrix: Mat Object: (fsa_fieldsplit_lambda_) 4 MPI processes type: schurcomplement rows=20, cols=20 Schur complement A11 - A10 inv(A00) A01 A11 Mat Object: (fsa_fieldsplit_lambda_) 4 MPI processes type: mpiaij rows=20, cols=20 total: nonzeros=20, allocated nonzeros=40 total number of mallocs used during MatSetValues calls =0 not using I-node (on process 0) routines A10 Mat Object: 4 MPI processes type: mpiaij rows=20, cols=3705 total: nonzeros=9501, allocated nonzeros=19002 total number of mallocs used during MatSetValues calls =0 using I-node (on process 0) routines: found 12 nodes, limit used is 5 KSP of A00 KSP Object: (fsa_fieldsplit_lambda_inner_) 4 MPI processes type: cg maximum iterations=10000, initial guess is zero tolerances: relative=1e-12, absolute=1e-50, divergence=10000 left preconditioning using PRECONDITIONED norm type for convergence test PC Object: (fsa_fieldsplit_lambda_inner_) 4 MPI processes type: gamg MG: type is MULTIPLICATIVE, levels=2 cycles=v Cycles per PCApply=1 Using Galerkin computed coarse grid matrices Coarse grid solver -- level ------------------------------- KSP Object: (fsa_fieldsplit_lambda_inner_mg_coarse_) 4 MPI processes type: preonly maximum iterations=1, initial guess is zero tolerances: relative=1e-05, absolute=1e-50, divergence=10000 left preconditioning using NONE norm type for convergence test PC Object: (fsa_fieldsplit_lambda_inner_mg_coarse_) 4 MPI processes type: bjacobi block Jacobi: number of blocks = 4 Local solve is same for all blocks, in the following KSP and PC objects: KSP Object: (fsa_fieldsplit_lambda_inner_mg_coarse_sub_) 1 MPI processes type: preonly maximum iterations=1, initial guess is zero tolerances: relative=1e-05, absolute=1e-50, divergence=10000 left preconditioning using NONE norm type for convergence test PC Object: (fsa_fieldsplit_lambda_inner_mg_coarse_sub_) 1 MPI processes type: lu LU: out-of-place factorization tolerance for zero pivot 2.22045e-14 using diagonal shift on blocks to prevent zero pivot [INBLOCKS] matrix ordering: nd factor fill ratio given 5, needed 4.18681 Factored matrix follows: Mat Object: 1 MPI processes type: seqaij rows=298, cols=298 package used to perform factorization: petsc total: nonzeros=13582, allocated nonzeros=13582 total number of mallocs used during MatSetValues calls =0 not using I-node routines linear system matrix = precond matrix: Mat Object: 1 MPI processes type: seqaij rows=298, cols=298 total: nonzeros=3244, allocated nonzeros=3244 total number of mallocs used during MatSetValues calls =0 not using I-node routines linear system matrix = precond matrix: Mat Object: 4 MPI processes type: mpiaij rows=298, cols=298 total: nonzeros=3244, allocated nonzeros=3244 total number of mallocs used during MatSetValues calls =0 not using I-node (on process 0) routines Down solver (pre-smoother) on level 1 ------------------------------- KSP Object: (fsa_fieldsplit_lambda_inner_mg_levels_1_) 4 MPI processes type: chebyshev Chebyshev: eigenvalue estimates: min = 0.147014, max = 3.0873 maximum iterations=2 tolerances: relative=1e-05, absolute=1e-50, divergence=10000 left preconditioning using nonzero initial guess using NONE norm type for convergence test PC Object: (fsa_fieldsplit_lambda_inner_mg_levels_1_) 4 MPI processes type: sor SOR: type = local_symmetric, iterations = 1, local iterations = 1, omega = 1 linear system matrix = precond matrix: Mat Object: (fsa_fieldsplit_phi_) 4 MPI processes type: mpiaij rows=3705, cols=3705 total: nonzeros=20597, allocated nonzeros=555750 total number of mallocs used during MatSetValues calls =0 not using I-node (on process 0) routines Up solver (post-smoother) same as down solver (pre-smoother) linear system matrix = precond matrix: Mat Object: (fsa_fieldsplit_phi_) 4 MPI processes type: mpiaij rows=3705, cols=3705 total: nonzeros=20597, allocated nonzeros=555750 total number of mallocs used during MatSetValues calls =0 not using I-node (on process 0) routines A01 Mat Object: 4 MPI processes type: mpiaij rows=3705, cols=20 total: nonzeros=9303, allocated nonzeros=9303 total number of mallocs used during MatSetValues calls =0 using I-node (on process 0) routines: found 190 nodes, limit used is 5 Mat Object: (fsa_fieldsplit_lambda_) 4 MPI processes type: mpiaij rows=20, cols=20 total: nonzeros=20, allocated nonzeros=40 total number of mallocs used during MatSetValues calls =0 not using I-node (on process 0) routines linear system matrix = precond matrix: Mat Object: 4 MPI processes type: nest rows=3725, cols=3725 Matrix object: type=nest, rows=2, cols=2 MatNest structure: (0,0) : prefix="fsa_fieldsplit_phi_", type=mpiaij, rows=3705, cols=3705 (0,1) : type=mpiaij, rows=3705, cols=20 (1,0) : type=mpiaij, rows=20, cols=3705 (1,1) : prefix="fsa_fieldsplit_lambda_", type=mpiaij, rows=20, cols=20 1 SNES Function norm 1.492618104103e-10 Nonlinear fsa_ solve converged due to CONVERGED_ITS iterations 1 |C phi - | / || = 7.4828E-14, ||= 6.0739E-02, ||/|phi|= 4.9042E-02 2d particle and heat flux is not ready yet for deltaf-f method 0 SNES Function norm 8.885599193467e+02 Residual norms for fsa_ solve. 0 KSP Residual norm 8.885599193467e+02 Linear fsa_fieldsplit_phi_ solve converged due to CONVERGED_RTOL iterations 19 Residual norms for fsa_fieldsplit_lambda_ solve. 0 KSP Residual norm 6.918218162709e-06 Linear fsa_fieldsplit_lambda_inner_ solve converged due to CONVERGED_RTOL iterations 15 1 KSP Residual norm 3.308907521504e-06 Linear fsa_fieldsplit_lambda_inner_ solve converged due to CONVERGED_RTOL iterations 15 2 KSP Residual norm 1.991826495532e-06 Linear fsa_fieldsplit_lambda_inner_ solve converged due to CONVERGED_RTOL iterations 15 3 KSP Residual norm 1.418506151180e-06 Linear fsa_fieldsplit_lambda_inner_ solve converged due to CONVERGED_RTOL iterations 15 4 KSP Residual norm 1.037197799935e-06 Linear fsa_fieldsplit_lambda_inner_ solve converged due to CONVERGED_RTOL iterations 15 5 KSP Residual norm 9.898224721600e-07 Linear fsa_fieldsplit_lambda_inner_ solve converged due to CONVERGED_RTOL iterations 15 6 KSP Residual norm 6.996816024385e-07 Linear fsa_fieldsplit_lambda_inner_ solve converged due to CONVERGED_RTOL iterations 15 7 KSP Residual norm 5.103453916167e-07 Linear fsa_fieldsplit_lambda_inner_ solve converged due to CONVERGED_RTOL iterations 15 8 KSP Residual norm 3.595059255268e-07 Linear fsa_fieldsplit_lambda_inner_ solve converged due to CONVERGED_RTOL iterations 15 9 KSP Residual norm 2.778444546775e-07 Linear fsa_fieldsplit_lambda_inner_ solve converged due to CONVERGED_RTOL iterations 15 10 KSP Residual norm 2.354855449662e-07 Linear fsa_fieldsplit_lambda_inner_ solve converged due to CONVERGED_RTOL iterations 15 11 KSP Residual norm 2.123161982793e-07 Linear fsa_fieldsplit_lambda_inner_ solve converged due to CONVERGED_RTOL iterations 15 12 KSP Residual norm 1.973410799160e-07 Linear fsa_fieldsplit_lambda_inner_ solve converged due to CONVERGED_RTOL iterations 15 13 KSP Residual norm 1.405394379926e-07 Linear fsa_fieldsplit_lambda_inner_ solve converged due to CONVERGED_RTOL iterations 15 14 KSP Residual norm 1.008579191914e-07 Linear fsa_fieldsplit_lambda_inner_ solve converged due to CONVERGED_RTOL iterations 15 15 KSP Residual norm 9.376877553657e-08 Linear fsa_fieldsplit_lambda_inner_ solve converged due to CONVERGED_RTOL iterations 15 16 KSP Residual norm 9.032230495680e-08 Linear fsa_fieldsplit_lambda_inner_ solve converged due to CONVERGED_RTOL iterations 15 17 KSP Residual norm 8.904015135877e-08 Linear fsa_fieldsplit_lambda_inner_ solve converged due to CONVERGED_RTOL iterations 15 18 KSP Residual norm 8.112651639997e-08 Linear fsa_fieldsplit_lambda_inner_ solve converged due to CONVERGED_RTOL iterations 15 19 KSP Residual norm 5.301105618499e-08 Linear fsa_fieldsplit_lambda_inner_ solve converged due to CONVERGED_RTOL iterations 15 20 KSP Residual norm 1.543158764519e-19 Linear fsa_fieldsplit_lambda_ solve converged due to CONVERGED_RTOL iterations 20 Linear fsa_fieldsplit_phi_ solve converged due to CONVERGED_RTOL iterations 19 1 KSP Residual norm 1.574912424389e-10 Linear fsa_ solve converged due to CONVERGED_RTOL iterations 1 KSP Object:(fsa_) 4 MPI processes type: fgmres GMRES: restart=50, using Classical (unmodified) Gram-Schmidt Orthogonalization with no iterative refinement GMRES: happy breakdown tolerance 1e-30 maximum iterations=3, initial guess is zero tolerances: relative=1e-12, absolute=1e-50, divergence=10000 right preconditioning using UNPRECONDITIONED norm type for convergence test PC Object:(fsa_) 4 MPI processes type: fieldsplit FieldSplit with Schur preconditioner, factorization FULL Preconditioner for the Schur complement formed from A11 Split info: Split number 0 Defined by IS Split number 1 Defined by IS KSP solver for A00 block KSP Object: (fsa_fieldsplit_phi_) 4 MPI processes type: cg maximum iterations=10000, initial guess is zero tolerances: relative=1e-12, absolute=1e-50, divergence=10000 left preconditioning using PRECONDITIONED norm type for convergence test PC Object: (fsa_fieldsplit_phi_) 4 MPI processes type: gamg MG: type is MULTIPLICATIVE, levels=2 cycles=v Cycles per PCApply=1 Using Galerkin computed coarse grid matrices Coarse grid solver -- level ------------------------------- KSP Object: (fsa_fieldsplit_phi_mg_coarse_) 4 MPI processes type: preonly maximum iterations=1, initial guess is zero tolerances: relative=1e-05, absolute=1e-50, divergence=10000 left preconditioning using NONE norm type for convergence test PC Object: (fsa_fieldsplit_phi_mg_coarse_) 4 MPI processes type: bjacobi block Jacobi: number of blocks = 4 Local solve is same for all blocks, in the following KSP and PC objects: KSP Object: (fsa_fieldsplit_phi_mg_coarse_sub_) 1 MPI processes type: preonly maximum iterations=1, initial guess is zero tolerances: relative=1e-05, absolute=1e-50, divergence=10000 left preconditioning using NONE norm type for convergence test PC Object: (fsa_fieldsplit_phi_mg_coarse_sub_) 1 MPI processes type: lu LU: out-of-place factorization tolerance for zero pivot 2.22045e-14 using diagonal shift on blocks to prevent zero pivot [INBLOCKS] matrix ordering: nd factor fill ratio given 5, needed 4.18681 Factored matrix follows: Mat Object: 1 MPI processes type: seqaij rows=298, cols=298 package used to perform factorization: petsc total: nonzeros=13582, allocated nonzeros=13582 total number of mallocs used during MatSetValues calls =0 not using I-node routines linear system matrix = precond matrix: Mat Object: 1 MPI processes type: seqaij rows=298, cols=298 total: nonzeros=3244, allocated nonzeros=3244 total number of mallocs used during MatSetValues calls =0 not using I-node routines linear system matrix = precond matrix: Mat Object: 4 MPI processes type: mpiaij rows=298, cols=298 total: nonzeros=3244, allocated nonzeros=3244 total number of mallocs used during MatSetValues calls =0 not using I-node (on process 0) routines Down solver (pre-smoother) on level 1 ------------------------------- KSP Object: (fsa_fieldsplit_phi_mg_levels_1_) 4 MPI processes type: richardson Richardson: damping factor=1 maximum iterations=2 tolerances: relative=1e-05, absolute=1e-50, divergence=10000 left preconditioning using nonzero initial guess using NONE norm type for convergence test PC Object: (fsa_fieldsplit_phi_mg_levels_1_) 4 MPI processes type: sor SOR: type = local_symmetric, iterations = 1, local iterations = 1, omega = 1 linear system matrix = precond matrix: Mat Object: (fsa_fieldsplit_phi_) 4 MPI processes type: mpiaij rows=3705, cols=3705 total: nonzeros=20597, allocated nonzeros=555750 total number of mallocs used during MatSetValues calls =0 not using I-node (on process 0) routines Up solver (post-smoother) same as down solver (pre-smoother) linear system matrix = precond matrix: Mat Object: (fsa_fieldsplit_phi_) 4 MPI processes type: mpiaij rows=3705, cols=3705 total: nonzeros=20597, allocated nonzeros=555750 total number of mallocs used during MatSetValues calls =0 not using I-node (on process 0) routines KSP solver for S = A11 - A10 inv(A00) A01 KSP Object: (fsa_fieldsplit_lambda_) 4 MPI processes type: gmres GMRES: restart=30, using Classical (unmodified) Gram-Schmidt Orthogonalization with no iterative refinement GMRES: happy breakdown tolerance 1e-30 maximum iterations=10000, initial guess is zero tolerances: relative=1e-12, absolute=1e-50, divergence=10000 left preconditioning using PRECONDITIONED norm type for convergence test PC Object: (fsa_fieldsplit_lambda_) 4 MPI processes type: bjacobi block Jacobi: number of blocks = 4 Local solve is same for all blocks, in the following KSP and PC objects: KSP Object: (fsa_fieldsplit_lambda_sub_) 1 MPI processes type: preonly maximum iterations=10000, initial guess is zero tolerances: relative=1e-05, absolute=1e-50, divergence=10000 left preconditioning using NONE norm type for convergence test PC Object: (fsa_fieldsplit_lambda_sub_) 1 MPI processes type: lu LU: out-of-place factorization tolerance for zero pivot 2.22045e-14 matrix ordering: nd factor fill ratio given 5, needed 1 Factored matrix follows: Mat Object: 1 MPI processes type: seqaij rows=20, cols=20 package used to perform factorization: petsc total: nonzeros=20, allocated nonzeros=20 total number of mallocs used during MatSetValues calls =0 not using I-node routines linear system matrix = precond matrix: Mat Object: (fsa_fieldsplit_lambda_) 1 MPI processes type: seqaij rows=20, cols=20 total: nonzeros=20, allocated nonzeros=20 total number of mallocs used during MatSetValues calls =0 not using I-node routines linear system matrix followed by preconditioner matrix: Mat Object: (fsa_fieldsplit_lambda_) 4 MPI processes type: schurcomplement rows=20, cols=20 Schur complement A11 - A10 inv(A00) A01 A11 Mat Object: (fsa_fieldsplit_lambda_) 4 MPI processes type: mpiaij rows=20, cols=20 total: nonzeros=20, allocated nonzeros=40 total number of mallocs used during MatSetValues calls =0 not using I-node (on process 0) routines A10 Mat Object: 4 MPI processes type: mpiaij rows=20, cols=3705 total: nonzeros=9501, allocated nonzeros=19002 total number of mallocs used during MatSetValues calls =0 using I-node (on process 0) routines: found 12 nodes, limit used is 5 KSP of A00 KSP Object: (fsa_fieldsplit_lambda_inner_) 4 MPI processes type: cg maximum iterations=10000, initial guess is zero tolerances: relative=1e-12, absolute=1e-50, divergence=10000 left preconditioning using PRECONDITIONED norm type for convergence test PC Object: (fsa_fieldsplit_lambda_inner_) 4 MPI processes type: gamg MG: type is MULTIPLICATIVE, levels=2 cycles=v Cycles per PCApply=1 Using Galerkin computed coarse grid matrices Coarse grid solver -- level ------------------------------- KSP Object: (fsa_fieldsplit_lambda_inner_mg_coarse_) 4 MPI processes type: preonly maximum iterations=1, initial guess is zero tolerances: relative=1e-05, absolute=1e-50, divergence=10000 left preconditioning using NONE norm type for convergence test PC Object: (fsa_fieldsplit_lambda_inner_mg_coarse_) 4 MPI processes type: bjacobi block Jacobi: number of blocks = 4 Local solve is same for all blocks, in the following KSP and PC objects: KSP Object: (fsa_fieldsplit_lambda_inner_mg_coarse_sub_) 1 MPI processes type: preonly maximum iterations=1, initial guess is zero tolerances: relative=1e-05, absolute=1e-50, divergence=10000 left preconditioning using NONE norm type for convergence test PC Object: (fsa_fieldsplit_lambda_inner_mg_coarse_sub_) 1 MPI processes type: lu LU: out-of-place factorization tolerance for zero pivot 2.22045e-14 using diagonal shift on blocks to prevent zero pivot [INBLOCKS] matrix ordering: nd factor fill ratio given 5, needed 4.18681 Factored matrix follows: Mat Object: 1 MPI processes type: seqaij rows=298, cols=298 package used to perform factorization: petsc total: nonzeros=13582, allocated nonzeros=13582 total number of mallocs used during MatSetValues calls =0 not using I-node routines linear system matrix = precond matrix: Mat Object: 1 MPI processes type: seqaij rows=298, cols=298 total: nonzeros=3244, allocated nonzeros=3244 total number of mallocs used during MatSetValues calls =0 not using I-node routines linear system matrix = precond matrix: Mat Object: 4 MPI processes type: mpiaij rows=298, cols=298 total: nonzeros=3244, allocated nonzeros=3244 total number of mallocs used during MatSetValues calls =0 not using I-node (on process 0) routines Down solver (pre-smoother) on level 1 ------------------------------- KSP Object: (fsa_fieldsplit_lambda_inner_mg_levels_1_) 4 MPI processes type: chebyshev Chebyshev: eigenvalue estimates: min = 0.147014, max = 3.0873 maximum iterations=2 tolerances: relative=1e-05, absolute=1e-50, divergence=10000 left preconditioning using nonzero initial guess using NONE norm type for convergence test PC Object: (fsa_fieldsplit_lambda_inner_mg_levels_1_) 4 MPI processes type: sor SOR: type = local_symmetric, iterations = 1, local iterations = 1, omega = 1 linear system matrix = precond matrix: Mat Object: (fsa_fieldsplit_phi_) 4 MPI processes type: mpiaij rows=3705, cols=3705 total: nonzeros=20597, allocated nonzeros=555750 total number of mallocs used during MatSetValues calls =0 not using I-node (on process 0) routines Up solver (post-smoother) same as down solver (pre-smoother) linear system matrix = precond matrix: Mat Object: (fsa_fieldsplit_phi_) 4 MPI processes type: mpiaij rows=3705, cols=3705 total: nonzeros=20597, allocated nonzeros=555750 total number of mallocs used during MatSetValues calls =0 not using I-node (on process 0) routines A01 Mat Object: 4 MPI processes type: mpiaij rows=3705, cols=20 total: nonzeros=9303, allocated nonzeros=9303 total number of mallocs used during MatSetValues calls =0 using I-node (on process 0) routines: found 190 nodes, limit used is 5 Mat Object: (fsa_fieldsplit_lambda_) 4 MPI processes type: mpiaij rows=20, cols=20 total: nonzeros=20, allocated nonzeros=40 total number of mallocs used during MatSetValues calls =0 not using I-node (on process 0) routines linear system matrix = precond matrix: Mat Object: 4 MPI processes type: nest rows=3725, cols=3725 Matrix object: type=nest, rows=2, cols=2 MatNest structure: (0,0) : prefix="fsa_fieldsplit_phi_", type=mpiaij, rows=3705, cols=3705 (0,1) : type=mpiaij, rows=3705, cols=20 (1,0) : type=mpiaij, rows=20, cols=3705 (1,1) : prefix="fsa_fieldsplit_lambda_", type=mpiaij, rows=20, cols=20 1 SNES Function norm 1.574897278577e-10 Nonlinear fsa_ solve converged due to CONVERGED_ITS iterations 1 |C phi - | / || = 1.4846E-13, ||= 3.5125E-02, ||/|phi|= 1.9588E-02 Application 6772851 resources: utime ~287s, stime ~11s, Rss ~28944, inblocks ~126888, outblocks ~231479