380   A->ops->mult                      = MatMult_ConstantDiagonal;  in MatCreate_ConstantDiagonal()
381   A->ops->multadd                   = MatMultAdd_ConstantDiagonal;  in MatCreate_ConstantDiagonal()
382   A->ops->multtranspose             = MatMult_ConstantDiagonal;  in MatCreate_ConstantDiagonal()
383   A->ops->multtransposeadd          = MatMultAdd_ConstantDiagonal;  in MatCreate_ConstantDiagonal()
384   A->ops->multhermitiantranspose    = MatMultHermitianTranspose_ConstantDiagonal;  in MatCreate_ConstantDiagonal()
385   A->ops->multhermitiantransposeadd = MatMultHermitianTransposeAdd_ConstantDiagonal;  in MatCreate_ConstantDiagonal()
386   A->ops->solve                     = MatSolve_ConstantDiagonal;  in MatCreate_ConstantDiagonal()
387   A->ops->solvetranspose            = MatSolve_ConstantDiagonal;  in MatCreate_ConstantDiagonal()
388   A->ops->norm                      = MatNorm_ConstantDiagonal;  in MatCreate_ConstantDiagonal()
389   A->ops->createsubmatrices         = MatCreateSubMatrices_ConstantDiagonal;  in MatCreate_ConstantDiagonal()
390   A->ops->duplicate                 = MatDuplicate_ConstantDiagonal;  in MatCreate_ConstantDiagonal()
391   A->ops->getrow                    = MatGetRow_ConstantDiagonal;  in MatCreate_ConstantDiagonal()
392   A->ops->restorerow                = MatRestoreRow_ConstantDiagonal;  in MatCreate_ConstantDiagonal()
393   A->ops->sor                       = MatSOR_ConstantDiagonal;  in MatCreate_ConstantDiagonal()
394   A->ops->shift                     = MatShift_ConstantDiagonal;  in MatCreate_ConstantDiagonal()
395   A->ops->scale                     = MatScale_ConstantDiagonal;  in MatCreate_ConstantDiagonal()
396   A->ops->getdiagonal               = MatGetDiagonal_ConstantDiagonal;  in MatCreate_ConstantDiagonal()
397   A->ops->view                      = MatView_ConstantDiagonal;  in MatCreate_ConstantDiagonal()
398   A->ops->zeroentries               = MatZeroEntries_ConstantDiagonal;  in MatCreate_ConstantDiagonal()
399   A->ops->destroy                   = MatDestroy_ConstantDiagonal;  in MatCreate_ConstantDiagonal()
400   A->ops->getinfo                   = MatGetInfo_ConstantDiagonal;  in MatCreate_ConstantDiagonal()
401   A->ops->equal                     = MatEqual_ConstantDiagonal;  in MatCreate_ConstantDiagonal()
402   A->ops->axpy                      = MatAXPY_ConstantDiagonal;  in MatCreate_ConstantDiagonal()
403   A->ops->setrandom                 = MatSetRandom_ConstantDiagonal;  in MatCreate_ConstantDiagonal()
404   A->ops->conjugate                 = MatConjugate_ConstantDiagonal;  in MatCreate_ConstantDiagonal()
405   A->ops->transpose                 = MatTranspose_ConstantDiagonal;  in MatCreate_ConstantDiagonal()
420   fact->ops->solve = MatMult_ConstantDiagonal;  in MatFactorNumeric_ConstantDiagonal()
427   fact->ops->lufactornumeric = MatFactorNumeric_ConstantDiagonal;  in MatFactorSymbolic_LU_ConstantDiagonal()
434   fact->ops->choleskyfactornumeric = MatFactorNumeric_ConstantDiagonal;  in MatFactorSymbolic_Cholesky_ConstantDiagonal()
446   (*B)->ops->ilufactorsymbolic      = MatFactorSymbolic_LU_ConstantDiagonal;  in MatGetFactor_constantdiagonal_petsc()
447   (*B)->ops->lufactorsymbolic       = MatFactorSymbolic_LU_ConstantDiagonal;  in MatGetFactor_constantdiagonal_petsc()
448   (*B)->ops->iccfactorsymbolic      = MatFactorSymbolic_Cholesky_ConstantDiagonal;  in MatGetFactor_constantdiagonal_petsc()
449   (*B)->ops->choleskyfactorsymbolic = MatFactorSymbolic_Cholesky_ConstantDiagonal;  in MatGetFactor_constantdiagonal_petsc()
451   (*B)->ops->shift       = NULL;  in MatGetFactor_constantdiagonal_petsc()
452   (*B)->ops->scale       = NULL;  in MatGetFactor_constantdiagonal_petsc()
453   (*B)->ops->mult        = NULL;  in MatGetFactor_constantdiagonal_petsc()
454   (*B)->ops->sor         = NULL;  in MatGetFactor_constantdiagonal_petsc()
455   (*B)->ops->zeroentries = NULL;  in MatGetFactor_constantdiagonal_petsc()