seq/kokkos/aijkok.kokkos.cxx

1832 // R and C are represented by rowperm and colperm in factors.
1833 // If R or C is identity (i.e, no reordering), then rowperm or colperm is empty.
1843   PetscIntKokkosView         &colperm      = factors->colperm;
1845   PetscBool                   col_identity = colperm.extent(0) ? PETSC_FALSE : PETSC_TRUE;
1873   if (!col_identity) PetscCallCXX(Kokkos::parallel_for(Kokkos::RangePolicy<>(exec, 0, m), KOKKOS_LAMBDA(const PetscInt i) { x(colperm(i)) = X(i); }));
1882 // R and C are represented by rowperm and colperm in factors.
1883 // If R or C is identity (i.e, no reordering), then rowperm or colperm is empty.
1894   PetscIntKokkosView         &colperm      = factors->colperm;
1896   PetscBool                   col_identity = colperm.extent(0) ? PETSC_FALSE : PETSC_TRUE;
1910     PetscCallCXX(Kokkos::parallel_for(Kokkos::RangePolicy<>(exec, 0, m), KOKKOS_LAMBDA(const PetscInt i) { B(i) = b(colperm(i)); }));
1995       IS        colperm = ((Mat_SeqAIJ *)B->data)->col;
1997       PetscCall(ISIdentity(colperm, &col_identity));
2001         PetscCall(ISGetIndices(colperm, &ip));
2002         factors->colperm = PetscIntKokkosView(NoInit("colperm"), n);
2003         PetscCallCXX(Kokkos::deep_copy(factors->colperm, PetscIntKokkosViewHost(const_cast<PetscInt *>(ip), n)));
2004         PetscCall(ISRestoreIndices(colperm, &ip));