<div dir="ltr">Oh, and I use:<div><br></div><div>-fieldsplit_e_pc_type lu<br>-fieldsplit_i1_pc_type lu<br></div><div><br></div><div>and </div><div><br></div><div>-fieldsplit_pc_type lu<br></div><div><br></div><div>does not seem to work.</div><div><br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Thu, Dec 31, 2020 at 4:25 PM Mark Adams <<a href="mailto:mfadams@lbl.gov">mfadams@lbl.gov</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div dir="ltr"><div>Ha! This seems to work. It Iooks good as far as I can tell. I don't think this is just a direct solver.</div><div>Thanks</div><div><br></div><div> ....</div><div> 3 SNES Function norm 1.557752448861e-11<br><b> 0 KSP Residual norm 1.557752448861e-11<br> 1 KSP Residual norm 7.987154898471e-27</b><br>KSP Object: 1 MPI processes<br> type: preonly<br> maximum iterations=10000, initial guess is zero<br> tolerances: relative=1e-05, absolute=1e-50, divergence=10000.<br> left preconditioning<br> using NONE norm type for convergence test<br>PC Object: 1 MPI processes<br> type: fieldsplit<br><b> FieldSplit with ADDITIVE composition: total splits = 2<br></b> Solver info for each split is in the following KSP objects:<br> Split number 0 Defined by IS<br> KSP Object: (fieldsplit_e_) 1 MPI processes<br> type: preonly<br> maximum iterations=10000, initial guess is zero<br> tolerances: relative=1e-05, absolute=1e-50, divergence=10000.<br> left preconditioning<br> using NONE norm type for convergence test<br> PC Object: (fieldsplit_e_) 1 MPI processes<br> type: lu<br> out-of-place factorization<br> tolerance for zero pivot 2.22045e-14<br> matrix ordering: nd<br> factor fill ratio given 5., needed 1.30805<br> Factored matrix follows:<br> Mat Object: 1 MPI processes<br> type: seqaij<br> rows=448, cols=448<br> package used to perform factorization: petsc<br> total: nonzeros=14038, allocated nonzeros=14038<br> using I-node routines: found 175 nodes, limit used is 5<br> linear system matrix = precond matrix:<br> Mat Object: (fieldsplit_e_) 1 MPI processes<br> type: seqaij<br> rows=448, cols=448<br> total: nonzeros=10732, allocated nonzeros=10732<br> total number of mallocs used during MatSetValues calls=0<br> using I-node routines: found 197 nodes, limit used is 5<br> Split number 1 Defined by IS<br> KSP Object: (fieldsplit_i1_) 1 MPI processes<br> type: preonly<br> maximum iterations=10000, initial guess is zero<br> tolerances: relative=1e-05, absolute=1e-50, divergence=10000.<br> left preconditioning<br> using NONE norm type for convergence test<br> PC Object: (fieldsplit_i1_) 1 MPI processes<br> type: lu<br> out-of-place factorization<br> tolerance for zero pivot 2.22045e-14<br> matrix ordering: nd<br> factor fill ratio given 5., needed 1.30805<br> Factored matrix follows:<br> Mat Object: 1 MPI processes<br> type: seqaij<br> rows=448, cols=448<br> package used to perform factorization: petsc<br> total: nonzeros=14038, allocated nonzeros=14038<br> using I-node routines: found 175 nodes, limit used is 5<br> linear system matrix = precond matrix:<br> Mat Object: (fieldsplit_i1_) 1 MPI processes<br> type: seqaij<br> rows=448, cols=448<br> total: nonzeros=10732, allocated nonzeros=10732<br> total number of mallocs used during MatSetValues calls=0<br> using I-node routines: found 197 nodes, limit used is 5<br> linear system matrix = precond matrix:<br> Mat Object: 1 MPI processes<br> type: seqaij<br> rows=896, cols=896<br> total: nonzeros=21464, allocated nonzeros=42928<br> total number of mallocs used during MatSetValues calls=0<br> using I-node routines: found 398 nodes, limit used is 5<br> 4 SNES Function norm 5.154807006139e-14<br> Nonlinear solve converged due to CONVERGED_SNORM_RELATIVE iterations 4<br> TSAdapt basic arkimex 0:1bee step 0 accepted t=0 + 1.000e-02 dt=1.100e-02 wlte=0.00161 wltea= -1 wlter= -1<br>1 TS dt 0.011 time 0.01<br> 1) species-0: charge density= -1.6022862392985e+01 z-momentum= 1.9463853500826e-19 energy= 9.6063874253791e+04<br> 1) species-1: charge density= 1.6029890912678e+01 z-momentum= -1.6132625562313e-18 energy= 9.6333617535820e+04<br> 1) Total: charge density= 7.0285196929696e-03, momentum= -1.4186240212230e-18, energy= 1.9239749178961e+05 (m_i[0]/m_e = 1835.47, 50 cells)<br>testSpitzer 1) time= 1.000e-02 n_e= 1.000e+00 E= 0.000e+00 J= 5.338e-08 J_re= -0.000e+00 -0 % Te_kev= 3.995e+00 Z_eff=1. E/J to eta ratio=0. (diff=-0.) constant E<br>[0] parallel consistency check OK<br>#PETSc Option Table entries:<br>-dm_landau_amr_levels_max 8<br>-dm_landau_amr_post_refine 0<br>-dm_landau_device_type cpu<br>-dm_landau_domain_radius 5<br>-dm_landau_Ez 0<br>-dm_landau_ion_charges 1<br>-dm_landau_ion_masses 1<br>-dm_landau_n 1,1<br>-dm_landau_thermal_temps 4,4<br>-dm_landau_type p4est<br>-dm_preallocate_only<br>-ex2_connor_e_field_units<br>-ex2_impurity_source_type pulse<br>-ex2_plot_dt 1<br>-ex2_pulse_rate 2e+0<br>-ex2_pulse_start_time 32<br>-ex2_pulse_width_time 15<br>-ex2_t_cold 1<br>-ex2_test_type spitzer<br><b>-fieldsplit_e_pc_type lu<br>-fieldsplit_i1_pc_type lu</b><br>-info :dm<br>-ksp_monitor<br>-ksp_type preonly<br>-ksp_view<br>-options_left<br><b>-pc_fieldsplit_type additive<br>-pc_type fieldsplit</b><br>-petscspace_degree 3<br>-snes_converged_reason<br>-snes_max_it 15<br>-snes_monitor<br>-snes_rtol 1.e-14<br>-snes_stol 1.e-14<br>-ts_adapt_clip .25,1.1<br>-ts_adapt_dt_max 1.<br>-ts_adapt_dt_min 1e-4<br>-ts_adapt_monitor<br>-ts_adapt_scale_solve_failed 0.75<br>-ts_adapt_time_step_increase_delay 5<br>-ts_arkimex_type 1bee<br>-ts_dt .1e-1<br>-ts_exact_final_time stepover<br>-ts_max_snes_failures -1<br>-ts_max_steps 1<br>-ts_max_time 100<br>-ts_monitor<br>-ts_rtol 1e-3<br>-ts_type arkimex<br>#End of PETSc Option Table entries<br>There are no unused options.<br>(base) 16:19 master= ~/Codes/petsc/src/ts/utils/dmplexlandau/tutorials$<br></div><div><br></div></div></div>
</blockquote></div>