207   PetscInt                     *Mi, *Mj, Mnz;  in MatSeqAIJCUSPARSEBuildFactoredMatrix_LU()  local
217       PetscCall(PetscMalloc1(Mnz, &Mj)); // Mj is temp  in MatSeqAIJCUSPARSEBuildFactoredMatrix_LU()
223 …  PetscCall(PetscArraycpy(Mj + Mi[i], Aj + Ai[i], llen));                           // entries of L  in MatSeqAIJCUSPARSEBuildFactoredMatrix_LU()
224 …Mj[Mi[i] + llen] = i;                                                             // diagonal entry  in MatSeqAIJCUSPARSEBuildFactoredMatrix_LU()
225 …PetscCall(PetscArraycpy(Mj + Mi[i] + llen + 1, Aj + adiag[i + 1] + 1, ulen - 1)); // entries of U …  in MatSeqAIJCUSPARSEBuildFactoredMatrix_LU()
233 …PetscCallCUDA(cudaMemcpy(fs->csrColIdx, Mj, sizeof(*fs->csrColIdx) * Mnz, cudaMemcpyHostToDevice));  in MatSeqAIJCUSPARSEBuildFactoredMatrix_LU()
272       PetscCall(PetscFree(Mj));  in MatSeqAIJCUSPARSEBuildFactoredMatrix_LU()
625   PetscInt                     *Mj, Mnz;  in MatSeqAIJCUSPARSEBuildFactoredMatrix_Cholesky()  local
636       PetscCall(PetscMalloc1(Mnz, &Mj)); // Mj[] is temp  in MatSeqAIJCUSPARSEBuildFactoredMatrix_Cholesky()
640         Mj[Ai[i]]     = i;                                              // diagonal entry  in MatSeqAIJCUSPARSEBuildFactoredMatrix_Cholesky()
641 …PetscCall(PetscArraycpy(Mj + Ai[i] + 1, Aj + Ai[i], ulen - 1)); // entries of U on the right of th…  in MatSeqAIJCUSPARSEBuildFactoredMatrix_Cholesky()
649       PetscCallCUDA(cudaMemcpy(fs->csrColIdx, Mj, sizeof(*Mj) * Mnz, cudaMemcpyHostToDevice));  in MatSeqAIJCUSPARSEBuildFactoredMatrix_Cholesky()
683       PetscCall(PetscFree(Mj));  in MatSeqAIJCUSPARSEBuildFactoredMatrix_Cholesky()