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 PetscErrorCode ierr; 112328674fSJunchao Zhang 122328674fSJunchao Zhang ierr = PetscInitialize(&argc, &args, (char*) 0, help);if (ierr) return ierr; 132328674fSJunchao Zhang comm = PETSC_COMM_WORLD; 14*5f80ce2aSJacob Faibussowitsch CHKERRQ(PetscOptionsGetInt(NULL,NULL, "-N", &N, NULL)); 15*5f80ce2aSJacob Faibussowitsch CHKERRQ(MatCreate(comm, &A)); 16*5f80ce2aSJacob Faibussowitsch CHKERRQ(MatSetSizes(A, PETSC_DECIDE, PETSC_DECIDE, N, N)); 17*5f80ce2aSJacob Faibussowitsch CHKERRQ(MatSetFromOptions(A)); 18*5f80ce2aSJacob Faibussowitsch CHKERRQ(MatSeqAIJSetPreallocation(A, 3, NULL)); 19*5f80ce2aSJacob Faibussowitsch CHKERRQ(MatMPIAIJSetPreallocation(A, 3, NULL, 2, NULL)); 20*5f80ce2aSJacob Faibussowitsch CHKERRQ(MatZeroEntries(A)); 21*5f80ce2aSJacob Faibussowitsch CHKERRQ(MatDestroy(&A)); 222328674fSJunchao Zhang ierr = PetscFinalize(); 232328674fSJunchao Zhang return ierr; 242328674fSJunchao Zhang } 252328674fSJunchao Zhang 262328674fSJunchao Zhang /*TEST 272328674fSJunchao Zhang test: 282328674fSJunchao Zhang requires: kokkos_kernels 292328674fSJunchao Zhang nsize: {{1 2}} 302328674fSJunchao Zhang output_file: output/ex252_1.out 312328674fSJunchao Zhang args: -mat_type aijkokkos 322328674fSJunchao Zhang TEST*/ 33