Lines Matching refs:slope
225 PetscScalar *x, *f, *slope; in FVRHSFunction_3WaySplit() local
243 …PetscCall(DMDAGetArray(da, PETSC_TRUE, &slope)); /* contains ghost points … in FVRHSFunction_3WaySplit()
282 slope[i * dof + j] = tmp; in FVRHSFunction_3WaySplit()
293 uL[j] = x[(i - 1) * dof + j] + slope[(i - 1) * dof + j] * hxs / 2; in FVRHSFunction_3WaySplit()
294 uR[j] = x[(i - 0) * dof + j] - slope[(i - 0) * dof + j] * hxs / 2; in FVRHSFunction_3WaySplit()
305 uL[j] = x[(i - 1) * dof + j] + slope[(i - 1) * dof + j] * hxs / 2; in FVRHSFunction_3WaySplit()
306 uR[j] = x[(i - 0) * dof + j] - slope[(i - 0) * dof + j] * hxm / 2; in FVRHSFunction_3WaySplit()
317 uL[j] = x[(i - 1) * dof + j] + slope[(i - 1) * dof + j] * hxm / 2; in FVRHSFunction_3WaySplit()
318 uR[j] = x[(i - 0) * dof + j] - slope[(i - 0) * dof + j] * hxs / 2; in FVRHSFunction_3WaySplit()
329 uL[j] = x[(i - 1) * dof + j] + slope[(i - 1) * dof + j] * hxm / 2; in FVRHSFunction_3WaySplit()
330 uR[j] = x[(i - 0) * dof + j] - slope[(i - 0) * dof + j] * hxm / 2; in FVRHSFunction_3WaySplit()
341 uL[j] = x[(i - 1) * dof + j] + slope[(i - 1) * dof + j] * hxm / 2; in FVRHSFunction_3WaySplit()
342 uR[j] = x[(i - 0) * dof + j] - slope[(i - 0) * dof + j] * hxf / 2; in FVRHSFunction_3WaySplit()
353 uL[j] = x[(i - 1) * dof + j] + slope[(i - 1) * dof + j] * hxf / 2; in FVRHSFunction_3WaySplit()
354 uR[j] = x[(i - 0) * dof + j] - slope[(i - 0) * dof + j] * hxm / 2; in FVRHSFunction_3WaySplit()
365 uL[j] = x[(i - 1) * dof + j] + slope[(i - 1) * dof + j] * hxf / 2; in FVRHSFunction_3WaySplit()
366 uR[j] = x[(i - 0) * dof + j] - slope[(i - 0) * dof + j] * hxf / 2; in FVRHSFunction_3WaySplit()
380 PetscCall(DMDARestoreArray(da, PETSC_TRUE, &slope)); in FVRHSFunction_3WaySplit()
399 PetscScalar *x, *f, *slope; in FVRHSFunctionslow_3WaySplit() local
415 PetscCall(DMDAGetArray(da, PETSC_TRUE, &slope)); in FVRHSFunctionslow_3WaySplit()
454 slope[i * dof + j] = tmp; in FVRHSFunctionslow_3WaySplit()
466 uL[j] = x[(i - 1) * dof + j] + slope[(i - 1) * dof + j] * hxs / 2; in FVRHSFunctionslow_3WaySplit()
467 uR[j] = x[(i - 0) * dof + j] - slope[(i - 0) * dof + j] * hxs / 2; in FVRHSFunctionslow_3WaySplit()
481 uL[j] = x[(i - 1) * dof + j] + slope[(i - 1) * dof + j] * hxs / 2; in FVRHSFunctionslow_3WaySplit()
482 uR[j] = x[(i - 0) * dof + j] - slope[(i - 0) * dof + j] * hxs / 2; in FVRHSFunctionslow_3WaySplit()
491 uL[j] = x[(i - 1) * dof + j] + slope[(i - 1) * dof + j] * hxs / 2; in FVRHSFunctionslow_3WaySplit()
492 uR[j] = x[(i - 0) * dof + j] - slope[(i - 0) * dof + j] * hxs / 2; in FVRHSFunctionslow_3WaySplit()
502 uL[j] = x[(i - 1) * dof + j] + slope[(i - 1) * dof + j] * hxs / 2; in FVRHSFunctionslow_3WaySplit()
503 uR[j] = x[(i - 0) * dof + j] - slope[(i - 0) * dof + j] * hxs / 2; in FVRHSFunctionslow_3WaySplit()
518 PetscCall(DMDARestoreArray(da, PETSC_TRUE, &slope)); in FVRHSFunctionslow_3WaySplit()
529 PetscScalar *x, *f, *slope; in FVRHSFunctionslowbuffer_3WaySplit() local
545 PetscCall(DMDAGetArray(da, PETSC_TRUE, &slope)); in FVRHSFunctionslowbuffer_3WaySplit()
584 slope[i * dof + j] = tmp; in FVRHSFunctionslowbuffer_3WaySplit()
596 uL[j] = x[(i - 1) * dof + j] + slope[(i - 1) * dof + j] * hxs / 2; in FVRHSFunctionslowbuffer_3WaySplit()
597 uR[j] = x[(i - 0) * dof + j] - slope[(i - 0) * dof + j] * hxs / 2; in FVRHSFunctionslowbuffer_3WaySplit()
607 uL[j] = x[(i - 1) * dof + j] + slope[(i - 1) * dof + j] * hxs / 2; in FVRHSFunctionslowbuffer_3WaySplit()
608 uR[j] = x[(i - 0) * dof + j] - slope[(i - 0) * dof + j] * hxs / 2; in FVRHSFunctionslowbuffer_3WaySplit()
621 uL[j] = x[(i - 1) * dof + j] + slope[(i - 1) * dof + j] * hxs / 2; in FVRHSFunctionslowbuffer_3WaySplit()
622 uR[j] = x[(i - 0) * dof + j] - slope[(i - 0) * dof + j] * hxm / 2; in FVRHSFunctionslowbuffer_3WaySplit()
631 uL[j] = x[(i - 1) * dof + j] + slope[(i - 1) * dof + j] * hxm / 2; in FVRHSFunctionslowbuffer_3WaySplit()
632 uR[j] = x[(i - 0) * dof + j] - slope[(i - 0) * dof + j] * hxs / 2; in FVRHSFunctionslowbuffer_3WaySplit()
642 uL[j] = x[(i - 1) * dof + j] + slope[(i - 1) * dof + j] * hxs / 2; in FVRHSFunctionslowbuffer_3WaySplit()
643 uR[j] = x[(i - 0) * dof + j] - slope[(i - 0) * dof + j] * hxs / 2; in FVRHSFunctionslowbuffer_3WaySplit()
656 uL[j] = x[(i - 1) * dof + j] + slope[(i - 1) * dof + j] * hxs / 2; in FVRHSFunctionslowbuffer_3WaySplit()
657 uR[j] = x[(i - 0) * dof + j] - slope[(i - 0) * dof + j] * hxs / 2; in FVRHSFunctionslowbuffer_3WaySplit()
667 PetscCall(DMDARestoreArray(da, PETSC_TRUE, &slope)); in FVRHSFunctionslowbuffer_3WaySplit()
678 PetscScalar *x, *f, *slope; in FVRHSFunctionmedium_3WaySplit() local
694 PetscCall(DMDAGetArray(da, PETSC_TRUE, &slope)); in FVRHSFunctionmedium_3WaySplit()
733 slope[i * dof + j] = tmp; in FVRHSFunctionmedium_3WaySplit()
745 uL[j] = x[(i - 1) * dof + j] + slope[(i - 1) * dof + j] * hxs / 2; in FVRHSFunctionmedium_3WaySplit()
746 uR[j] = x[(i - 0) * dof + j] - slope[(i - 0) * dof + j] * hxm / 2; in FVRHSFunctionmedium_3WaySplit()
756 uL[j] = x[(i - 1) * dof + j] + slope[(i - 1) * dof + j] * hxm / 2; in FVRHSFunctionmedium_3WaySplit()
757 uR[j] = x[(i - 0) * dof + j] - slope[(i - 0) * dof + j] * hxm / 2; in FVRHSFunctionmedium_3WaySplit()
770 uL[j] = x[(i - 1) * dof + j] + slope[(i - 1) * dof + j] * hxm / 2; in FVRHSFunctionmedium_3WaySplit()
771 uR[j] = x[(i - 0) * dof + j] - slope[(i - 0) * dof + j] * hxm / 2; in FVRHSFunctionmedium_3WaySplit()
780 uL[j] = x[(i - 1) * dof + j] + slope[(i - 1) * dof + j] * hxm / 2; in FVRHSFunctionmedium_3WaySplit()
781 uR[j] = x[(i - 0) * dof + j] - slope[(i - 0) * dof + j] * hxm / 2; in FVRHSFunctionmedium_3WaySplit()
791 uL[j] = x[(i - 1) * dof + j] + slope[(i - 1) * dof + j] * hxm / 2; in FVRHSFunctionmedium_3WaySplit()
792 uR[j] = x[(i - 0) * dof + j] - slope[(i - 0) * dof + j] * hxm / 2; in FVRHSFunctionmedium_3WaySplit()
805 uL[j] = x[(i - 1) * dof + j] + slope[(i - 1) * dof + j] * hxm / 2; in FVRHSFunctionmedium_3WaySplit()
806 uR[j] = x[(i - 0) * dof + j] - slope[(i - 0) * dof + j] * hxs / 2; in FVRHSFunctionmedium_3WaySplit()
816 PetscCall(DMDARestoreArray(da, PETSC_TRUE, &slope)); in FVRHSFunctionmedium_3WaySplit()
826 PetscScalar *x, *f, *slope; in FVRHSFunctionmediumbuffer_3WaySplit() local
842 PetscCall(DMDAGetArray(da, PETSC_TRUE, &slope)); in FVRHSFunctionmediumbuffer_3WaySplit()
881 slope[i * dof + j] = tmp; in FVRHSFunctionmediumbuffer_3WaySplit()
893 uL[j] = x[(i - 1) * dof + j] + slope[(i - 1) * dof + j] * hxm / 2; in FVRHSFunctionmediumbuffer_3WaySplit()
894 uR[j] = x[(i - 0) * dof + j] - slope[(i - 0) * dof + j] * hxm / 2; in FVRHSFunctionmediumbuffer_3WaySplit()
904 uL[j] = x[(i - 1) * dof + j] + slope[(i - 1) * dof + j] * hxm / 2; in FVRHSFunctionmediumbuffer_3WaySplit()
905 uR[j] = x[(i - 0) * dof + j] - slope[(i - 0) * dof + j] * hxm / 2; in FVRHSFunctionmediumbuffer_3WaySplit()
918 uL[j] = x[(i - 1) * dof + j] + slope[(i - 1) * dof + j] * hxm / 2; in FVRHSFunctionmediumbuffer_3WaySplit()
919 uR[j] = x[(i - 0) * dof + j] - slope[(i - 0) * dof + j] * hxf / 2; in FVRHSFunctionmediumbuffer_3WaySplit()
928 uL[j] = x[(i - 1) * dof + j] + slope[(i - 1) * dof + j] * hxf / 2; in FVRHSFunctionmediumbuffer_3WaySplit()
929 uR[j] = x[(i - 0) * dof + j] - slope[(i - 0) * dof + j] * hxm / 2; in FVRHSFunctionmediumbuffer_3WaySplit()
939 uL[j] = x[(i - 1) * dof + j] + slope[(i - 1) * dof + j] * hxm / 2; in FVRHSFunctionmediumbuffer_3WaySplit()
940 uR[j] = x[(i - 0) * dof + j] - slope[(i - 0) * dof + j] * hxm / 2; in FVRHSFunctionmediumbuffer_3WaySplit()
953 uL[j] = x[(i - 1) * dof + j] + slope[(i - 1) * dof + j] * hxm / 2; in FVRHSFunctionmediumbuffer_3WaySplit()
954 uR[j] = x[(i - 0) * dof + j] - slope[(i - 0) * dof + j] * hxm / 2; in FVRHSFunctionmediumbuffer_3WaySplit()
964 PetscCall(DMDARestoreArray(da, PETSC_TRUE, &slope)); in FVRHSFunctionmediumbuffer_3WaySplit()
975 PetscScalar *x, *f, *slope; in FVRHSFunctionfast_3WaySplit() local
991 PetscCall(DMDAGetArray(da, PETSC_TRUE, &slope)); in FVRHSFunctionfast_3WaySplit()
1028 slope[i * dof + j] = tmp; in FVRHSFunctionfast_3WaySplit()
1040 uL[j] = x[(i - 1) * dof + j] + slope[(i - 1) * dof + j] * hxm / 2; in FVRHSFunctionfast_3WaySplit()
1041 uR[j] = x[(i - 0) * dof + j] - slope[(i - 0) * dof + j] * hxf / 2; in FVRHSFunctionfast_3WaySplit()
1051 uL[j] = x[(i - 1) * dof + j] + slope[(i - 1) * dof + j] * hxf / 2; in FVRHSFunctionfast_3WaySplit()
1052 uR[j] = x[(i - 0) * dof + j] - slope[(i - 0) * dof + j] * hxf / 2; in FVRHSFunctionfast_3WaySplit()
1065 uL[j] = x[(i - 1) * dof + j] + slope[(i - 1) * dof + j] * hxf / 2; in FVRHSFunctionfast_3WaySplit()
1066 uR[j] = x[(i - 0) * dof + j] - slope[(i - 0) * dof + j] * hxm / 2; in FVRHSFunctionfast_3WaySplit()
1076 PetscCall(DMDARestoreArray(da, PETSC_TRUE, &slope)); in FVRHSFunctionfast_3WaySplit()