<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>