sandbox/bugs/zordering.c

    output_gfs() swap the component of tensors and vectors.

    To plot the components shown in the left side you have to set in gfsview the right side:

    • S_x_x = -S_x_y
    • S_x_y = S_x_x
    • S_y_x = -S_y_y
    • S_y_y = S_y_x
    • u_x = -u_y
    • u_y = u_x
    #include "run.h"

tensor S[];
vector u[];

int main() {
  init_grid(2);
  run();
}

event initialize (i = 0) {
  foreach() {
    S.x.x[] = 1;
    S.x.y[] = 2;
    S.y.x[] = 3;
    S.y.y[] = 4;
    u.x[] = 5;
    u.y[] = 6;
  }
}

event gfsview (i = 0) {
  FILE * fp = stdout;
  output_gfs (fp, t=t);
}