12328674fSJunchao Zhang static char help[] = "Test MatZeroEntries() on unassembled matrices \n\n"; 22328674fSJunchao Zhang 32328674fSJunchao Zhang #include <petscmat.h> 42328674fSJunchao Zhang 5d71ae5a4SJacob Faibussowitsch int main(int argc, char **args) 6d71ae5a4SJacob Faibussowitsch { 72328674fSJunchao Zhang Mat A; 82328674fSJunchao Zhang PetscInt N = 32; 92328674fSJunchao Zhang MPI_Comm comm; 102328674fSJunchao Zhang 11327415f7SBarry Smith PetscFunctionBeginUser; 12c8025a54SPierre Jolivet PetscCall(PetscInitialize(&argc, &args, NULL, help)); 132328674fSJunchao Zhang comm = PETSC_COMM_WORLD; 149566063dSJacob Faibussowitsch PetscCall(PetscOptionsGetInt(NULL, NULL, "-N", &N, NULL)); 159566063dSJacob Faibussowitsch PetscCall(MatCreate(comm, &A)); 169566063dSJacob Faibussowitsch PetscCall(MatSetSizes(A, PETSC_DECIDE, PETSC_DECIDE, N, N)); 179566063dSJacob Faibussowitsch PetscCall(MatSetFromOptions(A)); 189566063dSJacob Faibussowitsch PetscCall(MatSeqAIJSetPreallocation(A, 3, NULL)); 199566063dSJacob Faibussowitsch PetscCall(MatMPIAIJSetPreallocation(A, 3, NULL, 2, NULL)); 209566063dSJacob Faibussowitsch PetscCall(MatZeroEntries(A)); 219566063dSJacob Faibussowitsch PetscCall(MatDestroy(&A)); 229566063dSJacob Faibussowitsch PetscCall(PetscFinalize()); 23b122ec5aSJacob Faibussowitsch return 0; 242328674fSJunchao Zhang } 252328674fSJunchao Zhang 262328674fSJunchao Zhang /*TEST 272328674fSJunchao Zhang test: 282328674fSJunchao Zhang requires: kokkos_kernels 292328674fSJunchao Zhang nsize: {{1 2}} 30*3886731fSPierre Jolivet output_file: output/empty.out 312328674fSJunchao Zhang args: -mat_type aijkokkos 322328674fSJunchao Zhang TEST*/ 33