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*b122ec5aSJacob Faibussowitsch CHKERRQ(PetscInitialize(&argc, &args, (char*) 0, help)); 122328674fSJunchao Zhang comm = PETSC_COMM_WORLD; 135f80ce2aSJacob Faibussowitsch CHKERRQ(PetscOptionsGetInt(NULL,NULL, "-N", &N, NULL)); 145f80ce2aSJacob Faibussowitsch CHKERRQ(MatCreate(comm, &A)); 155f80ce2aSJacob Faibussowitsch CHKERRQ(MatSetSizes(A, PETSC_DECIDE, PETSC_DECIDE, N, N)); 165f80ce2aSJacob Faibussowitsch CHKERRQ(MatSetFromOptions(A)); 175f80ce2aSJacob Faibussowitsch CHKERRQ(MatSeqAIJSetPreallocation(A, 3, NULL)); 185f80ce2aSJacob Faibussowitsch CHKERRQ(MatMPIAIJSetPreallocation(A, 3, NULL, 2, NULL)); 195f80ce2aSJacob Faibussowitsch CHKERRQ(MatZeroEntries(A)); 205f80ce2aSJacob Faibussowitsch CHKERRQ(MatDestroy(&A)); 21*b122ec5aSJacob Faibussowitsch CHKERRQ(PetscFinalize()); 22*b122ec5aSJacob 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