261   PetscInt    i, idx[3], M, xm;  in RHSMatrixHeat()  local
283     idx[0] = M - 1;  in RHSMatrixHeat()
284     idx[1] = 0;  in RHSMatrixHeat()
285     idx[2] = 1;  in RHSMatrixHeat()
286     PetscCall(MatSetValues(A, 1, &mstart, 3, idx, v, INSERT_VALUES));  in RHSMatrixHeat()
292     idx[0] = M - 2;  in RHSMatrixHeat()
293     idx[1] = M - 1;  in RHSMatrixHeat()
294     idx[2] = 0;  in RHSMatrixHeat()
295     PetscCall(MatSetValues(A, 1, &mend, 3, idx, v, INSERT_VALUES));  in RHSMatrixHeat()
303     idx[0] = i - 1;  in RHSMatrixHeat()
304     idx[1] = i;  in RHSMatrixHeat()
305     idx[2] = i + 1;  in RHSMatrixHeat()
306     PetscCall(MatSetValues(A, 1, &i, 3, idx, v, INSERT_VALUES));  in RHSMatrixHeat()
318       idx[0] = M - 1;  in RHSMatrixHeat()
319       idx[1] = 1;  in RHSMatrixHeat()
320       PetscCall(MatSetValues(A, 1, &mstart, 2, idx, v, ADD_VALUES));  in RHSMatrixHeat()
326       idx[0] = M - 2;  in RHSMatrixHeat()
327       idx[1] = 0;  in RHSMatrixHeat()
328       PetscCall(MatSetValues(A, 1, &mend, 2, idx, v, ADD_VALUES));  in RHSMatrixHeat()
332       idx[0] = i - 1;  in RHSMatrixHeat()
333       idx[1] = i + 1;  in RHSMatrixHeat()
334       PetscCall(MatSetValues(A, 1, &i, 2, idx, v, ADD_VALUES));  in RHSMatrixHeat()
341       idx[0] = 0;  in RHSMatrixHeat()
342       idx[1] = 1;  in RHSMatrixHeat()
343       PetscCall(MatSetValues(A, 1, &mstart, 2, idx, v, ADD_VALUES));  in RHSMatrixHeat()
349       idx[0] = M - 1;  in RHSMatrixHeat()
350       idx[1] = 0;  in RHSMatrixHeat()
351       PetscCall(MatSetValues(A, 1, &mend, 2, idx, v, ADD_VALUES));  in RHSMatrixHeat()
355       idx[0] = i;  in RHSMatrixHeat()
356       idx[1] = i + 1;  in RHSMatrixHeat()
357       PetscCall(MatSetValues(A, 1, &i, 2, idx, v, ADD_VALUES));  in RHSMatrixHeat()