12328674fSJunchao Zhang static char help[] = "Test MatZeroEntries() on unassembled matrices \n\n"; 22328674fSJunchao Zhang 32328674fSJunchao Zhang #include <petscmat.h> 42328674fSJunchao Zhang 52328674fSJunchao Zhang int main(int argc, char **args) 62328674fSJunchao Zhang { 72328674fSJunchao Zhang Mat A; 82328674fSJunchao Zhang PetscInt N = 32; 92328674fSJunchao Zhang MPI_Comm comm; 102328674fSJunchao Zhang 11*9566063dSJacob Faibussowitsch PetscCall(PetscInitialize(&argc, &args, (char*) 0, help)); 122328674fSJunchao Zhang comm = PETSC_COMM_WORLD; 13*9566063dSJacob Faibussowitsch PetscCall(PetscOptionsGetInt(NULL,NULL, "-N", &N, NULL)); 14*9566063dSJacob Faibussowitsch PetscCall(MatCreate(comm, &A)); 15*9566063dSJacob Faibussowitsch PetscCall(MatSetSizes(A, PETSC_DECIDE, PETSC_DECIDE, N, N)); 16*9566063dSJacob Faibussowitsch PetscCall(MatSetFromOptions(A)); 17*9566063dSJacob Faibussowitsch PetscCall(MatSeqAIJSetPreallocation(A, 3, NULL)); 18*9566063dSJacob Faibussowitsch PetscCall(MatMPIAIJSetPreallocation(A, 3, NULL, 2, NULL)); 19*9566063dSJacob Faibussowitsch PetscCall(MatZeroEntries(A)); 20*9566063dSJacob Faibussowitsch PetscCall(MatDestroy(&A)); 21*9566063dSJacob Faibussowitsch PetscCall(PetscFinalize()); 22b122ec5aSJacob Faibussowitsch return 0; 232328674fSJunchao Zhang } 242328674fSJunchao Zhang 252328674fSJunchao Zhang /*TEST 262328674fSJunchao Zhang test: 272328674fSJunchao Zhang requires: kokkos_kernels 282328674fSJunchao Zhang nsize: {{1 2}} 292328674fSJunchao Zhang output_file: output/ex252_1.out 302328674fSJunchao Zhang args: -mat_type aijkokkos 312328674fSJunchao Zhang TEST*/ 32