Searched refs:len_slow1 (Results 1 – 1 of 1) sorted by relevance
396 PetscInt i, j, k, Mx, dof, xs, xm, len_slow1, len_slow2; in FVRHSFunctionslow2() local414 PetscCall(ISGetSize(ctx->iss, &len_slow1)); in FVRHSFunctionslow2()427 if (i < len_slow1 + len_slow2 + 1 && i > len_slow1 - 2) { in FVRHSFunctionslow2()457 if (i < len_slow1 + len_slow2 && i > len_slow1) { /* slow parts can be changed based on a */ in FVRHSFunctionslow2()467 for (j = 0; j < dof; j++) f[(i - len_slow1 - 1) * dof + j] -= ctx->flux[j] / hx; in FVRHSFunctionslow2()470 for (j = 0; j < dof; j++) f[(i - len_slow1) * dof + j] += ctx->flux[j] / hx; in FVRHSFunctionslow2()473 if (i == len_slow1 + len_slow2) { /* slow parts can be changed based on a */ in FVRHSFunctionslow2()483 for (j = 0; j < dof; j++) f[(i - len_slow1 - 1) * dof + j] -= ctx->flux[j] / hx; in FVRHSFunctionslow2()486 if (i == len_slow1) { /* slow parts can be changed based on a */ in FVRHSFunctionslow2()496 for (j = 0; j < dof; j++) f[(i - len_slow1) * dof + j] += ctx->flux[j] / hx; in FVRHSFunctionslow2()[all …]