* Note that the flow vectors must be scaled xy already (if x and y is shrinked
* before z)!
*/
-static int restrict(const darray_t *value, darray_t *nvalue, const darray_t *flow)
+static int restrict_array(const darray_t *value, darray_t *nvalue, const darray_t *flow)
{
int ks, wd, hd, kd;
darray_t *temp = NULL;
* substracted from the value. So the changes in the next level's value are
* prolonged to the value.
*/
-static int prolong(darray_t *value, const darray_t *nvalue, const darray_t *flow)
+static int prolong_array(darray_t *value, const darray_t *nvalue, const darray_t *flow)
{
int ws, hs, ks, wd, hd, kd;
int i, pixles;
curr_residual += residual(level, values[level], marks[level], init, nb_list[level]);
#endif
if (level < scalexyz)
- rc = restrict(values[level], values[level+1], flows[level]);
+ rc = restrict_array(values[level], values[level+1], flows[level]);
else
- rc = restrict(values[level], values[level+1], NULL);
+ rc = restrict_array(values[level], values[level+1], NULL);
if (rc < 0)
goto error;
}
/* go up */
for (level = nlevel-2; level >= 0; level--) {
if (level < scalexyz)
- rc = prolong(values[level], values[level+1], flows[level]);
+ rc = prolong_array(values[level], values[level+1], flows[level]);
else
- rc = prolong(values[level], values[level+1], NULL);
+ rc = prolong_array(values[level], values[level+1], NULL);
if (rc < 0)
goto error;
for (i = 0; i < iters; i++)