sandbox/bugs/segment.c
Segment traversal bad name
The global scalar field p used often in solvers collides with the coord p of foreach_segment. Proposal: use in foreach_segment r instead of p.
#include "utils.h"
scalar p[];
int main()
{
init_grid (8);
periodic (right);
periodic (left);
vector u[];
foreach() {
u.x[] = sin(2.*pi*x);
u.y[] = sin(6.*pi*x);
}
output_cells();
coord S[2][2] = {{{0.1, 0.1}, {0.76, 0.76}},
{{0.8, 0.1}, {0.2, 0.7}}};
for (int i = 0; i < 2; i++) {
fprintf (stderr, "%g %g\n%g %g\n\n",
S[i][0].x, S[i][0].y,
S[i][1].x, S[i][1].y);
foreach_segment (S[i], r)
for (int i = 0; i < 2; i++)
fprintf (stderr, "a %g %g %g %g %d\n", r[i].x, r[i].y, x, y, _n);
}
}