Lines Matching refs:shift
207 PetscInt shift; in MatGetRow_SeqSELL() local
212 shift = a->sliidx[row / a->sliceheight] + (row % a->sliceheight); in MatGetRow_SeqSELL()
216 for (j = 0; j < a->rlen[row]; j++) a->getrowcols[j] = a->colidx[shift + a->sliceheight * j]; in MatGetRow_SeqSELL()
221 for (j = 0; j < a->rlen[row]; j++) a->getrowvals[j] = a->val[shift + a->sliceheight * j]; in MatGetRow_SeqSELL()
781 PetscInt shift; in MatGetDiagonalMarkers_SeqSELL() local
788 … shift = a->sliidx[i / a->sliceheight] + i % a->sliceheight; /* starting index of the row i */ in MatGetDiagonalMarkers_SeqSELL()
790 if (a->colidx[shift + a->sliceheight * j] == i) { in MatGetDiagonalMarkers_SeqSELL()
791 a->diag[i] = shift + a->sliceheight * j; in MatGetDiagonalMarkers_SeqSELL()
809 … shift = a->sliidx[i / a->sliceheight] + i % a->sliceheight; /* starting index of the row i */ in MatGetDiagonalMarkers_SeqSELL()
812 if (a->colidx[shift + a->sliceheight * j] == i) { in MatGetDiagonalMarkers_SeqSELL()
813 a->diag[i] = shift + a->sliceheight * j; in MatGetDiagonalMarkers_SeqSELL()
963 PetscInt i, j, n, shift; in MatGetDiagonal_SeqSELL() local
983 shift = a->sliidx[i / a->sliceheight] + i % a->sliceheight; /* starting index of the row i */ in MatGetDiagonal_SeqSELL()
986 if (a->colidx[shift + a->sliceheight * j] == i) { in MatGetDiagonal_SeqSELL()
987 x[i] = a->val[shift + a->sliceheight * j]; in MatGetDiagonal_SeqSELL()
1047 PetscInt shift; in MatGetValues_SeqSELL() local
1055 … shift = a->sliidx[row / a->sliceheight] + (row % a->sliceheight); /* starting index of the row */ in MatGetValues_SeqSELL()
1056 cp = a->colidx + shift; /* pointer to the row */ in MatGetValues_SeqSELL()
1057 vp = a->val + shift; /* pointer to the row */ in MatGetValues_SeqSELL()
1086 PetscInt i, j, m = A->rmap->n, shift; in MatView_SeqSELL_ASCII() local
1108 shift = a->sliidx[i / a->sliceheight] + i % a->sliceheight; in MatView_SeqSELL_ASCII()
1111 … a->colidx[shift + a->sliceheight * j] + 1, (double)PetscRealPart(a->val[shift + a->sliceheight * … in MatView_SeqSELL_ASCII()
1113 …" %" PetscInt_FMT " %18.16e\n", i + 1, a->colidx[shift + a->sliceheight * j] + 1, (double)a->val[… in MatView_SeqSELL_ASCII()
1135 shift = a->sliidx[i / a->sliceheight] + i % a->sliceheight; in MatView_SeqSELL_ASCII()
1138 …if (PetscImaginaryPart(a->val[shift + a->sliceheight * j]) > 0.0 && PetscRealPart(a->val[shift + a… in MatView_SeqSELL_ASCII()
1139 …", a->colidx[shift + a->sliceheight * j], (double)PetscRealPart(a->val[shift + a->sliceheight * j]… in MatView_SeqSELL_ASCII()
1140 …} else if (PetscImaginaryPart(a->val[shift + a->sliceheight * j]) < 0.0 && PetscRealPart(a->val[sh… in MatView_SeqSELL_ASCII()
1141 …, a->colidx[shift + a->sliceheight * j], (double)PetscRealPart(a->val[shift + a->sliceheight * j])… in MatView_SeqSELL_ASCII()
1142 } else if (PetscRealPart(a->val[shift + a->sliceheight * j]) != 0.0) { in MatView_SeqSELL_ASCII()
1143 …ewer, " (%" PetscInt_FMT ", %g) ", a->colidx[shift + a->sliceheight * j], (double)PetscRealPart(a-… in MatView_SeqSELL_ASCII()
1146 …shift + a->sliceheight * j] != 0.0) PetscCall(PetscViewerASCIIPrintf(viewer, " (%" PetscInt_FMT ",… in MatView_SeqSELL_ASCII()
1168 shift = a->sliidx[i / a->sliceheight] + i % a->sliceheight; in MatView_SeqSELL_ASCII()
1170 if (jcnt < a->rlen[i] && j == a->colidx[shift + a->sliceheight * j]) { in MatView_SeqSELL_ASCII()
1199 shift = a->sliidx[i / a->sliceheight] + i % a->sliceheight; in MatView_SeqSELL_ASCII()
1202 …>colidx[shift + a->sliceheight * j] + fshift, (double)PetscRealPart(a->val[shift + a->sliceheight … in MatView_SeqSELL_ASCII()
1204 …%" PetscInt_FMT " %g\n", i + fshift, a->colidx[shift + a->sliceheight * j] + fshift, (double)a->va… in MatView_SeqSELL_ASCII()
1231 shift = a->sliidx[i / a->sliceheight] + i % a->sliceheight; in MatView_SeqSELL_ASCII()
1234 for (j = shift; j < a->diag[i]; j += a->sliceheight) { in MatView_SeqSELL_ASCII()
1236 if (PetscImaginaryPart(a->val[shift + a->sliceheight * j]) > 0.0) { in MatView_SeqSELL_ASCII()
1238 } else if (PetscImaginaryPart(a->val[shift + a->sliceheight * j]) < 0.0) { in MatView_SeqSELL_ASCII()
1262 for (j = a->diag[i] + 1; j < shift + a->sliceheight * a->rlen[i]; j += a->sliceheight) { in MatView_SeqSELL_ASCII()
1279 shift = a->sliidx[i / a->sliceheight] + i % a->sliceheight; in MatView_SeqSELL_ASCII()
1284 …", a->colidx[shift + a->sliceheight * j], (double)PetscRealPart(a->val[shift + a->sliceheight * j]… in MatView_SeqSELL_ASCII()
1286 …, a->colidx[shift + a->sliceheight * j], (double)PetscRealPart(a->val[shift + a->sliceheight * j])… in MatView_SeqSELL_ASCII()
1288 …ewer, " (%" PetscInt_FMT ", %g) ", a->colidx[shift + a->sliceheight * j], (double)PetscRealPart(a-… in MatView_SeqSELL_ASCII()
1291 …intf(viewer, " (%" PetscInt_FMT ", %g) ", a->colidx[shift + a->sliceheight * j], (double)a->val[sh… in MatView_SeqSELL_ASCII()
1308 PetscInt i, j, m = A->rmap->n, shift; in MatView_SeqSELL_Draw_Zoom() local
1326 shift = a->sliidx[i / a->sliceheight] + i % a->sliceheight; /* starting index of the row i */ in MatView_SeqSELL_Draw_Zoom()
1330 x_l = a->colidx[shift + a->sliceheight * j]; in MatView_SeqSELL_Draw_Zoom()
1332 if (PetscRealPart(a->val[shift + a->sliceheight * j]) >= 0.) continue; in MatView_SeqSELL_Draw_Zoom()
1338 shift = a->sliidx[i / a->sliceheight] + i % a->sliceheight; in MatView_SeqSELL_Draw_Zoom()
1342 x_l = a->colidx[shift + a->sliceheight * j]; in MatView_SeqSELL_Draw_Zoom()
1344 if (a->val[shift + a->sliceheight * j] != 0.) continue; in MatView_SeqSELL_Draw_Zoom()
1350 shift = a->sliidx[i / a->sliceheight] + i % a->sliceheight; in MatView_SeqSELL_Draw_Zoom()
1354 x_l = a->colidx[shift + a->sliceheight * j]; in MatView_SeqSELL_Draw_Zoom()
1356 if (PetscRealPart(a->val[shift + a->sliceheight * j]) <= 0.) continue; in MatView_SeqSELL_Draw_Zoom()
1376 shift = a->sliidx[i / a->sliceheight] + i % a->sliceheight; in MatView_SeqSELL_Draw_Zoom()
1380 x_l = a->colidx[shift + a->sliceheight * j]; in MatView_SeqSELL_Draw_Zoom()
1439 PetscInt i, shift, row_in_slice, row, nrow, *cp, lastcol, j, k; in MatAssemblyEnd_SeqSELL() local
1454 …shift = a->sliidx[i]; /* starting index of the s… in MatAssemblyEnd_SeqSELL()
1455 …cp = PetscSafePointerPlusOffset(a->colidx, shift); /* pointer to the column i… in MatAssemblyEnd_SeqSELL()
1456 …vp = PetscSafePointerPlusOffset(a->val, shift); /* pointer to the nonzero … in MatAssemblyEnd_SeqSELL()
1476 …if (a->sliidx[i + 1] != shift) lastcol = cp[row_in_slice - 1]; /* use the index from the previous … in MatAssemblyEnd_SeqSELL()
1479 for (k = nrow; k < (a->sliidx[i + 1] - shift) / a->sliceheight; ++k) { in MatAssemblyEnd_SeqSELL()
1536 PetscInt shift, i, k, l, low, high, t, ii, row, col, nrow; in MatSetValues_SeqSELL() local
1549 shift = a->sliidx[row / a->sliceheight] + row % a->sliceheight; /* starting index of the row */ in MatSetValues_SeqSELL()
1550 cp = a->colidx + shift; /* pointer to the row */ in MatSetValues_SeqSELL()
1551 vp = a->val + shift; /* pointer to the row */ in MatSetValues_SeqSELL()
1693 PetscInt n, m = A->rmap->n, i, j, shift; in MatSOR_SeqSELL() local
1715 … shift = a->sliidx[i / a->sliceheight] + i % a->sliceheight; /* starting index of the row i */ in MatSOR_SeqSELL()
1717 n = (diag[i] - shift) / a->sliceheight; in MatSOR_SeqSELL()
1718 …for (j = 0; j < n; j++) sum -= a->val[shift + a->sliceheight * j] * x[a->colidx[shift + a->slicehe… in MatSOR_SeqSELL()
1727 … shift = a->sliidx[i / a->sliceheight] + i % a->sliceheight; /* starting index of the row i */ in MatSOR_SeqSELL()
1729 n = a->rlen[i] - (diag[i] - shift) / a->sliceheight - 1; in MatSOR_SeqSELL()
1745 … shift = a->sliidx[i / a->sliceheight] + i % a->sliceheight; /* starting index of the row i */ in MatSOR_SeqSELL()
1747 n = (diag[i] - shift) / a->sliceheight; in MatSOR_SeqSELL()
1748 …for (j = 0; j < n; j++) sum -= a->val[shift + a->sliceheight * j] * x[a->colidx[shift + a->slicehe… in MatSOR_SeqSELL()
1751 n = a->rlen[i] - (diag[i] - shift) / a->sliceheight - 1; in MatSOR_SeqSELL()
1760 … shift = a->sliidx[i / a->sliceheight] + i % a->sliceheight; /* starting index of the row i */ in MatSOR_SeqSELL()
1765 …for (j = 0; j < n; j++) sum -= a->val[shift + a->sliceheight * j] * x[a->colidx[shift + a->slicehe… in MatSOR_SeqSELL()
1768 n = a->rlen[i] - (diag[i] - shift) / a->sliceheight - 1; in MatSOR_SeqSELL()