12328674fSJunchao Zhang static char help[] = "Test MatZeroEntries() on unassembled matrices \n\n";
22328674fSJunchao Zhang
32328674fSJunchao Zhang #include <petscmat.h>
42328674fSJunchao Zhang
main(int argc,char ** args)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