58   Mat          *jacob_mat, jacob_mat_coarse, *prolong_restr_mat;  in main()  local
384   PetscCall(PetscMalloc1(num_levels, &jacob_mat));  in main()
392 …e[level], U_l_size[level], U_g_size[level], U_g_size[level], jacob_ctx[level], &jacob_mat[level]));  in main()
393 … PetscCall(MatShellSetOperation(jacob_mat[level], MATOP_MULT, (void (*)(void))ApplyJacobian_Ceed));  in main()
394 …PetscCall(MatShellSetOperation(jacob_mat[level], MATOP_GET_DIAGONAL, (void (*)(void))GetDiag_Ceed)…  in main()
395     PetscCall(MatShellSetVecType(jacob_mat[level], vectype));  in main()
401   form_jacob_ctx->jacob_mat  = jacob_mat;  in main()
465 …PetscCall(SNESSetJacobian(snes, jacob_mat[fine_level], jacob_mat[fine_level], FormJacobian, form_j…  in main()
467 …PetscCall(SNESSetJacobian(snes, jacob_mat[0], jacob_mat_coarse, SNESComputeJacobianDefaultColor, N…  in main()
518         PetscCall(KSPSetOperators(ksp_smoother, jacob_mat[level], jacob_mat[level]));  in main()
810     PetscCall(MatDestroy(&jacob_mat[level]));  in main()
835   PetscCall(PetscFree(jacob_mat));  in main()