dmplexlandau/kokkos/landau.kokkos.cxx

44   ScalarType gg2[LANDAU_DIM];
45   ScalarType gg3[LANDAU_DIM][LANDAU_DIM];
50     for (int j = 0; j < LANDAU_DIM; j++) {
52       for (int k = 0; k < LANDAU_DIM; k++) gg3[j][k] = 0;
58     for (int j = 0; j < LANDAU_DIM; j++) {
60       for (int k = 0; k < LANDAU_DIM; k++) gg3[j][k] = rhs.gg3[j][k];
66     for (int j = 0; j < LANDAU_DIM; j++) {
68       for (int k = 0; k < LANDAU_DIM; k++) gg3[j][k] += src.gg3[j][k];
75     for (int j = 0; j < LANDAU_DIM; j++) {
77       for (int k = 0; k < LANDAU_DIM; k++) gg3[j][k] += src.gg3[j][k];
148   PetscCheck(LANDAU_DIM == dim, PETSC_COMM_WORLD, PETSC_ERR_PLIB, "dim %" PetscInt_FMT " != LANDAU_DIM %d", dim, LANDAU_DIM);
415   const PetscReal *d_z          = (LANDAU_DIM == 3) ? d_z_k->data() : NULL;
454   PetscCheck(LANDAU_DIM == dim, PETSC_COMM_WORLD, PETSC_ERR_PLIB, "dim %" PetscInt_FMT " != LANDAU_DIM %d", dim, LANDAU_DIM);
556                 PetscReal      refSpaceDer[LANDAU_DIM];
559                 for (d = 0; d < LANDAU_DIM; ++d) refSpaceDer[d] = 0.0;
615   #if LANDAU_DIM == 2
633   #if LANDAU_DIM == 3
640   #if LANDAU_DIM == 3
644   #if LANDAU_DIM == 2