Interfaces and simulations in SELALIB

Transcription

1 Interfaces and simulations in SELALIB M. Mehrenberger IRMA, Université de Strasbourg Garching, Selalib Day, November 2014 M. Mehrenberger (UDS) Interfaces and simulations in SELALIB Garching, November

2 Introduction Outline 0. Discussion on interfaces 1. Some examples of abstract interfaces Interpolation Characteristics Advection Poisson solvers Simulation 2. Other interfaces (not abstract) Initial function Remap 3. A specific example of abstract class 4. About simulations M. Mehrenberger (UDS) Interfaces and simulations in SELALIB Garching, November

3 Interfaces 0. Discussions on interfaces Use of Fortran 2003 capabilities is one of the most flexible way to define interfaces Define rules that are common Each people can focus on a specific implementation Main code remains unchanged How to propose a new abstract interface? Give mathematical description of the interface Give proposed implementation of the common interface (_base) Send the proposition to selalib work Finding the right interface takes time and can result from discussion Other interfaces (non fortran2003) are sometimes enough do not abuse with fortran2003 sometimes, first develop a type do then a wrapper to abstract, when everything is ready and works M. Mehrenberger (UDS) Interfaces and simulations in SELALIB Garching, November

4 1. Examples of abstract interfaces Examples Some examples are developed Other examples not developed here (but of interest!) Coordinate transformations Mesh Scalar fields Sparse matrix Multi-patch PIC M. Mehrenberger (UDS) Interfaces and simulations in SELALIB Garching, November

5 1. Examples of abstract interfaces INTERPOLATION Mathematical description Gives the value of f (η). f is function known on some grid points. η should be inside the domain. η can be a single value η R or several values η R N. Precomputations are done first, from known grid points values. Examples : splines, Hermite, Lagrange. Extension in 2D : η = (η 1, η 2 ) R N 1 N 2 and f (η) R N 1 N 2. M. Mehrenberger (UDS) Interfaces and simulations in SELALIB Garching, November

6 1. Examples of abstract interfaces INTERPOLATION Suggested implementation f u n c t i o n interpolate_value ( & interp, & eta ) & r e s u l t ( val ) sll_real64 : : val class ( sll_interpolator_1d_base ), i n t e n t ( i n ) : : interp sll_real64, i n t e n t ( i n ) : : eta & For info, usage is η R = eta, f (η) = val val= interp%interpolate_value ( eta ) M. Mehrenberger (UDS) Interfaces and simulations in SELALIB Garching, November

7 1. Examples of abstract interfaces INTERPOLATION Suggested implementation subroutine interpolate_array_values ( & interp, & num_pts, & eta, & val ) class ( sll_interpolator_1d_base ), & i n t e n t ( i n ) : : interp sll_int32, i n t e n t ( i n ) : : num_pts sll_real64, dimension ( : ), i n t e n t ( i n ) : : eta sll_real64, dimension ( : ), i n t e n t ( out ) : : val η R N = eta(1:num_pts), f (η) R N = val(1:num_pts) M. Mehrenberger (UDS) Interfaces and simulations in SELALIB Garching, November

8 1. Examples of abstract interfaces INTERPOLATION Suggested implementation f u n c t i o n compute_interpolants ( & interp, & input ) class ( sll_interpolator_1d_base ), & i n t e n t ( i n o u t ) : : interp sll_real64, dimension ( : ), i n t e n t ( i n ) : : input f at grid points = input M. Mehrenberger (UDS) Interfaces and simulations in SELALIB Garching, November

9 1. Examples of abstract interfaces INTERPOLATION Suggested implementation Some changes for extension in 2D class ( sll_interpolator_2d_base ) sll_real64, dimension ( :, : ), i n t e n t ( i n ) : : input sll_int32, i n t e n t ( i n ) : : num_pts1 sll_int32, i n t e n t ( i n ) : : num_pts2 sll_real64, dimension ( :, : ), i n t e n t ( i n ) : : eta1 sll_real64, dimension ( :, : ), i n t e n t ( i n ) : : eta2 sll_real64, dimension ( :, : ), i n t e n t ( out ) : : val η 1 R N 1 N 2 = eta1(1:num_pts1,1:num_pts2), η 2 R N 1 N 2 = eta2(1:num_pts1,1:num_pts2), f (η 1, η 2 ) R N 1 N 2 = val(1:num_pts1,1:num_pts2) M. Mehrenberger (UDS) Interfaces and simulations in SELALIB Garching, November

10 Exemple of use 1. Examples of abstract interfaces INTERPOLATION class ( sll_interpolator_2d_base ) sll_real64, dimension ( :, : ), i n t e n t ( i n ) : : input sll_int32, i n t e n t ( i n ) : : num_pts1 sll_int32, i n t e n t ( i n ) : : num_pts2 sll_real64, dimension ( :, : ), i n t e n t ( i n ) : : eta1 sll_real64, dimension ( :, : ), i n t e n t ( i n ) : : eta2 sll_real64, dimension ( :, : ), i n t e n t ( out ) : : val M. Mehrenberger (UDS) Interfaces and simulations in SELALIB Garching, November

11 1. Examples of abstract interfaces CHARACTERISTICS Characteristics We solve η (t) = A(η(t)) We suppose to know η( t) and want to compute η(0). Examples are : explicit Euler, implicit Euler, Verlet. In 2D, A = (A 1, A 2 ), with A 1 and A 2 known on a grid In 2D, η = (η 1, η 2 ). We suppose that η 1 ( t) and η 2 ( t) are on a 1D grid. M. Mehrenberger (UDS) Interfaces and simulations in SELALIB Garching, November

12 1. Examples of abstract interfaces CHARACTERISTICS compute_characteristics subroutine compute_characteristics ( & charac, & A, & dt, & input, & output ) class ( sll_characteristics_1d_base ) : : charac sll_real64, dimension ( : ), i n t e n t ( i n ) : : A sll_real64, i n t e n t ( i n ) : : dt sll_real64, dimension ( : ), i n t e n t ( i n ) : : input sll_real64, dimension ( : ), i n t e n t ( out ) : : output A = A, t = dt, η( t) = input, η(0) = output M. Mehrenberger (UDS) Interfaces and simulations in SELALIB Garching, November

13 Extension in 2D 1. Examples of abstract interfaces CHARACTERISTICS subroutine compute_characteristics (... ) class ( sll_characteristics_2d_base ) : : charac sll_real64, dimension ( :, : ), i n t e n t ( i n ) : : A1 sll_real64, dimension ( :, : ), i n t e n t ( i n ) : : A2 sll_real64, i n t e n t ( i n ) : : dt sll_real64, dimension ( : ), i n t e n t ( i n ) : : input1 sll_real64, dimension ( : ), i n t e n t ( i n ) : : input2 sll_real64, dimension ( :, : ), i n t e n t ( out ) : : output1 sll_real64, dimension ( :, : ), i n t e n t ( out ) : : output2 A 1 = A1, A 2 = A2, t = dt, η 1 ( t) = input1, η 2 ( t) = input2, η 1 (0) = output1, η 2 (0) = output2 M. Mehrenberger (UDS) Interfaces and simulations in SELALIB Garching, November

14 1. Examples of abstract interfaces ADVECTION Advection We solve t f + A x f = 0, over a time interval t. A does not depend on the variable t. We suppose to know f (t = 0), on a grid and then compute f (t = t) on the same grid. In 2D, A = (A 1, A 2 ) with A 1 and A 2 known on a grid. Advection can be initialized with characteristics and interpolation. Examples : BSL, CSL, FSL M. Mehrenberger (UDS) Interfaces and simulations in SELALIB Garching, November

15 1. Examples of abstract interfaces ADVECTION Constant advection : advect_1d_constant For constant advection, we use the following interface subroutine advect_1d_constant (& adv, & A, & dt, & input, & output ) class ( sll_advection_1d_base ) : : adv sll_real64, i n t e n t ( i n ) : : A sll_real64, i n t e n t ( i n ) : : dt sll_real64, dimension ( : ), i n t e n t ( i n ) : : input sll_real64, dimension ( : ), i n t e n t ( out ) : : output A R = A, t = dt, f (t = 0) = input, f (t = t) = output M. Mehrenberger (UDS) Interfaces and simulations in SELALIB Garching, November

16 advect 2d 1. Examples of abstract interfaces ADVECTION subroutine advect_2d (... ) class ( sll_advection_2d_base ) : : adv sll_real64, dimension ( :, : ), i n t e n t ( i n ) : : A1 sll_real64, dimension ( :, : ), i n t e n t ( i n ) : : A2 sll_real64, i n t e n t ( i n ) : : dt sll_real64, dimension ( :, : ), i n t e n t ( i n ) : : input sll_real64, dimension ( :, : ), i n t e n t ( out ) : : output A 1 = A1, A 2 = A2, t = dt, f (t = 0) = input, f (t = t) = output M. Mehrenberger (UDS) Interfaces and simulations in SELALIB Garching, November

17 1. Examples of abstract interfaces POISSON SOLVERS Poisson solvers compute potential Φ from charge density ρ or compute electric field E from charge density ρ in 2D, E = (E 1, E 2 ). ρ, Φ, E are on grid points. Examples : periodic, polar M. Mehrenberger (UDS) Interfaces and simulations in SELALIB Garching, November

18 In 1D 1. Examples of abstract interfaces POISSON SOLVERS subroutine compute_phi_from_rho (... ) class ( sll_poisson_1d_base ) : : poisson sll_real64, dimension ( : ), i n t e n t ( i n ) : : input sll_real64, dimension ( : ), i n t e n t ( out ) : : output ρ = input, Φ = output subroutine compute_e_from_rho (... ) class ( sll_poisson_1d_base ) : : poisson sll_real64, dimension ( : ), i n t e n t ( i n ) : : input sll_real64, dimension ( : ), i n t e n t ( out ) : : output ρ = input, E = output M. Mehrenberger (UDS) Interfaces and simulations in SELALIB Garching, November

19 In 2D 1. Examples of abstract interfaces POISSON SOLVERS subroutine compute_phi_from_rho (... ) class ( sll_poisson_2d_base ) : : poisson sll_real64, dimension ( :, : ), i n t e n t ( i n ) : : input sll_real64, dimension ( :, : ), i n t e n t ( out ) : : output ρ = input, Φ = output subroutine compute_e_from_rho (... ) class ( sll_poisson_2d_base ) : : poisson sll_real64, dimension ( :, : ), i n t e n t ( i n ) : : input sll_real64, dimension ( :, : ), i n t e n t ( out ) : : output1 sll_real64, dimension ( :, : ), i n t e n t ( out ) : : output2 ρ = input, E 1 = output1, E 2 = output2 M. Mehrenberger (UDS) Interfaces and simulations in SELALIB Garching, November

20 Simulation 1. Examples of abstract interfaces SIMULATION subroutine run ( sim ) class ( sll_simulation_base ), i n t e n t ( i n o u t ) : : sim & just run a simulation M. Mehrenberger (UDS) Interfaces and simulations in SELALIB Garching, November

21 2. Other interfaces Other interfaces Not all the interfaces need fortran 2003 abstract class Do not use fortran 2003, if not necessary try to choose the simplest solution (not always easy) Some examples developed Initial function Remap Examples not developed Point to point communication M. Mehrenberger (UDS) Interfaces and simulations in SELALIB Garching, November

22 Initial function 2. Other interfaces INITIAL FUNCTION f u n c t i o n s l l _ s c a l a r _ i n i t i a l i z e r _ 2 d ( & x1, & x2, & params ) & r e s u l t ( val ) sll_real64 : : val sll_real64, i n t e n t ( i n ) : : x1 sll_real64, i n t e n t ( i n ) : : x2 sll_real64, dimension ( : ), & i n t e n t ( i n ), o p t i o n a l : : params f (x 1, x 2 ) = val, x 1 = x1, x 2 = x2 params : storage of parameters describing f M. Mehrenberger (UDS) Interfaces and simulations in SELALIB Garching, November

23 2. Other interfaces REMAP Remap Handles parallelism with collective operations. permits to go from parallel in some direction(s) to parallel in other direction(s). Example of 4D array array of size N 1 N 2 N 3 N 4 dispatched in 32 processors f_parx1x2x4 is cut in directions x1, x2, x4 f_parx3 is cut in directions x3 remap changes from first to second configuration M. Mehrenberger (UDS) Interfaces and simulations in SELALIB Garching, November

24 Remap 2. Other interfaces REMAP par4d_x1x2x4 => new_layout_4d ( sll_world_collective ) c a l l i n i t i a l i z e _ l a y o u t _ w i t h _ d i s t r i b u t e d _ 4 D _ a r r a y ( & 129, & 257, & 33, & 129, & 4, & 4, & 1, & 2, & par4d_x1x2x4 ) N 1 = 129, N 2 = 257, N 3 = 33, N 4 = 129. for f_parx1x2x4 of size M. Mehrenberger (UDS) Interfaces and simulations in SELALIB Garching, November

25 Remap 2. Other interfaces REMAP par4d_x3 => new_layout_4d ( sll_world_collective ) c a l l i n i t i a l i z e _ l a y o u t _ w i t h _ d i s t r i b u t e d _ 4 D _ a r r a y ( & 129, & 257, & 33, & 129, & 1, & 1, & 32, & 1, & par4d_x3 ) N 1 = 129, N 2 = 257, N 3 = 33, N 4 = 129. for f_parx3 of size M. Mehrenberger (UDS) Interfaces and simulations in SELALIB Garching, November

26 Remap 2. Other interfaces REMAP par4d_x1x2x4_to_x3 => new_remap_plan ( & par4d_x1x2x4, & par4d_x3, & f_parx1x2x4 ) c a l l apply_remap_4d ( & par4d_x1x2x4_to_x3, & f_parx1x2x4, & f_parx3 ) M. Mehrenberger (UDS) Interfaces and simulations in SELALIB Garching, November

27 3. A specific example of abstract class Specific example of abstract class An abstract class permit to give some common rules An abstract class is initialized to a specific implementation through new initialize Initializationnew are not common to the interface M. Mehrenberger (UDS) Interfaces and simulations in SELALIB Garching, November

28 3. A specific example of abstract class INTERPOLATOR An example for cubic splines interp => new_cubic_spline_interpolator_1d ( & 129, & 0. _f64, & 1. _f64, & SLL_PERIODIC ) val = interp%interpolate_value ( eta ) M. Mehrenberger (UDS) Interfaces and simulations in SELALIB Garching, November

29 3. A specific example of abstract class INTERPOLATOR Initialization of cubic splines interpolator with new f u n c t i o n new_cubic_spline_interpolator_1d (... )& r e s u l t ( interp ) type ( sll_cubic_spline_interpolator_1d ), & p o i n t e r : : interp sll_int32, i n t e n t ( i n ) : : num_points sll_real64, i n t e n t ( i n ) : : xmin sll_real64, i n t e n t ( i n ) : : xmax sll_int32, i n t e n t ( i n ) : : bc_type sll_real64, i n t e n t ( i n ), o p t i o n a l : : s l o p e _ l e f t sll_real64, i n t e n t ( i n ), o p t i o n a l : : s l o p e _ r i g h t M. Mehrenberger (UDS) Interfaces and simulations in SELALIB Garching, November

30 3. A specific example of abstract class INTERPOLATOR Initialization of cubic splines interpolator with initialize subroutine initialize_cubic_spline_interpolator_1d (.. ) class ( sll_cubic_spline_interpolator_1d ), & i n t e n t ( i n o u t ) : : interp sll_int32, i n t e n t ( i n ) : : num_points sll_real64, i n t e n t ( i n ) : : xmin sll_real64, i n t e n t ( i n ) : : xmax sll_int32, i n t e n t ( i n ) : : bc_type sll_real64, i n t e n t ( i n ), o p t i o n a l : : s l o p e _ l e f t sll_real64, i n t e n t ( i n ), o p t i o n a l : : s l o p e _ r i g h t M. Mehrenberger (UDS) Interfaces and simulations in SELALIB Garching, November

31 3. A specific example of abstract class INTERPOLATOR Comments new has always the same structure (no need to think!) allocate the pointer call to initialize new and initialize have a specific name and with this name, we can choose the specific method M. Mehrenberger (UDS) Interfaces and simulations in SELALIB Garching, November

32 4. About simulations Simulation Simulations can be done, using Selalib as external library Simulations can be done, as programs inside Selalib (but not a module) Simulations can be done as a module (top level) Each simulation has his own flexibility Different simulations can be developed independently M. Mehrenberger (UDS) Interfaces and simulations in SELALIB Garching, November

33 4. About simulations Initialization of a simulation Incremental degree of flexibility Initialization with default parameters Initialization with filename : Use of filenames.nml Use of keywords "SLL_CUBIC_SPLINES" translated into select case at initialization Initialization by hand (most flexible feature) Abstract class is a parameter of the simulation Specific abstract class implementation coded outside the module. M. Mehrenberger (UDS) Interfaces and simulations in SELALIB Garching, November

34 4. About simulations Test of a simulation ctest selalib-result.sh run small example with optimized parameters give the results, with a simple diagnostic M. Mehrenberger (UDS) Interfaces and simulations in SELALIB Garching, November

35 Swirling deformation flow 4. About simulations EXAMPLES &geometry mesh_case_x1= "SLL_LOGICAL_MESH" num_ cells_x1 = 64 x1_min = x1_max = mesh_case_x2= "SLL_LOGICAL_MESH" num_ cells_x2 = 64 x2_min = x2_max = M. Mehrenberger (UDS) Interfaces and simulations in SELALIB Garching, November

36 Swirling deformation flow 4. About simulations EXAMPLES & i n i t i a l _ f u n c t i o n i n i t i a l _ f u n c t i o n _ c a s e = "SLL_COS_BELL" xc_1 = 1. xc_2 = 0.2 &t i m e _ i t e r a t i o n s dt = number_iterations = 10 freq_diag = 1 freq_diag_time = 1 time_loop_case = "SLL_PREDICTOR_CORRECTOR" M. Mehrenberger (UDS) Interfaces and simulations in SELALIB Garching, November

37 Swirling deformation flow 4. About simulations EXAMPLES &advector advect2d_case = " SLL_BSL" f_interp2d_case = "SLL_CUBIC_SPLINES" phi_interp2d_case = "SLL_CUBIC_SPLINES" charac2d_case = " SLL_VERLET" A_interp_case = "SLL_CUBIC_SPLINES" a d v e c t i o n _ f i e l d _ c a s e & = "SLL_SWIRLING_DEFORMATION_FLOW" time_period = 1.5 M. Mehrenberger (UDS) Interfaces and simulations in SELALIB Garching, November

38 Swirling deformation flow (0.1s) M. Mehrenberger (UDS) Interfaces and simulations in SELALIB Garching, November

39 guiding center simulation cartesian &geometry mesh_case_x1= "SLL_LANDAU_MESH" num_ cells_x1 = 32 x1_min = 0.0 nbox_x1 = 1 mesh_case_x2= "SLL_LANDAU_MESH" num_ cells_x2 = 32 x2_min = 0.0 nbox_x2 = 1 M. Mehrenberger (UDS) Interfaces and simulations in SELALIB Garching, November

40 guiding center simulation cartesian & i n i t i a l _ f u n c t i o n i n i t i a l _ f u n c t i o n _ c a s e = "SLL_KHP1" kmode_x1 = 0.5 kmode_x2 = 1. eps = 1e 6! &t i m e _ i t e r a t i o n s dt = 0.5 number_iterations = 120 freq_diag = 100 freq_diag_time = 1 time_ loop_ case = " SLL_LEAP_FROG" M. Mehrenberger (UDS) Interfaces and simulations in SELALIB Garching, November

41 guiding center simulation cartesian &advector advect2d_case = " SLL_BSL" charac2d_case = " SLL_VERLET" f_interp2d_case = "SLL_CUBIC_SPLINES" phi_interp2d_case = "SLL_CUBIC_SPLINES" A_interp_case = "SLL_CUBIC_SPLINES" &poisson poisson_case = "SLL_PHI_FROM_RHO" poisson_solver = "SLL_POISSON_FFT" M. Mehrenberger (UDS) Interfaces and simulations in SELALIB Garching, November

42 Guiding center simulation cartesian (0.23s) M. Mehrenberger (UDS) Interfaces and simulations in SELALIB Garching, November

43 guiding center simulation polar &geometry mesh_case= "SLL_POLAR_MESH" num_ cells_x1 = 64 r_min = 1. r_max = 10. num_ cells_x2 = 32 & i n i t i a l _ f u n c t i o n i n i t i a l _ f u n c t i o n _ c a s e = "SLL_DIOCOTRON" r_minus = 4. r_plus = 5. kmode_x2 = 3. eps = 1. e 6 M. Mehrenberger (UDS) Interfaces and simulations in SELALIB Garching, November

44 guiding center simulation polar &t i m e _ i t e r a t i o n s dt = 0.5 number_iterations = 160 freq_diag = 80 freq_diag_time = 1 time_loop_case = "SLL_PREDICTOR_CORRECTOR" &poisson poisson_case = "SLL_PHI_FROM_RHO" poisson_solver = "SLL_POLAR_FFT" M. Mehrenberger (UDS) Interfaces and simulations in SELALIB Garching, November

45 guiding center simulation polar &advector advect2d_case = " SLL_BSL" f_interp2d_case = "SLL_HERMITE" hermite_degree_eta1 = 6 hermite_degree_eta2 = 6 phi_interp2d_case = "SLL_CUBIC_SPLINES" charac2d_case = " SLL_VERLET" A_interp_case = "SLL_CUBIC_SPLINES" M. Mehrenberger (UDS) Interfaces and simulations in SELALIB Garching, November

46 Guiding center simulation polar (0.9s) M. Mehrenberger (UDS) Interfaces and simulations in SELALIB Garching, November

47 vlasov poisson cartesian (landau) &geometry mesh_case_x1 = "SLL_LANDAU_MESH" num_ cells_x1 = 32 x1_min = 0.0 nbox_x1 = 1 mesh_case_x2 = "SLL_LOGICAL_MESH" num_ cells_x2 = 64 x2_min = 6.0 x2_max = 6.0 M. Mehrenberger (UDS) Interfaces and simulations in SELALIB Garching, November

48 vlasov poisson cartesian (landau) & i n i t i a l _ f u n c t i o n i n i t i a l _ f u n c t i o n _ c a s e kmode = 0.5 eps = = "SLL_LANDAU" &t i m e _ i t e r a t i o n s dt = 0.1 number_iterations = 600 freq_diag = 100 freq_diag_time = 1 nb_mode = 20 t i m e _ i n i t = 0. s p l i t _ c a s e = "SLL_STRANG_VTV" M. Mehrenberger (UDS) Interfaces and simulations in SELALIB Garching, November

49 vlasov poisson cartesian (landau) &advector advector_x1 = " SLL_SPLINES " order_x1 = 4 advector_x2 = " SLL_SPLINES " order_x2 = 4 &poisson poisson_solver = " SLL_FFT " &d r i v e d r i v e _ t y p e = "SLL_NO_DRIVE" M. Mehrenberger (UDS) Interfaces and simulations in SELALIB Garching, November

50 vlasov poisson cartesian (landau) (0.6s) M. Mehrenberger (UDS) Interfaces and simulations in SELALIB Garching, November

51 vlasov poisson cartesian (bump on tail) &geometry mesh_case_x1 = "SLL_LANDAU_MESH" num_ cells_x1 = 256 x1_min = 0.0 nbox_x1 = 3 mesh_case_x2 = "SLL_LOGICAL_MESH" num_ cells_x2 = 256 x2_min = 8.0 x2_max = 8.0 & i n i t i a l _ f u n c t i o n i n i t i a l _ f u n c t i o n _ c a s e = "SLL_BUMP_ON_TAIL" kmode = 0.3 eps = M. Mehrenberger (UDS) Interfaces and simulations in SELALIB Garching, November

52 vlasov poisson cartesian (bump on tail) &t i m e _ i t e r a t i o n s dt = 2.5 number_iterations = 40 freq_diag = 10 freq_diag_time = 1 s p l i t _ c a s e = "SLL_ORDER6VPnew1_VTV" &advector advector_x1 = "SLL_LAGRANGE" order_x1 = 10 advector_x2 = "SLL_LAGRANGE" order_x2 = 10 &poisson poisson_solver = " SLL_FFT " M. Mehrenberger (UDS) Interfaces and simulations in SELALIB Garching, November

53 vlasov poisson cartesian (bump on tail) &advector advector_x1 = "SLL_LAGRANGE" order_x1 = 10 advector_x2 = "SLL_LAGRANGE" order_x2 = 10 &poisson poisson_solver = " SLL_FFT " &d r i v e d r i v e _ t y p e = "SLL_NO_DRIVE" M. Mehrenberger (UDS) Interfaces and simulations in SELALIB Garching, November

54 vlasov poisson cartesian (bump on tail) (2s) M. Mehrenberger (UDS) Interfaces and simulations in SELALIB Garching, November

55 vlasov poisson cartesian (two stream instability) &geometry mesh_case_x1 = "SLL_LANDAU_MESH" num_ cells_x1 = 256 x1_min = 0.0 nbox_x1 = 1 mesh_case_x2 = "SLL_LOGICAL_MESH" num_ cells_x2 = 1024 x2_min = 6.0 x2_max = 6.0 & i n i t i a l _ f u n c t i o n i n i t i a l _ f u n c t i o n _ c a s e = "SLL_TWO_STREAM_INSTABILITY" kmode = 0.5 eps = M. Mehrenberger (UDS) Interfaces and simulations in SELALIB Garching, November

56 vlasov poisson cartesian (two stream instability) &t i m e _ i t e r a t i o n s dt = 1.75! 0.5 number_iterations = 20 freq_diag = 5 freq_diag_time = 1 s p l i t _ c a s e = "SLL_ORDER6VPnew1_VTV" M. Mehrenberger (UDS) Interfaces and simulations in SELALIB Garching, November

57 vlasov poisson cartesian (two stream instability) &poisson poisson_solver = " SLL_FFT " &d r i v e d r i v e _ t y p e = "SLL_NO_DRIVE" M. Mehrenberger (UDS) Interfaces and simulations in SELALIB Garching, November

58 vlasov poisson cartesian (two stream instability) (3s) M. Mehrenberger (UDS) Interfaces and simulations in SELALIB Garching, November

59 vlasov poisson cartesian (KEEN wave) &geometry mesh_case_x1 = "SLL_LANDAU_MESH" num_ cells_x1 = 128 x1_min = 0.0 nbox_x1 = 1 mesh_case_x2 = "SLL_LOGICAL_MESH" num_ cells_x2 = 174 x2_min = 6.0 x2_max = 6.0 mesh_case_x2 = "SLL_TWO_GRID_MESH" x2_fine_min = 0.375! 0.36 x2_fine_max = 2. 25! density_x2_min_to_x2_fine_min = 1 density_x2_fine_min_to_x2_fine_max = 12 density_ x2_ fine_ max_ to_ x2_ max = 1 M. Mehrenberger (UDS) Interfaces and simulations in SELALIB Garching, November

60 vlasov poisson cartesian (KEEN wave) & i n i t i a l _ f u n c t i o n i n i t i a l _ f u n c t i o n _ c a s e = "SLL_LANDAU" kmode = eps = 0.0 &t i m e _ i t e r a t i o n s dt = 2. number_iterations = 100 freq_diag = 20 f r e q _ d i a g _ r e s t a r t = 100 freq_diag_time = 1 nb_mode = 20 t i m e _ i n i t = 0. s p l i t _ c a s e = "SLL_ORDER6VPnew1_VTV" M. Mehrenberger (UDS) Interfaces and simulations in SELALIB Garching, November

61 vlasov poisson cartesian (KEEN wave) &advector advector_x1 = "SLL_LAGRANGE" order_x1 = 18 advector_x2 = "SLL_NON_UNIFORM_CUBIC_SPLINES" order_x2 = 18 advection_ form_ x2 = " SLL_CONSERVATIVE" i n t e g r a t i o n _ c a s e = "SLL_CONSERVATIVE" &poisson poisson_solver = " SLL_FFT " M. Mehrenberger (UDS) Interfaces and simulations in SELALIB Garching, November

62 vlasov poisson cartesian (KEEN wave) &d r i v e d r i v e _ t y p e = "SLL_KEEN_DRIVE" keen_t0 = 0. keen_tl = 69. keen_tr = 307. keen_twl = 20. keen_twr = 20. k e e n _ t f l a t = 100. k e e n _ t u r n _ d r i v e _ o f f =. t r u e. keen_edrmax = 0.2 keen_omegadr = M. Mehrenberger (UDS) Interfaces and simulations in SELALIB Garching, November

63 vlasov poisson cartesian (KEEN wave) (1.5s) M. Mehrenberger (UDS) Interfaces and simulations in SELALIB Garching, November

64 vlasov poisson cartesian (beam) &geometry mesh_case_x1 = "SLL_LOGICAL_MESH" num_ cells_x1 = 128 x1_min = 4.0 x1_max = 4.0 mesh_case_x2 = "SLL_LOGICAL_MESH" num_ cells_x2 = 128 x2_min = 4.0 x2_max = 4.0 M. Mehrenberger (UDS) Interfaces and simulations in SELALIB Garching, November

65 vlasov poisson cartesian (beam) & i n i t i a l _ f u n c t i o n i n i t i a l _ f u n c t i o n _ c a s e = "SLL_LANDAU" kmode = 0.5 eps = i n i t i a l _ f u n c t i o n _ c a s e = "SLL_BEAM" alpha_gaussian = 0. 2 &t i m e _ i t e r a t i o n s dt = 0.5 number_iterations = 41 freq_diag = 10 freq_diag_time = 1 s p l i t _ c a s e = "SLL_ORDER6VPnew1_VTV" M. Mehrenberger (UDS) Interfaces and simulations in SELALIB Garching, November

66 vlasov poisson cartesian (beam) &advector advector_x1 = " SLL_SPLINES " order_x1 = 4 advector_x2 = " SLL_SPLINES " order_x2 = 4 f a c t o r _ x 1 = factor_x2_1 = factor_x2_rho = 1. &poisson poisson_solver = "SLL_POLAR" &d r i v e d r i v e _ t y p e = "SLL_NO_DRIVE" M. Mehrenberger (UDS) Interfaces and simulations in SELALIB Garching, November

67 vlasov poisson cartesian (beam) (0.54s) M. Mehrenberger (UDS) Interfaces and simulations in SELALIB Garching, November

68 vlasov poisson cartesian no splitting (beam) (0.4s) &geometry mesh_case_x1 = "SLL_LOGICAL_MESH" num_ cells_x1 = 64 x1_min = 4.0 x1_max = 4.0 mesh_case_x2 = "SLL_LOGICAL_MESH" num_ cells_x2 = 64 x2_min = 4.0 x2_max = 4.0 M. Mehrenberger (UDS) Interfaces and simulations in SELALIB Garching, November

69 vlasov poisson cartesian no splitting (beam) (0.4s) & i n i t i a l _ f u n c t i o n! i n i t i a l _ f u n c t i o n _ c a s e = "SLL_LANDAU" kmode = 0.5 eps = i n i t i a l _ f u n c t i o n _ c a s e = "SLL_BEAM" alpha_gaussian = 0. 2 &t i m e _ i t e r a t i o n s dt = 0.25 number_iterations = 80 freq_diag = 20 freq_diag_time = 1 time_loop_case = "SLL_PREDICTOR_CORRECTOR" M. Mehrenberger (UDS) Interfaces and simulations in SELALIB Garching, November

70 vlasov poisson cartesian no splitting (beam) (0.4s) &advector advect2d_case = " SLL_BSL" f_interp2d_case = "SLL_CUBIC_SPLINES" phi_interp2d_case = "SLL_CUBIC_SPLINES"! charac2d_case = "SLL_EULER" charac2d_case = " SLL_VERLET" A_interp_case = "SLL_CUBIC_SPLINES" f a c t o r _ x 1 = ! factor_x2_1 = ! factor_x2_rho = 1. &poisson poisson_solver = "SLL_POLAR" &d r i v e M. Mehrenberger (UDS) Interfaces and simulations in SELALIB Garching, November

71 vlasov poisson cartesian no splitting (beam) (0.4s) M. Mehrenberger (UDS) Interfaces and simulations in SELALIB Garching, November

72 vlasov poisson cartesian 4d (landau) &geometry mesh_case_x1= "SLL_LANDAU_MESH" num_ cells_x1 = 16 x1_min = 0.0 nbox_x1 = 1 mesh_case_x2= "SLL_LANDAU_MESH" num_ cells_x2 = 16 x2_min = 0.0 nbox_x2 = 1 mesh_case_x3= "SLL_LOGICAL_MESH" num_ cells_x3 = 32 x3_min = 6. x3_max = 6. mesh_case_x4= "SLL_LOGICAL_MESH" num_ cells_x4 = 32 x4_min = 6. x4_max = 6. M. Mehrenberger (UDS) Interfaces and simulations in SELALIB Garching, November

73 vlasov poisson cartesian 4d (landau) & i n i t i a l _ f u n c t i o n i n i t i a l _ f u n c t i o n _ c a s e = "SLL_LANDAU" kmode_x1 = 0.5 kmode_x2 = 0.5 eps = 1e 3 &t i m e _ i t e r a t i o n s dt = 2 number_iterations = 5 freq_diag = 20 freq_diag_time = 1 s p l i t _ c a s e = "SLL_ORDER6VPnew1_VTV" M. Mehrenberger (UDS) Interfaces and simulations in SELALIB Garching, November

74 vlasov poisson cartesian 4d (landau) &advector advector_x1 = "SLL_LAGRANGE" order_x1 = 4 advector_x2 = "SLL_LAGRANGE" order_x2 = 4 advector_x3 = "SLL_LAGRANGE" order_x3 = 4 advector_x4 = "SLL_LAGRANGE" order_x4 = 4 M. Mehrenberger (UDS) Interfaces and simulations in SELALIB Garching, November

75 vlasov poisson cartesian 4d (landau) (2.3s) M. Mehrenberger (UDS) Interfaces and simulations in SELALIB Garching, November

76 drift kinetic 4d 4. About simulations EXAMPLES &mesh num_ cells_x1 = 32 num_ cells_x2 = 32 num_ cells_x3 = 8 num_ cells_x4 = 32 r_min = 0.1 r_max = z_min = 0. z_max = v_min = 7.32 v_max = M. Mehrenberger (UDS) Interfaces and simulations in SELALIB Garching, November

77 drift kinetic 4d 4. About simulations EXAMPLES &e q u i l i b r i u m tau0 = 1. rho_peak = 0.5 kappan = 13.2! d e l t a r n = 0.2! 2.9 kappati = 66.! d e l t a r T i = 0.1! 1.45 kappate = 66.! d e l t a r T e = 0.1! 1.45 QN_case = "SLL_QUASI_NEUTRAL_WITHOUT_ZONAL_FLOW" poisson2d_bc_rmin = "SLL_NEUMANN_MODE_0" poisson2d_bc_rmax = " SLL_DIRICHLET " M. Mehrenberger (UDS) Interfaces and simulations in SELALIB Garching, November

78 drift kinetic 4d 4. About simulations EXAMPLES &p e r t u r b a t i o n perturb_choice = 1 mmode = 5 nmode = 1 eps_perturb = 1.0e 6 M. Mehrenberger (UDS) Interfaces and simulations in SELALIB Garching, November

79 drift kinetic 4d 4. About simulations EXAMPLES &sim_params dt = 16 number_iterations = 100 freq_diag_time = 1 freq_diag = 100 time_loop_case = "SLL_TIME_LOOP_PREDICTOR_CORRECTOR" charac2d_case = "SLL_CHARAC_EULER" charac2d_tol = 1e 12 charac2d_maxiter = 1000 advect2d_case = " SLL_BSL" interp_x1x2 = "SLL_CUBIC_SPLINES" phi_interp_x1x2 = "SLL_CUBIC_SPLINES" p h i _ i n t e r p _ x 3 = "SLL_CUBIC_SPLINES" advector_x3 = " SLL_SPLINES " order_x3 = 4 advector_x4 = " SLL_SPLINES " order_x4 = 4 poisson2d_case = " POLAR_FFT" M. Mehrenberger (UDS) Interfaces and simulations in SELALIB Garching, November

80 drift kinetic 4d (21s) 4. About simulations EXAMPLES M. Mehrenberger (UDS) Interfaces and simulations in SELALIB Garching, November

81 drift kinetic 4d one mu (22s) &mesh num_ cells_x1 = 32 num_ cells_x2 = 32 num_ cells_x3 = 8 num_ cells_x4 = 32 r_min = 0.1 r_max = z_min = 0. z_max = v_min = 7.32 v_max = M. Mehrenberger (UDS) Interfaces and simulations in SELALIB Garching, November

82 drift kinetic 4d one mu (22s) &e q u i l i b r i u m tau0 = 1. rho_peak = 0.5 kappan = 13.2! d e l t a r n = 0.2! 2.9 kappati = 66.! d e l t a r T i = 0.1! 1.45 kappate = 66.! d e l t a r T e = 0.1! 1.45 QN_case = "SLL_QUASI_NEUTRAL_WITHOUT_ZONAL_FLOW" poisson2d_bc_rmin = "SLL_NEUMANN_MODE_0" poisson2d_bc_rmax = " SLL_DIRICHLET " M. Mehrenberger (UDS) Interfaces and simulations in SELALIB Garching, November

83 drift kinetic 4d one mu (22s) &p e r t u r b a t i o n perturb_choice = 1 mmode = 5 nmode = 1 eps_perturb = 1.0e 6 M. Mehrenberger (UDS) Interfaces and simulations in SELALIB Garching, November

84 drift kinetic 4d one mu (22s) &sim_params dt = 32 number_iterations = 50! number_iterations = 1 freq_diag_time = 1 freq_diag = 100! time_loop_case = "SLL_TIME_LOOP_EULER" time_loop_case = "SLL_TIME_LOOP_PREDICTOR_CORRECTOR" charac2d_case = "SLL_CHARAC_VERLET" charac2d_tol = 1e 12 charac2d_maxiter = 1000 advect2d_case = " SLL_BSL" interp_x1x2 = "SLL_CUBIC_SPLINES" interp_x1x2 = "SLL_HERMITE" hermite_degree_eta1 = 6 hermite_degree_eta2 = 6 M. Mehrenberger (UDS) Interfaces and simulations in SELALIB Garching, November

85 drift kinetic 4d one mu (22s) phi_interp_x1x2 = "SLL_CUBIC_SPLINES" p h i _ i n t e r p _ x 3 = "SLL_CUBIC_SPLINES" advector_x3 = " SLL_SPLINES " order_x3 = 4 advector_x4 = " SLL_SPLINES " order_x4 = 4 poisson2d_case = " POLAR_FFT" gyroaverage_case = "HERMITE_C1_PRECOMPUTE" mu = 2.2e 0 gyroaverage_n_points = 1024 gyroaverage_interp_degree_x1 = 3 gyroaverage_interp_degree_x2 = 3 delta_f_method = 1 M. Mehrenberger (UDS) Interfaces and simulations in SELALIB Garching, November

86 drift kinetic 4d one mu (22s) M. Mehrenberger (UDS) Interfaces and simulations in SELALIB Garching, November

87 drift kinetic 4d field aligned &mesh num_ cells_x1 = 32 num_ cells_x2 = 32 num_ cells_x3 = 8 num_ cells_x4 = 32 r_min = 0.1 r_max = z_min = 0. z_max = v_min = 7.32 v_max = M. Mehrenberger (UDS) Interfaces and simulations in SELALIB Garching, November

88 drift kinetic 4d field aligned &e q u i l i b r i u m tau0 = 1. rho_peak = 0.5 kappan = 13.2! d e l t a r n = 0.2! 2.9 kappati = 66.! d e l t a r T i = 0.1! 1.45 kappate = 66.! d e l t a r T e = 0.1! 1.45 QN_case = "SLL_QUASI_NEUTRAL_WITHOUT_ZONAL_FLOW" poisson2d_bc_rmin = "SLL_NEUMANN_MODE_0" poisson2d_bc_rmax = " SLL_DIRICHLET " i o t a 0 = 0.8 Dr_iota0 = 0 B0 = 1. M. Mehrenberger (UDS) Interfaces and simulations in SELALIB Garching, November

89 drift kinetic 4d field aligned &p e r t u r b a t i o n perturb_choice = 1 mmode = 5 nmode = 3 eps_perturb = 1.0e 6 M. Mehrenberger (UDS) Interfaces and simulations in SELALIB Garching, November

90 drift kinetic 4d field aligned &sim_params dt = 16 number_iterations = 100! 8000! number_iterations = 1 freq_diag_time = 1 freq_diag = 50! time_loop_case = "SLL_TIME_LOOP_EULER" time_loop_case = "SLL_TIME_LOOP_PREDICTOR_CORRECTOR" charac2d_case = "SLL_CHARAC_VERLET" charac2d_tol = 1e 12 charac2d_maxiter = 2! 1000 advect2d_case = " SLL_BSL" interp_x1x2 = "SLL_CUBIC_SPLINES" phi_interp_x1x2 = "SLL_CUBIC_SPLINES" p h i _ i n t e r p _ x 3 = "SLL_CUBIC_SPLINES" M. Mehrenberger (UDS) Interfaces and simulations in SELALIB Garching, November

91 drift kinetic 4d field aligned advector_x2 = " SLL_SPLINES " order_x2 = 4 advector_x3 = " SLL_SPLINES " order_x3 = 4 advector_x4 = " SLL_SPLINES " order_x4 = 4 poisson2d_case = " POLAR_FFT" l a g r a n g e _ s t e n c i l _ l e f t = 2 l a g r a n g e _ s t e n c i l _ r i g h t = 3 d e r i v _ s t e n c i l _ l e f t = 3 d e r i v _ s t e n c i l _ r i g h t = 3 u s e _ f i e l d _ a l i g n e d _ i n t e r p o l a t i o n = T u s e _ f i e l d _ a l i g n e d _ d e r i v a t i v e = T M. Mehrenberger (UDS) Interfaces and simulations in SELALIB Garching, November

92 drift kinetic 4d field aligned (60s) M. Mehrenberger (UDS) Interfaces and simulations in SELALIB Garching, November