Lines Matching refs:icol
564 const PetscInt *lidx, *icol, *gidx; in PCGAMGProlongator_Classical_Standard() local
632 PetscCall(MatGetRow(lA, i, &ncols, &icol, NULL)); in PCGAMGProlongator_Classical_Standard()
642 if (lcid[icol[j]] >= 0) { in PCGAMGProlongator_Classical_Standard()
643 pcontrib[icol[j]] = 1.; in PCGAMGProlongator_Classical_Standard()
645 ci = icol[j]; in PCGAMGProlongator_Classical_Standard()
646 PetscCall(MatRestoreRow(lA, i, &ncols, &icol, NULL)); in PCGAMGProlongator_Classical_Standard()
647 PetscCall(MatGetRow(lA, ci, &ncols, &icol, NULL)); in PCGAMGProlongator_Classical_Standard()
649 if (lcid[icol[k]] >= 0) pcontrib[icol[k]] = 1.; in PCGAMGProlongator_Classical_Standard()
651 PetscCall(MatRestoreRow(lA, ci, &ncols, &icol, NULL)); in PCGAMGProlongator_Classical_Standard()
652 PetscCall(MatGetRow(lA, i, &ncols, &icol, NULL)); in PCGAMGProlongator_Classical_Standard()
656 if (lcid[icol[j]] >= 0 && pcontrib[icol[j]] != 0.) { in PCGAMGProlongator_Classical_Standard()
657 lni = lcid[icol[j]]; in PCGAMGProlongator_Classical_Standard()
663 pcontrib[icol[j]] = 0.; in PCGAMGProlongator_Classical_Standard()
665 ci = icol[j]; in PCGAMGProlongator_Classical_Standard()
666 PetscCall(MatRestoreRow(lA, i, &ncols, &icol, NULL)); in PCGAMGProlongator_Classical_Standard()
667 PetscCall(MatGetRow(lA, ci, &ncols, &icol, NULL)); in PCGAMGProlongator_Classical_Standard()
669 if (lcid[icol[k]] >= 0 && pcontrib[icol[k]] != 0.) { in PCGAMGProlongator_Classical_Standard()
670 lni = lcid[icol[k]]; in PCGAMGProlongator_Classical_Standard()
676 pcontrib[icol[k]] = 0.; in PCGAMGProlongator_Classical_Standard()
679 PetscCall(MatRestoreRow(lA, ci, &ncols, &icol, NULL)); in PCGAMGProlongator_Classical_Standard()
680 PetscCall(MatGetRow(lA, i, &ncols, &icol, NULL)); in PCGAMGProlongator_Classical_Standard()
686 PetscCall(MatRestoreRow(lA, i, &ncols, &icol, &vcol)); in PCGAMGProlongator_Classical_Standard()
705 PetscCall(MatGetRow(lA, i, &ncols, &icol, &vcol)); in PCGAMGProlongator_Classical_Standard()
708 if (lcid[icol[j]] >= 0) { in PCGAMGProlongator_Classical_Standard()
710 pcontrib[icol[j]] += pentry; in PCGAMGProlongator_Classical_Standard()
711 } else if (icol[j] != i) { in PCGAMGProlongator_Classical_Standard()
713 ci = icol[j]; in PCGAMGProlongator_Classical_Standard()
715 PetscCall(MatRestoreRow(lA, i, &ncols, &icol, &vcol)); in PCGAMGProlongator_Classical_Standard()
716 PetscCall(MatGetRow(lA, ci, &ncols, &icol, &vcol)); in PCGAMGProlongator_Classical_Standard()
720 if (ci == icol[k]) jdiag = PetscRealPart(vcol[k]); in PCGAMGProlongator_Classical_Standard()
723 if (lcid[icol[k]] >= 0 && jdiag * PetscRealPart(vcol[k]) < 0.) { in PCGAMGProlongator_Classical_Standard()
731 if (lcid[icol[k]] >= 0 && jdiag * PetscRealPart(vcol[k]) < 0.) { in PCGAMGProlongator_Classical_Standard()
733 pcontrib[icol[k]] += pjentry; in PCGAMGProlongator_Classical_Standard()
739 PetscCall(MatRestoreRow(lA, ci, &ncols, &icol, &vcol)); in PCGAMGProlongator_Classical_Standard()
740 PetscCall(MatGetRow(lA, i, &ncols, &icol, &vcol)); in PCGAMGProlongator_Classical_Standard()
748 if (lcid[icol[j]] >= 0 && pcontrib[icol[j]] != 0.) { in PCGAMGProlongator_Classical_Standard()
750 if (PetscAbsScalar(pcontrib[icol[j]]) > 0.0) { in PCGAMGProlongator_Classical_Standard()
751 lni = lcid[icol[j]]; in PCGAMGProlongator_Classical_Standard()
752 pvcol[pncols] = -pcontrib[icol[j]] * diag; in PCGAMGProlongator_Classical_Standard()
756 pcontrib[icol[j]] = 0.; in PCGAMGProlongator_Classical_Standard()
759 ci = icol[j]; in PCGAMGProlongator_Classical_Standard()
760 PetscCall(MatRestoreRow(lA, i, &ncols, &icol, &vcol)); in PCGAMGProlongator_Classical_Standard()
761 PetscCall(MatGetRow(lA, ci, &ncols, &icol, &vcol)); in PCGAMGProlongator_Classical_Standard()
763 if (lcid[icol[k]] >= 0 && pcontrib[icol[k]] != 0.) { in PCGAMGProlongator_Classical_Standard()
764 if (PetscAbsScalar(pcontrib[icol[k]]) > 0.0) { in PCGAMGProlongator_Classical_Standard()
765 lni = lcid[icol[k]]; in PCGAMGProlongator_Classical_Standard()
766 pvcol[pncols] = -pcontrib[icol[k]] * diag; in PCGAMGProlongator_Classical_Standard()
770 pcontrib[icol[k]] = 0.; in PCGAMGProlongator_Classical_Standard()
773 PetscCall(MatRestoreRow(lA, ci, &ncols, &icol, &vcol)); in PCGAMGProlongator_Classical_Standard()
774 PetscCall(MatGetRow(lA, i, &ncols, &icol, &vcol)); in PCGAMGProlongator_Classical_Standard()
776 pcontrib[icol[j]] = 0.; in PCGAMGProlongator_Classical_Standard()
778 PetscCall(MatRestoreRow(lA, i, &ncols, &icol, &vcol)); in PCGAMGProlongator_Classical_Standard()