Searched refs:uglyrmapd (Results 1 – 1 of 1) sorted by relevance
196 static PetscInt *uglyrmapd = NULL, *uglyrmapo = NULL; /* mapping from the local ordering to the "di… variable218 PetscCall(PetscMalloc1(n, &uglyrmapd)); in MatMPIBAIJDiagonalScaleLocalSetUp()221 for (j = 0; j < bs; j++) uglyrmapd[(r_rmapd[i] - 1) * bs + j - cstart] = i * bs + j; in MatMPIBAIJDiagonalScaleLocalSetUp()259 if (!uglyrmapd) PetscCall(MatMPIBAIJDiagonalScaleLocalSetUp(A, scale)); in MatDiagonalScaleLocal_MPIBAIJ()265 …for (i = 0; i < n; i++) d[i] = s[uglyrmapd[i]]; /* copy "diagonal" (true local) portion of scale i… in MatDiagonalScaleLocal_MPIBAIJ()277 PetscCall(PetscFree(uglyrmapd)); in MatDiagonalScaleLocal_MPIBAIJ()