Searched refs:dxm (Results 1 – 2 of 2) sorted by relevance
| /petsc/src/dm/impls/da/ |
| H A D | da.c | 887 const PetscInt *dxm = (const PetscInt *)idxm; in DMDAMapMatStencilToGlobal() local
915 … dxm += 3 - dim; /* Input is {k,j,i,c}; move the pointer to the first used index, e.g., j in 2D */ in DMDAMapMatStencilToGlobal()
918 …if (tmp < 0 || dxm[j] < starts[j] || dxm[j] >= (starts[j] + dims[j])) tmp = -1; /* Beyond the ghos… in DMDAMapMatStencilToGlobal()
919 else tmp = tmp * dims[j] + (dxm[j] - starts[j]); in DMDAMapMatStencilToGlobal()
923 if (dof > 1) dxm += sdim; /* c is already counted in sdim */ in DMDAMapMatStencilToGlobal()
924 else dxm += sdim + 1; /* skip the unused c */ in DMDAMapMatStencilToGlobal()
|
| /petsc/src/mat/interface/ |
| H A D | matrix.c | 1792 …PetscInt *starts = mat->stencil.starts, *dxm = (PetscInt *)idxm, *dxn = (PetscInt *)idxn, sdim = d… in MatSetValuesStencil() local
1810 for (j = 0; j < 3 - sdim; j++) dxm++; in MatSetValuesStencil()
1811 tmp = *dxm++ - starts[0]; in MatSetValuesStencil()
1813 if ((*dxm++ - starts[j + 1]) < 0 || tmp < 0) tmp = -1; in MatSetValuesStencil()
1814 else tmp = tmp * dims[j] + *(dxm - 1) - starts[j + 1]; in MatSetValuesStencil()
1816 if (mat->stencil.noc) dxm++; in MatSetValuesStencil()
1907 …PetscInt *starts = mat->stencil.starts, *dxm = (PetscInt *)idxm, *dxn = (PetscInt *)idxn, sdim = d… in MatSetValuesBlockedStencil() local
1926 for (j = 0; j < 3 - sdim; j++) dxm++; in MatSetValuesBlockedStencil()
1927 tmp = *dxm++ - starts[0]; in MatSetValuesBlockedStencil()
1929 if ((*dxm++ - starts[j + 1]) < 0 || tmp < 0) tmp = -1; in MatSetValuesBlockedStencil()
[all …]
|