Lines Matching refs:jj
357 PetscInt arrsz, kk, ii, jj, nloc, ndatarows, ndf; in PCSetCoordinates_AGG() local
390 for (jj = 0; jj < ndatarows; jj++) { in PCSetCoordinates_AGG()
391 if (ii == jj) data[ii * M + jj] = 1.0; in PCSetCoordinates_AGG()
392 else data[ii * M + jj] = 0.0; in PCSetCoordinates_AGG()
513 …PetscInt Istart, my0, Iend, nloc, clid, flid = 0, aggID, kk, jj, ii, mm, nSelected, minsz, ng… in formProl0() local
539 PetscCall(MatGetOwnershipRangeColumn(a_Prol, &ii, &jj)); in formProl0()
541 …jj - ii) / nSAvec, PETSC_COMM_SELF, PETSC_ERR_PLIB, "nSelected %" PetscInt_FMT " != (jj %" PetscIn… in formProl0()
553 PetscCall(PetscCDCountAt(agg_llists, mm, &jj)); in formProl0()
554 if (jj > 0) { in formProl0()
563 PetscCall(PetscBLASIntCast(jj, &asz)); in formProl0()
591 for (jj = 0; jj < N; jj++) { in formProl0()
592 PetscReal d = data[jj * data_stride + ii]; in formProl0()
594 qqc[jj * Mdata + aggID * bs + ii] = d; in formProl0()
604 for (jj = 0; jj < N; jj++, kk++) qqc[jj * Mdata + ii] = .0; in formProl0()
616 for (jj = 0; jj < nSAvec; jj++) { in formProl0()
618 …PetscCheck(data[jj * out_data_stride + ii] == PETSC_MAX_REAL, PETSC_COMM_SELF, PETSC_ERR_PLIB, "da… in formProl0()
619 if (ii <= jj) data[jj * out_data_stride + ii] = PetscRealPart(qqc[jj * Mdata + ii]); in formProl0()
620 else data[jj * out_data_stride + ii] = 0.; in formProl0()
630 for (jj = 0; jj < N; jj++) qqr[N * ii + jj] = qqc[jj * Mdata + ii]; in formProl0()
783 for (PetscInt row = 0, grow = Istart, ncol_row, jj; row < nloc; row++, grow++) { in PCGAMGCreateGraph_AGG() local
785 for (ncol_row = jj = 0; jj < ncols; jj++, nnz0++) { in PCGAMGCreateGraph_AGG()
786 PetscScalar sv = PetscAbs(PetscRealPart(vals[jj])); in PCGAMGCreateGraph_AGG()
788 PetscInt cid = idx[jj] + Istart; //diag in PCGAMGCreateGraph_AGG()
791 if (c != a) cid = garray[idx[jj]]; in PCGAMGCreateGraph_AGG()
792 AA[ncol_row] = vals[jj]; in PCGAMGCreateGraph_AGG()
1282 PetscInt Istart, Iend, nloc, ii, jj, kk, my0, nLocalSelected, bs; in PCGAMGConstructProlongator_AGG() local
1348 for (jj = 0; jj < col_bs; jj++) { in PCGAMGConstructProlongator_AGG()
1351 const PetscReal *tp = PetscSafePointerPlusOffset(pc_gamg->data, jj * bs * nloc + kk); in PCGAMGConstructProlongator_AGG()
1357 …if (!jj && !kk) { /* now I know how many total nodes - allocate TODO: move below and do in one 'co… in PCGAMGConstructProlongator_AGG()
1361 tp2 = PetscSafePointerPlusOffset(data_w_ghost, jj * bs * stride + kk); in PCGAMGConstructProlongator_AGG()
1428 PetscInt jj; in PCGAMGOptimizeProlongator_AGG() local
1510 for (jj = 0; jj < pc_gamg_agg->nsmooths; jj++) { in PCGAMGOptimizeProlongator_AGG()