55   if (ptap->AP_loc) { /* used by alg_rap */  in MatProductCtxDestroy_MPIAIJ_PtAP()
56     Mat_SeqAIJ *ap = (Mat_SeqAIJ *)ptap->AP_loc->data;  in MatProductCtxDestroy_MPIAIJ_PtAP()
59     PetscCall(MatDestroy(&ptap->AP_loc));  in MatProductCtxDestroy_MPIAIJ_PtAP()
102   Mat                  AP_loc, C_loc, C_oth;  in MatPtAPNumeric_MPIAIJ_MPIAIJ_scalable()  local
112 …PetscCheck(ptap->AP_loc, PetscObjectComm((PetscObject)C), PETSC_ERR_ARG_WRONGSTATE, "PtAP cannot b…  in MatPtAPNumeric_MPIAIJ_MPIAIJ_scalable()
123   AP_loc = ptap->AP_loc;  in MatPtAPNumeric_MPIAIJ_MPIAIJ_scalable()
124   ap     = (Mat_SeqAIJ *)AP_loc->data;  in MatPtAPNumeric_MPIAIJ_MPIAIJ_scalable()
140   PetscCall(ISLocalToGlobalMappingApply(ptap->ltog, api[AP_loc->rmap->n], apj, apj));  in MatPtAPNumeric_MPIAIJ_MPIAIJ_scalable()
148 …PetscCall(ISGlobalToLocalMappingApply(ptap->ltog, IS_GTOLM_DROP, api[AP_loc->rmap->n], apj, &nout,…  in MatPtAPNumeric_MPIAIJ_MPIAIJ_scalable()
149 …Check(api[AP_loc->rmap->n] == nout, PETSC_COMM_SELF, PETSC_ERR_ARG_INCOMP, "Incorrect mapping %" P…  in MatPtAPNumeric_MPIAIJ_MPIAIJ_scalable()
153   PetscCall(MatMatMultNumeric_SeqAIJ_SeqAIJ_Scalable(ptap->Rd, AP_loc, ptap->C_loc));  in MatPtAPNumeric_MPIAIJ_MPIAIJ_scalable()
154   PetscCall(MatMatMultNumeric_SeqAIJ_SeqAIJ_Scalable(ptap->Ro, AP_loc, ptap->C_oth));  in MatPtAPNumeric_MPIAIJ_MPIAIJ_scalable()
348   PetscCall(MatCreateSeqAIJWithArrays(PETSC_COMM_SELF, am, pN, api, apj, apv, &ptap->AP_loc));  in MatPtAPSymbolic_MPIAIJ_MPIAIJ_scalable()
349   PetscCall(MatSeqAIJCompactOutExtraColumns_SeqAIJ(ptap->AP_loc, &ptap->ltog));  in MatPtAPSymbolic_MPIAIJ_MPIAIJ_scalable()
357   ptap->AP_loc->info.mallocs           = nspacedouble;  in MatPtAPSymbolic_MPIAIJ_MPIAIJ_scalable()
358   ptap->AP_loc->info.fill_ratio_given  = fill;  in MatPtAPSymbolic_MPIAIJ_MPIAIJ_scalable()
359   ptap->AP_loc->info.fill_ratio_needed = apfill;  in MatPtAPSymbolic_MPIAIJ_MPIAIJ_scalable()
362 …PetscCall(PetscInfo(ptap->AP_loc, "Scalable algorithm, AP_loc reallocs %" PetscInt_FMT "; Fill rat…  in MatPtAPSymbolic_MPIAIJ_MPIAIJ_scalable()
363 …PetscCall(PetscInfo(ptap->AP_loc, "Use MatPtAP(A,B,MatReuse,%g,&C) for best AP_loc performance.;\n…  in MatPtAPSymbolic_MPIAIJ_MPIAIJ_scalable()
365     PetscCall(PetscInfo(ptap->AP_loc, "Scalable algorithm, AP_loc is empty \n"));  in MatPtAPSymbolic_MPIAIJ_MPIAIJ_scalable()
370   PetscCall(MatProductCreate(ptap->Ro, ptap->AP_loc, NULL, &ptap->C_oth));  in MatPtAPSymbolic_MPIAIJ_MPIAIJ_scalable()
434   PetscCall(MatProductCreate(ptap->Rd, ptap->AP_loc, NULL, &ptap->C_loc));  in MatPtAPSymbolic_MPIAIJ_MPIAIJ_scalable()
1637   PetscCall(MatCreateSeqAIJWithArrays(PETSC_COMM_SELF, am, pN, api, apj, apv, &ptap->AP_loc));  in MatPtAPSymbolic_MPIAIJ_MPIAIJ()
1638   PetscCall(MatSetType(ptap->AP_loc, ((PetscObject)p->A)->type_name));  in MatPtAPSymbolic_MPIAIJ_MPIAIJ()
1645   ptap->AP_loc->info.mallocs           = nspacedouble;  in MatPtAPSymbolic_MPIAIJ_MPIAIJ()
1646   ptap->AP_loc->info.fill_ratio_given  = fill;  in MatPtAPSymbolic_MPIAIJ_MPIAIJ()
1647   ptap->AP_loc->info.fill_ratio_needed = apfill;  in MatPtAPSymbolic_MPIAIJ_MPIAIJ()
1650 …PetscCall(PetscInfo(ptap->AP_loc, "Nonscalable algorithm, AP_loc reallocs %" PetscInt_FMT "; Fill …  in MatPtAPSymbolic_MPIAIJ_MPIAIJ()
1651 …PetscCall(PetscInfo(ptap->AP_loc, "Use MatPtAP(A,B,MatReuse,%g,&C) for best AP_loc performance.;\n…  in MatPtAPSymbolic_MPIAIJ_MPIAIJ()
1653     PetscCall(PetscInfo(ptap->AP_loc, "Nonscalable algorithm, AP_loc is empty \n"));  in MatPtAPSymbolic_MPIAIJ_MPIAIJ()
1661   PetscCall(MatProductCreate(ptap->Ro, ptap->AP_loc, NULL, &ptap->C_oth));  in MatPtAPSymbolic_MPIAIJ_MPIAIJ()
1725   PetscCall(MatProductCreate(ptap->Rd, ptap->AP_loc, NULL, &ptap->C_loc));  in MatPtAPSymbolic_MPIAIJ_MPIAIJ()
1863   Mat                  AP_loc, C_loc, C_oth;  in MatPtAPNumeric_MPIAIJ_MPIAIJ()  local
1874 …PetscCheck(ptap->AP_loc, PetscObjectComm((PetscObject)C), PETSC_ERR_ARG_WRONGSTATE, "PtAP cannot b…  in MatPtAPNumeric_MPIAIJ_MPIAIJ()
1884   AP_loc = ptap->AP_loc;  in MatPtAPNumeric_MPIAIJ_MPIAIJ()
1885   ap     = (Mat_SeqAIJ *)AP_loc->data;  in MatPtAPNumeric_MPIAIJ_MPIAIJ()
1902   PetscCall(MatSeqAIJGetArrayWrite(AP_loc, &apa));  in MatPtAPNumeric_MPIAIJ_MPIAIJ()
1917   PetscCall(MatSeqAIJRestoreArrayWrite(AP_loc, &apa));  in MatPtAPNumeric_MPIAIJ_MPIAIJ()
1924   PetscCall(PetscObjectStateIncrease((PetscObject)AP_loc));  in MatPtAPNumeric_MPIAIJ_MPIAIJ()