aboutsummaryrefslogtreecommitdiff
path: root/residual_calc.h
diff options
context:
space:
mode:
authorAnton Khirnov <anton@khirnov.net>2019-06-14 15:29:02 +0200
committerAnton Khirnov <anton@khirnov.net>2019-06-14 16:10:14 +0200
commita1f10a1338e77dd147e7c94f3f401d2703181bca (patch)
tree6ee70770635391904b4146a4ed0a9f6c691609cd /residual_calc.h
parentd61fcbfbdb1517c39d6d894f687f087a8c728b27 (diff)
rescalc: improve reflection boundary conditions
Make parameter names more clear/consistent, document them, implement missing 1U boundary.
Diffstat (limited to 'residual_calc.h')
-rw-r--r--residual_calc.h9
1 files changed, 8 insertions, 1 deletions
diff --git a/residual_calc.h b/residual_calc.h
index 8c9b628..1498bc9 100644
--- a/residual_calc.h
+++ b/residual_calc.h
@@ -67,6 +67,12 @@ int mg2di_residual_calc_init(ResidualCalcContext *ctx);
/**
* Calculate the residual.
+ *
+ * @param reflect Flags indicating if dst should be extended by reflection.
+ * If the 'MG2D_BOUNDARY_x'th bit is set in reflect, then dst should be
+ * reflected reflect_dist points beyond its bounds in the corresponding
+ * direction.
+ * @param reflect_dist number of points to fill by reflection
*/
int mg2di_residual_calc(ResidualCalcContext *ctx, size_t size[2],
double *residual_max,
@@ -75,6 +81,7 @@ int mg2di_residual_calc(ResidualCalcContext *ctx, size_t size[2],
const double *rhs, ptrdiff_t rhs_stride,
const double * const diff_coeffs[MG2D_DIFF_COEFF_NB],
ptrdiff_t diff_coeffs_stride,
- double u_mult, double res_mult, int mirror_line);
+ double u_mult, double res_mult,
+ int reflect, size_t reflect_dist);
#endif // MG2D_RESIDUAL_CALC_H