23 PetscErrorCode landau_field_print_access_callback(DM dm, Vec x, PetscInt local_field, PetscInt grid…  in landau_field_print_access_callback()  argument
31   species = ctx->species_offset[grid] + local_field;  in landau_field_print_access_callback()
32   val     = (PetscScalar)(LAND_PACK_IDX(b_id, grid) + (species + 1) * 10);  in landau_field_print_access_callback()
34 …_FMT " and local field %" PetscInt_FMT " with %" PetscInt_FMT " grids\n", grid, b_id, local_field,…  in landau_field_print_access_callback()
93 …DM dm, Vec X, PetscReal time, PetscReal temps[], PetscReal ns[], PetscInt grid, PetscReal shifts[]…  in SetMaxwellians()  argument
100 …for (PetscInt ii = ctx->species_offset[grid], i0 = 0; ii < ctx->species_offset[grid + 1]; ii++, i0…  in SetMaxwellians()
110 …ign(ctx->charges[ctx->species_offset[grid]]) * ctx->electronShift * ctx->m_0 / ctx->masses[ctx->sp…  in SetMaxwellians()
114 …data[0].shift = ctx->electronShift * shifts[grid] * PetscSqrtReal(data[0].kT_m) / ctx->v_0; // shi…  in SetMaxwellians()
256       for (PetscInt grid = 0; grid < ctx->num_grids; grid++) {  in Monitor()  local
258         DM          dm      = ctx->plex[grid];  in Monitor()
260 …      PetscReal   vz_0 = 0, n, energy, e_perp, e_par, m_s = ctx->masses[ctx->species_offset[grid]];  in Monitor()
261         Vec         Xloc = XsubArray[LAND_PACK_IDX(ctx->batch_view_idx, grid)];  in Monitor()
289 …if (grid == 0) PetscCall(PetscPrintf(PETSC_COMM_WORLD, "step %4d) time= %e temperature (eV): ", (i…  in Monitor()
290 …PetscCall(PetscPrintf(PETSC_COMM_WORLD, "%s T= %9.4e T_par= %9.4e T_perp= %9.4e ", (grid == 0) ? "…  in Monitor()
291         if (n_cm3[grid] == 0) n_cm3[grid] = ctx->n_0 * n * 1e-6; // does not change m^3 --> cm^3  in Monitor()
374     for (PetscInt grid = 0; grid < ctx->num_grids; grid++) {  in main()  local
376 …PetscCall(SetMaxwellians(ctx->plex[grid], XsubArray[LAND_PACK_IDX(b_id, grid)], 0.0, ctx->thermal_…  in main()